Assessing Land Dynamics and Sustainability on the Paciﬁc Coast of Nicaragua: A Method Based on Comprehensive Land Units

: In the coastal zones, varied uses converge, some of them of priority interest. In this study, an integrated method for the planning and management of the territory is proposed, which includes the evaluation of sustainability. A total of 15 different land-use classes were estimated in 80 sampling units distributed regularly along the Paciﬁc coastline of Nicaragua and classiﬁed to determine land management sectors. For each of the identiﬁed sectors, the ecological, economic, social, and productive dimensions were evaluated independently, handling a total of 53 variables from different databases, by means of ordination multivariate factor analysis. Subsequently, the four dimensions were integrated into a model and the results were evaluated based on their similarity with theoretical development scenarios, assessed by discriminant analysis. Among these, the scenarios considered as a goal for sustainability in the studied area were present. On the Paciﬁc coast of Nicaragua, productive and economic activities are currently prioritized, without having an integrated planning scheme for the entire territory, which includes nature conservation. The main contribution has been to provide a method for evaluating the land in an integrative and multidimensional way, while at the same time qualifying the different territorial sectors from a sustainable development. Even under a context of relative scarcity of information for some relevant aspects, the dimension-values assessment is largely solved by ordering the territorial sectors with a multivariate strategy, so that they are classiﬁed in relative and not absolute terms, which allows the strategy to be very useful for countries lacking some databases and cartography. This holistic and comprehensive vision of the entire territory facilitates social participation and contributes to decision-making aimed at advancing toward sustainability.


Introduction
Coastal zones can be defined as the areas where natural terrestrial and marine processes overlap and interact. They constitute unique areas of high biodiversity where varied uses converge, some of them essential for human activity, such as commerce, tourism, industry (refinery, ports), residential areas, and communities dependent on fishery. They also provide aesthetic values and areas for sports and recreation, and play a significant role in our cultural heritage. Yet, while activities such as fishery and tourism are beneficial to the economy and human well-being, they can also modify the natural stability of these areas, and place significant pressures, including loss of and damage to biodiversity, contamination, and introduction of nonindigenous species [1][2][3]. The confluence of such different uses makes coastal areas highly sensitive to the environment and triggers a conflict over how space is occupied [4,5]. With a constantly growing strategic value from different perspectives-physical-natural, economic-productive and legal-coadministrative [6]-a mounting stress on such valuable ecosystems stands out. Given the recognized importance of coastal resources and the ecosystem services they provide to humans, sustainable management of their natural resources becomes essential and constitutes a growing body of scientific literature. However, the subsequent spatial planning policies, which should be harmoniously combined with the needs arising from the management of their natural resources, remain complex and far behind.
Land planning and management to approach such an overall view are also complex. Reasons are due to the interplay of ecosystems of different nature, the necessary coherence between the broad (national or regional) and the local planning scales, and the need of an integrated approach that reconciles all demands avoiding conflicts [7][8][9]. Therefore, understanding the interaction between expectations of use in a given territory, while recognizing its multifunctional nature, is key to successful planning of humanized landscapes and the values generated by their own idiosyncrasies [10][11][12], and policies must be effective in addressing these challenges [13].
As it has recently been highlighted, there is a gap in contemporary land-planning theory with regard to the survey methodology of land uses, and rethinking the methodologies for planning and management is needed [14,15]. Reviews of the evolution of urban land-use survey methodologies during the last century reveal that methodologies have been seriously trivialized, with simplistic and empirical approaches to land-use survey methodology [16]. However, in the context of multicriteria decision analysis, numerous methods are applicable in different decision-making situations, although advice for selecting the most appropriate for the specific field of application and problem has not been thoroughly investigated [16].
Land-use planning has been the usual tool for an efficient management of spatially distributed resources, including the adaptation of land-use paradigms to changing environmental or socioeconomic circumstances [17]. However, for confronting the new challenges imposed by global change, improved tools able to plan sustainability are needed, and should include the sustainable development goals established by the United Nations [18,19]. Territorial management must ultimately facilitate adoption of measures capable of anticipating new territorial conflicts and trade-offs, including those derived from global change. Moreover, they have to be available to decision-makers, in order for them to propose responsible policies on sustainability. Furthermore, this is especially important in developing countries, as in the case of the region considered here.
In countries with greater access to information, cartographic sources, well-established planning tools, and complex methodologies have been used for land management and evaluation [20][21][22][23][24], but until now they have not evaluated the sustainability of the current and future possible actions. This is worrying, as the United Nations has recognized the 21st century as the century of sustainability. Therefore, now is the time in which it is essential to face the main social and environmental challenges arising from global change.
The aim of this study is to propose a methodology that points directly to contexts in which available territorial information to start with is limited, and there is a need for research and resource analysis. This is the case of the coastal areas in many developing countries, where complexity is high and available data are often scarce. Social and institutional stakeholders could use this information to design programs and public policies, and to evaluate their effects on the territory's capacity to provide different ecosystem services [5,11]. Moreover, our objective is to integrate territorial analysis with evaluation of sustainability, thus adding to a resulting territorial planning the possibility to evaluate how the adopted decisions fit with previous established sustainability goals. We also consider handling different scenarios of sustainability in the analysis, thus opening a final discussion on the different degrees of model fitting obtained.
To achieve these goals, we go through a two-step process. In the first phase, we establish Land-Management Sectors based on the varied types of land use, to carefully identify possible conflicts among uses. In the second phase, we evaluate the Sustainability status corresponding to the different sectors, considering the ecological, economic, social, and productive dimensions of sustainability, and contrast them with the previously constructed theoretical development scenarios, which include the scenarios that could be considered as a target for most of the studied sectors.

