<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article
  PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1 20151215//EN" "https://jats.nlm.nih.gov/publishing/1.1/JATS-journalpublishing1.dtd">
<article article-type="research-article" dtd-version="1.1" specific-use="sps-1.8" xml:lang="es" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
	<front>
		<journal-meta>
			<journal-id journal-id-type="publisher-id">Avances</journal-id>
			<journal-title-group>
				<journal-title>Avances: Investigación en Ingeniería</journal-title>
				<abbrev-journal-title abbrev-type="publisher">Avances</abbrev-journal-title>
			</journal-title-group>
			<issn publication-format="print">1794-4953</issn>
            <issn publication-format="electronic">2619-6581</issn>
			<publisher>
				<publisher-name>Universidad Libre</publisher-name>
			</publisher>
		</journal-meta>
		<article-meta>
			<article-id pub-id-type="doi">10.18041/1794-4953/avances.1.5727</article-id>
			<article-categories>
				<subj-group subj-group-type="heading">
					<subject>Artículos</subject>
				</subj-group>
			</article-categories>
			<title-group>
				<article-title>Solución analítica de la ecuación de Richards en suelos no saturados basada en el modelo de Gardner. Caso de difusividad constant</article-title>
				<trans-title-group xml:lang="en">
					<trans-title>Analytical Solution of Richards' Equation in Unsaturated Soils based on Gardner's Model. Constant Diffusivity case study</trans-title>
				</trans-title-group>
			</title-group>
            <contrib-group>
				<contrib contrib-type="author" corresp="no">
					<name>
						<surname>Sánchez Quispe</surname>
						<given-names>Hugo Rolando</given-names>
					</name>
                    <xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
					<email>hugo.sanchez@espoch.edu.ec</email>
				</contrib>
				<contrib contrib-type="author" corresp="no">
                    <name>
						<surname>Chipantiza</surname>
						<given-names>Juan Gabriel</given-names>
					</name>
                    <xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
					<email>juan.chipantiza@espoch.edu.ec</email>
				</contrib>
                <contrib contrib-type="author" corresp="no">
					<name>
						<surname>Bonilla</surname>
						<given-names>Amanda</given-names>
					</name>
                    <xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
					<email>amanda.bonilla@espoch.edu.ec</email>
				</contrib>
                <contrib contrib-type="author" corresp="no">
					<name>
						<surname>Mendoza</surname>
						<given-names>Julio César</given-names>
					</name>
                    <xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
					<email>cmendoza@ucla.edu.ve</email>
				</contrib>
			</contrib-group>
            <aff id="aff1">
                <label>1</label>
                <institution content-type="original">Ingeniero automotriz, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch</institution>
                <institution content-type="normalized">Ingeniero automotriz, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch </institution>
                <institution content-type="orgname">Ingeniero automotriz, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch </institution>
                <addr-line>
                <city>Riobamba </city>
                </addr-line>
                <country country="EC">Ecuador</country>
            </aff>
             <aff id="aff2">
                <label>2</label>
                <institution content-type="original">Ingeniero agrícola, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch </institution>
                <institution content-type="normalized">Ingeniero agrícola, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch </institution>
                <institution content-type="orgname">Ingeniero agrícola, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch </institution>
                <addr-line>
                <city>Riobamba</city>
				</addr-line>
                <country country="EC">Ecuador</country>
             </aff>
             <aff id="aff3">
                <label>3</label>
                <institution content-type="original">Ingeniero agropecuaria, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch.</institution>
                <institution content-type="normalized">Ingeniero agropecuaria, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch.</institution>
                <institution content-type="orgname">Ingeniero agropecuaria, MSc, Profesor Investigador, Escuela Superior Politécnica del Chimborazo -Espoch.</institution>
                <addr-line>
                <city>Riobamba</city>
				</addr-line>
                <country country="EC">Ecuador</country>
             </aff>
             <aff id="aff4">
                <label>4</label>
                <institution content-type="original">Ingeniero Mecánico, MSc, Profesor Investigador Universidad Centroccidental Lisandro Alvarado</institution>
                <institution content-type="normalized">Ingeniero Mecánico, MSc, Profesor Investigador Universidad Centroccidental Lisandro Alvarado</institution>
                <institution content-type="orgname">Ingeniero Mecánico, MSc, Profesor Investigador Universidad Centroccidental Lisandro Alvarado</institution>
                <addr-line>
                <city>Barquesimeto</city>
				</addr-line>
                <country country="VE">Venezuela</country>
             </aff>
			 <pub-date pub-type="collection">
				<season>Jan-Jun</season>
				<year>2020</year>
			</pub-date>
			<volume>17</volume>
			<issue>01</issue>
			<fpage>1</fpage>
			<lpage>19</lpage>
			<history>
				<date date-type="received">
					<day>28</day>
					<month>05</month>
					<year>2019</year>
				</date>
				<date date-type="accepted">
					<day>10</day>
					<month>03</month>
					<year>2020</year>
				</date>
			</history>
			<permissions>
				<license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by-nc-sa/4.0" xml:lang="es">
					<license-p>Esta obra está bajo licencia internacional Creative Commons Atribución-NoComercial-CompartirIgual 4.0 Internacional (CC BY-NC-SA 4.0)</license-p>
				</license>
			</permissions>
			<abstract>
			<title>RESUMEN</title>
				<p>El artículo se centra en encontrar la solución analítica de la ecuación de Richards, usada para modelar el movimiento del agua en un medio poroso no saturado y para restringir el estudio a la circulación vertical descendente del fluido, con fuentes de humedad nulas y difusividad constante en el suelo. Se empleó la relación constitutiva de Gardner. Los resultados arrojados muestran que empleando la relación de Gardner se consigue una solución para el contenido de humedad en el suelo de tipo exponencial, tanto en la componente espacial como en la temporal y dependiente del autovalor γ, el cual restringe las soluciones particulares del modelo de Gardner a un conjunto determinado de valores.</p>
			</abstract>
			<trans-abstract xml:lang="en">
			<title>ABSTRACT</title>
			    <p>The main objective of this research was to find the analytical solution of the Richards equation, which is used to model the movement of water in an unsaturated porous medium, restricting the study to the downward vertical circulation of the fluid, with sources of moisture null and constant diffusivity in the soil.  The results indicate that using  Gardner's relationship, a solution for the moisture content in the soil of an exponential type is achieved, both in the spatial and temporal component,  reliant on the self-value, which restricts the particular solutions from  Gardner's model to a certain set of values. The main objective of this work was to find the analytical solution of the Richards equation, which is used to model the movement of water in an unsaturated porous medium, restricting the study to the downward vertical circulation of the fluid, with sources of moisture null and constant diffusivity in the soil. The constitutive relationship of Gardner was used. The results show that, using the Gardner relationship, a solution for the moisture content in the soil of an exponential type is achieved, both in the spatial and temporal component and dependent on the self-value γ, which restricts the particular solutions from the Gardner model to a certain set of values.</p>
			</trans-abstract>
            <kwd-group xml:lang="es">
			<title>Palabras clave</title>
				<kwd>infiltración</kwd>
				<kwd>movimiento del agua</kwd>
				<kwd>medio poroso</kwd>
                <kwd>permeabilidad</kwd>
			</kwd-group>
			<kwd-group xml:lang="en">
				<title>Keywords</title>
				<kwd>infiltration</kwd>
				<kwd>permeability</kwd>
				<kwd>porous medium</kwd>
				<kwd>water movement</kwd>
			</kwd-group>
            <counts>
				<fig-count count="25"/>
				<table-count count="1"/>
				<equation-count count="0"/>
				<ref-count count="27"/>
				<page-count count="8"/>
			</counts>
		</article-meta>
	</front>
    <body>
	    <sec sec-type="intro">
		    <title>Introduction</title>
            <p>Al flujo de agua en el suelo lo regulan una serie de factores físicos, entre los cuales se destacan la gravedad, la permeabilidad del suelo, la humedad, la densidad y la viscosidad del agua. Determinar el movimiento del agua en el suelo es importante, dado que así se establece con precisión la capacidad de agua disponible para los cultivos [<xref  ref-type="bibr" rid="r1">1</xref>], la capacidad de recarga de un acuífero [<xref  ref-type="bibr" rid="r2">2</xref>] y los riesgos de lixiviación e inundación [<xref  ref-type="bibr" rid="r3">3</xref>]. Ello permitirá una adecuada planificación y diseño de las obras de riego y mitigar los problemas de contaminación de acuíferos. </p>
            <p>El movimiento de agua y otros fluidos en el suelo es regido por la Ley de Darcy [<xref  ref-type="bibr" rid="r4">4</xref>], y partir de ella se han desarrollado ecuaciones como la de Richards [<xref  ref-type="bibr" rid="r5">5</xref>], que se puede resolver matemáticamente mediante modelos como el de Gardner [<xref  ref-type="bibr" rid="r6">6</xref>] y Brooks-Corey [<xref  ref-type="bibr" rid="r7">7</xref>], cuyo propósito es que a través de la modelización matemática se realice una estimación óptima de los parámetros hidrológicos del suelo. No obstante, las soluciones analíticas no han sido precisas, por lo que los valores obtenidos no describen adecuadamente el flujo de agua en el suelo. La necesidad de encontrar modelos matemáticos para estimar los parámetros hidrológicos del suelo radica es que los métodos de campo (p. ej., los infiltrómetros) o los de laboratorios (p. ej., los permeámetros) presentan una alta variabilidad, y no son confiables los resultados.</p>
            <p>El suelo y el agua son dos recursos fundamentales en las ciencias agrícolas y en otras muchas especialidades [<xref  ref-type="bibr" rid="r8">8</xref>]. Estructuralmente, los suelos constituyen un medio particulado y poroso, que se define como un material conformado por una parte sólida o matriz y unos espacios o poros, los cuales pueden ser ocupados por una o varias fases de fluido, ya sean líquidas o gaseosas [<xref  ref-type="bibr" rid="r9">9</xref>].</p>
            <p>Al estar en la superficie del suelo, el agua infiltra hacia su interior y ocupa las cavidades intersticiales obligadas por fuerzas de distintas naturaleza, intensidad y dirección. El movimiento del agua en el suelo ocurre en el caso saturado y en el no saturado. El flujo de agua en un suelo no saturado es el que ocurre sujeto a cualquier condición de humedad menor a la de saturación. </p>
            <p>En la fase líquida de un suelo no saturado, tanto el agua como el aire se movilizan por la influencia de los potenciales presentes en este medio. Se admite que el flujo no saturado lo describe la Ley de Darcy, con la diferencia de que la conductividad hidráulica no saturada es inconstante y que varía, predominantemente, en función de la humedad y, en consecuencia, de la succión matricial [<xref  ref-type="bibr" rid="r10">10</xref>].</p>
            <p>La conductividad hidráulica (Ks) representa la mayor o la menor facilidad con que un medio poroso permite el paso de un fluido a través de este por unidad de área transversal a la dirección del flujo. Puede expresarse bajo condiciones saturadas y no saturadas. Para el caso saturado, la conductividad hidráulica se representa por un solo valor que se mide en condiciones de saturación máxima del suelo. </p>
            <p>Para el caso no saturado, la conductividad hidráulica varía con el contenido de humedad. No depende de la presión de succión, sino del contenido de humedad [<xref  ref-type="bibr" rid="r11">11</xref>], y es también función de las características del fluido (viscosidad y peso específico) y de las propiedades geométricas del medio poroso [<xref  ref-type="bibr" rid="r12">12</xref>]. Esta función es altamente no lineal con el contenido de humedad, por lo que, junto a las curvas de succión, introduce grandes complicaciones a la solución de la ecuación de Richards [<xref  ref-type="bibr" rid="r13">13</xref>]. </p>
            <p>La conductividad hidráulica de un suelo no saturado varía con su contenido de humedad. A medida que se infiltra el agua en el perfil del suelo, se llega a condiciones de saturación en las cuales estas adquieren un valor constante [<xref  ref-type="bibr" rid="r14">14</xref>]. </p>
            <p>Matemáticamente, la conductividad hidráulica se expresa como:</p>
			<p>
				<fig id="ec1">
					<label>Ecuación</label>
					<caption>
						<title>1 y 2</title>
					</caption>
					<graphic xlink:href="art5727_ec1y2.png" />
				</fig>
			</p>
            <p>Dónde: k es la permeabilidad intrínseca del suelo; η es la viscosidad del fluido; ρ es la densidad; g es aceleración de la gravedad y θ es la humedad presente en el suelo.</p>
            <p>Entre tanto, la difusividad se interpreta como la transferencia o desplazamiento de moléculas individuales a través de un fluido por medio de los desplazamientos individuales y desordenados de las moléculas, debido a una diferencia de concentraciones. La difusión ocurre en sistemas de fluidos estancados o en fluidos que se están moviendo. Es un coeficiente relativo al comportamiento de un suelo en relación con su capacidad para transmitir agua, la cual a su vez depende las condiciones texturales del suelo y a su porosidad [<xref  ref-type="bibr" rid="r15">15</xref>].</p>
            <p>La conductividad (K) y la difusividad hidráulica (D) para un medio poroso no saturado se relacionan por medio de la ecuación (3):</p>
			<p>
				<fig id="ec3">
					<label>Ecuación</label>
					<caption>
						<title>3</title>
					</caption>
					<graphic xlink:href="art5727_ec3.png" />
				</fig>
			</p>
            <p>El objetivo de investigación fue desarrollar una solución analítica, a través de la ecuación de Richards, mediante el modelo matemático de Gardner, a fin de desarrollar una propuesta que describa de modo adecuado el flujo de agua en suelos no saturados y estime de manera precisa los parámetros hidrológicos.</p>
        </sec>
        <sec> 
            <title>1. Metodología</title>
            <p>Para resolver la ecuación de Richards se empleó el modelo de Gardner, y previo a la demostración matemática y a la resolución de la ecuación se describen cada uno de estos modelos</p>
            <sec>
                <title>1.1. Ecuación de Richards y modelo de Gardner</title>
				<sec>
                	<title>1.1.1. Ecuación de Richards</title>
                	<p>La ecuación fue formulada por Richards [<xref  ref-type="bibr" rid="r16">16</xref>], quien consolidó los aportes realizados por un grupo importante de investigadores ligados al estudio de la física de suelos, formulando una teoría macroscópica general relacionada con el movimiento del agua en un medio poroso no saturado.</p>
                	<p>La formulación matemática de la ecuación obtenida por Richards deriva de la combinación de un balance de masa sustentado en la ecuación de continuidad y el principio físico del momentum expresado en la Ley de Darcy. </p>
                	<p>La teoría de Richards se puede formular en el marco de la moderna “teoría del continuo de las mezclas”, a condición de que se reconozca desde el principio la existencia de la separación del medio sólido, del líquido y de la fase gaseosa [<xref  ref-type="bibr" rid="r17">17</xref>]. Esta teoría rige la dinámica de un fluido incomprensible que fluye en un medio poroso isotérmico, tridimensional, parcialmente saturado, anisotrópico y heterogéneo.</p>
                	<p>El modelo matemático que dedujo Richards es una ecuación diferencial en derivadas parciales, parabólicas y no lineales [<xref  ref-type="bibr" rid="r18">18</xref>] que se expresa: </p>
					<p>
						<fig id="ec4">
							<label>Ecuación</label>
							<caption>
								<title>4</title>
							</caption>
							<graphic xlink:href="art5727_ec4.png" />
						</fig>
					</p>	
                	<p>Dónde: K es la conductividad hidráulica; h es el potencial del agua en el suelo; t es el tiempo; θ es la humedad presente en el suelo; TA es un tensor adimensional de segundo orden que representa la naturaleza anisótropa del medio, y S_i es una fuente exterior.</p>
            	</sec>
				<sec>
					<title>1.1.2. Modelo de Gardner</title>
					<p>Originalmente, la ecuación de Gardner [<xref  ref-type="bibr" rid="r6">6</xref>] se propuso para evaluar el coeficiente de permeabilidad en un suelo no saturado y también fue una de las primeras empleadas para representar la curva característica suelo-agua. Sin embargo, la forma matemática propuesta para la ecuación de conductividad se ha deducido para inferir el contenido de humedad en función de la succión del suelo [<xref  ref-type="bibr" rid="r18">18</xref>]</p>
					<p>
						<fig id="ec5">
							<label>Ecuación</label>
							<caption>
								<title>5 y 6</title>
							</caption>
							<graphic xlink:href="art5727_ec5y6.png" />
						</fig>
					</p>
					<p>Dónde: Ks es la conductividad hidráulica para el caso saturado; θr es el contenido de agua residual presente en el suelo; θs es la razón entre el volumen de agua en el interior de los poros del suelo y el volumen total de vacíos. </p>
					<p>h_c es el valor relacionado con la presión de entrada de aire al suelo, y α_g es un valor que caracteriza el tipo de suelo</p>
				</sec>
			</sec>	
			<sec>
				<title>1.2. Formulación matemática</title>
				<p>En la <xref  ref-type="fig" rid="f1">figura 1</xref> se muestran las condiciones utilizadas en el estudio.</p>
				<p>
					<fig id="f1">
						<label>Figura 1</label>
						<caption>
							<title>Esquema de las condiciones utilizadas en el estudio</title>
						</caption>
						<graphic xlink:href="art5727_f1.png" />
					</fig>
				</p>
				<p>Sea Ω ⊂ R3 el espacio donde se desarrolla la infiltración del fluido y ∂Ω el entorno de Ω. El contenido de humedad en el suelo se define como θ(x_i,t), donde x_i∈Ω, con i = 1, 2, 3 y t > 0.</p>
				<p>El proceso de infiltración del agua en un medio poroso no saturado se rige por la ecuación de Richards, la cual, como se dijo, deriva de la combinación de la ley de conservación de la masa y del campo de velocidades calculado conforme a la Ley de Darcy [<xref  ref-type="bibr" rid="r19">19</xref>].</p>
				<p>En términos matemáticos, la ley de conservación de la masa se escribe como:</p>
				<p>
					<fig id="ec7">
						<label>Ecuación</label>
						<caption>
							<title>7</title>
						</caption>
						<graphic xlink:href="art5727_ec7.png" />
					</fig>
				</p>
				<p>Mientras que la Ley de Darcy expone que:</p>
				<p>
					<fig id="ec8">
						<label>Ecuación</label>
						<caption>
							<title>8</title>
						</caption>
						<graphic xlink:href="art5727_ec8.png" />
					</fig>
				</p>
				<p>Al sustituir la ecuación (6) en la (5) y desarrollando convenientemente, se obtiene la ecuación:</p>
				<p>
					<fig id="ec81">
						<label>Ecuación</label>
						<caption>
							<title>8 (1)</title>
						</caption>
						<graphic xlink:href="art5727_ec81.png" />
					</fig>
				</p>
				<p>Al expandir el primer miembro de esta ecuación y teniendo en cuenta que el valor del potencial en la zona no saturada se reduce a sus componentes matricial y gravitatoria, se consigue una expresión para la ecuación de Richards de la forma:</p>
				<p>
					<fig id="ec9">
						<label>Ecuación</label>
						<caption>
							<title>9</title>
						</caption>
						<graphic xlink:href="art5727_ec9.png" />
					</fig>
				</p>
				<p>Al restringir el estudio a un movimiento del agua exclusivamente vertical (xi = xj) = z y sin fuentes o sumideros (S = 0), se logra reducir el modelo de tres tridimensional a una dimensión, tal como se observa al escribir la ecuación (7):</p>
				<p>
					<fig id="ec91">
						<label>Ecuación</label>
						<caption>
							<title>9 (1)</title>
						</caption>
						<graphic xlink:href="art5727_ec91.png" />
					</fig>
				</p>
				<p>Desarrollando la derivada parcial de z se consigue:</p>
				<p>
					<fig id="ec92">
						<label>Ecuación</label>
						<caption>
							<title>9 (2)</title>
						</caption>
						<graphic xlink:href="art5727_ec92.png" />
					</fig>
				</p>
				<p>Ahora:</p>
				<p>
					<fig id="ec93">
						<label>Ecuación</label>
						<caption>
							<title>9 (3)</title>
						</caption>
						<graphic xlink:href="art5727_ec93.png" />
					</fig>
				</p>
				<p>Entonces:</p>
				<p>
					<fig id="ec10">
						<label>Ecuación</label>
						<caption>
							<title>10</title>
						</caption>
						<graphic xlink:href="art5727_ec10.png" />
					</fig>
				</p>
			</sec>	
			<sec>
				<title>1.3. Solución analítica</title>
				<p>La ecuación (10) es altamente no lineal, lo cual impide encontrarle una solución exacta; entonces, para resolver y encontrar un resultado analítico, se debe considerar el caso hipotético en el cual la difusividad D sea del orden C0 respecto a z, lo cual se basa en los bajos valores de difusividad que son de 0,005 m2 s-1 en suelos arcillosos y de 0,008 m2 s-1 en suelos arenosos [<xref  ref-type="bibr" rid="r20">20</xref>]. Dado que los valores son extremadamente pequeños, se pueden asumir los valores de D como constantes y así eliminar la no linealidad de la ecuación como lo planteó [<xref  ref-type="bibr" rid="r21">21</xref>]. Así, tenemos que:</p>
				<p>
					<fig id="ec101">
						<label>Ecuación</label>
						<caption>
							<title>10 (1)</title>
						</caption>
						<graphic xlink:href="art5727_ec101.png" />
					</fig>
				</p>
				<p>Con lo que ecuación (10) se reduce a:</p>
				<p>
					<fig id="ec11">
						<label>Ecuación</label>
						<caption>
							<title>11</title>
						</caption>
						<graphic xlink:href="art5727_ec11.png" />
					</fig>
				</p>
				<p>De las ecuaciones (5) y (6) se consigue que:</p>
				<p>
					<fig id="ec111">
						<label>Ecuación</label>
						<caption>
							<title>11 (1)</title>
						</caption>
						<graphic xlink:href="art5727_ec111.png" />
					</fig>
				</p>
				<p>y</p>
				<p>
					<fig id="ec12">
						<label>Ecuación</label>
						<caption>
							<title>12</title>
						</caption>
						<graphic xlink:href="art5727_ec12.png" />
					</fig>
				</p>
				<p>Al sustituir (12) en (11) se obtiene:</p>
				<p>
					<fig id="ec13">
						<label>Ecuación</label>
						<caption>
							<title>13</title>
						</caption>
						<graphic xlink:href="art5727_ec13.png" />
					</fig>
				</p>
				<p>Sujeta a la condición inicial:</p>
				<p>
					<fig id="ec131">
						<label>Ecuación</label>
						<caption>
							<title>13 (1)</title>
						</caption>
						<graphic xlink:href="art5727_ec131.png" />
					</fig>
				</p>
				<p>y a las condiciones de contorno de Dirichlet:</p>
				<p>
					<fig id="ec132">
						<label>Ecuación</label>
						<caption>
							<title>13 (2)</title>
						</caption>
						<graphic xlink:href="art5727_ec132.png" />
					</fig>
				</p>
				<p>La solución exacta de la ecuación diferencial (13) tiene la forma:</p>
				<p>
					<fig id="ec14">
						<label>Ecuación</label>
						<caption>
							<title>14</title>
						</caption>
						<graphic xlink:href="art5727_ec14.png" />
					</fig>
				</p>
				<p>Desde el punto de vista práctico, con la ecuación (14) se modela el flujo de agua u otro fluido en el suelo, que tiene un comportamiento exponencial, siendo un flujo rápido en principio, hasta alcanzar valores contantes cuando el suelo se ha saturado. </p>
				<p>El modelado matemático permite la estimación adecuada de parámetros hidrológicos como infiltración, riesgos de lixiviación, y con ello se calcula la demanda hídrica de los cultivos [<xref  ref-type="bibr" rid="r22">22</xref>], el tiempo y frecuencia de riego [<xref  ref-type="bibr" rid="r23">23</xref>], el riesgo de escurrimiento y erosión [<xref  ref-type="bibr" rid="r24">24</xref>], así como los riesgos de lixiviación de agroquímicos o hidrocarburos [<xref  ref-type="bibr" rid="r25">25</xref>], en estudios de impacto ambiental</p>
			</sec>	
		</sec>
		<sec>
			<title>2. Resultados y discusión</title>
			<p>Para ejemplificar cómo se aplica la solución analítica, se usó el modelo propuesto para evaluar la ecuación (14) en un suelo tipo franco-limoso, el cual presenta las características mostradas en la <xref  ref-type="table" rid="t1">tabla 1</xref>.</p>
			<p>
				<table-wrap id="t1">
					<label>Tabla 1</label>
					<caption>
						<title>Características hidrodinámicas del suelo tipo franco-limoso G.E.3 [<xref  ref-type="bibr" rid="r16">16</xref>]</title>
					</caption>
					<table>
						<tr>
							<td>0s</td>
							<td>0r</td>
							<td>Ag</td>
							<td>Ks</td>
							<td>he</td>
							<td>D</td>
							<td>L</td>
						</tr>
						<tr>
							<td>0,396</td>
							<td>0,131</td>
							<td>0,0007</td>
							<td>0,207</td>
							<td>0</td>
							<td>5</td>
							<td>100</td>
						</tr>
					</table>
				</table-wrap>
			</p>
			<p>Las unidades de cada uno de estos elementos son:</p>
			<p>
				<fig id="ec15">
					<label>Ecuación</label>
					<caption>
						<title>15</title>
					</caption>
					<graphic xlink:href="art5727_ec15.png" />
				</fig>
			</p>
			<p>En el caso de la difusividad constante y el empleo de la relación de Gardner, la solución de la ecuación de Richards viene dada por la ecuación (14). Considerando los posibles valores de γ, se presentan los siguientes casos evaluados a continuación: </p>
			<p>Caso I: γ = 0</p>
			<p>Con esta condición se consigue que:</p>
			<p>
				<fig id="ec16">
					<label>Ecuación</label>
					<caption>
						<title>16</title>
					</caption>
					<graphic xlink:href="art5727_ec16.png" />
				</fig>
			</p>
			<p>Esta solución es físicamente inaceptable, ya que se puede interpretar como que no hubo infiltración del agua ubicada en la capa superficial del suelo, y se debe recordar que la condición de contorno en la referencia z = 0 toma el valor θ(z = 0) = θs. También se puede interpretar este resultado como que si la capa superficial del suelo tuviese un valor nulo de θ, lo cual sería una contradicción a la condición de partida impuesta para θ en esta referencia.</p>
			<p>Caso II: γ &lt; 0</p>
			<p>Tomando en consideración que los valores de γ sean negativos, tendríamos una solución que crece de modo exponencial. Ahora bien, como el caso estudiado se trata de un problema de infiltración de un líquido (en este caso agua) en el suelo, se espera que con el transcurrir del tiempo el problema transitorio se convierta, a lo sumo, en uno estacionario en que el almacenamiento del agua en las cavidades intersticiales del suelo alcance el valor límite de humedecimiento saturado y no que este crezca de manera exponencial hacia el infinito.</p>
			<p>Caso III: γ &gt; 0</p>
			<p>Para valores positivos de γ, el término exp(-γt) experimenta un decaimiento en su valor, lo cual es indicativo de estar en presencia de una solución física aceptable, pues describe el comportamiento de un suelo no saturado, que se va humedeciendo rápidamente. Así, los valores de infiltración y permeabilidad en los primeros minutos y en los primeros centímetros son altos, pero cuando el espacio poroso pasa a una condición saturada, estos valores se hacen constantes y los potenciales de presión en el suelo tienden a 0.</p>
			<p>La <xref  ref-type="fig" rid="f2">figura 2</xref> muestra el comportamiento del proceso de humedecimiento del suelo a medida que el agua se va posicionando de los poros al interior del suelo. La zona superior de la figura, que corresponde a la parte superficial del suelo y a sus estratos más vecinos, se conoce como zona de saturación y transmisión, la cual alcanza una profundidad de 80 cm, y a partir de esta profundidad se inicia la admisión del agua hacia las partes más interiores del suelo, lo que lleva a la saturación de todo el espacio poros; mientras que la tensión se hace prácticamente nula. </p>
			<p>
				<fig id="f2">
					<label>Figura 2</label>
					<caption>
						<title>Evolución de la humedad en función de la profundidad para el modelo de Gardner cuando γ = 1</title>
					</caption>
					<graphic xlink:href="art5727_f2.png" />
				</fig>
			</p>
			<p>La otra zona corresponde a los últimos 20 cm (<xref  ref-type="fig" rid="f3">figura 3</xref>), donde el contenido de humedad es menor que el observado en la primera zona; pero que tiende a aumentar. Debido al proceso de infiltración, el tiempo considerado para la ocurrencia del proceso de humedecimiento en ambas dos zonas fue de aproximadamente una hora.</p>
			<p>
				<fig id="f3">
					<label>Figura 3</label>
					<caption>
						<title>Evolución de la tensión en función de la humedad para el modelo de Gardner cuando γ = 0</title>
					</caption>
					<graphic xlink:href="art5727_f3.png" />
				</fig>
			</p>
			<p>La importancia de la determinación mediante un modelo matemático se debe a que resultan costosos los métodos por parámetros hidrológicos en campo [<xref  ref-type="bibr" rid="r15">15</xref>]. Entre tanto, por medio de la determinación de la difusividad en condiciones controladas y de la formulación de modelos matemáticos basados en la ecuación de Richards, se han hecho determinaciones precisas de infiltración [<xref  ref-type="bibr" rid="r26">26</xref>], conductividad hidráulica [<xref  ref-type="bibr" rid="r27">27</xref>] y sortividad [<xref  ref-type="bibr" rid="r15">15</xref>], cuyos valores se usan para el cálculo con fines de sistemas de riego, movimiento de fluidos en el área de hidrocarburos y riesgos de contaminación por lixiviación con agroquímicos.</p>
		</sec>
		<sec>
			<title>Conclusiones</title>
			<p>La ecuación de Richards empleada para estimar la conductividad hidráulica en el suelo es una función altamente no lineal, por lo que su solución analítica solo puede conseguirse si se reduce su no linealidad, imponiendo ciertas restricciones que idealicen, entre otras, las propiedades físicas y mecánicas del suelo infiltrado por el agua. Para ello se consideró que la difusividad sea constante, lo cual permitió eliminar la no linealidad de la ecuación de Richards y posibilitó conseguirle una solución exacta, con lo que se estimarían de manera precisa los parámetros hidrológicos del suelo. </p>
			<p>Al evaluar la ecuación de Richards con el modelo para la conductividad hidráulica propuesto por Gardner, y considerando un suelo con ciertas condiciones iniciales y de contorno, se obtuvo la solución descrita por medio de la ecuación (14).</p>
			<p>Una solución particular conseguida para la ecuación (14) muestra que la respuesta para el contenido de humedad en el suelo es de tipo exponencial, tanto en la componente espacial como en la temporal.</p>
		</sec>
	</body>
	<back>
		<ref-list>
			<title>Referencias</title>
			<ref id="r1">
				<mixed-citation>E. G. Carnero Carnero y G. G. Carnero Guzmán, “Análisis numérico de infiltración de agua en un suelo residual no saturado utilizando la curva de retención de agua”, Rev. Investig. Altoandinas, vol. 20, n.º 4, pp. 439-450, 2018. <ext-link xlink:href="http://dx.doi.org/10.18271/ria.2018.421">http://dx.doi.org/10.18271/ria.2018.421</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Carnero Carnero</surname>
							<given-names>E. G.</given-names>
						</name>
						<name>
							<surname>Carnero Guzmán</surname>
							<given-names>G. G.</given-names>
						</name>
					</person-group>
					<year>2018</year>
					<article-title>Análisis numérico de infiltración de agua en un suelo residual no saturado utilizando la curva de retención de agua</article-title>
					<source>Rev. Investig. Altoandinas</source>
					<volume>20</volume>
					<issue>4</issue>
					<fpage>430</fpage>
					<lpage>450</lpage>
					<ext-link ext-link-type="uri">http://dx.doi.org/10.18271/ria.2018.421</ext-link>
					<pub-id pub-id-type="doi">10.18271/ria.2018.421</pub-id>
				</element-citation>
			</ref>
            <ref id="r2">
               <mixed-citation>R. Aguilar-García y M. A. Ortega-Guerrero, “Análisis de la dinámica del agua en la zona no saturada en un suelo sujeto a prácticas de conservación: implicaciones en la gestión de acuíferos y adaptación al cambio climático”, Rev. Mex. Cienc. Geol., vol. 34, n.º 2, pp. 91-104, 2017</mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Aguilar-García</surname>
							<given-names>R</given-names>
						</name>
						<name>
							<surname>Ortega-Guerrero</surname>
							<given-names>M .A.</given-names>
						</name>
                    </person-group>
					<year>2017</year>
					<article-title>Análisis de la dinámica del agua en la zona no saturada en un suelo sujeto a prácticas de conservación: implicaciones en la gestión de acuíferos y adaptación al cambio climático</article-title>
					<source>Rev. Mex. Cienc. Geol.</source>
					<volume>34</volume>
					<issue>2</issue>
					<fpage>91</fpage>
					<lpage>104</lpage>
				</element-citation>
			</ref>
            <ref id="r3">
                <mixed-citation>F. Zamora, N. Rodríguez, D. Torres y H. Yendis, “Efecto del riego con aguas residuales sobre propiedades químicas de suelos de la planicie de Coro, Estado Falcón”, Bioagro, vol. 20, n.º 3, pp. 193-199, 2008</mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Zamora</surname>
							<given-names>F</given-names>
						</name>
						<name>
							<surname>Rodríguez</surname>
							<given-names>N.</given-names>
						</name>
						<name>
							<surname>Torres</surname>
							<given-names>D.</given-names>
						</name>
						<name>
							<surname>Yendis</surname>
							<given-names>H.</given-names>
						</name>
                    </person-group>
					<year>2008</year>
					<article-title>Efecto del riego con aguas residuales sobre propiedades químicas de suelos de la planicie de Coro, Estado Falcón</article-title>
					<source>Bioagro</source>
					<volume>20</volume>
					<issue>3</issue>
					<fpage>193</fpage>
					<lpage>199</lpage>
				</element-citation>
            </ref>
            <ref id="r4">
                <mixed-citation>C. Fuentes, B. de León, J. Y. Parlange y A. C. Antonino, “El sistema de ecuaciones de Saint-Venant y Richards del riego por gravedad: 1. La ley potencial de resistencia hidráulica”, Tecnol. Cienc. Agua, vol. 19, n.º 2, pp.65-75, 2015</mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Fuentes</surname>
							<given-names>C</given-names>
						</name>
						<name>
							<surname>de León</surname>
							<given-names>B</given-names>
						</name>
						<name>
							<surname>Parlange</surname>
							<given-names>J. Y.</given-names>
						</name>
						<name>
							<surname>Antonino</surname>
							<given-names>A. C.</given-names>
						</name>
                    </person-group>
					<year>2015</year>
					<article-title>El sistema de ecuaciones de Saint-Venant y Richards del riego por gravedad: 1. La ley potencial de resistencia hidráulica</article-title>
					<source>Tecnol. Cienc. Agua</source>
					<volume>19</volume>
					<issue>2</issue>
                    <fpage>65</fpage>
                    <lpage>75</lpage>
				</element-citation>
			</ref>
            <ref id="r5">
                <mixed-citation>O. Seguel, I. Fuentes, J. Dörner y J. Cuevas, J. “Propiedades físicas e hidráulicas de suelo y su relación con el movimiento de contaminantes”, Agro Sur, vol. 43, n.º 1, pp. 41-50, 2015. <ext-link xlink:href="https://doi.org/10.4206/agrosur.2015.v43n1-06">https://doi.org/10.4206/agrosur.2015.v43n1-06</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Seguel</surname>
							<given-names>O</given-names>
						</name>
						<name>
							<surname>Fuentes</surname>
							<given-names>I</given-names>
						</name>
						<name>
							<surname>Dörner</surname>
							<given-names>J</given-names>
						</name>
						<name>
							<surname>Cuevas</surname>
							<given-names>J</given-names>
						</name>
                    </person-group>
					<year>2015</year>
					<article-title>Propiedades físicas e hidráulicas de suelo y su relación con el movimiento de contaminantes</article-title>
					<source>Agro Sur</source>
					<volume>43</volume>
					<issue>1</issue>
                    <fpage>41</fpage>
                    <lpage>50</lpage>
					<pub-id pub-id-type="doi">10.4206/agrosur.2015.v43n1-06</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.4206/agrosur.2015.v43n1-06</ext-link>
				</element-citation>
            </ref>
            <ref id="r6">
                <mixed-citation>W. R. Gardner, “Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table”, Soil Sci., vol. 85, n.º 4, pp. 228-232, 1958</mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Gardner</surname>
							<given-names>W. R.</given-names>
						</name>
                    </person-group>
					<year>1958</year>
					<article-title>Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table</article-title>
					<source>Soil Sci</source>
					<volume>85</volume>
					<issue>4</issue>
                    <fpage>228</fpage>
                    <lpage>232</lpage>
				</element-citation>
            </ref>
            <ref id="r7">
               <mixed-citation>D. J. Timlin, L. R. Ahuja, Y. Pachepsky, R. D. Williams, D. Gimenez y W. Rawls, “Use of Brooks-Corey parameters to improve estimates of saturated conductivity from effective porosity”, Soil Sci., vol. 63, n.º 5, pp.1086-1092, 1999. <ext-link xlink:href="https://doi.org/10.2136/sssaj1999.6351086x">https://doi.org/10.2136/sssaj1999.6351086x</ext-link></mixed-citation>
			   <element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Timlin</surname>
							<given-names>D. J.</given-names>
						</name>
						<name>
							<surname>Ahuja</surname>
							<given-names>L. R.</given-names>
						</name>
						<name>
							<surname>Pachepsky</surname>
							<given-names>Y</given-names>
						</name>
						<name>
							<surname>Williams</surname>
							<given-names>R. D.</given-names>
						</name>
						<name>
							<surname>Gimenez</surname>
							<given-names>D.</given-names>
						</name>
						<name>
							<surname>Rawls</surname>
							<given-names>W.</given-names>
						</name>
                    </person-group>
					<year>1999</year>
					<article-title>Use of Brooks-Corey parameters to improve estimates of saturated conductivity from effective porosity</article-title>
					<source>Soil Sci.</source>
					<volume>63</volume>
					<issue>5</issue>
                    <fpage>1086</fpage>
                    <lpage>1092</lpage>
					<pub-id pub-id-type="doi">10.2136/sssaj1999.6351086x</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.2136/sssaj1999.6351086x</ext-link>
				</element-citation>
            </ref>
            <ref id="r8">
                <mixed-citation>A. Y. Hoekstra, “Water footprint assessment: evolvement of a new research field”, Water Resour. Manag, vol. 31, n.º 10, pp. 3061-3081, 2017. <ext-link xlink:href="https://doi.org/10.1007/s11269-017-1618-5">https://doi.org/10.1007/s11269-017-1618-5</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Hoekstra</surname>
							<given-names>A. Y.</given-names>
						</name>
                    </person-group>
					<year>2017</year>
					<article-title>Water footprint assessment: evolvement of a new research field</article-title>
					<source>Water Resour. Manag.</source>
					<volume>31</volume>
					<issue>10</issue>
                    <fpage>3061</fpage>
                    <lpage>3081</lpage>
					<pub-id pub-id-type="doi">10.1007/s11269-017-1618-5</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.1007/s11269-017-1618-5</ext-link>
				</element-citation>
            </ref>  
            <ref id="r9">
                 <mixed-citation>J. A. Blanco, “Bosques, suelo y agua: explorando sus interacciones”, Rev. Ecosistemas, vol. 26, n.º 2, pp. 1-9, 2017. <ext-link xlink:href="https://doi.org/10.7818/ECOS.2017.26-2.01">https://doi.org/10.7818/ECOS.2017.26-2.01</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Blanco</surname>
							<given-names>J. A.</given-names>
						</name>
                    </person-group>
					<year>2017</year>
					<article-title>Bosques, suelo y agua: explorando sus interacciones</article-title>
					<source>Rev. Ecosistemas</source>
					<volume>26</volume>
					<issue>2</issue>
                    <fpage>1</fpage>
                    <lpage>9</lpage>
					<pub-id pub-id-type="doi">10.7818/ECOS.2017.26-2.01</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.7818/ECOS.2017.26-2.01</ext-link>
				</element-citation>
			</ref>  
            <ref id="r10">
                <mixed-citation>M. Durán-Pinzón, J. Páez-Arango y P. J. García-Guarín, “Numerical modeling and experimental analysis for flows in a homogeneous porous medium through soils”, Iteckne, vol. 15, n.º 1, pp. 24-33, 2018. <ext-link xlink:href="http://dx.doi.org/10.15332/iteckne.v15i1.1961">http://dx.doi.org/10.15332/iteckne.v15i1.1961</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Durán-Pinzón</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Páez-Arango</surname>
							<given-names>J.</given-names>
						</name>
						<name>
							<surname>García-Guarín</surname>
							<given-names>P. J.</given-names>
						</name>
                    </person-group>
					<year>2018</year>
					<article-title>Numerical modeling and experimental analysis for flows in a homogeneous porous medium through soils</article-title>
					<source>Iteckne</source>
					<volume>15</volume>
					<issue>1</issue>
                    <fpage>24</fpage>
                    <lpage>33</lpage>
					<pub-id pub-id-type="doi">10.15332/iteckne.v15i1.1961</pub-id>
					<ext-link ext-link-type="uri">http://dx.doi.org/10.15332/iteckne.v15i1.1961</ext-link>
				</element-citation>
			</ref>  
			<ref id="r11">
				<mixed-citation>R. Angulo-Jaramillo, V. Bagarello, M. Iovino y L. Lassabatere, “Saturated soil hydraulic conductivity”, en Infiltration Measurements for Soil Hydraulic Characterization. Springer, 2016, pp. 43-180. <ext-link xlink:href="https://doi.org/10.1007/978-3-319-31788-5_2">https://doi.org/10.1007/978-3-319-31788-5_2</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>Angulo-Jaramillo</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Bagarello</surname>
							<given-names>V.</given-names>
						</name>
						<name>
							<surname>Iovino</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Lassabatere</surname>
							<given-names>L.</given-names>
						</name>
                    </person-group>
					<year>2018</year>
					<article-title>“Saturated soil hydraulic conductivity”</article-title>
					<source>Infiltration Measurements for Soil Hydraulic Characterization. Springer</source>
                    <fpage>43</fpage>
                    <lpage>180</lpage>
					<pub-id pub-id-type="doi">10.1007/978-3-319-31788-5_2</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.1007/978-3-319-31788-5_2</ext-link>
				</element-citation>
            </ref>
			<ref id="r12">
				<mixed-citation>M. Zavala, H. Saucedo, C. Fuentes y C. Bautista. “Modelo de conductividad hidráulica dual para el movimiento del agua en suelos Macroporosos”, Agrociencia, vol. 46, n.º 3, pp. 205-220, 2012 [en línea]. Disponible: <ext-link xlink:href="http://www.scielo.org.mx/scielo.php?script=sci_arttext&amp;pid=S1405-31952012000300001">http://www.scielo.org.mx/scielo.php?script=sci_arttext&amp;pid=S1405-31952012000300001</ext-link></mixed-citation>
				<element-citation publication-type="webpage">
                    <person-group person-group-type="author">
						<name>
							<surname>Zavala</surname>
							<given-names>M.</given-names>
						</name>
						<name>
							<surname>Saucedo</surname>
							<given-names>H.</given-names>
						</name>
						<name>
							<surname>Fuentes</surname>
							<given-names>C.</given-names>
						</name>
						<name>
							<surname>Bautista</surname>
							<given-names>C.</given-names>
						</name>
                    </person-group>
					<year>2012</year>
					<article-title>Modelo de conductividad hidráulica dual para el movimiento del agua en suelos Macroporosos</article-title>
					<source>Agrociencia</source>
					<volume>46</volume>
					<issue>3</issue>
                    <fpage>205</fpage>
                    <lpage>220</lpage>
					<ext-link ext-link-type="uri">http://www.scielo.org.mx/scielo.php?script=sci_arttext&amp;pid=S1405-31952012000300001</ext-link>
				</element-citation>
            </ref>
			<ref id="r13">
				<mixed-citation>A. A. Torrealba Almarza, “Aplicación de un modelo numérico para estimar curvas de evaporación a partir de información de granulometría básica”, Tesis de Magíster en Ciencias de la Ingeniería, Pontificia Universidad Católica de Chile, 2008</mixed-citation>
				<element-citation publication-type="book">
                    <person-group person-group-type="author">
						<name>
							<surname>Torrealba Almarza</surname>
							<given-names>A. A.</given-names>
						</name>
                    </person-group>
					<year>2008</year>
					<source>Aplicación de un modelo numérico para estimar curvas de evaporación a partir de información de granulometría básica</source>
					<publisher-name>Pontificia Universidad Católica de Chile</publisher-name>
				</element-citation>
            </ref>
			<ref id="r14">
				<mixed-citation>E. F. García-Aristizábal, J. C. Guzmán-Martínez, A. N. Gallego-Hernández y C. A. Vega-Posada, “Análisis numérico unidimensional del proceso de infiltración en suelos parcialmente saturados reforzados con geosintéticos permeables”, Scientia et Technica, vol. 24, n.º 1, pp. 46-56, 2019. <ext-link xlink:href="http://dx.doi.org/10.22517/23447214.18611">http://dx.doi.org/10.22517/23447214.18611</ext-link></mixed-citation>
				<element-citation publication-type="journal">
                    <person-group person-group-type="author">
						<name>
							<surname>García-Aristizábal</surname>
							<given-names>E. F.</given-names>
						</name>
						<name>
							<surname>Guzmán-Martínez</surname>
							<given-names>J. C.</given-names>
						</name>
						<name>
							<surname>Gallego-Hernández</surname>
							<given-names>A. N.</given-names>
						</name>
						<name>
							<surname>Vega-Posada</surname>
							<given-names>C. A.</given-names>
						</name>
                    </person-group>
					<year>2019</year>
					<article-title>Análisis numérico unidimensional del proceso de infiltración en suelos parcialmente saturados reforzados con geosintéticos permeables</article-title>
					<source>Scientia et Technica</source>
					<volume>24</volume>
					<issue>1</issue>
                    <fpage>46</fpage>
                    <lpage>56</lpage>
					<pub-id pub-id-type="doi">10.22517/23447214.18611</pub-id>
					<ext-link ext-link-type="uri">http://dx.doi.org/10.22517/23447214.18611</ext-link>
				</element-citation>
            </ref>
			<ref id="r15">
				<mixed-citation>L. A. Lozano, G. Soracco, R. Villarreal, R. Filgueira y G. Oliverio, “Difusividad del agua en el suelo: un método simple de laboratorio para su determinación”, Cienc. Suelo, vol. 35, n.º 1, pp. 11-19, 2017 [en línea]. Disponible: <ext-link xlink:href="http://hdl.handle.net/11336/56701">http://hdl.handle.net/11336/56701</ext-link></mixed-citation>
				<element-citation publication-type="webpage">
                    <person-group person-group-type="author">
						<name>
							<surname>Lozano</surname>
							<given-names>L. A.</given-names>
						</name>
						<name>
							<surname>Soracco</surname>
							<given-names>G.</given-names>
						</name>
						<name>
							<surname>Villarreal</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Filgueira</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Oliverio</surname>
							<given-names>G.</given-names>
						</name>
                    </person-group>
					<year>2017</year>
					<article-title>Difusividad del agua en el suelo: un método simple de laboratorio para su determinación</article-title>
					<source>Cienc. Suelo</source>
					<volume>35</volume>
					<issue>1</issue>
                    <fpage>11</fpage>
                    <lpage>19</lpage>
					<ext-link ext-link-type="uri">http://hdl.handle.net/11336/56701</ext-link>
				</element-citation>
			</ref>
			<ref id="r16">
				<mixed-citation>M. Th. van Genuchten, “A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1”, Soil Sci., vol. 44, n.º 5, pp. 892-898, 1980</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>van Genuchten</surname>
							<given-names>M. Th.</given-names>
						</name>
                    </person-group>
					<year>1980</year>
                    <article-title>A closed-form equation for predicting the hydraulic conductivity of unsaturated soils 1</article-title>
					<source>Soil Sci</source>
                    <volume>44</volume>
                    <issue>5</issue>
                    <fpage>892</fpage>
                    <lpage>898</lpage>
                </element-citation>
			</ref>
			<ref id="r17">
				<mixed-citation>C. van Duijn y G. Pieters, Steady Flows in Unsatured soils are Stable. Vol. 16. Eindhoven: Eindhoven University of Technology, 2003</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>van Duijn</surname>
							<given-names>C.</given-names>
						</name>
						<name>
							<surname>Pieters</surname>
							<given-names>G.</given-names>
						</name>
                    </person-group>
					<year>2003</year>
					<source>Steady Flows in Unsatured soils are Stable</source>
					<publisher-name>Eindhoven: Eindhoven University of Technology</publisher-name>
                    <volume>44</volume>
                </element-citation>
            </ref>
			<ref id="r18">
				<mixed-citation>V. Arrarás. “Hydrofem: un modelo tridimensional para la simulación del flujo de agua en medios porosos de saturación variable. Estudio de la zona no saturada del suelo”, en Estudios de la zona no saturada de suelo, vol. VI. Instituto Tecnológico Agrario de Castilla y León, 2003, pp. 271-278 [en línea]. Disponible: <ext-link xlink:href="https://abe.ufl.edu/faculty/carpena/files/pdf/zona_no_saturada/estudios_de_la_zona_v6/p271-278.pdf">https://abe.ufl.edu/faculty/carpena/files/pdf/zona_no_saturada/estudios_de_la_zona_v6/p271-278.pdf</ext-link></mixed-citation>
         	    <element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Arrarás</surname>
							<given-names>V.</given-names>
						</name>
                	</person-group>
					<year>2003</year>
					<source>Hydrofem: un modelo tridimensional para la simulación del flujo de agua en medios porosos de saturación variable. Estudio de la zona no saturada del suelo</source>
					<publisher-name>Instituto Tecnológico Agrario de Castilla y León</publisher-name>
                    <fpage>271</fpage>
                    <lpage>278</lpage>
					<ext-link ext-link-type="uri">https://abe.ufl.edu/faculty/carpena/files/pdf/zona_no_saturada/estudios_de_la_zona_v6/p271-278.pdf</ext-link>
                </element-citation>
			</ref>
			<ref id="r19">
				<mixed-citation>J. Samper. Apuntes de hidrología del subsuelo. La Coruña: ETS Ingenieros de Caminos, Canales y Puertos-Universidad de La Coruña. Curso 2013-2014</mixed-citation>
				<element-citation publication-type="book">
					<person-group person-group-type="author">
						<name>
							<surname>Samper</surname>
							<given-names>J.</given-names>
						</name>
                	</person-group>
					<year>2014</year>
					<source>Apuntes de hidrología del subsuelo. La Coruña: ETS Ingenieros de Caminos, Canales y Puertos-Universidad de La Coruña</source>
                </element-citation>
            </ref>
			<ref id="r20">
				<mixed-citation>M. van Genuchten, J. Simunet, F. Leij, y M. Sejna, “RETC version 6: Code for quantifying the hydraulic functions of unsaturated soils”, US Salinity Laboratory, USDA. ARS, USA, 2000</mixed-citation>
         	       	<element-citation publication-type="book">
                    <person-group person-group-type="author">
						<name>
							<surname>van Genuchten</surname>
							<given-names>M.</given-names>
						</name>
                        <name>
							<surname>Simunet</surname>
							<given-names>J.</given-names>
						</name>
						<name>
							<surname>Leij</surname>
							<given-names>F.</given-names>
						</name>
						<name>
							<surname>Sejna</surname>
							<given-names>M.</given-names>
						</name>
                	</person-group>
					<year>2000</year>
					<source>RETC version 6: Code for quantifying the hydraulic functions of unsaturated soils</source>
					<publisher-loc>USA</publisher-loc>
                    <publisher-name>US Salinity Laboratory</publisher-name>
				</element-citation>
			</ref>
			<ref id="r21">
				<mixed-citation>G. López-Canteñs, J. A. Ostos-Santos, A. Rojano-Aguilar, A. Vázquez-Alarcón y E. Arteaga-Tovar, “Aproximación de la ecuación de difusión en el suelo con serie de Fourier”, Terra Latinoamericana, vol. 24, n.º 4, pp. 443-450, 2006</mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>López-Canteñs</surname>
							<given-names>G.</given-names>
						</name>
						<name>
							<surname>Ostos-Santos</surname>
							<given-names>J. A.</given-names>
						</name>
						<name>
							<surname>Rojano-Aguilar</surname>
							<given-names>A.</given-names>
						</name>
						<name>
							<surname>Vázquez-Alarcón</surname>
							<given-names>A.</given-names>
						</name>
						<name>
							<surname>Arteaga-Tovar</surname>
							<given-names>E.</given-names>
						</name>
                    </person-group>
					<year>2006</year>
                    <article-title>Aproximación de la ecuación de difusión en el suelo con serie de Fourier</article-title>
					<source>Terra Latinoamericana</source>
                    <volume>24</volume>
                    <issue>4</issue>
                    <fpage>443</fpage>
                    <lpage>450</lpage>
                </element-citation>
            </ref>
			<ref id="r22">
				<mixed-citation>E. O. L. Robles, I. C. Silva, M. I. Y. Díaz, H. G. Rodríguez, J. M. Monsiváis y S. J. B. Pulido, “Ajuste de modelos empíricos de infiltración en un Umbrisol bajo diferentes tratamientos silvícolas”, Rev. Mex. Cienc. Forestales, vol. 11, n.º 57, pp. 132-152, 2020. <ext-link xlink:href="https://doi.org/10.29298/rmcf.v11i57.643">https://doi.org/10.29298/rmcf.v11i57.643</ext-link></mixed-citation>
         	        <element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Robles</surname>
							<given-names>E. O. L.</given-names>
						</name>
                        <name>
							<surname>Silva</surname>
							<given-names>I. C.</given-names>
						</name>
                        <name>
							<surname>Díaz</surname>
							<given-names>M. I. Y.</given-names>
						</name>
						<name>
							<surname>Rodríguez</surname>
							<given-names>H. G.</given-names>
						</name>
						<name>
							<surname>Monsiváis</surname>
							<given-names>J. M.</given-names>
						</name>
						<name>
							<surname>Pulido</surname>
							<given-names>S. J. B.</given-names>
						</name>
                	</person-group>
					<year>2020</year>
                    <article-title>Ajuste de modelos empíricos de infiltración en un Umbrisol bajo diferentes tratamientos silvícolas</article-title>
					<source>Rev. Mex. Cienc. Forestales</source>
                    <volume>11</volume>
					<issue>57</issue>
                    <fpage>132</fpage>
					<lpage>152</lpage>
                    <pub-id pub-id-type="doi">10.29298/rmcf.v11i57.643</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.29298/rmcf.v11i57.643</ext-link>
					</element-citation>
				</ref>
			<ref id="r23">
				<mixed-citation>M. Servin Palestina, R. A. Sánchez Gutiérrez, O. Ramírez Valle, M. A. Galindo Reyes y H. Gutiérrez Bañuelos, “Modelos para programación y optimización de agua de riego en avena forrajera”, Rev. Mex. Cienc. Pecu., vol. 9, n.º 4, pp. 667-684, 2018. <ext-link xlink:href="http://dx.doi.org/10.22319/rmcp.v9i4.4404">http://dx.doi.org/10.22319/rmcp.v9i4.4404</ext-link></mixed-citation>
         	    <element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Servin Palestina</surname>
							<given-names>M.</given-names>
						</name>
                        <name>
							<surname>Sánchez Gutiérrez</surname>
							<given-names>R. A.</given-names>
						</name>
                        <name>
							<surname>Ramírez Valle</surname>
							<given-names>O.</given-names>
						</name>
						<name>
							<surname>Galindo Reyes</surname>
							<given-names>M. A.</given-names>
						</name>
						<name>
							<surname>Gutiérrez Bañuelos</surname>
							<given-names>H.</given-names>
						</name>
                	</person-group>
					<year>2018</year>
                    <article-title>Modelos para programación y optimización de agua de riego en avena forrajera</article-title>
					<source>Rev. Mex. Cienc. Pecu.</source>
                    <volume>9</volume>
					<issue>4</issue>
                    <fpage>667</fpage>
					<lpage>684</lpage>
                    <pub-id pub-id-type="doi">10.22319/rmcp.v9i4.4404</pub-id>
					<ext-link ext-link-type="uri">http://dx.doi.org/10.22319/rmcp.v9i4.4404</ext-link>
				</element-citation>
			</ref>
			<ref id="r24">
				<mixed-citation>J. L. González-Barrios et al., “Caracterización de la porosidad edáfica como indicador de la calidad física del suelo”, Terra Latinoamericana, vol. 29, n.º 4, pp. 369-377, 2011. <ext-link xlink:href="https://doi.org/10.13140/2.1.1262.4002">https://doi.org/10.13140/2.1.1262.4002</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>González-Barrios</surname>
							<given-names>J. L.</given-names>
						</name>
                	</person-group>
					<year>2011</year>
                    <article-title>Caracterización de la porosidad edáfica como indicador de la calidad física del suelo</article-title>
					<source>Terra Latinoamericana</source>
                    <volume>29</volume>
					<issue>4</issue>
                    <fpage>369</fpage>
					<lpage>377</lpage>
                    <pub-id pub-id-type="doi">10.13140/2.1.1262.4002</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.13140/2.1.1262.4002</ext-link>
				</element-citation>
            </ref>
			<ref id="r25">
				<mixed-citation>M. R. Campos, G. C. Lazo, F. G. Robaina y T. L Seijas, “Curvas tensión humedad y función conductividad hidráulica en suelos pardos”, Revista Ingeniería Agrícola, vol. 10, n.º 2, 2020 [en línea]. Disponible: <ext-link xlink:href="https://revistas.unah.edu.cu/index.php/IAgric/article/view/1241">https://revistas.unah.edu.cu/index.php/IAgric/article/view/1241</ext-link></mixed-citation>
				<element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Campos</surname>
							<given-names>M. R.</given-names>
						</name>
						<name>
							<surname>Lazo</surname>
							<given-names>G. C.</given-names>
						</name>
						<name>
							<surname>Robaina</surname>
							<given-names>F. C.</given-names>
						</name>
						<name>
							<surname>Seijas</surname>
							<given-names>T. L.</given-names>
						</name>
                	</person-group>
					<year>2020</year>
                    <article-title>Curvas tensión humedad y función conductividad hidráulica en suelos pardos</article-title>
					<source>Revista Ingeniería Agrícola</source>
                    <volume>10</volume>
					<issue>2</issue>
					<ext-link ext-link-type="uri">https://revistas.unah.edu.cu/index.php/IAgric/article/view/1241</ext-link>
				</element-citation>
            </ref>
			<ref id="r26">
				<mixed-citation>H. Zhu, T. Liu, B. Xue y G. Wang, “Modified Richards’ Equation to improve estimates of soil moisture in two-layered soils after infiltration”, Water, vol. 10, n.º 9, pp. 1174, 2018. <ext-link xlink:href="https://doi.org/10.3390/w10091174">https://doi.org/10.3390/w10091174</ext-link></mixed-citation>
                <element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Zhu</surname>
							<given-names>H.</given-names>
						</name>
						<name>
							<surname>Liu</surname>
							<given-names>T.</given-names>
						</name>
						<name>
							<surname>Xue</surname>
							<given-names>B.</given-names>
						</name>
						<name>
							<surname>Wang</surname>
							<given-names>G.</given-names>
						</name>
                	</person-group>
					<year>2018</year>
                    <article-title>Modified Richards’ Equation to improve estimates of soil moisture in two-layered soils after infiltration</article-title>
					<source>Water</source>
                    <volume>10</volume>
					<issue>9</issue>
					<fpage>1174</fpage>
					<lpage>1174</lpage>
					<pub-id pub-id-type="doi">10.3390/w10091174</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.3390/w10091174</ext-link>
				</element-citation>
			</ref>
			<ref id="r27">
				<mixed-citation>R. Villarreal, L. A. Lozano, E. M. Melani, M. P. Salazar, M. F. Otero Y C. G. Soracco, “Diffusivity and sorptivity determination at different soil water contents from horizontal infiltration”, Geoderma, vol. 338, pp. 88-96, 2019. <ext-link xlink:href="https://doi.org/10.1016/j.geoderma.2018.11.045">https://doi.org/10.1016/j.geoderma.2018.11.045</ext-link></mixed-citation>
         	    <element-citation publication-type="journal">
					<person-group person-group-type="author">
						<name>
							<surname>Villarreal</surname>
							<given-names>R.</given-names>
						</name>
						<name>
							<surname>Lozano</surname>
							<given-names>L. A.</given-names>
						</name>
						<name>
							<surname>Melani</surname>
							<given-names>E. M.</given-names>
						</name>
						<name>
							<surname>Salazar</surname>
							<given-names>M. P.</given-names>
						</name>
						<name>
							<surname>Otero</surname>
							<given-names>M. F.</given-names>
						</name>
						<name>
							<surname>Soracco</surname>
							<given-names>C. G.</given-names>
						</name>
                	</person-group>
					<year>2019</year>
                    <article-title>Diffusivity and sorptivity determination at different soil water contents from horizontal infiltration</article-title>
					<source>Geoderma</source>
                    <volume>338</volume>
					<fpage>88</fpage>
					<lpage>96</lpage>
					<pub-id pub-id-type="doi">10.1016/j.geoderma.2018.11.045</pub-id>
					<ext-link ext-link-type="uri">https://doi.org/10.1016/j.geoderma.2018.11.045</ext-link>
				</element-citation>
			</ref>
		</ref-list>
	</back>
</article>