Abstract. - The properties of grazed grassland soils are frequently degraded in tropical areas. Theobjective of this work was to review the literature on the theory and results of sampling and analysis ofthe spatial variability of grassland soil properties worldwide and in Ecuador. The methodology, results,and discussion of 15 articles were reviewed, available in the databases: Scopus, SciELO, ScienceDirect,Scimago Journal & Country Rank, Dialnet, and SpringerLink, and the respective theory of 5 books. It wasfound that worldwide there are soil sampling works that evaluate the spatial variability of their properties,while at the national level, the samples do not consider such variability. Therefore, it is crucial to conductresearch in Ecuador on the spatial variability of soil properties with grasslands, allowing reliable decisionsto be made for their management and conservation.Keywords: Semivariogram, kriging, soil spatial variability, grassland.Resumen. - Las propiedades de los suelos con pastizales de pastoreo frecuentemente se degradan en laszonas tropicales. El objetivo de este trabajo fue revisar la literatura de la teoría y resultados del muestreo yanálisis de la variabilidad espacial de las propiedades del suelo de pastizales a nivel mundial y en Ecuador. Serevisó la metodología, resultados y discusión de 15 artículos del tema, disponibles en las bases de datos:Scopus, SciELO, ScienceDirect, Scimago Journal & Country Rank, Dialnet y SpingerLink; y la teoría respectivade 5 libros. Se encontró que a nivel mundial existen trabajos de muestreo de suelos que evalúan lavariabilidad espacial de sus propiedades; mientras que a nivel nacional los muestreos no consideran dichavariabilidad. Es importante llevar a cabo en Ecuador investigaciones de la variabilidad espacial depropiedades del suelo de pastizales, que permitan tomar decisiones fiables para su manejo y conservación.Palabras clave: Semivariograma, krigeado, variabilidad espacial del suelo, pastizal.ISSN-E: 2697-3650Minerva JournalSaquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesSaquicela Rojas Rodrigo Alberto https://orcid.org/0000-0003-2608-3150rodrigo.saquicela@ute.edu.ecMaestrante del Instituto de Posgradode la Universidad Técnica de ManabíPortoviejo-EcuadorUniversidad UTESanto Domingo-Ecuador Sampling the spatial variability of grassland soil properties138Recibido (11/07/2022), Aceptado (17/09/2022)Vol.3, Spec. Issue N°1, (pp. 138-149) Muestreo de la variabilidad espacial depropiedades del suelo de pastizalesUseche Castro Lellyhttps://orcid.org/0000-0002-4294-9009lelly.useche@utm.edu.ecUniversidad Técnica de ManabíPortoviejo-EcuadorGonzález Pedraza Ana Francisca https://orcid.org/0000-0002-4392-3724ana.gonzalez2@unipamplona.edu.coUniversidad de PamplonaPamplona-Colombiahttps://doi.org/10.47460/minerva.v1iSpecial.89
I. INTRODUCCIÓN El suelo es un conjunto dinámico formado por partes físicas, químicas y biológicas que básicamentemantienen la productividad biológica y conservan la calidad del ambiente [8]. Dichas partes son propiedadesque varían espacialmente en el suelo debido al manejo, clima y biota [19]. El pastoreo en un manejo queinfluye en sus propiedades, como la compactación que está relacionada con el pisoteo del ganado [6]. Así, anivel mundial hay investigaciones de la variabilidad espacial de propiedades del suelo de agricultura ypastizales [7], [12], [13], [17]. En Ecuador se desconoce la situación actual del muestreo de las propiedades de los suelos con pastizalesde pastoreo, con fines de planificar diseños de muestreos que evalúen la variabilidad espacial y predigan conprecisión el valor esperado de estas propiedades. El objetivo de esta revisión de literatura fue conocer el estado actual del diseño y análisis de muestreo de lavariabilidad espacial de propiedades de suelos de pastizales y el comportamiento de la variación de estas.Esta revisión se basa en el hecho de que las propiedades del suelo varían principalmente por el manejoagronómico, clima, biota y tipo de suelo [19]. La medición de propiedades del suelo con predicciones precisaspermite tomar decisiones fiables para su manejo y conservación. De esta forma, existen herramientas paramedir y analizar la variabilidad espacial de propiedades del suelo, como el muestreo basado en el modelo,semivariograma y krigeado [2]. Se analizaron diversas publicaciones científicas relacionadas con la medición de la variabilidad espacial delas propiedades del suelo de pastizales a nivel mundial y en Ecuador. También se examinaron libros acercadel muestreo y análisis de la variabilidad espacial del suelo. Se resumió la información encontrada. El presente artículo tiene la siguiente estructura: resumen, traducción del resumen, introducción, desarrollode la revisión de la literatura, metodología, resultados, conclusiones y referencias.ISSN-E: 2697-3650Minerva Journal139II. DESARROLLO A. Propiedades del suelo El suelo es una unidad dinámica con componentes físicos, químicos y biológicos que interactúan de formacompleja y definen su funcionamiento para variados fines. El funcionamiento se refiere a la calidad del suelo,que es su capacidad para impulsar la salud de animales y plantas, conservar la calidad ambiental y mantenerla productividad biológica [8]. Los componentes del suelo son materia con cualidades y atributos denominados propiedades, que sirvenpara distinguir entre muestras de suelo [5], [8]. El suelo tiene propiedades físicas, químicas y biológicas. Lasfísicas no cambian la composición de la muestra mientras se mide u observa la propiedad [8]. Las químicasson los cambios en la composición de una muestra en condiciones de reacción química [5]. Finalmente, lasbiológicas son cambios causados por los organismos en las propiedades físicas y químicas [8]. Entre las propiedades físicas del suelo están la conductividad hidráulica, retención de agua, porosidad ypropiedades mecánicas. Entre las químicas se encuentran el potencial de hidrógeno (pH), estado redox,retención de iones en el suelo, salinidad y sodicidad. Por último, las biológicas se consideran lamineralización potencial de nitrógeno o carbono, tasa de respiración, lombrices de tierra, biomasabacteriana, diversidad bacteriana, presencia de patógenos [8]. Por otro lado, las propiedades del suelo tienen variación y correlación espacial debido al manejoagronómico, biota, tipo de suelo y clima [19]. Las observaciones más cercanas entre son frecuentementemás parecidas que las más lejanas [20]. Se ha reportado investigaciones de variabilidad espacial depropiedades físicas, químicas y biológicas del suelo [1], [6]. Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
140ISSN-E: 2697-3650Minerva Journal Medir con precisión las propiedades del suelo ayuda a su diagnóstico y toma de decisiones para sumanejo y conservación. Para tal efecto, se han desarrollado herramientas para el diseño del muestreode la variación espacial de las propiedades del suelo, tales como el muestreo basado en el diseño, queno considera la variación espacial; y el basado en el modelo, que considera la variación espacial [20]. B. Muestreo del suelo basado en el diseño Existen dos enfoques para el muestreo de propiedades del suelo: los basados en el diseño y losbasados en el modelo. El muestreo basado en el diseño sirve para muestrear con objetividad y sinsesgos. La inferencia se basa completamente en el diseño del muestreo y la independencia obtenidacon la selección aleatoria de las muestras. Cada punto de muestreo en el espacio y tiempo se considerafijo, así como el valor de la variable que se obtenga en dicho punto [19], [18]. Con respecto a la inferencia, se infiere parámetros poblacionales (como media y varianza) únicos(¿cuánto es su valor?) para todo el dominio espacial, siempre que el tamaño de la muestra seaadecuado, con los supuestos de independencia espacial (que no exista correlación espacial de lasvariables) y que la población (de datos) esté distribuida idénticamente en el dominio espacial. Losestimadores de los parámetros son insesgados [18], [19], [20]. De esta forma, el muestreo basado en el diseño, para poblaciones idénticamente distribuidas eindependientes, utiliza los diseños de muestreos aleatorios clásicos: simple, sistemático, estratificado,conglomerados, además del muestreo aleatorio en dos pasos. El último consiste en la elección aleatoriade estratos y elección aleatoria de muestras dentro de los estratos elegidos. En los muestreos simple,sistemático y conglomerados se supone la independencia y población idéntica dentro del dominioespacial; mientras que en el estratificado y de dos pasos la suposición es dentro de los estratos. Noobstante, la variación y correlación espacial de las variables está presente en el dominio espacial [20],lo que no permitiría muestreos fiables basados en el diseño. C. Muestreo del suelo basado en el modelo El muestreo basado en el modelo permite muestrear sin el supuesto de independencia espacial, portanto, no es necesaria la aleatorización para elegir cada punto de muestreo [19]. Los valores de lavariable en cada punto de muestreo se consideran aleatorios y representa un muestreo único, lapoblación, de infinitos muestreos posibles a un universo representado por un modelo estocástico (convariables aleatorias que cambian entre instantes de muestreo), la superpoblación. De esa manera, laaleatoriedad del diseño se genera por un conjunto de modelos estocásticos [18]. Además, se asume laindependencia espacial al muestrear variables aleatorias de un universo estocástico. Un modeloestocástico usado es el krigeado [19], [20]. En lo referente a la inferencia, las variables se suponen aleatorias en los puntos de muestreo de undiseño basado en el modelo, así, la verdadera media espacial también es aleatoria. Por tanto, sepredice (¿dónde están los valores?) dicha media en el dominio espacial, no se estima con un único valor[18]. De esta manera se puede predecir valores de la variable en puntos no muestreados, elaborarmapas de la variación de las propiedades del suelo y estimar parámetros del modelo estocástico [20]. Por otro lado, existen los siguientes objetivos a cumplir cuando se muestrea con base en el diseño:minimizar la estimación de la varianza del error, cubrir por igual el dominio espacial y cubrir por igualen el espacio de las características [20]. Para minimizar la estimación de la varianza del error se usa elkrigeado, que estima imparcialmente los valores de una variable en puntos no muestreados con unmínimo error de estimación cuadrático medio. También existen las medias de superficie con nohomogeneidad [20].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
Por otro lado, para cubrir por igual el dominio espacial se asignan los puntos de muestreo de talforma que cubran el dominio lo más uniformemente posible por medio de una cuadrícula finita, segúnel criterio de la mínima distancia promedio entre los puntos adyacentes de la cuadrícula aplicada a sucentro [20]. De esta forma, la cuadrícula triangular equilátera cumple esta condición para un dominiosin restricciones en asignar puntos de muestreo [11]. También existen otras opciones para la distanciamínima, como la media ponderada de las distancias más cortas y la distancia media cuadrática basadaen polígonos de Thiessen. Otras formas de cuadrícula usadas en la práctica son: muestreo encuadrícula, muestreo en transecto, muestreo secuencial y muestreo anidado [20]. Para cubrir por igual el espacio de las características del dominio espacial se distribuye las unidadesde muestreo de una forma tal que permita obtener la mejor distribución posible de la población. Parael efecto, se aplican dos criterios principales: el muestreo de hipercubo latino y el de Warrick-Myers. Encaso de desconocer distribuciones de la población previas al muestreo se utiliza la experiencia en casossimilares o datos auxiliares [20]. En lo referente al tamaño de la cuadrícula, para el muestreo de la variabilidad de propiedadesquímicas del suelo, varía entre 1 m a 100 m por lado. Los puntos de muestreo se ubican principalmenteen el centroide o en las esquinas de las cuadrículas. También pueden estar distribuidos por toda lacelda [19]. Por otro lado, la cantidad de puntos de muestreo para estimar semivarianzas (base paraestimar el modelo estocástico) dentro de límites de confianza aceptables se sugiere al menos entre 100a 150 puntos de muestreo [11]. Así mismo, se han desarrollado técnicas estadísticas para el análisisespacial de datos, como el semivariograma, para analizar la estructura de la variación espacial depropiedades del suelo; y el krigeado, que predice las variables espaciales en función delsemivariograma [2]. D. Análisis estructural El análisis estructural consiste en encontrar una función que explique la estructura de la dependenciao correlación espacial de las variables medidas. La función es la base para la predicción espacial de lasvariables y se denominan funciones de covarianza (covariogramas) y semivariogramas. La función decovarianza empírica se indica en (1), donde es la covarianza empírica para la clase de distancia deintervalo; es el número de pares de distancias separadas por la distancia h; y y sonlos valores de la muestra medidos en la posición i e i+h, respectivamente; es la verdadera media de lapoblación y se estima con el promedio de ; N= es el número total de datos [10].141ISSN-E: 2697-3650Minerva Journal La función del semivariograma empírico usualmente se calcula con el estimador de momentos deMatheron (2), donde es el semivariograma empírico para la clase de distancia de intervalo [2].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
142ISSN-E: 2697-3650Minerva Journal No obstante, cuando la varianza residual es pequeña, el semivariograma se puede calcular de formamenos sesgada con el estimador de Cressie-Hawkins [2]. El semivariograma se completa al modelar elsemivariograma en función de la distancia , con los modelos exponencial (3), esférico (4), gausiano (5);donde es la ordenada de origen del semivariograma (pepita o nugget) que representa la varianza totalde los errores aleatorios o sin estructura espacial y los errores por no medir con menor distancia deintervalo ; es la varianza estructural o meseta del semivariograma (sill); es el rango e indica la distancia desde el origen hasta cuando se llega a la meseta. El rango práctico es la distancia desde el origen hastacuando la semivarianza el 95 % de la varianzal total; siempre que el semivariograma llegueasintóticamente a la meseta; indica la varianza umbral y a partir de ahí las observaciones sonindependientes [10]. Además, hay otros modelos como el exponencial-coseno, de onda, De Wijsian,cúbico, pentaesférico, pepita puro [10]. Los métodos de ajuste de los datos al modelo son por mínimoscuadrados y modelos lineales mixtos [2].No obstante, cuando la varianza residual es pequeña, elsemivariograma se puede calcular de forma menos sesgada con el estimador de Cressie-Hawkins [2]. Elsemivariograma se completa al modelar el semivariograma en función de la distancia , con los modelosexponencial (3), esférico (4), gausiano (5); donde es la ordenada de origen del semivariograma (pepita onugget) que representa la varianza total de los errores aleatorios o sin estructura espacial y los errorespor no medir con menor distancia de intervalo ; es la varianza estructural o meseta del semivariograma(sill); es el rango e indica la distancia desde el origen hasta cuando se llega a la meseta. El rangopráctico es la distancia desde el origen hasta cuando la semivarianza el 95 % de la varianzal total;siempre que el semivariograma llegue asintóticamente a la meseta; indica la varianza umbral y a partirde ahí las observaciones son independientes [10]. Además, hay otros modelos como el exponencial-coseno, de onda, De Wijsian, cúbico, pentaesférico, pepita puro [10]. Los métodos de ajuste de los datosal modelo son por mínimos cuadrados y modelos lineales mixtos [2]. La varianza estructural relativa es la razón entre la varianza estructural y la varianza umbral e indica elgrado de estructuración espacial. Mientras más alto sea el valor de la varianza estructural relativa, máseficientes serán las predicciones estadísticas. Se clasifica en bajo, < 25 %; medio, 25 % a 75 %; alto, > 75 %[2]. E. Predicción espacial La predicción de la variable en el dominio espacial (campo, de terreno, donde se muestrea) se denominapredicción espacial. El método de predicción común es el krigeado (kriging) que tiene su base en elsemivariograma empírico. El krigeado genera el valor esperado con el mejor estimador lineal insesgado y elerror de estimación (varianza krigeada). Existe el krigeado puntual y en bloques. El puntual se usa parapredecir la variable en lugares (puntos) no muestreados mucho menores que las distancias espaciales;mientas que el krigeado en bloques es para predecir el promedio de la variable en una superficie mayor queel puntual, por medio de estimaciones puntuales en dicha superficie [2].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal Entre el krigeado puntal existen el krigeado ordinario (6), simple (7) y universal (8) para la predicciónespacial, donde es la variable predicha en cualquier nuevo punto; son los pesos asignados a cadaobservación con suma de uno; es función de las coordenadas espaciales. El ordinario y simple es paravariables con distribución normal y el universal para variables con media que depende de la varianza. Elkrigeado ordinario se utiliza para estimar localmente la media de la variable; el krigeado simple, cuandose conoce la media poblacional de la variable; el universal, para estimar la influencia espacial de los datoscon las coordenadas espaciales [2].143III. METODOLOGÍA La información de muestreo de variabilidad espacial de las propiedades del suelo de pastizales se hizo conel motor de búsqueda Google Académico. Las palabras clave usadas fueron las siguientes: “spatial variabilitysoil physical”, “kriging”, “gaussian secuantial simulation”, “soil degradation”, “soil properties spatial variationgrass”, “soil sampling”, “grid soil sampling” y “variabilidade espacial do solo pastagem”. Se eligieron 20documentos científicos, de los cuales 15 eran artículos y 5 eran libros. Los artículos que generó la búsquedaestuvieron en bases de datos indexadas como Scopus, SciELO, ScienceDirect, Scimago Journal & CountryRank, Dialnet, SpingerLink. Los libros consultados fueron publicados por editoras como Elsiever, Springer,Pearson, Brujas, Mundi-Prensa y Wiley. Los artículos elegidos se revisaron principalmente las metodologías,resultados y conclusiones, y se escribieron los respectivos resúmenes. En los libros se revisaron conceptosde propiedades del suelo y de muestreo de la variabilidad espacial de dichas propiedades. También existe el krigeado indicador para variables binarias o variables continuas transformadas abinarias; el krigeado disyuntivo (y el krigeado indicador) para interpolación no lineal, el krigeadointrínseco para variables heterogéneas con varias medias; kriging recursivo con el filtro Kalman ykrigeado con la máxima entropía bayesiana que tiene menos restricciones e incorporar másconocimiento previo [10]. En cuanto al krigeado en bloques el más común es el krigeado ordinario enbloques, donde es la variable predicha en el dominio centrado en el punto de interés [2].IV. RESULTADOS Se presentan los resultados de la revisión de literatura de diseño de muestreo y análisis de la variabilidadespacial de propiedades del suelo de pastizales. A. Diseños de muestreo de la variabilidad espacial de propiedades del suelo a nivel mundial En Portugal se diseñó un muestreo para medir la variabilidad espacial de la textura, humedad, pH, materiaorgánica, nitrógeno, fósforo y potasio del suelo de un pastizal cultivado y con pastoreo de ovejas y ganadobovino. Se tomaron 76 muestras compuestas en un terreno de 6 ha dividido en cuadrículas de 28 m por 28m. Las muestras se tomaron de 0-30 cm de profundidad del suelo. La muestra compuesta, para obtener unamuestra por cuadrícula, consistió en cinco submuestras, cuatro muestreadas en las esquinas de lacuadrícula y una de su centro. Los datos se analizaron con estadística descriptiva, correlación y mapas devariabilidad espacial con un interpolador de distancia al cuadrado inverso [17].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal144 En otra investigación llevada a cabo en Brasil en un pasto cultivado y con pastoreo de bovinos se muestreópropiedades del suelo para medir la variabilidad espacial y elaborar mapas de fertilidad del suelo. Semuestreó el pH, contenidos de carbono orgánico y arcilla, fósforo, calcio y magnesio disponibles, potasioaluminio e hidrógeno intercambiables, capacidad de intercambio catiónico efectiva, saturación de bases yconductividad eléctrica. El pastizal tenía 5,3 ha divido en 270 potreros donde se tomaron 73 muestras desuelo, excepto para la conductividad eléctrica, que fueron 4 794. Se utilizó el muestreo por zonas de manejo,que consiste en delimitar zonas de similares rendimientos del pasto y presencia de malezas, dentro de lascuales se tomó una muestra compuesta de al menos 10 submuestras. No se hizo el muestro por cuadrículas,no se indica las dimensiones de las zonas de muestreo. La profundidad de muestreo fue de 0-20 cm deprofundidad. Se usó el semivariograma ajustado a modelos empíricos y programas de sistemas deinformación geográfica para elaborar mapas de fertilidad del suelo. No se indica la cantidad total demuestras para la conductividad eléctrica [7]. Así mismo, en Brasil, se muestreó un pastizal cultivado (degradado) para pastoreo las propiedades del suelopara evaluar su variabilidad espacial. Se muestreó el contenido de arena, limo, arcilla, de 0-20 cm deprofundidad, la densidad aparente de 0-10 cm y de 10-20 cm, resistencia a la penetración cada 5 cm deprofundidad de 0-30 cm, permeabilidad e infiltración tridimensional de 0-20 cm. El pH, potasio, calcio,magnesio, fósforo, aluminio, materia orgánica, suma de bases, acidez potencial, capacidad de intercambiocatiónico, saturación de bases y de aluminio, se muestrearon de 0-10 cm y de 10-20 cm. El pastizal tenía 3,67ha y se tomaron 50 muestras de suelo cada 30 m aproximadamente, con espaciados irregulares (nocuadrículas) entre puntos de muestreo (definidos con sistema de posicionamiento global, GPS, denavegación). Los datos se analizaron con estadísticas descriptiva, semivariograma ajustado a modelos ykrigeado ordinario [6]. También en Brasil, se investigó las propiedades del suelo de un pasto cultivado y con pastoreo para evaluarla variabilidad espacial. Las propiedades muestreadas fueron: pH, calcio, magnesio y aluminiointercambiables, fósforo, acidez potencial, capacidad de intercambio catiónico, suma de bases, saturación debases, saturación de aluminio y carbono orgánico. el muestreo se hizo en un área de 56 m por 80 m encuadrículas de 8 m con espaciado regular. Los puntos de muestreo fueron las intersecciones de lascuadrículas, 88 en total. En cada punto de muestreo se tomaron muestras de 0-5 cm, 5-10 cm y 10-20 cm deprofundidad. Los datos se analizaron con la estadística descriptiva, semivariograma ajustado al modeloesférico y krigeado [4]. En Brasil se muestreó un pastizal de pastoreo para medir la variabilidad espacial de flujo de dióxido decarbono, temperatura y humedad del suelo. El lugar de muestreo fue de 70 m por 70 m, en cuadrículasregulares de 10 m de espaciamiento. Los puntos de muestreo fueron las intersecciones de las cuadrículas yla profundidad de muestreo fue de 0-20 cm. Los datos fueron analizados con estadística descriptiva,semivariograma ajustados a modelos y krigeado [1]. Por otro lado, en Venezuela se muestrearon las propiedades del suelo de un cultivo agrícola: arena, arcilla,limo, pH y materia orgánica. Una finca de 286 ha se dividió en cinco unidades de suelo según la variabilidadanalizada en una fotografía aérea (escala 1:25 000). Se muestreó en el primer horizonte genésico (no seindica la profundidad) por puntos de muestreo separados sistemáticamente entre de 50-100 m y seobtuvo un total de 67 muestras que cubría toda la finca. Los datos se analizaron con la estadísticadescriptiva, semivariograma ajustado a modelos teóricos y kigreado puntual [13]. En otro lugar, Arabia Saudita, se muestreó un suelo con pasto cultivado para medir la variabilidad espacialde sus propiedades: pH, conductividad eléctrica, textura y compactación. Se muestreó 16 ha de suelo pormedio de cuadrículas regulares de 40 m por 40 m de separación. Los puntos de muestreo fueron 86, laprofundidad de muestreo fue de 0-20 cm, excepto para la compactación que fue de 0-15 m. Los datosfueron analizados con semivariograma y krigeado ordinario [3].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal145 Finalmente, en España se diseñó un muestreo para medir la variabilidad de la conductividad eléctrica,resistencia a la penetración y humedad gravimétrica del suelo de un cultivo de maíz, que antes tenía pastode corte. El diseño consistió en 40 puntos de muestreo cada 10 cm hasta los 90 cm de profundidad,distribuidos en 6 ha, con muestreo no aleatorio dentro de los transectos donde se midió la conductividadeléctrica. La mínima distancia entre dos puntos adyacentes fue de 12 m y máxima de 113 m. La mínimadistancia entre transectos fue de 4 m y la máxima de 30 m. Se usó la conductividad eléctrica del suelo paraelegir los puntos de muestreo, por estudios previos que indican una relación de la conductividad eléctricacon la textura, humedad y salinidad del suelo. Los datos se analizaron con la estadística descriptiva,modelado del variograma y krigeado ordinario y universal [12]. B. Diseños de muestreo de propiedades del suelo en Ecuador En una investigación realizada en El Oro se muestreó las propiedades del suelo de un pastizal cultivadopara pastoreo y bosque: textura, pH, fósforo, potasio, calcio, hierro, magnesio, carbono orgánico, flora total ybacterias. Se muestreó 3 ha de 2 parcelas por medio de 2 transectos de 30 m cada uno, separados entre cada 50 m de distancia. Los puntos de muestreo fueron a los 0 m, 15 m y 30 m dentro de cada transecto pormedio de calicatas escalonadas. Una de las parcelas tenía forma de L y se muestreó a los 0 m y 25 m endirección vertical y de 25-50 m en dirección horizontal. La profundidad de muestreo fue de 0-10 cm, 10-20cm y 20-30 cm. También se muestrearon 2 ha de bosque nativo en 4 transectos ubicados a lo largo de lospastizales. El análisis estadístico se hizo con análisis de varianza y correlación. No se reporta análisis de lavariabilidad espacial del suelo [9]. En otra investigación ejecutada en Quevedo se muestreó propiedades del suelo en cinco sistemas de usodel suelo: bosque primario, pasto cultivado para pastoreo, cacao, palma aceitera y maíz. Las propiedadesevaluadas fueron: conductividad hidráulica, densidad aparente, densidad real, porosidad total, porosidad deaireación, textura, arcilla dispersa en agua, grado de floculación, materia orgánica e índice de materiaorgánica. En cada sistema se tomaron muestras hasta 60 cm de profundidad cada 10 cm, por triplicado y seobtuvo 18 muestras en cada sistema. No se reporta las distancias de muestreo ni el área de muestreo. Seutilizó el análisis de varianza, correlación y análisis de componentes principales para el analizar los datos. Nose indica análisis de variabilidad espacial del suelo [14]. Por otro lado, en Machala se investigó las propiedades del suelo de cinco agrosistemas: banano (8,32 ha),ciclo corto (1,74 ha), pastos cultivados para pastoreo (4,89 ha), cacao (5,98 ha) y bosque secundario (2,43 ha).El área total de los agrosistemas es de 23,36 ha. Las propiedades medidas fueron: densidad real, textura, pH,materia orgánica, conductividad eléctrica y capacidad de intercambio catiónico. Dentro de cadaagroecosistema se fijaron aleatoriamente tres puntos de muestreo y las muestras se tomaron cada 15 cmhasta los 30 cm de profundidad por medio de calicatas. En total se obtuvieron 6 muestras poragroecosistema. No se indica las distancias entre puntos de muestreo. El análisis estadístico empleado fue elanálisis de varianza. No se reportan análisis de la variabilidad espacial del suelo [15]. Finalmente, en Pastaza se muestrearon propiedades de suelos con varios usos que fueron 3 chacrasdistintas (terrenos con varios cultivos), 3 pastos distintos y 1 boque natural: densidad aparente,conductividad hidráulica saturada, porosidad total y de aireación, pH, calcio, magnesio, potasio, fósforo, zinc,aluminio, cobre, acidez intercambiable, carbono orgánico total, número de lombrices. En cada suelo sedefinió un transecto que cubrió el terreno y se muestreó sistemáticamente en 5 parcelas de 10 m por 10 m,donde se tomó una muestra compuesta de 5 submuestras, a profundidades de 0-10 cm y de 10-30 cm. Lasmuestras para el análisis físico y químico del suelo se tomaron separadas, que fueron 20 en total porparcela. No se indica en área de muestreo ni la distancia entre parcelas de muestreo. El análisis utilizado fueel análisis de varianza, análisis de componentes principales y correlación. No se indica análisis de variabilidadespacial del suelo [16].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal146 C. Variabilidad espacial de las propiedades del suelo a nivel mundial Una de las investigaciones de muestreo de la variabilidad de propiedades del suelo cultivado con pastos depastoreo que fue hecha en Brasil, se observa que el semivariograma indica independencia de propiedadesde fertilidad del suelo entre 62 m para el Mg a 10 000 m para el pH. Los modelos que se ajustan alsemivariograma fueron el esférico, gausiano y exponencial. Para el suelo en particular muestreado serecomienda distancias de 62 m para el muestreo de la variabilidad espacial de las propiedades de fertilidaddel suelo [7]. En otra investigación, asimismo en Brasil, la independencia espacial de propiedades físicas de un suelo conpastizal cultivado con pastoreo está entre 35 m para la resistencia a la penetración y 110 m para lainfiltración. En cuanto a las propiedades químicas de fertilidad del suelo se observa la distancia deindependencia espacial entre 50 m para la capacidad de intercambio catiónico y 140 m de distancia para elpotasio. Se ajustan los modelos: exponencial, gausiano y, en su mayoría, el esférico. Con respecto a larelación entre variables, la compactación del suelo no tiene relación con la textura del suelo, si tiene con elpisoteo del ganado. En general hay dependencia entre los mapas de pH, capacidad de intercambio catiónicoy saturación de bases con los de materia orgánica, calcio y magnesio, las propiedades físicas del suelopermiten detectar áreas degradadas del pastizal y la producción del pasto está más influida por laspropiedades físicas que por las químicas [6]. Por otro lado, en España, en un suelo cultivado con maíz que anteriormente tenía pasto de corte, seobtiene un semivariograma de la resistencia a la penetración hasta 50 cm de profundidad ajustado almodelo semiesférico; sin embargo, de 50-90 cm solo hay efecto pepita puro debido a que se tomaronfrecuentemente datos poco útiles. Además, la conductividad eléctrica aparente (medida directamente encampo) representa eficientemente la variación espacial de la resistencia a la penetración, porque ambasvariables están relacionadas con la humedad, además de mejorar ligeramente el krigeado al ser usada comocokrigeado universal [12]. Finalmente, en Venezuela se muestreó las propiedades de un suelo de cultivo agrícola y se reporta que laindependencia espacial está entre 347,3 m para el limo y arena y 469,7 m para el pH. La dependenciaespacial es fuerte para la materia orgánica, pH y limo; moderada para la arena; débil para la arcilla. El modeloque se ajusta a la semivarianza fue el esférico. El krigeado muestra que el terreno se puede dividir enunidades de suelo con alta coincidencia a la división con base en la fotografía aérea, que inicialmente se usópara dividir el terreno en dichas unidades y muestrear [13]. D. Variabilidad de las propiedades de suelos en Ecuador En una investigación realizada en los Ríos se muestreó propiedades físicas del suelo de sistemas de uso delsuelo. El suelo con pasto cultivado para pastoreo presentó cambios de las propiedades principalmente enlas capas superficiales del suelo, al igual que los monocultivos, al ser comparados con el suelo de bosqueprimario. Los cambios de las propiedades perjudicaron el movimiento vertical del agua en el suelo. Según elanálisis de componentes principales el suelo del pastizal tiene los mayores contenidos de humedadvolumétrica, limo, conductividad hidráulica y relación materia orgánica/(limo más arcilla). Asimismo, en todoslos sistemas de uso de suelo la densidad aparente y materia orgánica, la densidad real y arcilla dispersa enagua, la conductividad hidráulica y materia orgánica/(limo más arcilla) tienen las mayores correlacionespositivas; mientras que la humedad volumétrica con la densidad real y arcilla dispersa en agua, la arena conla densidad aparente y materia orgánica tuvieron las mayores correlaciones negativas [14]. Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal147 Por otro lado, en El Oro se muestreó en suelos destinados a varios agrosistemas y se encuentra que ladensidad real, materia orgánica, del suelo con pastos cultivados para pastoreo es significativamente menorque el suelo de bosque, banano, cacao y maíz; mientras que la capacidad de intercambio catiónico yconductividad eléctrica del suelo de pastos son significativamente mayores que los demás cultivosnombrados [15]. Por último, en Pastaza se observó que el suelo destinado al cultivo de pastos tuvo diferencias significativasde propiedades físicas, químicas y biológicas al ser comparadas con el suelo de bosque y chacras. Tambiénse evaluó los suelos con un índice de calidad y reportan que los destinados a los pastos y chacras en generaltienen una calidad levemente superior; sin embargo, las chacras y bosque tienen mayor calidad en la capasuperficial del suelo. Los suelos de 10-30 cm de profundidad de los pastos tienen los mayores valores dedensidad aparente, porosidad de retención y pH; mientras que a profundidades de 0-10 cm el aluminio,potasio, magnesio y calcio tienen los mayores valores. Por otro lado, la conductividad hidráulica saturada y laporosidad de aireación, el potasio, calcio y magnesio, la densidad aparente y porosidad de retención tienenlas mayores correlaciones positivas; mientras que la densidad aparente con la conductividad hidráulicasaturada y porosidad de aireación, y la porosidad de retención con la materia orgánica tienen las mayorescorrelaciones negativas [16].CONCLUSIONES Se ha realizado una revisión preliminar de literatura acerca del muestreo de la variabilidad espacial depropiedades del suelo de pastizales. A nivel mundial existen diseños de muestreo de suelos con pastizalesque evalúan la variabilidad espacial del suelo; mientras que a nivel nacional (Ecuador) los muestreosrevisados no estudian la variabilidad espacial del suelo (relación de las propiedades con sus coordenadasdonde fue medida), generalmente se muestrea dentro de parcelas de un terreno para estimar parámetrosdescriptivos de las propiedades del suelo, comparaciones entre diferentes usos del suelo, correlaciones orelaciones entre propiedades. Los diseños del muestreo de la variabilidad espacial del suelo empleados son en cuadrículas regulares,transectos o sin un patrón regular de muestreo. Los muestreos en cuadrículas regulares tienen de 8 m hasta40 m de lado de cuadrícula y los puntos de muestreo son generalmente en las intersecciones de lacuadrícula. Mientras que los muestreos en transectos o sin un patrón regular tienen distancias de 4 m a 50m entre sí, con puntos de muestreo espaciados de 12 m a 113 m. El número de puntos de muestreo va de40 a 88 distribuidos en áreas de 0,45 ha hasta 700 ha. La profundidad de muestreo va de 0-90 cm, en capasde al menos 5 cm de espesor. Las técnicas usadas para estudiar la variabilidad espacial del suelo son el semivariograma ajustado a losmodelos esférico, gausiano y exponencial, y el krigeado ordinario y universal. En general, las propiedades delsuelo empiezan a tener independencia espacial a partir de los 35 m de distancia. Finalmente, es importante llevar a cabo en Ecuador investigaciones de la variabilidad espacial depropiedades del suelo de pastizales, que permitan tomar decisiones fiables para su manejo y conservación.REFERENCIAS[1] D. Silva, M. Costa, B. Mantovanelli, L. Coutrim, M. Rodrigues, J. Cunha, “Variabilidade espacial da emissãode CO2, temperatura e umidade do solo em área de pastagem na região Amazônica, Brasil”, Revista deCiências Agroveterinárias, vol. 18, no. 1, pp. 119-126, noviembre 2019, [En línea],http://revistas.udesc.br/index.php/agroveterinaria/article/view/10155/pdf.[2] M. Córdoba, P. Paccioretti, F. Giannini, C. Bruno, M. Balzarini. “Guía para el análisis de datos espaciales enagricultura”, ed., Serie Estadística Aplicada, Córdoba: Brujas, 2019, [En línea],https://ri.conicet.gov.ar/bitstream/handle/11336/128391/CONICET_Digital_Nro.ª5eab968-1409-4511-b6d6-07134144803b_A.pdf?sequence=2&isAllowed=y.Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal148[3] E. Tola, K. Gaadi, C. Biradar, R. Madugundu, A. Zeyada, A. Kayad, “Characterization of spatial variability ofsoil physicochemical properties and its impact on Rhodes grass productivity”, Saudi Journal of BiologicalSciences, vol. 24, pp. 421-429, abril 2017, [En línea], https://doi.org/10.1016/j.sjbs.2016.04.013[4] A. Lima, M. Campos, J. Cunha, L. Silva, F. Oliveira, B. Mantovanelli, E. de Brito, R. Gomes, 2019, “Spatialvariability and sampling density of chemical attributes in archaeological black earths under pasture insouthern Amazonas, Brazil, SOIL Discuss, junio 2019, [Preimpreso], [En línea], https://doi.org/10.5194/soil-2019-26[5] R. Petrucci, F. Herring, J. Madura, C. Bissonnette, Química general, Madrid: Pearson, 2011.[6] C. Grego, C. Gonçalves, S. Nogueira, F. Andrade, A. de Oliveira, C. Ferrer, A. dos Santos, J. de Abreu,“Variabilidade espacial do solo e da biomassa epígea de pastagem, identificada por meio de geostatística”,Pesquisa agropecuaria brasileira, vol. 47, np. 9, pp. 1404-1412, septiembre 2012, [En línea],https://doi.org/10.1590/S0100-204X2012000900026[7] A. Bernardi, K. Santos, G. Bettiol, L. Rabello, R. Ferreira, R. Inamasu, “Spatial variability of soil propertiesand yield of a grazed alfalfa pasture in Brazil”, Precision Agriculture, vol. 17, pp. 737-752, abril 2016, [En línea],https://doi.org/10.1007/s11119-016-9446-9.[8] A. Delgado, J. Gómez, “The Soil. Physical, Chemical and Biological Properties”, en Principles of Agronomyfor Sustainable, Springer, 2016. [En línea], https://link.springer.com/chapter/10.1007/978-3-319-46116-8_2[9] S. Barrezueta, K. Velepucha, L. Hurtado, E. Jaramillo, “Soil properties and storage of organic carbon in theland use pasture and forest”, Revista de Ciencias Agrícolas, vol. 36, no. 2, pp. 31-45, abril 2019, [En línea],https://doi.org/10.22267/rcia.193602.116[10] J. Wu, J. He, G. Christakos, Quantitative analysis and modeling of earth and environmental data. Space-time and spacetime data considerations, India: Elsevier, 2022, [En línea], https://doi.org/10.1016/C2018-0-00055-9[11] M. Oliver, R. Webster, “Basic Steps in Geostatistics: The Variogram and Kriging”, Springer, 2015, [En línea],https://link.springer.com/book/10.1007/978-3-319-15865-5?error=cookies_not_supported&code=76d5fdcb-b77b-4658-a544-6bf0d847312f[12] G. Machado, J. Dafonte, J. Bueno, M. Valcacer, E. Farias. “Using Soil Apparent Electrical Conductivity toOptimize Sampling of Soil Penetration Resistance and to Improve the Estimations of Spatial Patterns of SoilCompaction”, The Scientific World Journal, pp. 1-12, diciembre 2014, [En línea],http://dx.doi.org/10.1155/2014/269480[13] A. González, J. Rey, J. Atencio,” Variabilidad espacial de los suelos de la unidad experimental la glorietacon fines agropecuarios”, Revista de la Facultad de Agronomía (LUZ), suplemento 1, pp. 539-553, junio 2014,[En linea], https://www.revfacagronluz.org.ve/PDF/suplemento_2014/ing/ingsupl12014539553.pdf[14] Novillo I.; Carrillo, M.; Cargua, J.; Moreira, V.; Albán, K.; Morales, K, “Propiedades físicas del suelo endiferentes sistemas agrícolas en la provincia de Los Ríos, Ecuador”. Temas Agrarios, vol. 23, no. 2, pp. 177 -187, mayo 2018,https://www.researchgate.net/publication/329046726_Propiedades_fisicas_del_suelo_en_diferentes_sistemas_agricolas_en_la_provincia_de_Los_Rios_Ecuador[15] I. Rodríguez, H. Pérez, R. García, A. Quezada. “Efecto del manejo agrícola en propiedades físicas yquímicas del suelo en diferentes agroecosistemas”, Universidad y Sociedad, vol. 12, no. 5, pp. 389-398,septiembre 2020, [En línea], http://scielo.sld.cu/pdf/rus/v12n5/2218-3620-rus-12-05-389.pdf[16] C. Bravo, F. Goyes, Y. Arteaga, Y. Garcia, D. Changoluisa. “A soil quality index for seven productivelandscapes in the Andean-Amazonian foothills of Ecuador”, Land Degradation & Development, vol. 32, pp.2226-2241, enero 2021, [En línea], https://doi.org/10.1002/ldr.3897[17] J. Serrano, S. Shahidian, J. Silva, “Spatial variability and temporal stability of apparent soilelectrical conductivity in a Mediterranean pasture”, Precision Agriculture, vol. 18, pp. 245-263, julio 2017, [Enlínea], https://doi.org/10.1007/s11119-016-9460-y[18] D. Brus, J. de Gruijter, “Random sampling or geostatistical modelling? Choosing between design-basedand model-based sampling strategies for soil (with discussion)”, Geoderma, vol. 80, pp. 1-44, octubre 1997,[En línea], https://doi.org/10.1016/S0016-7061(97)00072-4.Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)
ISSN-E: 2697-3650Minerva Journal149[19] P. Lawrence, W. Roper, T. Morris, K. Guillard, “Guiding soil sampling strategies using classical and spatialstatistics: A review”, Agronomy Journal, vol. 112, pp. 493-510, enero 2020. [En línea],https://doi.org/10.1002/agj2.20048.[20] J. Wang, A. Stein, B. Gao, Y. Ge, “A review of spatial sampling”, Spatial Statistics, vol. 2, pp. 1-14, agosto2012, [En línea], https://doi.org/10.1016/j.spasta.2012.08.001.Rodrigo Saquicela, ingeniero agropecuario, especializado en nutrición vegetal.Maestrante del Instituto Posgrado de la Universidad Técnica de Manabí, UTM,Portoviejo, Ecuador. Docente de la Universidad UTE.Lelly María Useche Castro, ingeniero industrial, PhD. es estadística. Directora dedel Grupo de Análisis Multivariante y Estocástico (G.A.M.E). Instituto de CienciasBásicas, Departamento de Matemáticas y Estadística. Universidad Técnica deManabí, Ecuador.Ana Francisca González Pedraza, Ingeniera Agrónoma, Magister en manejo de losrecursos agua y suelo. Doctora en Ciencias, mención Ecología. Docente de laUniversidad de Pamplona, Pamplona. Norte de Santander, Colombia.Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3, Spec. Issue N°1, (pp. 138-149)