Materials and Methods
The study area corresponds to the coastal line of the Pacific of Nicaragua (450 km), administratively divided into five departments and 12 municipalities. It is part of the Pacific Coastal Plain, delimited by the Los Maribios volcanic mountain range ( Figure 1).
Land 2021, 10, x FOR PEER REVIEW identify possible conflicts among uses. In the second phase, we evaluate the Sustain ity status corresponding to the different sectors, considering the ecological, econ social, and productive dimensions of sustainability, and contrast them with the p ously constructed theoretical development scenarios, which include the scenarios could be considered as a target for most of the studied sectors.

Materials and Methods
The study area corresponds to the coastal line of the Pacific of Nicaragua (450 administratively divided into five departments and 12 municipalities. It is part o Pacific Coastal Plain, delimited by the Los Maribios volcanic mountain range (Figur The dominant climate is subtropical, with "Subtropical Humid and Hot Fore the Chinandega and León provinces, "Subtropical Dry and Hot" at the north and sou Lake Xolotlán, and "Tropical Dry Forest" on the southern coast [25,26]. Six prot natural areas from the Habitat/Species Management categories of the International U for Conservation of Nature (IUCN) are present in the area: four Nature Reserves and Wildlife Refuges [27].

Identification of Land-Management Sectors
We established 80 sampling units (SUs) regularly distributed every 5 km alon coast, generated from the sea-land boundary. The SUs had a circular area of 2.5 k dius (1.963 ha), a surface that allows the detection of the different land uses at a de scale.
Based on aerial orthophotos, we defined 15 types of land use: Beach, Wetland tuary, Aquaculture, Mangrove Forest, Saltworks, Artificial Lagoons, Dry Tropical F Gallery Forest, Scrub, Forest Plantations, Agriculture, Livestock, Roads, and Infras tures [26]. Land uses were expressed as percentage of the SU area covered. Orthop from 2013 were corrected by the Nicaraguan Institute of Territorial Studies (INE The dominant climate is subtropical, with "Subtropical Humid and Hot Forest" in the Chinandega and León provinces, "Subtropical Dry and Hot" at the north and south of Lake Xolotlán, and "Tropical Dry Forest" on the southern coast [25,26]. Six protected natural areas from the Habitat/Species Management categories of the International Union for Conservation of Nature (IUCN) are present in the area: four Nature Reserves and two Wildlife Refuges [27].

Identification of Land-Management Sectors
We established 80 sampling units (SUs) regularly distributed every 5 km along the coast, generated from the sea-land boundary. The SUs had a circular area of 2.5 km radius (1.963 ha), a surface that allows the detection of the different land uses at a detailed scale.
Based on aerial orthophotos, we defined 15 types of land use: Beach, Wetland, Estuary, Aquaculture, Mangrove Forest, Saltworks, Artificial Lagoons, Dry Tropical Forest, Gallery Forest, Scrub, Forest Plantations, Agriculture, Livestock, Roads, and Infrastructures [26]. Land uses were expressed as percentage of the SU area covered. Orthophotos from 2013 were corrected by the Nicaraguan Institute of Territorial Studies (INETER), with the support of Google Earth Images (free version). The size of each pixel was 1 m, and the digitization scale of the polygons was 1:1000. The percentage of sea occupation at every SU was also noted. Combinations of land uses among the coastal sampling unites were detected, which were explored by principal component analysis (PCA) [28] to interpret those main combinations found. We also performed a hierarchical classification analysis (Ward's statistic) to obtain homogeneous linear segments (consolidated coastal sections formed by SUs belonging to the same class) that we identified as components of the Land Management Sectors (LMSs)-each of the classes identified in the cluster analyses, connected by transition zones (sections formed by SUs of different classes). The statistical analyses were performed with the SPSS software and Community Analysis Package (version 3.1) for the multivariate classification and PCA, and ArcGis 10 software was used for the territorial analyses.

Multidimensional Analyses for the Evaluation of Sustainability
The sustainability analysis applied the model proposed by Gómez-Sal et al. [29]. This methodology is based on identification of relevant dimensions to assess sustainability of the land uses, and posterior comparison to reference development scenarios [10,29]. For this study, four relevant development dimensions were considered: Ecological, Economical, Productive, and Social, which were assessed through the data from 53 variables (Table 1).
An important amount of the social, economic, and productive data was obtained through surveys directed to the leaders of 30 communities located on the Pacific coastline of Nicaragua. Surveys were carried out from October to December 2013. Other relevant information came from the National Fisheries survey [30], which includes data on vessels, total catches, and number of fishermen, among others. The information for the maritime ports came from the official Report of the Nicaraguan State for 2013 [31]. Finally, the data on the ecological variables were obtained from the analysis of land uses in each of the 80 sampling units. A matrix was elaborated with this set of variables according to their role for the evaluation of the four dimensions considered, where a variable may be involved in the evaluation of more than one dimension. The continuous variables were transformed into discrete variables (0 = absent, 1 = present), developing a matrix of numerical data for each dimension in which the LMS and the set of 53 variables selected were listed. The LMSs were ordered for each dimension using detrended correspondence analysis (DCA) [32].
The data of the new variables (the coordinates of the different LMSs in the first axis of the DCA) for each of the analyzed dimensions of development were subject to a reassignment of values to facilitate their integration into the models that represent different development scenarios [29]. The new values varied between 0 and 100, so that the lowest coordinate value (including negative values) corresponded to 0, and the highest to 100. In this way, a gradient was established from low to high values in the ordination of LMSs for each of the evaluative dimensions. The results of the reallocation allowed a comparable approach (between 0 and 100) for the different dimensions. To support the consistency of the model, the independent nature of the dimensions considered was analyzed by means of a correlation between the first axes of the respective DCA. A summary scheme of the procedure can be seen in Figure 2. The continuous variables were transformed into discrete variables (0 = absent, 1 = present), developing a matrix of numerical data for each dimension in which the LMS and the set of 53 variables selected were listed. The LMSs were ordered for each dimension using detrended correspondence analysis (DCA) [32].
The data of the new variables (the coordinates of the different LMSs in the first axis of the DCA) for each of the analyzed dimensions of development were subject to a reassignment of values to facilitate their integration into the models that represent different development scenarios [29]. The new values varied between 0 and 100, so that the lowest coordinate value (including negative values) corresponded to 0, and the highest to 100. In this way, a gradient was established from low to high values in the ordination of LMSs for each of the evaluative dimensions. The results of the reallocation allowed a comparable approach (between 0 and 100) for the different dimensions. To support the consistency of the model, the independent nature of the dimensions considered was analyzed by means of a correlation between the first axes of the respective DCA. A summary scheme of the procedure can be seen in Figure 2.

Integration of the Dimensions and Comparison with Theoretical Scenarios of Ssustainable Development
In addition to the independent estimation of each dimension of the respective sections of the coastline, the objective was to obtain an integrative evaluation of the four dimensions on each section and the LMS of which each was a part. This integrated evaluation was compared to the theoretical scenarios previously established [29,33] (Table 2, Figure 3). The assignment of the sections to one or another scenario was carried out through discriminant analysis.
The statistical analyses were carried out with the statistical program PAST v3.12 [34] (free version).

Integration of the Dimensions and Comparison with Theoretical Scenarios of Ssustainable Development
In addition to the independent estimation of each dimension of the respective sections of the coastline, the objective was to obtain an integrative evaluation of the four dimensions on each section and the LMS of which each was a part. This integrated evaluation was compared to the theoretical scenarios previously established [29,33] (Table 2, Figure 3). The assignment of the sections to one or another scenario was carried out through discriminant analysis.
The statistical analyses were carried out with the statistical program PAST v3.12 [34] (free version). Modified from Gómez-Sal et al. [29], with the addition of the social dimension and adapted to the development context of the present research.  [29], with the addition of the social dimension and adapted to the development context of the present research. Table 2. Theoretical scenarios considered to estimate the degree of similarity with the models obtained at the sectors of study, and their descriptions. Based on the scenarios proposed by Gómez-Sal [29,35] and Gómez-Sal and González García [10], expanded with a new scenario (Close conservation).

Theoretical Scenarios Description
Traditional sustainable Production of a variety of products is compatible with a high natural value. Attention and work are required to be sustainable Traditional unsustainable Forced production by degrading the ecosystem Industrial unsustainable Production dependent on external inputs and subject to market cycles, implying a deep impoverishment of the ecosystem Impoverished Production leads to a degraded ecosystem, with little chance of sustaining viable and profitable activities Supported or funded Maintenance of valuable natural features which attract profitable economic activities; e.g., natural parks, protected landscapes, botanical gardens, etc.

Sound sustainable
Sensible or rational use that preserves ecosystem processes and increases sustaining capacity and natural capital. The economic value is enhanced through an increase of quality in its products: including valuable landscape, sustainable tourism, ecological agriculture, etc.
Close conservation Maintenance of natural features with a minimum of economic profitability, productivity, and social activities; for example, natural parks, nature reserves, etc.
Ideal or utopian Can only exist under special conditions; e.g., in ecosystems that benefit from the inflow of external resources (fertility, water, temporary valuable species, products of high economic value, etc.)

Land-Use Spatial Patterns
The Pacific coastal zone of Nicaragua presents a high diversity of land uses covering both terrestrial and marine-transition ecosystems. Tropical Dry Forest is present throughout the entire coastal zone (29.3% coverage in the entire area), with the highest coverage in the southern sector (with more than 40.4% for each sampling unit) (Figure 4(Aa)) and important extensions in the north (the Cosigüina Volcano protected area). Gallery Forest represents 2% of the total sampled area and is mainly distributed in the south (Figure 4(Ab)). The Scrub land-use class covers 3.3% of the total area and is located in the central part of the coastline (Figure 4(Ac)). The largest Infrastructures area is located in the South and is tourist-oriented (the municipality of San Juan del Sur; Figure 4(Ad)); therefore, it appears linked to the highest density of roads in this area. The Infrastructures in the coastal zone are mainly associated with urbanized areas. The highest coverage of housing and industry is located in the southern zone due to the touristic activity. On the other hand, in the north there is a restriction of infrastructures due to the greater presence of protected areas, mangroves, and the scarcity of road accesses. The communities that occupy more than 10 ha of the coastal zone correspond to cities with more than 1000 inhabitants. San Juan del Sur and Corinto, the country's main port, represent the highest population densities on the coast, with more than 15,000 inhabitants.
Agriculture is the second category with the highest coverage. It is distributed in a higher percentage in the center and north (Figure 4(Ae)). The areas dedicated to Livestock farming are located along the entire coast, but the greatest coverage is found in the central zone (Figure 4(Af)). Productive activities include wood and energy, mainly in the south and north (Figure 4(Ag)). In 80% of the sampling units, roads were present (Figure 4(Ah)).
Regarding the transition environments (zone of influence of the sea), Mangrove Forest is distributed in the north (Figure 4(Ba)), as well as in the Estero Real (a relevant wetland basin) and surrounding areas. Wetlands were identified in seven sampling units (1.0%; Figure 4(Bb)). Beaches are distributed in 96% of the coastline, being generally longer in the north (Figure 4(Bc)). Aquaculture, with shrimp production (Litopennaeus vannamei and L. stylirostri) is linked to the presence of mangroves, so the coverage related to this activity predominates in the northern zone, located between the municipalities of León and El Viejo (Figure 4(Bd)). Aquaculture areas that are also present in the central part of the coast are linked to breeding areas for broodstock and larval production.  The land-use analysis allowed the identification of similar territorial combinations of land uses, which were explored by PCA ( Table 3). The first component of the PCA was strongly influenced by the presence of the Mangrove ecosystem and Aquaculture activities, categories that depend on the presence of estuaries. In the second component, the categories associated with environments with a greater presence of the sea (Beach and Estuaries) predominated, accompanied by productive activities such as Aquaculture and Livestock. Finally, the third component was dominated by the Tropical Dry Forest (natural systems little intervened) compared to productive activities.  The land-use analysis allowed the identification of similar territorial combinations of land uses, which were explored by PCA ( Table 3). The first component of the PCA was strongly influenced by the presence of the Mangrove ecosystem and Aquaculture activities, categories that depend on the presence of estuaries. In the second component, the categories associated with environments with a greater presence of the sea (Beach and Estuaries) predominated, accompanied by productive activities such as Aquaculture and Livestock. Finally, the third component was dominated by the Tropical Dry Forest (natural systems little intervened) compared to productive activities. From the hierarchical classification analysis, we obtained two large groups of units related to terrestrial and marine/transitional environments ( Figure 5). The main dividing variable was the presence or absence of Mangroves. From the hierarchical classification analysis, we obtained two large groups of units related to terrestrial and marine/transitional environments ( Figure 5). The main dividing variable was the presence or absence of Mangroves. Of the 80 sample units analyzed, 54 corresponded to the terrestrial environments group, while 26 corresponded to the marine-coastal environments group. For each of the two large groups, another level of classification identified eight more detailed groups, whose main features are detailed in Table 4. Table 4. Description of the groups of land-use units. The names reflect the dominant land-use coverage. Total number of sampling units represented is also indicated.

Land Use Units Description Number
Natural (little intervened) The largest ground cover is represented by the tropical dry forest, followed by coverage related to infrastructures (housing, industry, ports), agricultural and livestock activities. This category presents beaches in all the sampled units.

16
Diversified productive It is a heterogeneous unit, with predominance of areas dedicated to livestock combined with the presence of the tropical dry forest. Agriculture and infrastructures accompanied by scrub cover and gallery forest are also important. Of the 80 sample units analyzed, 54 corresponded to the terrestrial environments group, while 26 corresponded to the marine-coastal environments group. For each of the two large groups, another level of classification identified eight more detailed groups, whose main features are detailed in Table 4. Table 4. Description of the groups of land-use units. The names reflect the dominant land-use coverage. Total number of sampling units represented is also indicated.

Land Use Units Description Number
Natural (little intervened) The largest ground cover is represented by the tropical dry forest, followed by coverage related to infrastructures (housing, industry, ports), agricultural and livestock activities. This category presents beaches in all the sampled units.

16
Diversified productive It is a heterogeneous unit, with predominance of areas dedicated to livestock combined with the presence of the tropical dry forest. Agriculture and infrastructures accompanied by scrub cover and gallery forest are also important.

Land Use Units Description Number
Mangrove ecosystems (little intervened) Dominated by mangrove ecosystem, with wide estuaries and low proportion of soil dedicated to aquaculture. 10 Mangrove ecosystem with productive activities Mosaic of types of cover, integrating the presence of mangrove and dry forest ecosystems with different productive activities such as livestock, agriculture, saltwork, and aquaculture. The presence of infrastructure is also relevant.

Estuary
Presence of the mangrove forest with the largest estuaries of all categories. Absence of wetlands, so that immediately after the mangrove forest and beaches there are fragments of dry forest, with small areas dedicated to crops. 5 The spatial distribution of the land use units along the coast is shown in Figure 6. For each of the groups, its canonical unit was identified as the one showing most similar values to the average of the group (Figure 6). Mosaic of types of cover, integrating the presence of mangrove and dry forest ecosystems with different productive activities such as livestock, agriculture, saltwork, and aquaculture.
The presence of infrastructure is also relevant.

Estuary
Presence of the mangrove forest with the largest estuaries of all categories. Absence of wetlands, so that immediately after the mangrove forest and beaches there are fragments of dry forest, with small areas dedicated to crops. 5 The spatial distribution of the land use units along the coast is shown in Figure 6. For each of the groups, its canonical unit was identified as the one showing most similar values to the average of the group (Figure 6).

Land-Management Sectors
The transference of the different land-use units to the territory allowed the identification of 21 homogeneous sections that contained land units of the same group. The length of these sections varied between 10 and 25 km, and they could include two or more land units. They were not always contiguous; in those cases, they were connected by heterogeneous zones considered transitional (Figure 7).

Land-Management Sectors
The transference of the different land-use units to the territory allowed the identification of 21 homogeneous sections that contained land units of the same group. The length of these sections varied between 10 and 25 km, and they could include two or more land units. They were not always contiguous; in those cases, they were connected by heterogeneous zones considered transitional (Figure 7).
In order to obtain homogeneous land sectors covering the whole coastline, the land units from those transitional zones were assigned to one of the contiguous sections, according to their similarity, based on the results of the ordination analysis (PCA) ( Table 3) performed as a complement to the classification for practical management. Finally, 8 sectors covered the coast, each of them dominated by one type of the land-use units. This proposal of Land-Management Sectors (LMS) (Figure 8) made it possible to simplify the complexity of the territory, creating relatively homogeneous spaces that facilitated its understanding as well as the analysis and management of the processes that converge in it, endowing the territory with different functions and meaning. In order to obtain homogeneous land sectors covering the whole coastline, the land units from those transitional zones were assigned to one of the contiguous sections, according to their similarity, based on the results of the ordination analysis (PCA) ( Table 3) performed as a complement to the classification for practical management. Finally, 8 sectors covered the coast, each of them dominated by one type of the land-use units. This proposal of Land-Management Sectors (LMS) (Figure 8) made it possible to simplify the complexity of the territory, creating relatively homogeneous spaces that facilitated its understanding as well as the analysis and management of the processes that converge in it, endowing the territory with different functions and meaning.  In order to obtain homogeneous land sectors covering the whole coastline, the land units from those transitional zones were assigned to one of the contiguous sections, according to their similarity, based on the results of the ordination analysis (PCA) ( Table 3) performed as a complement to the classification for practical management. Finally, 8 sectors covered the coast, each of them dominated by one type of the land-use units. This proposal of Land-Management Sectors (LMS) (Figure 8) made it possible to simplify the complexity of the territory, creating relatively homogeneous spaces that facilitated its understanding as well as the analysis and management of the processes that converge in it, endowing the territory with different functions and meaning.

Multicriteria Evaluation of the Land-Management Sectors
From the results for each of the evaluative dimensions, we were able to obtain different readings for the sustainability for each LMS, taking as a reference the coordinates of the 21 sections in the respective first axes of the DCAs (Figure 9). With this information, it was possible to appreciate the different importance of the processes (ecological, productive, social, and economic) considered of interest for evaluating the development sustainability in the different sectors of the coastal strip.
The highest values (blue) of the ecological dimension were presented in the north, and the lowest (red) toward the center and south of the coastline. The productive dimension showed varied high average values toward the southern half of the coast, while the highest values of the economic dimension were located in the central zone and in the extreme south. Finally, the high values of the social dimension appeared disjointed in different areas of the coastline.

Multicriteria Evaluation of the Land-Management Sectors
From the results for each of the evaluative dimensions, we were able to obtain different readings for the sustainability for each LMS, taking as a reference the coordinates of the 21 sections in the respective first axes of the DCAs (Figure 9). With this information, it was possible to appreciate the different importance of the processes (ecological, productive, social, and economic) considered of interest for evaluating the development sustainability in the different sectors of the coastal strip.
The highest values (blue) of the ecological dimension were presented in the north, and the lowest (red) toward the center and south of the coastline. The productive dimension showed varied high average values toward the southern half of the coast, while the highest values of the economic dimension were located in the central zone and in the extreme south. Finally, the high values of the social dimension appeared disjointed in different areas of the coastline.

Theoretical Scenarios of Sustainable Development
Correlations between the first axes of the evaluative dimensions were carried out to analyze their possible redundancy. Table 5 shows a low correlation between the dimensions, indicating that they could be interpreted as independent and represented as Cartesian axes of a hypothetical factorial space. The result of the discriminant analysis assigned the 21 coastal sections to five of the eight theoretical scenarios of sustainability (Table 6). Of the analyzed sections, 50% corresponded to scenarios that reflected degradation of ecosystems and loss of natural value. Likewise, the predominance of the sections corresponding to the "traditional sustainable" and "traditional unsustainable" scenarios was observed, with the latter mainly associated with intensive agriculture. The theoretical scenarios "impoverished", "utopian", and "close conservation" were not represented. The largest proportion of the coastline was classified into the traditional unsustainable scenario, in which the excesses of the production system caused the degradation of the ecosystems, compromising future well-being of local communities. In general, these were sections characterized by the intensive production of sugar cane and peanuts.
The spatial pattern distribution of the five identified development scenarios on the Pacific coastal zone of Nicaragua was heterogeneous ( Figure 10).
The largest number of sections was related to unsustainable traditional and industrial scenarios in the south, as well as to the sustainable traditional and sound scenarios in the north. The industrial unsustainability scenario was located in San Juan del Sur, Villa El Carmen, San Rafael del Sur, and Corinto.
The scenarios characterized by nature conservation and environmental services appear divided between the supported or funded and sound sustainable scenarios, with three sections each. Five of the sections corresponded to protected areas, classified as Natural Reserve and Wildlife Refuge. The only protected area that was classified into the traditional sustainable scenario was the Chacocente Wildlife Refuge (Santa Teresa).
Attending to the distribution of scenarios per municipalities, the coastal zone of El Viejo was revealed as highly complex in terms of land management due to the presence of four different scenarios (Table 7). León and San Juan del Sur contained three scenarios, The largest number of sections was related to unsustainable traditional and industrial scenarios in the south, as well as to the sustainable traditional and sound scenarios in the north. The industrial unsustainability scenario was located in San Juan del Sur, Villa El Carmen, San Rafael del Sur, and Corinto.
The scenarios characterized by nature conservation and environmental services appear divided between the supported or funded and sound sustainable scenarios, with three sections each. Five of the sections corresponded to protected areas, classified as Natural Reserve and Wildlife Refuge. The only protected area that was classified into the traditional sustainable scenario was the Chacocente Wildlife Refuge (Santa Teresa).
Attending to the distribution of scenarios per municipalities, the coastal zone of El Viejo was revealed as highly complex in terms of land management due to the presence of four different scenarios (Table 7). León and San Juan del Sur contained three scenarios, each of them with a section in the sound sustainable scenario. The rest of the municipalities were made up of one or two scenarios for their coastal lines.

General Land-Use Analysis
The Mangrove and Tropical Dry Forest were the natural ecosystems predominating in the coastal strip studied. These ecosystems are under constant pressure due to the increase in productive activities such as aquaculture, livestock, agriculture, and infrastructures [36][37][38][39], so their ecological integrity and capacity to provide ecosystem services is threatened. Even so, those landscapes managed through traditional productive activities (diversification and local ecological knowledge) maintain valuable characteristics for the formulation of more sustainable development models [29]. The main impacts leading to mangrove loss are associated with the expansion of shrimp farms, deforestation (for firewood), advancing of the agricultural frontier, urbanization, sedimentation, and water pollution [40]. The most fragile ecosystems are threatened by social and economic pressures from both the coastal zone itself, and from adjacent areas [41]. Infrastructure plays an important role in urbanized landscapes. The largest coverage of land in homes, industries, and roads is located in the southern zone due to the increase in tourist activity. On the other hand, in the north, this expansion is more limited, due to the creation of protected areas, the presence of the mangrove forest, and the scarcity of road infrastructure, also due to less accessibility from the most populated areas of the country.

Land-Management Sectors
The analysis based on a regular sampling of land use units was demonstrated to be effective in identifying land-management sectors, some of them with a multifunctional nature and others with different degrees of specialization [10]. These LMSs can also be analyzed according to their ability to provide different ecosystem services [12].
This result is an important contribution to understanding the territory, since it characterizes different sections and LMSs from a multicriteria perspective. In previous works referring to the classification and planning of the Pacific coastal zones of Nicaragua, only ecological and conservation aspects were considered, forgetting their interaction with other uses that can be relevant as social and economic drivers of change [41,42]. The identified land-management sectors allow characterizing land complexity through sections that facilitate management of the processes that overlap in the territory, and endow it with different functions and meaning. The proposal of LMS can serve as a tool for management [43], with a perspective of natural resources conservation and sustainability building. These classification-based analyses offer a disciplinary "neutral" tool, which can be used to improve communication of data and concepts, social participation, and active collaboration between multiple disciplines [44].
The characterization of coastal land by means of the regularly distributed sampling units provided detailed and accurate information, allowing management of natural and anthropic variables together [23]. Coastal areas in developing countries, with outstanding natural values, are characterized by high complexity [6,8,45] and often poor data. For this reason, we considered that the approach adopted in this work was suitable for organization, research, and analysis of data in similar contexts [44]. Heterogeneous landscapes require the formulation of comprehensive policies aimed at strengthening sustainable management with well-established goals in systems that are multifunctional [11].
The coastal zone classification through LMSs can help decision-makers to define management plans with the active participation of key stakeholders, taking into account local knowledge and the prosperity aspirations of communities. In similar coastal environments, there are experiences in which actions carried out at the local level have contributed to progress in integrated management [46][47][48].

Comprehensive Evaluation of Sustainable Development
Sustainable development does not exist as a static equilibrium with fixed limits [49]. Technological innovation and changes in social organization make sustainable development a dynamic process. An important part of the process includes learning to understand how change-drivers influence the behavior of the concerned systems (ecological, social, cultural, etc.) [50]. Transitions between sustainability scenarios [51] are not rigid boundaries. Figure 11 represents the transitions between the different sustainability scenarios.
Considering the traditional sustainable scenario as a starting reference, it is possible to move to a sound sustainable scenario, with relevant ecological and social values, or to others in which the degradation of ecosystems and their services predominate. Another possible option would be to move toward close conservation of nature, which also reduces the flow of ecosystem services to humans. This figure shows a dynamic vision of the possible options [52]. Drivers of change can come from public policies or from natural or social events. The desirable goal is to combine social progress and the prudent use of natural resources with the economic system oriented by social and ecological objectives, in such a way that they can ensure favorable living conditions for current and future generations, reaching the so-called sound sustainable scenario.

Comprehensive Evaluation of Sustainable Development
Sustainable development does not exist as a static equilibrium with fixed limits [49]. Technological innovation and changes in social organization make sustainable development a dynamic process. An important part of the process includes learning to understand how change-drivers influence the behavior of the concerned systems (ecological, social, cultural, etc.) [50]. Transitions between sustainability scenarios [51] are not rigid boundaries. Figure 11 represents the transitions between the different sustainability scenarios. Figure 11. Description of hypothetical transitions between scenarios, starting from the Traditional sustainable one.
Considering the traditional sustainable scenario as a starting reference, it is possible to move to a sound sustainable scenario, with relevant ecological and social values, or to others in which the degradation of ecosystems and their services predominate. Another possible option would be to move toward close conservation of nature, which also reduces the flow of ecosystem services to humans. This figure shows a dynamic vision of the possible options [52]. Drivers of change can come from public policies or from natural or social events. The desirable goal is to combine social progress and the prudent use of natural resources with the economic system oriented by social and ecological objectives, in such a way that they can ensure favorable living conditions for current and future generations, reaching the so-called sound sustainable scenario.
Sustainable development analyses are often carried out on a sectorial basis or with separate dimensions, without taking into account a broader spectrum of perspectives that require the inclusion of further dimensions concerning sustainability [53]. In view of some of the results of these policies, it should be pointed out as a mistake to define land-management programs taking into account partial views, which lead to the unsustainability of the human uses or to ecosystem degradation [10,29,48]. Sectorial methodologies are limited and sometimes lead to biased and short-term management actions and regulation measures. Each dimension provides different readings that, separately, can lead to errors in the characterization of development, and therefore to an inappropriate intervention proposal. The methodology proposed here integrates several dimensions relevant for land management and facilitates the adoption of policies closer to the regional and local reality [54].
The low correlation between the development dimensions considered reflects a high potential of the applied method, which depends on an correct selection of the variables used for dimension assessment [29]. Through the integration of the values of the four dimensions considered, we obtained an accurate characterization of the 21 coastal sections' sustainability, according to the similarity of the results to five of the eight possible theoretical scenarios. Assigning a certain theoretical scenario to each section of the land-management sectors allows subsequent monitoring of progress in development Sustainable development analyses are often carried out on a sectorial basis or with separate dimensions, without taking into account a broader spectrum of perspectives that require the inclusion of further dimensions concerning sustainability [53]. In view of some of the results of these policies, it should be pointed out as a mistake to define landmanagement programs taking into account partial views, which lead to the unsustainability of the human uses or to ecosystem degradation [10,29,48]. Sectorial methodologies are limited and sometimes lead to biased and short-term management actions and regulation measures. Each dimension provides different readings that, separately, can lead to errors in the characterization of development, and therefore to an inappropriate intervention proposal. The methodology proposed here integrates several dimensions relevant for land management and facilitates the adoption of policies closer to the regional and local reality [54].
The low correlation between the development dimensions considered reflects a high potential of the applied method, which depends on an correct selection of the variables used for dimension assessment [29]. Through the integration of the values of the four dimensions considered, we obtained an accurate characterization of the 21 coastal sections' sustainability, according to the similarity of the results to five of the eight possible theoretical scenarios. Assigning a certain theoretical scenario to each section of the land-management sectors allows subsequent monitoring of progress in development according to the degree of over time approximation, with respect to an established goal scenario.
A key feature of the evaluation model is to consider the production system as a separate dimension. The concept of production includes technology, management, and planning, and ultimately informs us about how resources are obtained and managed, beyond their economic value. Sustainability is not a feasible objective if one starts from the dominant economic development models [55]. It is necessary to promote new approaches to resource management, the role of the market, and institutional policies, which consider the relevance of how these activities are carried out; in short, the production and governance systems.
According to Wehbe et al. [56], each of the dimensions demands actions that are its own and others that are projected toward the rest of the dimensions, to improve human well-being. In fact, the dimensions, although theoretically independent, interact in practice, due to their hierarchical and functional relationships [29,[57][58][59]. It is essential to combine them, understanding their trade-offs and synergies, in order to advance toward human well-being, managing a comprehensive approach to sustainability [60].
The sound sustainable scenario can be a goal to be achieved for most of the analyzed coastal sections, so its continuity in space must be guaranteed; otherwise, it can lead to some sections being isolated, fragile, and exposed to degradation. This scenario is adequate to prioritize the conservation of resources, together with an increase in the economic value of marketable products, while guaranteeing people's quality of life, which has been called a win-win situation, by combining high values for different dimensions [61].
The traditional sustainable scenario is distributed mainly in the north of the coastal zone, associated with the communities of artisanal fishermen. These are sectors where natural resources are exploited in a rational way by local communities, achieving a positive effect on resource maintenance and conservation [62]. In this regard, Palacios et al. [63] suggested the need to consider development scenarios at different scales, from a larger territorial scale to that of local communities.
The classification of the coastal sections in the different development scenarios also gives a reference to the status of each of the municipalities in the coastal zone. Our proposal of eight sectors with their representation in the municipalities can serve decisionmakers when establishing and executing local strategies. The sound sustainable scenario is represented in several municipalities, and could serve as support and reference to establish or expand sustainability goals for the rest of the sections in each municipality.
The application of the proposed method to coastal zone management requires development objectives that must take into account both the zone as a whole and the particularities of each sector. The definition of these objectives by the decision-makers is an opportunity to facilitate the participation of local populations and stakeholders. All information related to the sectors should be used for this purpose. Once the objectives have been established; i.e., a sound sustainable scenario as a target, progress toward it can be evaluated periodically by analyzing each sector's progress in a multidimensional way.

Conclusions
The Pacific coast of Nicaragua is characterized by its outstanding nature conservation values. However, in recent years, production and economic activities have been prioritized with little planning. Its heterogeneity adds difficulties when establishing management plans. The results indicated important differences along the coast, with conservation predominating in the northern half, although with serious threats of degradation by aquaculture, and a greater impact of infrastructure in the southern half of the coastal fringe.
The classification of land-management sectors characterized from some selected development/sustainability dimensions is an important contribution to the elaboration of action plans directed to each sector. The results are useful both for the municipalities and for the elaboration of strategies that include the whole coastal territory.
The identification of management sectors based on four dimensions (ecological, productive, economic, and social) allows for a comprehensive approach of coastal land reality. The methodology used includes the characterization of different land sections according to their similarities with respect to different possible development scenarios. This classification helps us to understand the development processes occurring in the coastal fringe both on a local scale (communities or municipalities) and on a regional scale (the whole studied area). This sustainability status characterization can serve as an important tool for decision-makers and their interaction with local communities.
The main contribution of this study has been to provide a method for evaluating the land in an integrative and multidimensional way, while at the same time qualifying the different territorial sectors from the perspective of sustainability. Even in a context of relative scarcity of information for some relevant aspects, this was largely solved by ordering the territorial section with a multivariate strategy, so that they were classified in relative and not absolute terms, which allows the strategy to be very useful for countries lacking some databases and cartography.
According to the results, the applied development evaluation model is a valid and suitable tool in developing countries, where the lack of information in public records is frequent. In any case, the potentialities of the tool imply in-depth knowledge of the current land situation and recommending actions toward the recovery or strengthening of those dimensions with lower values, in order to approach an established goal scenario. Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.