I.INTRODUCCIÓNElsueloesunconjuntodinámicoformadoporpartesfísicas,químicasybiológicasquebásicamentemantienenlaproductividadbiológicayconservanlacalidaddelambiente[8].Dichaspartessonpropiedadesquevaríanespacialmenteenelsuelodebidoalmanejo,climaybiota[19].Elpastoreoenunmanejoqueinfluyeensuspropiedades,comolacompactaciónqueestárelacionadaconelpisoteodelganado[6].Así,anivelmundialhayinvestigacionesdelavariabilidadespacialdepropiedadesdelsuelodeagriculturaypastizales[7],[12],[13],[17].EnEcuadorsedesconocelasituaciónactualdelmuestreodelaspropiedadesdelossuelosconpastizalesdepastoreo,confinesdeplanificardiseñosdemuestreosqueevalúenlavariabilidadespacialyprediganconprecisiónelvaloresperadodeestaspropiedades.Elobjetivodeestarevisióndeliteraturafueconocerelestadoactualdeldiseñoyanálisisdemuestreodelavariabilidadespacialdepropiedadesdesuelosdepastizalesyelcomportamientodelavariacióndeestas.Estarevisiónsebasaenelhechodequelaspropiedadesdelsuelovaríanprincipalmenteporelmanejoagronómico,clima,biotaytipodesuelo[19].Lamedicióndepropiedadesdelsueloconprediccionesprecisaspermitetomardecisionesfiablesparasumanejoyconservación.Deestaforma,existenherramientasparamediryanalizarlavariabilidadespacialdepropiedadesdelsuelo,comoelmuestreobasadoenelmodelo,semivariogramaykrigeado[2].SeanalizarondiversaspublicacionescientíficasrelacionadasconlamedicióndelavariabilidadespacialdelaspropiedadesdelsuelodepastizalesanivelmundialyenEcuador.Tambiénseexaminaronlibrosacercadelmuestreoyanálisisdelavariabilidadespacialdelsuelo.Seresumiólainformaciónencontrada.Elpresenteartículotienelasiguienteestructura:resumen,traduccióndelresumen,introducción,desarrollodelarevisióndelaliteratura,metodología,resultados,conclusionesyreferencias.ISSN-E:2697-3650MinervaJournal139II.DESARROLLOA.PropiedadesdelsueloElsueloesunaunidaddinámicaconcomponentesfísicos,químicosybiológicosqueinteractúandeformacomplejaydefinensufuncionamientoparavariadosfines.Elfuncionamientoserefierealacalidaddelsuelo,queessucapacidadparaimpulsarlasaluddeanimalesyplantas,conservarlacalidadambientalymantenerlaproductividadbiológica[8].Loscomponentesdelsuelosonmateriaconcualidadesyatributosdenominadospropiedades,quesirvenparadistinguirentremuestrasdesuelo[5],[8].Elsuelotienepropiedadesfísicas,químicasybiológicas.Lasfísicasnocambianlacomposicióndelamuestramientrassemideuobservalapropiedad[8].Lasquímicassonloscambiosenlacomposicióndeunamuestraencondicionesdereacciónquímica[5].Finalmente,lasbiológicassoncambioscausadosporlosorganismosenlaspropiedadesfísicasyquímicas[8].Entrelaspropiedadesfísicasdelsueloestánlaconductividadhidráulica,retencióndeagua,porosidadypropiedadesmecánicas.Entrelasquímicasseencuentranelpotencialdehidrógeno(pH),estadoredox,retencióndeionesenelsuelo,salinidadysodicidad.Porúltimo,lasbiológicasseconsideranlamineralizaciónpotencialdenitrógenoocarbono,tasaderespiración,lombricesdetierra,biomasabacteriana,diversidadbacteriana,presenciadepatógenos[8].Porotrolado,laspropiedadesdelsuelotienenvariaciónycorrelaciónespacialdebidoalmanejoagronómico,biota,tipodesueloyclima[19].Lasobservacionesmáscercanasentresísonfrecuentementemásparecidasquelasmáslejanas[20].Sehareportadoinvestigacionesdevariabilidadespacialdepropiedadesfísicas,químicasybiológicasdelsuelo[1],[6].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) 140ISSN-E:2697-3650MinervaJournalMedirconprecisiónlaspropiedadesdelsueloayudaasudiagnósticoytomadedecisionesparasumanejoyconservación.Paratalefecto,sehandesarrolladoherramientasparaeldiseñodelmuestreodelavariaciónespacialdelaspropiedadesdelsuelo,talescomoelmuestreobasadoeneldiseño,quenoconsideralavariaciónespacial;yelbasadoenelmodelo,queconsideralavariaciónespacial[20].B.MuestreodelsuelobasadoeneldiseñoExistendosenfoquesparaelmuestreodepropiedadesdelsuelo:losbasadoseneldiseñoylosbasadosenelmodelo.Elmuestreobasadoeneldiseñosirveparamuestrearconobjetividadysinsesgos.Lainferenciasebasacompletamenteeneldiseñodelmuestreoylaindependenciaobtenidaconlaselecciónaleatoriadelasmuestras.Cadapuntodemuestreoenelespacioytiemposeconsiderafijo,asícomoelvalordelavariablequeseobtengaendichopunto[19],[18].Conrespectoalainferencia,seinfiereparámetrospoblacionales(comomediayvarianza)únicos(¿cuántoessuvalor?)paratodoeldominioespacial,siemprequeeltamañodelamuestraseaadecuado,conlossupuestosdeindependenciaespacial(quenoexistacorrelaciónespacialdelasvariables)yquelapoblación(dedatos)estédistribuidaidénticamenteeneldominioespacial.Losestimadoresdelosparámetrossoninsesgados[18],[19],[20].Deestaforma,elmuestreobasadoeneldiseño,parapoblacionesidénticamentedistribuidaseindependientes,utilizalosdiseñosdemuestreosaleatoriosclásicos:simple,sistemático,estratificado,conglomerados,ademásdelmuestreoaleatorioendospasos.Elúltimoconsisteenlaelecciónaleatoriadeestratosyelecciónaleatoriademuestrasdentrodelosestratoselegidos.Enlosmuestreossimple,sistemáticoyconglomeradossesuponelaindependenciaypoblaciónidénticadentrodeldominioespacial;mientrasqueenelestratificadoydedospasoslasuposiciónesdentrodelosestratos.Noobstante,lavariaciónycorrelaciónespacialdelasvariablesestápresenteeneldominioespacial[20],loquenopermitiríamuestreosfiablesbasadoseneldiseño.C.MuestreodelsuelobasadoenelmodeloElmuestreobasadoenelmodelopermitemuestrearsinelsupuestodeindependenciaespacial,portanto,noesnecesarialaaleatorizaciónparaelegircadapuntodemuestreo[19].Losvaloresdelavariableencadapuntodemuestreoseconsideranaleatoriosyrepresentaunmuestreoúnico,lapoblación,deinfinitosmuestreosposiblesaununiversorepresentadoporunmodeloestocástico(convariablesaleatoriasquecambianentreinstantesdemuestreo),lasuperpoblación.Deesamanera,laaleatoriedaddeldiseñosegeneraporunconjuntodemodelosestocásticos[18].Además,seasumelaindependenciaespacialalmuestrearvariablesaleatoriasdeununiversoestocástico.Unmodeloestocásticousadoeselkrigeado[19],[20].Enloreferentealainferencia,lasvariablessesuponenaleatoriasenlospuntosdemuestreodeundiseñobasadoenelmodelo,así,laverdaderamediaespacialtambiénesaleatoria.Portanto,sepredice(¿dóndeestánlosvalores?)dichamediaeneldominioespacial,noseestimaconunúnicovalor[18].Deestamanerasepuedepredecirvaloresdelavariableenpuntosnomuestreados,elaborarmapasdelavariacióndelaspropiedadesdelsueloyestimarparámetrosdelmodeloestocástico[20].Porotrolado,existenlossiguientesobjetivosacumplircuandosemuestreaconbaseeneldiseño:minimizarlaestimacióndelavarianzadelerror,cubrirporigualeldominioespacialycubrirporigualenelespaciodelascaracterísticas[20].Paraminimizarlaestimacióndelavarianzadelerrorseusaelkrigeado,queestimaimparcialmentelosvaloresdeunavariableenpuntosnomuestreadosconunmínimoerrordeestimacióncuadráticomedio.Tambiénexistenlasmediasdesuperficieconnohomogeneidad[20].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) Porotrolado,paracubrirporigualeldominioespacialseasignanlospuntosdemuestreodetalformaquecubraneldominiolomásuniformementeposiblepormediodeunacuadrículafinita,segúnelcriteriodelamínimadistanciapromedioentrelospuntosadyacentesdelacuadrículaaplicadaasucentro[20].Deestaforma,lacuadrículatriangularequiláteracumpleestacondiciónparaundominiosinrestriccionesenasignarpuntosdemuestreo[11].Tambiénexistenotrasopcionesparaladistanciamínima,comolamediaponderadadelasdistanciasmáscortasyladistanciamediacuadráticabasadaenpolígonosdeThiessen.Otrasformasdecuadrículausadasenlaprácticason:muestreoencuadrícula,muestreoentransecto,muestreosecuencialymuestreoanidado[20].Paracubrirporigualelespaciodelascaracterísticasdeldominioespacialsedistribuyelasunidadesdemuestreodeunaformatalquepermitaobtenerlamejordistribuciónposibledelapoblación.Paraelefecto,seaplicandoscriteriosprincipales:elmuestreodehipercubolatinoyeldeWarrick-Myers.Encasodedesconocerdistribucionesdelapoblaciónpreviasalmuestreoseutilizalaexperienciaencasossimilaresodatosauxiliares[20].Enloreferentealtamañodelacuadrícula,paraelmuestreodelavariabilidaddepropiedadesquímicasdelsuelo,varíaentre1ma100mporlado.Lospuntosdemuestreoseubicanprincipalmenteenelcentroideoenlasesquinasdelascuadrículas.Tambiénpuedenestardistribuidosportodalacelda[19].Porotrolado,lacantidaddepuntosdemuestreoparaestimarsemivarianzas(baseparaestimarelmodeloestocástico)dentrodelímitesdeconfianzaaceptablessesugierealmenosentre100a150puntosdemuestreo[11].Asímismo,sehandesarrolladotécnicasestadísticasparaelanálisisespacialdedatos,comoelsemivariograma,paraanalizarlaestructuradelavariaciónespacialdepropiedadesdelsuelo;yelkrigeado,quepredicelasvariablesespacialesenfuncióndelsemivariograma[2].D.AnálisisestructuralElanálisisestructuralconsisteenencontrarunafunciónqueexpliquelaestructuradeladependenciaocorrelaciónespacialdelasvariablesmedidas.Lafuncióneslabaseparalapredicciónespacialdelasvariablesysedenominanfuncionesdecovarianza(covariogramas)ysemivariogramas.Lafuncióndecovarianzaempíricaseindicaen(1),dondeeslacovarianzaempíricaparalaclasededistanciadeintervalo;eselnúmerodeparesdedistanciasseparadasporladistanciah;yysonlosvaloresdelamuestramedidosenlaposicióniei+h,respectivamente;eslaverdaderamediadelapoblaciónyseestimaconelpromediode;N=eselnúmerototaldedatos[10].141ISSN-E:2697-3650MinervaJournalLafuncióndelsemivariogramaempíricousualmentesecalculaconelestimadordemomentosdeMatheron(2),dondeeselsemivariogramaempíricoparalaclasededistanciadeintervalo[2].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) 142ISSN-E:2697-3650MinervaJournalNoobstante,cuandolavarianzaresidualespequeña,elsemivariogramasepuedecalculardeformamenossesgadaconelestimadordeCressie-Hawkins[2].Elsemivariogramasecompletaalmodelarelsemivariogramaenfuncióndeladistancia,conlosmodelosexponencial(3),esférico(4),gausiano(5);dondeeslaordenadadeorigendelsemivariograma(pepitaonugget)querepresentalavarianzatotaldeloserroresaleatoriososinestructuraespacialyloserrorespornomedirconmenordistanciadeintervalo;eslavarianzaestructuralomesetadelsemivariograma(sill);eselrangoeindicaladistanciadesdeelorigenhastacuandosellegaalameseta.Elrangoprácticoesladistanciadesdeelorigenhastacuandolasemivarianzael95%delavarianzaltotal;siemprequeelsemivariogramallegueasintóticamentealameseta;indicalavarianzaumbralyapartirdeahílasobservacionessonindependientes[10].Además,hayotrosmodeloscomoelexponencial-coseno,deonda,DeWijsian,cúbico,pentaesférico,pepitapuro[10].Losmétodosdeajustedelosdatosalmodelosonpormínimoscuadradosymodeloslinealesmixtos[2].Noobstante,cuandolavarianzaresidualespequeña,elsemivariogramasepuedecalculardeformamenossesgadaconelestimadordeCressie-Hawkins[2].Elsemivariogramasecompletaalmodelarelsemivariogramaenfuncióndeladistancia,conlosmodelosexponencial(3),esférico(4),gausiano(5);dondeeslaordenadadeorigendelsemivariograma(pepitaonugget)querepresentalavarianzatotaldeloserroresaleatoriososinestructuraespacialyloserrorespornomedirconmenordistanciadeintervalo;eslavarianzaestructuralomesetadelsemivariograma(sill);eselrangoeindicaladistanciadesdeelorigenhastacuandosellegaalameseta.Elrangoprácticoesladistanciadesdeelorigenhastacuandolasemivarianzael95%delavarianzaltotal;siemprequeelsemivariogramallegueasintóticamentealameseta;indicalavarianzaumbralyapartirdeahílasobservacionessonindependientes[10].Además,hayotrosmodeloscomoelexponencial-coseno,deonda,DeWijsian,cúbico,pentaesférico,pepitapuro[10].Losmétodosdeajustedelosdatosalmodelosonpormínimoscuadradosymodeloslinealesmixtos[2].Lavarianzaestructuralrelativaeslarazónentrelavarianzaestructuralylavarianzaumbraleindicaelgradodeestructuraciónespacial.Mientrasmásaltoseaelvalordelavarianzaestructuralrelativa,máseficientesseránlasprediccionesestadísticas.Seclasificaenbajo,<25%;medio,25%a75%;alto,>75%[2].E.PredicciónespacialLaprediccióndelavariableeneldominioespacial(campo,deterreno,dondesemuestrea)sedenominapredicciónespacial.Elmétododeprediccióncomúneselkrigeado(kriging)quetienesubaseenelsemivariogramaempírico.Elkrigeadogeneraelvaloresperadoconelmejorestimadorlinealinsesgadoyelerrordeestimación(varianzakrigeada).Existeelkrigeadopuntualyenbloques.Elpuntualseusaparapredecirlavariableenlugares(puntos)nomuestreadosmuchomenoresquelasdistanciasespaciales;mientasqueelkrigeadoenbloquesesparapredecirelpromediodelavariableenunasuperficiemayorqueelpuntual,pormediodeestimacionespuntualesendichasuperficie[2].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournalEntreelkrigeadopuntalexistenelkrigeadoordinario(6),simple(7)yuniversal(8)paralapredicciónespacial,dondeeslavariablepredichaencualquiernuevopunto;sonlospesosasignadosacadaobservaciónconsumadeuno;esfuncióndelascoordenadasespaciales.Elordinarioysimpleesparavariablescondistribuciónnormalyeluniversalparavariablesconmediaquedependedelavarianza.Elkrigeadoordinarioseutilizaparaestimarlocalmentelamediadelavariable;elkrigeadosimple,cuandoseconocelamediapoblacionaldelavariable;eluniversal,paraestimarlainfluenciaespacialdelosdatosconlascoordenadasespaciales[2].143III.METODOLOGÍALainformacióndemuestreodevariabilidadespacialdelaspropiedadesdelsuelodepastizalessehizoconelmotordebúsquedaGoogleAcadémico.Laspalabrasclaveusadasfueronlassiguientes:“spatialvariabilitysoilphysical”,“kriging”,“gaussiansecuantialsimulation”,“soildegradation”,“soilpropertiesspatialvariationgrass”,“soilsampling”,“gridsoilsampling”y“variabilidadeespacialdosolopastagem”.Seeligieron20documentoscientíficos,deloscuales15eranartículosy5eranlibros.LosartículosquegenerólabúsquedaestuvieronenbasesdedatosindexadascomoScopus,SciELO,ScienceDirect,ScimagoJournal&CountryRank,Dialnet,SpingerLink.LoslibrosconsultadosfueronpublicadosporeditorascomoElsiever,Springer,Pearson,Brujas,Mundi-PrensayWiley.Losartículoselegidosserevisaronprincipalmentelasmetodologías,resultadosyconclusiones,yseescribieronlosrespectivosresúmenes.Enloslibrosserevisaronconceptosdepropiedadesdelsueloydemuestreodelavariabilidadespacialdedichaspropiedades.Tambiénexisteelkrigeadoindicadorparavariablesbinariasovariablescontinuastransformadasabinarias;elkrigeadodisyuntivo(yelkrigeadoindicador)parainterpolaciónnolineal,elkrigeadointrínsecoparavariablesheterogéneasconvariasmedias;krigingrecursivoconelfiltroKalmanykrigeadoconlamáximaentropíabayesianaquetienemenosrestriccioneseincorporarmásconocimientoprevio[10].Encuantoalkrigeadoenbloqueselmáscomúneselkrigeadoordinarioenbloques,dondeeslavariablepredichaeneldominiocentradoenelpuntodeinterés[2].IV.RESULTADOSSepresentanlosresultadosdelarevisióndeliteraturadediseñodemuestreoyanálisisdelavariabilidadespacialdepropiedadesdelsuelodepastizales.A.DiseñosdemuestreodelavariabilidadespacialdepropiedadesdelsueloanivelmundialEnPortugalsediseñóunmuestreoparamedirlavariabilidadespacialdelatextura,humedad,pH,materiaorgánica,nitrógeno,fósforoypotasiodelsuelodeunpastizalcultivadoyconpastoreodeovejasyganadobovino.Setomaron76muestrascompuestasenunterrenode6hadivididoencuadrículasde28mpor28m.Lasmuestrassetomaronde0-30cmdeprofundidaddelsuelo.Lamuestracompuesta,paraobtenerunamuestraporcuadrícula,consistióencincosubmuestras,cuatromuestreadasenlasesquinasdelacuadrículayunadesucentro.Losdatosseanalizaronconestadísticadescriptiva,correlaciónymapasdevariabilidadespacialconuninterpoladordedistanciaalcuadradoinverso[17].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal144EnotrainvestigaciónllevadaacaboenBrasilenunpastocultivadoyconpastoreodebovinossemuestreópropiedadesdelsueloparamedirlavariabilidadespacialyelaborarmapasdefertilidaddelsuelo.SemuestreóelpH,contenidosdecarbonoorgánicoyarcilla,fósforo,calcioymagnesiodisponibles,potasioaluminioehidrógenointercambiables,capacidaddeintercambiocatiónicoefectiva,saturacióndebasesyconductividadeléctrica.Elpastizaltenía5,3hadividoen270potrerosdondesetomaron73muestrasdesuelo,exceptoparalaconductividadeléctrica,quefueron4794.Seutilizóelmuestreoporzonasdemanejo,queconsisteendelimitarzonasdesimilaresrendimientosdelpastoypresenciademalezas,dentrodelascualessetomóunamuestracompuestadealmenos10submuestras.Nosehizoelmuestroporcuadrículas,noseindicalasdimensionesdelaszonasdemuestreo.Laprofundidaddemuestreofuede0-20cmdeprofundidad.Seusóelsemivariogramaajustadoamodelosempíricosyprogramasdesistemasdeinformacióngeográficaparaelaborarmapasdefertilidaddelsuelo.Noseindicalacantidadtotaldemuestrasparalaconductividadeléctrica[7].Asímismo,enBrasil,semuestreóunpastizalcultivado(degradado)parapastoreolaspropiedadesdelsueloparaevaluarsuvariabilidadespacial.Semuestreóelcontenidodearena,limo,arcilla,de0-20cmdeprofundidad,ladensidadaparentede0-10cmyde10-20cm,resistenciaalapenetracióncada5cmdeprofundidadde0-30cm,permeabilidadeinfiltracióntridimensionalde0-20cm.ElpH,potasio,calcio,magnesio,fósforo,aluminio,materiaorgánica,sumadebases,acidezpotencial,capacidaddeintercambiocatiónico,saturacióndebasesydealuminio,semuestrearonde0-10cmyde10-20cm.Elpastizaltenía3,67haysetomaron50muestrasdesuelocada30maproximadamente,conespaciadosirregulares(nocuadrículas)entrepuntosdemuestreo(definidosconsistemadeposicionamientoglobal,GPS,denavegación).Losdatosseanalizaronconestadísticasdescriptiva,semivariogramaajustadoamodelosykrigeadoordinario[6].TambiénenBrasil,seinvestigólaspropiedadesdelsuelodeunpastocultivadoyconpastoreoparaevaluarlavariabilidadespacial.Laspropiedadesmuestreadasfueron:pH,calcio,magnesioyaluminiointercambiables,fósforo,acidezpotencial,capacidaddeintercambiocatiónico,sumadebases,saturacióndebases,saturacióndealuminioycarbonoorgánico.elmuestreosehizoenunáreade56mpor80mencuadrículasde8mconespaciadoregular.Lospuntosdemuestreofueronlasinterseccionesdelascuadrículas,88entotal.Encadapuntodemuestreosetomaronmuestrasde0-5cm,5-10cmy10-20cmdeprofundidad.Losdatosseanalizaronconlaestadísticadescriptiva,semivariogramaajustadoalmodeloesféricoykrigeado[4].EnBrasilsemuestreóunpastizaldepastoreoparamedirlavariabilidadespacialdeflujodedióxidodecarbono,temperaturayhumedaddelsuelo.Ellugardemuestreofuede70mpor70m,encuadrículasregularesde10mdeespaciamiento.Lospuntosdemuestreofueronlasinterseccionesdelascuadrículasylaprofundidaddemuestreofuede0-20cm.Losdatosfueronanalizadosconestadísticadescriptiva,semivariogramaajustadosamodelosykrigeado[1].Porotrolado,enVenezuelasemuestrearonlaspropiedadesdelsuelodeuncultivoagrícola:arena,arcilla,limo,pHymateriaorgánica.Unafincade286hasedividióencincounidadesdesuelosegúnlavariabilidadanalizadaenunafotografíaaérea(escala1:25000).Semuestreóenelprimerhorizontegenésico(noseindicalaprofundidad)porpuntosdemuestreoseparadossistemáticamenteentresíde50-100myseobtuvountotalde67muestrasquecubríatodalafinca.Losdatosseanalizaronconlaestadísticadescriptiva,semivariogramaajustadoamodelosteóricosykigreadopuntual[13].Enotrolugar,ArabiaSaudita,semuestreóunsueloconpastocultivadoparamedirlavariabilidadespacialdesuspropiedades:pH,conductividadeléctrica,texturaycompactación.Semuestreó16hadesuelopormediodecuadrículasregularesde40mpor40mdeseparación.Lospuntosdemuestreofueron86,laprofundidaddemuestreofuede0-20cm,exceptoparalacompactaciónquefuede0-15m.Losdatosfueronanalizadosconsemivariogramaykrigeadoordinario[3].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal145Finalmente,enEspañasediseñóunmuestreoparamedirlavariabilidaddelaconductividadeléctrica,resistenciaalapenetraciónyhumedadgravimétricadelsuelodeuncultivodemaíz,queantesteníapastodecorte.Eldiseñoconsistióen40puntosdemuestreocada10cmhastalos90cmdeprofundidad,distribuidosen6ha,conmuestreonoaleatoriodentrodelostransectosdondesemidiólaconductividadeléctrica.Lamínimadistanciaentredospuntosadyacentesfuede12mymáximade113m.Lamínimadistanciaentretransectosfuede4mylamáximade30m.Seusólaconductividadeléctricadelsueloparaelegirlospuntosdemuestreo,porestudiospreviosqueindicanunarelacióndelaconductividadeléctricaconlatextura,humedadysalinidaddelsuelo.Losdatosseanalizaronconlaestadísticadescriptiva,modeladodelvariogramaykrigeadoordinarioyuniversal[12].B.DiseñosdemuestreodepropiedadesdelsueloenEcuadorEnunainvestigaciónrealizadaenElOrosemuestreólaspropiedadesdelsuelodeunpastizalcultivadoparapastoreoybosque:textura,pH,fósforo,potasio,calcio,hierro,magnesio,carbonoorgánico,floratotalybacterias.Semuestreó3hade2parcelaspormediode2transectosde30mcadauno,separadosentresícada50mdedistancia.Lospuntosdemuestreofueronalos0m,15my30mdentrodecadatransectopormediodecalicatasescalonadas.UnadelasparcelasteníaformadeLysemuestreóalos0my25mendirecciónverticalyde25-50mendirecciónhorizontal.Laprofundidaddemuestreofuede0-10cm,10-20cmy20-30cm.Tambiénsemuestrearon2hadebosquenativoen4transectosubicadosalolargodelospastizales.Elanálisisestadísticosehizoconanálisisdevarianzaycorrelación.Nosereportaanálisisdelavariabilidadespacialdelsuelo[9].EnotrainvestigaciónejecutadaenQuevedosemuestreópropiedadesdelsueloencincosistemasdeusodelsuelo:bosqueprimario,pastocultivadoparapastoreo,cacao,palmaaceiteraymaíz.Laspropiedadesevaluadasfueron:conductividadhidráulica,densidadaparente,densidadreal,porosidadtotal,porosidaddeaireación,textura,arcilladispersaenagua,gradodefloculación,materiaorgánicaeíndicedemateriaorgánica.Encadasistemasetomaronmuestrashasta60cmdeprofundidadcada10cm,portriplicadoyseobtuvo18muestrasencadasistema.Nosereportalasdistanciasdemuestreonieláreademuestreo.Seutilizóelanálisisdevarianza,correlaciónyanálisisdecomponentesprincipalesparaelanalizarlosdatos.Noseindicaanálisisdevariabilidadespacialdelsuelo[14].Porotrolado,enMachalaseinvestigólaspropiedadesdelsuelodecincoagrosistemas:banano(8,32ha),ciclocorto(1,74ha),pastoscultivadosparapastoreo(4,89ha),cacao(5,98ha)ybosquesecundario(2,43ha).Eláreatotaldelosagrosistemasesde23,36ha.Laspropiedadesmedidasfueron:densidadreal,textura,pH,materiaorgánica,conductividadeléctricaycapacidaddeintercambiocatiónico.Dentrodecadaagroecosistemasefijaronaleatoriamentetrespuntosdemuestreoylasmuestrassetomaroncada15cmhastalos30cmdeprofundidadpormediodecalicatas.Entotalseobtuvieron6muestrasporagroecosistema.Noseindicalasdistanciasentrepuntosdemuestreo.Elanálisisestadísticoempleadofueelanálisisdevarianza.Nosereportananálisisdelavariabilidadespacialdelsuelo[15].Finalmente,enPastazasemuestrearonpropiedadesdesuelosconvariosusosquefueron3chacrasdistintas(terrenosconvarioscultivos),3pastosdistintosy1boquenatural:densidadaparente,conductividadhidráulicasaturada,porosidadtotalydeaireación,pH,calcio,magnesio,potasio,fósforo,zinc,aluminio,cobre,acidezintercambiable,carbonoorgánicototal,númerodelombrices.Encadasuelosedefinióuntransectoquecubrióelterrenoysemuestreósistemáticamenteen5parcelasde10mpor10m,dondesetomóunamuestracompuestade5submuestras,aprofundidadesde0-10cmyde10-30cm.Lasmuestrasparaelanálisisfísicoyquímicodelsuelosetomaronseparadas,quefueron20entotalporparcela.Noseindicaenáreademuestreoniladistanciaentreparcelasdemuestreo.Elanálisisutilizadofueelanálisisdevarianza,análisisdecomponentesprincipalesycorrelación.Noseindicaanálisisdevariabilidadespacialdelsuelo[16].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal146C.VariabilidadespacialdelaspropiedadesdelsueloanivelmundialUnadelasinvestigacionesdemuestreodelavariabilidaddepropiedadesdelsuelocultivadoconpastosdepastoreoquefuehechaenBrasil,seobservaqueelsemivariogramaindicaindependenciadepropiedadesdefertilidaddelsueloentre62mparaelMga10000mparaelpH.Losmodelosqueseajustanalsemivariogramafueronelesférico,gausianoyexponencial.Paraelsueloenparticularmuestreadoserecomiendadistanciasde62mparaelmuestreodelavariabilidadespacialdelaspropiedadesdefertilidaddelsuelo[7].Enotrainvestigación,asimismoenBrasil,laindependenciaespacialdepropiedadesfísicasdeunsueloconpastizalcultivadoconpastoreoestáentre35mparalaresistenciaalapenetracióny110mparalainfiltración.Encuantoalaspropiedadesquímicasdefertilidaddelsueloseobservaladistanciadeindependenciaespacialentre50mparalacapacidaddeintercambiocatiónicoy140mdedistanciaparaelpotasio.Seajustanlosmodelos:exponencial,gausianoy,ensumayoría,elesférico.Conrespectoalarelaciónentrevariables,lacompactacióndelsuelonotienerelaciónconlatexturadelsuelo,sitieneconelpisoteodelganado.EngeneralhaydependenciaentrelosmapasdepH,capacidaddeintercambiocatiónicoysaturacióndebasesconlosdemateriaorgánica,calcioymagnesio,laspropiedadesfísicasdelsuelopermitendetectaráreasdegradadasdelpastizalylaproduccióndelpastoestámásinfluidaporlaspropiedadesfísicasqueporlasquímicas[6].Porotrolado,enEspaña,enunsuelocultivadoconmaízqueanteriormenteteníapastodecorte,seobtieneunsemivariogramadelaresistenciaalapenetraciónhasta50cmdeprofundidadajustadoalmodelosemiesférico;sinembargo,de50-90cmsolohayefectopepitapurodebidoaquesetomaronfrecuentementedatospocoútiles.Además,laconductividadeléctricaaparente(medidadirectamenteencampo)representaeficientementelavariaciónespacialdelaresistenciaalapenetración,porqueambasvariablesestánrelacionadasconlahumedad,ademásdemejorarligeramenteelkrigeadoalserusadacomocokrigeadouniversal[12].Finalmente,enVenezuelasemuestreólaspropiedadesdeunsuelodecultivoagrícolaysereportaquelaindependenciaespacialestáentre347,3mparaellimoyarenay469,7mparaelpH.Ladependenciaespacialesfuerteparalamateriaorgánica,pHylimo;moderadaparalaarena;débilparalaarcilla.Elmodeloqueseajustaalasemivarianzafueelesférico.Elkrigeadomuestraqueelterrenosepuededividirenunidadesdesueloconaltacoincidenciaaladivisiónconbaseenlafotografíaaérea,queinicialmenteseusóparadividirelterrenoendichasunidadesymuestrear[13].D.VariabilidaddelaspropiedadesdesuelosenEcuadorEnunainvestigaciónrealizadaenlosRíossemuestreópropiedadesfísicasdelsuelodesistemasdeusodelsuelo.Elsueloconpastocultivadoparapastoreopresentócambiosdelaspropiedadesprincipalmenteenlascapassuperficialesdelsuelo,aligualquelosmonocultivos,alsercomparadosconelsuelodebosqueprimario.Loscambiosdelaspropiedadesperjudicaronelmovimientoverticaldelaguaenelsuelo.Segúnelanálisisdecomponentesprincipaleselsuelodelpastizaltienelosmayorescontenidosdehumedadvolumétrica,limo,conductividadhidráulicayrelaciónmateriaorgánica/(limomásarcilla).Asimismo,entodoslossistemasdeusodesueloladensidadaparenteymateriaorgánica,ladensidadrealyarcilladispersaenagua,laconductividadhidráulicaymateriaorgánica/(limomásarcilla)tienenlasmayorescorrelacionespositivas;mientrasquelahumedadvolumétricaconladensidadrealyarcilladispersaenagua,laarenaconladensidadaparenteymateriaorgánicatuvieronlasmayorescorrelacionesnegativas[14].Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal147Porotrolado,enElOrosemuestreóensuelosdestinadosavariosagrosistemasyseencuentraqueladensidadreal,materiaorgánica,delsueloconpastoscultivadosparapastoreoessignificativamentemenorqueelsuelodebosque,banano,cacaoymaíz;mientrasquelacapacidaddeintercambiocatiónicoyconductividadeléctricadelsuelodepastossonsignificativamentemayoresquelosdemáscultivosnombrados[15].Porúltimo,enPastazaseobservóqueelsuelodestinadoalcultivodepastostuvodiferenciassignificativasdepropiedadesfísicas,químicasybiológicasalsercomparadasconelsuelodebosqueychacras.Tambiénseevaluólossuelosconuníndicedecalidadyreportanquelosdestinadosalospastosychacrasengeneraltienenunacalidadlevementesuperior;sinembargo,laschacrasybosquetienenmayorcalidadenlacapasuperficialdelsuelo.Lossuelosde10-30cmdeprofundidaddelospastostienenlosmayoresvaloresdedensidadaparente,porosidadderetenciónypH;mientrasqueaprofundidadesde0-10cmelaluminio,potasio,magnesioycalciotienenlosmayoresvalores.Porotrolado,laconductividadhidráulicasaturadaylaporosidaddeaireación,elpotasio,calcioymagnesio,ladensidadaparenteyporosidadderetencióntienenlasmayorescorrelacionespositivas;mientrasqueladensidadaparenteconlaconductividadhidráulicasaturadayporosidaddeaireación,ylaporosidadderetenciónconlamateriaorgánicatienenlasmayorescorrelacionesnegativas[16].CONCLUSIONESSeharealizadounarevisiónpreliminardeliteraturaacercadelmuestreodelavariabilidadespacialdepropiedadesdelsuelodepastizales.Anivelmundialexistendiseñosdemuestreodesuelosconpastizalesqueevalúanlavariabilidadespacialdelsuelo;mientrasqueanivelnacional(Ecuador)losmuestreosrevisadosnoestudianlavariabilidadespacialdelsuelo(relacióndelaspropiedadesconsuscoordenadasdondefuemedida),generalmentesemuestreadentrodeparcelasdeunterrenoparaestimarparámetrosdescriptivosdelaspropiedadesdelsuelo,comparacionesentrediferentesusosdelsuelo,correlacionesorelacionesentrepropiedades.Losdiseñosdelmuestreodelavariabilidadespacialdelsueloempleadossonencuadrículasregulares,transectososinunpatrónregulardemuestreo.Losmuestreosencuadrículasregularestienende8mhasta40mdeladodecuadrículaylospuntosdemuestreosongeneralmenteenlasinterseccionesdelacuadrícula.Mientrasquelosmuestreosentransectososinunpatrónregulartienendistanciasde4ma50mentresí,conpuntosdemuestreoespaciadosde12ma113m.Elnúmerodepuntosdemuestreovade40a88distribuidosenáreasde0,45hahasta700ha.Laprofundidaddemuestreovade0-90cm,encapasdealmenos5cmdeespesor.Lastécnicasusadasparaestudiarlavariabilidadespacialdelsuelosonelsemivariogramaajustadoalosmodelosesférico,gausianoyexponencial,yelkrigeadoordinarioyuniversal.Engeneral,laspropiedadesdelsueloempiezanatenerindependenciaespacialapartirdelos35mdedistancia.Finalmente,esimportantellevaracaboenEcuadorinvestigacionesdelavariabilidadespacialdepropiedadesdelsuelodepastizales,quepermitantomardecisionesfiablesparasumanejoyconservación.REFERENCIAS[1]D.Silva,M.Costa,B.Mantovanelli,L.Coutrim,M.Rodrigues,J.Cunha,“VariabilidadeespacialdaemissãodeCO2,temperaturaeumidadedosoloemáreadepastagemnaregiãoAmazônica,Brasil”,RevistadeCiênciasAgroveterinárias,vol.18,no.1,pp.119-126,noviembre2019,[Enlí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íaparaelanálisisdedatosespacialesenagricultura”,1ªed.,SerieEstadísticaAplicada,Córdoba:Brujas,2019,[Enlí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.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal148[3]E.Tola,K.Gaadi,C.Biradar,R.Madugundu,A.Zeyada,A.Kayad,“CharacterizationofspatialvariabilityofsoilphysicochemicalpropertiesanditsimpactonRhodesgrassproductivity”,SaudiJournalofBiologicalSciences,vol.24,pp.421-429,abril2017,[Enlí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.deBrito,R.Gomes,2019,“SpatialvariabilityandsamplingdensityofchemicalattributesinarchaeologicalblackearthsunderpastureinsouthernAmazonas,Brazil”,SOILDiscuss,junio2019,[Preimpreso],[Enlínea],https://doi.org/10.5194/soil-2019-26[5]R.Petrucci,F.Herring,J.Madura,C.Bissonnette,“Químicageneral”,Madrid:Pearson,2011.[6]C.Grego,C.Gonçalves,S.Nogueira,F.Andrade,A.deOliveira,C.Ferrer,A.dosSantos,J.deAbreu,“Variabilidadeespacialdosoloedabiomassaepígeadepastagem,identificadapormeiodegeostatística”,Pesquisaagropecuariabrasileira,vol.47,np.9,pp.1404-1412,septiembre2012,[Enlínea],https://doi.org/10.1590/S0100-204X2012000900026[7]A.Bernardi,K.Santos,G.Bettiol,L.Rabello,R.Ferreira,R.Inamasu,“SpatialvariabilityofsoilpropertiesandyieldofagrazedalfalfapastureinBrazil”,PrecisionAgriculture,vol.17,pp.737-752,abril2016,[Enlínea],https://doi.org/10.1007/s11119-016-9446-9.[8]A.Delgado,J.Gómez,“TheSoil.Physical,ChemicalandBiologicalProperties”,enPrinciplesofAgronomyforSustainable,Springer,2016.[Enlínea],https://link.springer.com/chapter/10.1007/978-3-319-46116-8_2[9]S.Barrezueta,K.Velepucha,L.Hurtado,E.Jaramillo,“Soilpropertiesandstorageoforganiccarboninthelandusepastureandforest”,RevistadeCienciasAgrícolas,vol.36,no.2,pp.31-45,abril2019,[Enlínea],https://doi.org/10.22267/rcia.193602.116[10]J.Wu,J.He,G.Christakos,“Quantitativeanalysisandmodelingofearthandenvironmentaldata.Space-timeandspacetimedataconsiderations”,India:Elsevier,2022,[Enlínea],https://doi.org/10.1016/C2018-0-00055-9[11]M.Oliver,R.Webster,“BasicStepsinGeostatistics:TheVariogramandKriging”,Springer,2015,[Enlí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.“UsingSoilApparentElectricalConductivitytoOptimizeSamplingofSoilPenetrationResistanceandtoImprovetheEstimationsofSpatialPatternsofSoilCompaction”,TheScientificWorldJournal,pp.1-12,diciembre2014,[Enlínea],http://dx.doi.org/10.1155/2014/269480[13]A.González,J.Rey,J.Atencio,”Variabilidadespacialdelossuelosdelaunidadexperimentallaglorietaconfinesagropecuarios”,RevistadelaFacultaddeAgronomía(LUZ),suplemento1,pp.539-553,junio2014,[Enlinea],https://www.revfacagronluz.org.ve/PDF/suplemento_2014/ing/ingsupl12014539553.pdf[14]NovilloI.;Carrillo,M.;Cargua,J.;Moreira,V.;Albán,K.;Morales,K,“PropiedadesfísicasdelsueloendiferentessistemasagrícolasenlaprovinciadeLosRíos,Ecuador”.TemasAgrarios,vol.23,no.2,pp.177-187,mayo2018,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.“Efectodelmanejoagrícolaenpropiedadesfísicasyquímicasdelsueloendiferentesagroecosistemas”,UniversidadySociedad,vol.12,no.5,pp.389-398,septiembre2020,[Enlí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.“AsoilqualityindexforsevenproductivelandscapesintheAndean-AmazonianfoothillsofEcuador”,LandDegradation&Development,vol.32,pp.2226-2241,enero2021,[Enlínea],https://doi.org/10.1002/ldr.3897[17]J.Serrano,S.Shahidian,J.Silva,“SpatialvariabilityandtemporalstabilityofapparentsoilelectricalconductivityinaMediterraneanpasture”,PrecisionAgriculture,vol.18,pp.245-263,julio2017,[Enlínea],https://doi.org/10.1007/s11119-016-9460-y[18]D.Brus,J.deGruijter,“Randomsamplingorgeostatisticalmodelling?Choosingbetweendesign-basedandmodel-basedsamplingstrategiesforsoil(withdiscussion)”,Geoderma,vol.80,pp.1-44,octubre1997,[Enlí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.IssueN°1,(pp.138-149) ISSN-E:2697-3650MinervaJournal149[19]P.Lawrence,W.Roper,T.Morris,K.Guillard,“Guidingsoilsamplingstrategiesusingclassicalandspatialstatistics:Areview”,AgronomyJournal,vol.112,pp.493-510,enero2020.[Enlínea],https://doi.org/10.1002/agj2.20048.[20]J.Wang,A.Stein,B.Gao,Y.Ge,“Areviewofspatialsampling”,SpatialStatistics,vol.2,pp.1-14,agosto2012,[Enlínea],https://doi.org/10.1016/j.spasta.2012.08.001.RodrigoSaquicela,ingenieroagropecuario,especializadoennutriciónvegetal.MaestrantedelInstitutoPosgradodelaUniversidadTécnicadeManabí,UTM,Portoviejo,Ecuador.DocentedelaUniversidadUTE.LellyMaríaUsecheCastro,ingenieroindustrial,PhD.esestadística.DirectoradedelGrupodeAnálisisMultivarianteyEstocástico(G.A.M.E).InstitutodeCienciasBásicas,DepartamentodeMatemáticasyEstadística.UniversidadTécnicadeManabí,Ecuador.AnaFranciscaGonzálezPedraza,IngenieraAgrónoma,Magisterenmanejodelosrecursosaguaysuelo.DoctoraenCiencias,menciónEcología.DocentedelaUniversidaddePamplona,Pamplona.NortedeSantander,Colombia.Saquicela et al. Muestreo de la variabilidad espacial de propiedades del suelo de pastizalesVol.3,Spec.IssueN°1,(pp.138-149)