Earth Sciences Research Journal

Estimación de atributos estructurales en un fragmento de bosque tropical natural utilizando imágenes de muy alta resolución obtenidas con sistemas aéreos no tripulados

ABSTRACT Structural attributes are fundamental biophysical parameters of forests, useful for environmental monitoring and planning. Canopy height (CH) is an important input for estimating several biophysical parameters such as aboveground biomass and carbon stocks, and can be associated with forest degradation, deforestation, emission reduction. Thus, an accurate CH estimation is a crucial issue in climate change, helping to increase biomass estimation accuracy, and support REDD+ initiatives. Very-high-resolution (VHR) imagery from unmanned aircraft systems (UAS's) have been studied as a low cost means for CH estimation at local scales, however, estimation the accuracy is a factor that determines its effectiveness. We evaluated the ability of VHR imagery from UAS's to derive structural attributes, specifically tree-crown area and height, in a tropical forest fragment located in the foothills of the Andes, in the humid tropical forests of the region known as Biogeographic Chocó in Colombia South America. We used a structure from motion (SfM) approach to derive the forest fragment's CH, and we applied mean-shift algorithms to identify single tree crowns. We performed accuracy assessment using tree height derived from field campaigns and visual interpretation of VHR imagery. Results showed a RMSE of 3.6 m of the canopy height model (CHM) with a R2 = 0.75; the total accuracy for delineating tree crowns was 73.9%. We found that using VHR imagery collected by UASs, specific trees and canopy gaps can be identified in forest fragments, which is an important step to determine forest structure.

RESUMEN Los atributos estructurales son parámetros biofísicos fundamentales de un bosque, útiles para el monitoreo y la planificación ecológica y ambiental. La altura del dosel de un árbol es un insumo importante para la estimación de varios parámetros biofísicos como la biomasa aérea y las reservas de carbono, y puede relacionarse con la degradación forestal, la deforestación y la reducción de emisiones. Por lo tanto, una estimación precisa de la altura del dosel es un tema crucial en los estudios de cambio climático y las iniciativas REDD+. Las imágenes VHR de sistemas de aeronaves no tripuladas se han estudiado como un medio de bajo costo para la estimación de la altura del dosel a escalas locales, pero la precisión en su estimación es un factor que determina su utilidad. Se evaluó la capacidad de las imágenes VHR de los sistemas de aeronaves no tripuladas para derivar atributos estructurales, específicamente el área y la altura de las copas de los árboles, en un fragmento de bosque tropical natural ubicado en las estribaciones de la Cordillera de los Andes, en los bosques tropicales húmedos de la región conocida como Chocó Biogeográfico, Sudamérica. La región es una de las áreas con mayor biodiversidad del mundo y tiene un alto nivel de endemismo, pero también tiene el mayor riesgo de pérdida de recursos naturales. Se usó un enfoque de estructura a partir de movimiento para derivar modelos de altura del dosel del fragmento de bosque, y se aplicaron algoritmos de desplazamiento medio para identificar las copas de árbol. La evaluación de la precisión se realizó utilizando datos de referencia derivados de campañas de campo e interpretación visual de imágenes VHR. El error cuadrático medio estimado de la población de errores verticales para el modelo de altura del dosel fue de 3,6 m. La precisión total para delinear las copas de los árboles fue del 73,9%. Se encuentra que utilizando imágenes VHR, se pueden identificar y monitorear fácilmente árboles específicos y brechas, lo cual es un paso importante en los programas de conservación. También se discute la utilidad de estos hallazgos en el contexto de bosques fragmentados y las compensaciones entre el precio de un sistema LIDAR y la precisión de este enfoque.

Ensayo experimental para detectar la altura de reflexión efectiva con información de la relación señal/ruido en GPS

ABSTRACT: This study aims to estimate effective reflector heights and height differences using the basic geometrical principle of multipath theory by controlling the signal quality for estimations. The geometry of the reflecting signal allows computing the effective reflector height, which is extracted from where the signal reflects on the ground and arrives at the GPS antenna phase center. To achieve that, an experimental case with two stations was conducted in the snow-free environment and GPS receivers were mounted on reflectors, which allowed to measure daily in-situ reflector heights and artificial decrement variations. The reflections from the roof surface were tracked with the first-Fresnel zones. To validate the estimated reflector heights in a controlled test environment, twelve different combinations within four simulated scenarios as a combination of decrement values have been implemented and accuracy analysis was performed. Here, a vertical shift procedure on reflectors was applied. Meanwhile, the vertical shift amount was tracked in each computation to determine, which reflected signal could be able to use for assigning reflector height as effective. Comparisons of the estimated heights and in-situ measurements show congruency with ±1.2 cm to ±8 cm accuracy. The best overall accuracy of the model among the four scenarios is computed as ±2.2 cm. When the vertical shift decrements are considered, the RMSE values are estimated within ±2.92 cm to ±3.96 cm. Although the RMSEs of the differences show a good agreement with estimated reflector heights, it is found that some reflector height estimations are statistically insignificant.

RESUMEN: Este estudio se enfoca en calcular la altura de reflexión efectiva y las diferencias de altura a través de principios geométricos de la teoría de multitrayectos al controlar la calidad de la señal para realizar los cálculos. La geometría de la señal reflejada permite computar la altura de reflexión efectiva, la cual se toma entre la señal que se refleja en el piso y la que llega a la antena GPS. Para lograr esto se realizó un ensayo experimental con dos estaciones en un ambiente sin nieve y con los receptores GPS instalados en los puntos de reflexión, lo que permite medir diariamente las alturas reflejadas in-situ y las variaciones en la reducción artificial. Las reflexiones tomadas en la superficie del techo se monitorizaron con las primeras zonas de Fresnel. Para validar las alturas de reflexión estimadas en un ambiente controlado se implementaron doce combinaciones en cuatro escenarios simulados como valores de reducción y se realizaron análisis de exactitud. En este punto, se aplicó un procedimiento de migración vertical en los puntos de reflexión. Además, la dimensión de la migración vertical se monitorizó en cada cómputo para determinar cual señal reflejada es efectiva para ser asignada como altura de reflexión. La comparación de las alturas estimadas y las medidas in-situ muestra congruencia con un rango de exactitud de ±1.2 cm a ±8. El mejor promedio de exactitud del modelo entre los cuatro escenarios se computó con un margen de ±2.2 cm. Cuando se considera la reducción en la migración vertical, los valores RMSE se estiman en el rango de ±2.92 cm a ±3.96 cm. A pesar de las diferencias de los valores RMSE muestran coincidencia con la altura de reflexión estimada se encontró que algunas estimaciones de alturas reflejadas son estadísticamente insignificantes.

Identificación de fallas con potencial sísmico de riesgo con base en análisis morfotectónico en la ciudad de Kerman (sudeste asiático)

ABSTRACT The morphotectonic analysis is a useful time-saving and cost-effective method to assess tectonic activity, especially in large regions. In this paper, the morphotectonic study was carried out in Kerman city which is in the southeast of Iran and is a part of the Central Iran structural zone. Despite the occurrence of numerous historical and instrumental earthquakes around Kerman city, there is a gap in earthquake records of this city that has obscured its seismicity status. The aim of this study is to identify active faults and the associated seismic hazard to Kerman city. For this purpose, seven geomorphic indices namely stream length-gradient index (SL), drainage basin asymmetry (Af), hypsometric integral (Hi), ratio of valley-floor width to valley height (Vf), index of drainage basin shape (Bs), index of mountain front sinuosity (Smf) and Transverse Topographic Symmetry Factor (T) were examined for 51 basins of the study area, using digital elevation model (DEM) and geological maps in a GIS environment. Then, two kinds of relative tectonic activity indices were calculated for each basin, i.e., Iat and Iat_T, from the combination of these indices. The Iat_T was chosen because its results show better agreement with the structural geology and seismic records of this region. Finally, the study area was divided into three regions according to the Iat_T values. The morphotectonic analyses indicate that the Kuhbanan fault system, especially its southern splays, has the potential to produce serious seismic hazards to Kerman city in the future; the Mahan-Jupar fault-related folds, the southern tip of the Jorjafk fault, the probable fault of Zangi-Abad, the Rafsanjan-Zarand fault system, and Gowk fault are considered as other threats to Kerman city. Furthermore, this study reveals that morphotectonic analysis is a reliable tool to evaluate fault capability and to determine the fault types, and therefore to estimate seismotectonic hazard.

RESUMEN El análisis morfotectónico es un método útil, económico y ahorrador de tiempo para evaluar la actividad tectónica, especialmente en regiones grandes. En este artículo, el estudio morfotectónico se realizó en la ciudad de Kerman, en el sudeste de Irán, y que es parte de la zona estructural de Centro Irán. Si bien han ocurrido numerosos terremotos históricos o de detección instrumental en Kerman, hay una brecha en los registros de estos movimientos que oscurecen el estatus sísmico. El objetivo de este estudio es identificar las fallas activas y el riesgo sísmico asociado para la ciudad de Kerman. Para este propósito se examinaron siete índices geomórficos (índice del gradiente del perfil longitudinal, SL; asimetría de cuenca, Af; integral hipsométrica, Hi; relación de la parte baja del valle con la altura del valle, Vf; índice de la forma de cuenca, Bs; índice de sinuosidad de frentes montañosos, Smf, y factor de simetría topográfica transversal, T) en 51 cuencas en el área de estudio con Modelos de Elevación Digital (DEM) y mapas geológicos en un ambiente de Sistemas de Información Geográfica. Luego se calcularon dos clases de índices de actividad tectónica relativa para cada cuenca, por ejemplo Iat e Iat_T, desde la combinación de estos índices. El índice Iat_T fue elegido porque sus resultados muestran concordancia con la geología estructural y los registros sísmicos de la región. Finalmente, el área de estudio se dividió en tres regiones, de acuerdo con los valores Iat_T. Los análisis morfotectónicos indican que el sistema de fallas de Kuhbanan, especialmente en su extensión al sur, tiene el potencial de producir riesgos sísmicos serios a la ciudad de Kerman en el futuro; los pliegues relacionados con la falla Mahan-Jupar, la punta sur de la falla Jorjafk, la posible falla de Zangi-Abad, el sistema de fallas de Rafsanjan-Zarand, y la falla de Gowk se consideran como otras amenazas a Kerman. Además, este estudio muestra que el análisis morfotectónico es una herramienta confiable para evaluar la capacidad de las fallas y para determinar los tipos de falla y, además, para estimar el riesgo sismotectónico.

Efectos de la altura de la vegetación y de la inmersión relativa en vegetación rígida en la corriente de un canal abierto

ABSTRACT: Vegetation commonly existing in natural rivers interact with water flow in a strongly nonlinear manner. In this research, vegetation height and a concept of relative submergence are adopted to experimentally investigate the effects of rigid submerged vegetation on flow structure. The results reveal that the relative submergences decrease along the submerged vegetation section. From the perspective of case-average, the linear relationship between water depth and vegetation height demonstrates a decrease trend of average relative submergence with the increase of vegetation height. The differences between the maximum and minimum water depths of submerged cases can be described by a positive linear function with vegetation height and a negative power function with relative submergence, by which, the raised water depths induced by vegetation block can be predicted under similar hydraulic conditions. The findings of the research, especially the introduction of relative submergence and its functional relationship with the water depth alternation induced by submerged vegetation, may improve the comprehensive understandings of the hydraulic behaviors in vegetated open channel and give some guidance to the practices of river ecological restoration.

RESUMEN La vegetación en las corrientes de un río natural interactúa con el flujo de agua de una manera no lineal. En esta investigación se estudia la altura de la vegetación y un concepto de inmersión relativa para investigar experimentalmente los efectos de la vegetación rígida sumergida en una estructura de flujo. Los resultados revelan que las submersiones relativas decrecen a lo largo de la sección de vegetación sumergida. Desde la perspectiva del caso promedio, la relación lineal entre la profundidad del agua y la altura de la vegetación demuestra una tendencia a la baja del promedio de la inmersión relativa a medida que se incrementa la altura de la vegetación. Las diferencias entre la profundidad mínima y máxima del agua en los casos de inmersión puede describirse a través de una función lineal positiva con la altura de la vegetación y una función de poder negativa con inmersión relativa, por lo que el incremento del nivel del agua ocasionado por el bloqueo de la vegetación se puede predecir con condiciones hidráulicas similares. Los hallazgos de la investigación, especialmente la introducción de la inmersión relativa y su relación funcional con la alternación de la profundidad del agua inducida por la vegetación sumergida, pueden mejorar los amplios conocimientos de los comportamientos hidráulicos en un canal abierto con vegetación y pueden ofrecer información para la restauración ecológica de ríos.

Controles estructurales de mena y determinaciones geoestadísticas de zonas de alta ley en depósitos vetiformes hospedados en zonas de cizalla, El Bagre-Antioquia, Colombia

ABSTRACT Shear zone hosted lode gold type deposits at the El Bagre district, Colombia, are the only known Permian gold mineralization in Colombia Andes. Multiple fluid flows are inferred in the formation of these veins, probably synchronous to several shearing episodes. Nonetheless, orientation and control of high-grade gold concentrations (i.e., ore-shoots) are poorly known. In the El Bagre district, three main lode systems are studied with N-NNW-S-SSE orientation and variable dips (40 to 80° to the La Ye and 30 to 45° to Mangos and Cordero). Geostatistical analyses with ellipses determined by variogram modelling, reveal the anisotropy of gold mineralization and the orientation of ore-shoots. The long axes of the anisotropy ellipses trend SSE, certainly related to sinistral shearing and NNW-SSE faulting during D3 deformation event. This study evidences the use of a simple and effective tool for lode gold deposits exploration.

RESUMEN Los depósitos vetiformes de oro hospedados en zonas de cizalla en El Bagre, al noreste del departamento de Antioquia, son la única evidencia de mineralización de oro Pérmica hasta la fecha en los Andes Colombianos. Se han inferido múltiples flujos de fluidos en la formación de estas vetas, probablemente sincrónicos a varios episodios de cizallamiento. No obstante, la orientación y el control de los altos tenores (ore-shoots) de oro son pobremente conocidas. En el distrito El Bagre, tres sistemas de veta principales fueron estudiados con una orientación N-NNW-S-SSE y con buzamientos variables (40 a 80° para La Ye y de 30 a 45° para Mangos y Cordero). Los análisis geoestadísticos con elipses determinadas por el modelamiento variográfico, revelan la anisotropia de la mineralización de oro y la orientación de los ore-shoots. Los ejes principales de los elipsoides de anisotropia tienden al SSE, ciertamente relacionados con el cizallamiento si-nestral y el fallamiento NNW durante el evento de deformación D3. Este estudio evidencia el uso de una herramienta simple y efectiva para la exploración de depósitos de oro vetiformes.

Mineralogía de arcillas y señales geoquímicas como indicadores rápidos de procesos de meteorización química diferencial en una litosecuencia silícea de perfiles de suelos

ABSTRACT This paper presents detailed mineralogical results together with a geochemical characterization for a sequence of six natural soil profiles. Bedrock samples (R series) and overlying soil samples (S series) were characterized. The soil profiles are distributed in a series of Paleozoic lithological units from lower Ordovician to upper Carboniferous in age (Iberian Massif, NW Iberia). The lithological influence on mineral properties and geochemical composition and, how different weathering may be occurring under very similar temperate and acidic conditions, have been studied. Field observations together with laboratory analyses were indicative of differential weathering. So, a series of selected chemical indices and relations were applied to clarify this assumption. The mineralogy was analysed by Scanning electron microscopy (SEMEDS), X-ray diffraction (XRD) of rock powder and soil oriented aggregates. X-ray fluorescence spectrometry (XRF) and inductively coupled plasma mass spectroscopy (ICP-MS) were applied to analyse chemical composition. The first results showed how major elements, SiO2, Al2O3 and Fe2O3, slightly enriched in the soil profiles, are consistent with the dominant mineralogy: quartz, chlorite, muscovite and/or illite, together with kaolinite and albite. The bases K2O, Na2O, CaO and MgO are also coherent with mineral composition and experience little variation, but are gradually removed in the profiles. The mobility of major elements leads to a general loss of bases and, in general, a slight enrichment in silica and sesquioxides. SiO2 is enriched, firstly accumulated in soils and partially depleted by dissolution as colloidal form. Al2O3 in some soils is slightly less than in former rocks, so other physical processes are expected to take place, involving clay removal with consequent aluminium depletion too. A special emphasis has been given to albite coexisting with kaolinite, firstly supposed to be directly inherited from parent substrates when present, but finally the chemical index PIA shows it was mainly due to mineral alteration of plagioclases. The best correlations to explain the albite alteration and kaolinitization progress were obtained with chemical indices PIA, CIW, CIA and Al2O3/Na2O ratio. This together with mineralogical signatures, suggest that kaolinite is the result of gradual dissolution due to the acid hydrolysis of albite in such acidic environments, which may also be attributed to the organic matter influence

RESUMEN Este trabajo presenta resultados mineralógicos detallados junto con una caracterización geoquímica para una secuencia de seis perfiles de suelos naturales. Se caracterizaron muestras de rocas (serie R) y muestras de suelos suprayacentes (serie S). Los perfiles de suelo se distribuyen en una serie de unidades litológicas paleozoicas de edad desde el Ordovícico inferior hasta el Carbonífero superior (Macizo Ibérico, NW de Iberia). Se ha estudiado la influencia litológica en las propiedades mineralógicas y en la composición geoquímica y cómo puede ocurrir una alteración diferente en los perfiles bajo condiciones templadas y ácidas muy similares. Las observaciones de campo junto con los análisis de laboratorio fueron indicativos de alteración diferencial. Por lo tanto, se aplicaron una serie de índices químicos de meteorización y algunas relaciones seleccionadas para aclarar esta suposición. La mineralogía se analizó mediante microscopía electrónica de barrido (SEM-EDS), difracción de rayos X (XRD) de muestras de polvo de roca y agregados orientados del suelo. Se aplicaron espectrometría de fluorescencia de rayos X (XRF) y espectroscopia de masas de plasma acoplado inductivamente (ICP-MS) para analizar la composición química. Los primeros resultados mostraron cómo los elementos mayoritarios, SiO2, Al2O3 y Fe2O3, ligeramente enriquecidos en los perfiles de suelo, son consistentes con la mineralogía dominante: cuarzo, clorita, moscovita y/o ilita, junto con caolinita y albita. Las bases K2O, Na2O, CaO y MgO también son coherentes con la composición mineral y experimentan poca variación, pero se eliminan gradualmente en los perfiles. La movilidad de los elementos mayoritarios lleva a una pérdida general de bases y, en general, un ligero enriquecimiento en sílice y sesquióxidos. El SiO2, primero se acumula en los suelos y se pierde parcialmente por disolución en forma coloidal. El Al2O3, en algunos suelos es ligeramente menor que en las rocas, por lo que se espera que parte se pierda también en forma coloidal o tengan lugar otros procesos físicos que involucran la remoción de la arcilla, con la consiguiente pérdida del aluminio también. Se ha dado un énfasis especial a la presencia de albita que coexiste con la caolinita; en primer lugar, se supone que se heredó directamente de los materiales originales en los que está presente, pero finalmente el índice PIA muestra que se debe principalmente a la alteración de las plagioclasas. Las mejores correlaciones para explicar esta alteración de la albita y el progreso de la caolinita, se obtuvieron con los índices químicos PIA, CIW, CIA y la relación Al2O3/Na2O, que junto con firmas mineralógicas, sugieren que la caolinita es el resultado de la disolución gradual debido a la hidrólisis ácida de la albita en dichos ambientes ácidos, que también puede ser atribuido a la influencia de la materia orgánica.

Efectos de la acidez del agua y las propiedades ingenieriles en la durabilidad de la roca

ABSTRACT Fifteen sedimentary rock samples from four rock types, including sandstone, limestone, travertine, and conglomerate, Keywords: Rock; Physical properties; Mechanical were collected from northern Iran's Damghan area. Mineralogical studies and laboratory experiments were performed to assess the rocks' physical properties, Schmidt hardness, P-wave velocity, slake-durability index, uniaxial compressive, point load, Brazilian tensile, and block punch strengths. The studied rock samples are mainly composed of quartz, calcite, and dolomite with different textures. The slake-durability test was carried out up to ten cycles in fluids with different pHs. The utilized test fluids were natural water with a pH of 7 and sulfuric acid solutions with pH of 5.5 and 4. Based on the results, the slake-durability index is affected by the pH of the test fluids. Also, the different rock types had distinctive slaking behaviors. The decreasing rate of the slake-durability index at initial cycles is higher than the end cycles. Regression analyses showed that the evaluated parameters are correlated to the slake-durability index. In other words, the slake-durability index of the studied rock samples is greatly affected by their mineral composition, texture, physical properties such as porosity, and it is closely related to the strength parameters of the rocks. These preliminary results led to the extraction of empirical equations for determining essential characteristics of the rocks from the slake-durability index.

RESUMEN Para este trabajo se recolectaron quince muestras tipo sedimentarias de cuatro clases de rocas, donde se incluyen areniscas, calizas, travertinos y conglomerados, en el área Damghan, en el norte de Irán. Estudios mineralógicos y experimentos de laboratorio se realizaron para medir las propiedades físicas de las rocas, el índice de dureza de Schmidt, la onda de presión, el índice de durabilidad, el de compresión uniaxial, el punto de carga, la prueba brasileña de tracción y la resistencia de golpe de bloque. Las muestras de roca estudiadas se componen principalmente de cuarzo, calcita y dolomita con diferentes texturas. El índice de durabilidad se realizó en diez ciclos con fluidos de diferente alcalinidad. Los fluidos utilizados para el índice fueron agua natural con pH de siete y soluciones de ácido sulfúrico con pH 5.5 y 4. Con base en los resultados, el índice de durabilidad se afecta por el pH de los fluidos del ensayo. También los tipos de las rocas tienen diferente comportamiento en la durabilidad. La caida en el índice de la durabilidad de los ciclos iniciales es mayor que en los ciclos finales. Los análisis de regresión muestran que los parámetros evaluados están correlacionados con el índice de durabilidad. En otras palabras, el índice de durabilidad de las muestras de rocas estudiadas está ampliamente afectado por su composición mineral, textura, y propiedades físicas como la porosidad, y está relacionado de cerca con los parámetros de resistencia de las rocas. Estos resultados preliminares condujeron al diseño de ecuaciones empíricas para determinar las características esenciales de las rocas en el índice de durabilidad.

Análisis amplio de las características de la "línea roja" cultivable en Donggang, China

ABSTRACT: In the face of endless idle and even abandoned farmland, it is trying to protect farmland resources and ensure food supply by determining the bottom line (red line) of arable land that cannot be changed in use. The red line has landscape characteristics and functions, but there remain limited understandings of how to judge the dominant landscape characteristics of cultivated land. Therefore, this study divides landscape characteristics into aesthetic value and stability according to landscape pattern. A comprehensive evaluation system was developed to explore the dominant landscape features of the cultivated resources based on a multi-factor comprehensive evaluation. On this basis, this study conducts an empirical study in Donggang, a coastal city in China. The results show that cultivated land resources with high aesthetic value are distributed near cities, major roads, and important natural resources. The cultivated land resources with greater continuity and compatibility with adjacent plots are relatively stable, mainly paddy fields. The red line of cultivated land should focus on the characteristics of cultivated land stability. At the same time, the aesthetic value of cultivated land resources is also an important factor (related degree is 0.852). According to this evaluation system, the dominant area of arable land protection in Donggang City is 88,027.21hm2, which is the preferred area for the red line of arable land.

RESUMEN: Con el fin de resolver el alto índice de vacancia o incluso el abandono de tierras de cultivo, en este trabajo se busca proteger los recursos de tierras cultivables y asegurar el suministro de comida al determinar el margen inferior (línea roja) del terreno apto para sembrar al cual no se le debe cambiar el uso. Para efectos de este estudio, se estableción la línea roja en 1800 millones de mu (cada mu corresponde a 1/15 de hectárea). La línea roja tiene funciones y características paisajísticas pero deben estar definidos los criterios de cómo juzgar las características dominantes paisajísticas de las tierras cultivadas. Además, en este estudio se dividen las características paisajísticas entre valores estéticos y de estabilidad de acuerdo con los patrones del paisaje. Se desarrolló un sistema amplio de evaluación para explorar las características dominantes del paisaje en los terrenos cultivados con base en la evaluación amplia multifactorial. Con base en este enunciado, este trabajo consiste en un estudio empírico en la ciudad costera de Donggang, China. Los resultados empíricos muestran que las tierras cultivadas cerca de las ciudades, carreteras principales y reservas naturales tienen un valor mayor estético, cerca de 85440 hm2, y que corresponden a un 81.80 % de las tierras cultivadas. Las tierras adyacentes cultivadas y las tierras altamente compatibles para el cultivo en las planicies tienen una alta estabilidad, y el tipo de paisaje dominante es el arrozal, con una área de unos 91682 hm2, que corresponde al 87.77 %. Los análisis de las características del paisaje con las ventajas compuestas de las tierras cultivadas muestra que el área de los recursos cultivados tanto con las ventajas estéticas y de estabilidad es de cerca de 82.653hm2, que corresponde al 79,13%. Este es el límite sugerido para establecer la línea roja de protección a la tierra cultivable.