Uso actual del suelo de los pastizales templados de Argentina : impacto del manejo sobre la provisión de servicios ecosistémicos intermedios relacionados con la dinámica del carbono y el nitrógeno

Las transformaciones realizadas por el ser humano alteran la composición y estructura del ecosistema para dirigir la energía fijada a la obtención de bienes con valor de mercado. Sin embargo, dichas alteraciones en el ambiente también reducen la provisión de servicios ecosistémicos. En esta tesis se analiza el impacto de las transformaciones agrícolas sobre la provisión de servicios ecosistémicos a nivel regional (Región Pampeana), en particular el carbono orgánico del suelo (COS) y las emisiones de N2O. Para ello se evalúa, en la principal zona agrícola del país, la dinámica temporal de las rotaciones agrícola-ganaderas y la importancia relativa de distintos tipos de labranza. La aproximación metodológica incluyó el uso de técnicas de teledetección, ensayos a campo y modelos de simulación. Los resultados obtenidos señalan que en la Pampa Ondulada, la agricultura continua fue el uso dominante y se encontró que un 28 por ciento de la superficie bajo este uso es monocultivo de soja. En cuanto a los sistemas de labranza la siembra directa representa el 70 por ciento del área bajo cultivos de verano. Mediante un modelo de simulación se determinó que de las secuencias de cultivos relevadas en la región la soja1ra/maíz bajo labranza convencional y sin fertilización fue la de mayor pérdida de COS para un período de 60 años (37 por ciento). Por el contrario, rotaciones de soja/trigo-soja2da (6 años) pastura (4 años) bajo siembra directa y altos niveles de fertilización presentaron un aumento del 10 por ciento en el mismo período. La información espacial muestra menores pérdidas de COS hacia el este y mayores hacia el oeste de la zona de estudio. En la pampa interior se comparó la emisión de N2O entre un cultivo de soja, solo y en doble cultivo, con el pastizal. Los resultados no mostraron diferencias significativas entre las diferentes coberturas

Transferencia de nitrógeno de ecosistemas agrícolas recientemente intensificados a acuíferos: efecto de la variación en las precipitaciones

A partir del monitoreo a campo de cultivos, suelo y agua subterránea, modelos de simulación y un experimento en lisímetros se estudió el intercambio de nitrógeno (N) inorgánico entre cultivos y el agua subterránea en sistemas agrícolas recientemente intensificados con suelos ricos en materia orgánica (MO). El monitoreo (1998-2010) en tres posiciones del paisaje de la Pampa Ondulada permitió captar períodos plurianuales húmedos y secos (1506-444 mm año-1) disparadores de cambios intensos en las profundidades freáticas (mayor 1 a menor a 6 m) que afectaron la función del agua subterránea y de los cultivos como moduladores de los flujos de agua y N. Con niveles freáticos superficiales, se redujo la capacidad de almacenaje del agua drenada hacia la zona no saturada, produciéndose un flujo lateral subsuperficial que redistribuyó el N lixiviado desde posiciones altas hacia los bajos y contribuyó a la recarga del acuífero aguas abajo. Este flujo, favorecido por fuertes gradientes hidráulicos, y la concentración local de solutos disparada por el consumo de agua subterránea por los cultivos en la posición de pie de loma, produjeron picos de concentración de cloruro (menor a 500 mg Cl. l-1) y nitrato (menor a 45 mg N-NO3. l-1) en el agua freática superficial. Durante períodos normales a secos los drenajes de diferentes eventos de lluvia se alojan en la zona no saturada, disipando y retrasando la respuesta del nivel freático a eventos individuales. El agua de lluvia fluye verticalmente arrastrando N del suelo hacia una superficie freática profunda (con mayor desfasaje temporal entre la concentración de nitrato del drenaje y la del agua freática a medida que el nivel se profundiza) y luego descarga en el arroyo. La aplicación de un fertilizante enriquecido con 15N a un cultivo de maíz en lisímetros permitió demostrar que la fertilización representa un aporte insignificante (1 por ciento) al flujo de N lixiviado y que el N derivado de la mineralización de la MO del suelo constituye una fuente importante de N lixiviable durante períodos con balance hídrico positivo y/o baja demanda de los cultivos. A partir de los resultados obtenidos se propone un modelo conceptual de los flujos hidrológicos y de N en paisajes ondulados de uso agrícola para situaciones climáticas contrastantes

Worksheet for astronomical variables (VARAST 1.0)

p.171-177

Worksheet for astronomical variables (VARAST 1.0)

p.171-177

Efectos de la temperatura sobre el peso seco, concentración de aceite y la proporción de ácidos grasos en el aceite de frutos de olivo (Olea europaea L. var Arauco)

Las respuestas de los componentes del rendimiento y de la calidad del aceite a la temperatura aún no han sido estudiadas en olivo ni en otras oleaginosas que acumulan aceite principalmente en mesocarpo (cerca de 95 por ciento del total en el caso de olivo). La información disponible en olivo se basa en correlaciones ligadas a variaciones en ubicación geográfica, altitud y años en los que pudieron covariar con la temperatura otros factores (radiación, disponibilidad hídrica, nutrición, etc). El objetivo de esta tesis fue evaluar el rol de la temperatura durante el crecimiento de frutos de olivo sobre el peso seco de los mismos, su contenido y concentración de aceite y las proporciones de ácidos grasos del aceite. Se realizaron experimentos manipulativos aplicados a ramas fructíferas creciendo a campo, utilizando cámaras transparentes con control y registro de temperatura. Los rangos de temperaturas medias alcanzados (entre 15-32 ºC) fueron mucho mayores a los rangos de 1 y 4 ºC logrados en estudios correlativos y permitieron generar relaciones funcionales para el peso seco del fruto, su concentración de aceite, y la proporción de ácidos grasos en el aceite. El contraste entre aplicaciones de temperaturas mayores o menores al ambiente durante períodos largos (114 días) y cortos (30 días) demostraron que tanto la proporción de aceite del fruto como la de ácido oleico en el aceite de fruto entero (mesocarpo + semilla) disminuye al aumentar la temperatura. Si bien en ambas estructuras del fruto la temperatura tuvo un efecto negativo sobre la concentración de aceite, las respuestas a la temperatura de la proporción de ácidos grasos de la semilla y el mesocarpo, medidos en forma separada, difirieron entre sí. Los patrones de respuesta de las proporciones de ácidos grasos en el aceite de semilla mostraron similitudes con las conocidas para especies oleaginosas de semilla como girasol. Por contraste, en el aceite de mesocarpo de olivo la proporción de ácido oleico disminuyó con el aumento de la temperatura, patrón opuesto a lo manifestado en semillas de olivo y en las de oleaginosas anuales. En esta tesis también se evaluó el efecto de la temperatura mínima nocturna sobre el crecimiento de los frutos y la calidad del aceite. Los resultados sugieren que la temperatura mínima nocturna y la amplitud térmica diarias son las dimensiones del régimen térmico diario con los que mejor se asociaron los cambios en el peso seco del fruto, la proporción de aceite y las proporciones de ácidos grasos en el mismo. Los resultados de esta tesis sirven para guiar la selección de sitios de plantación de nuevos olivares basados en los registros térmicos zonales, y para explicar los desajustes de la proporción de ácidos grasos del aceite de oliva producido en el NOA y la normativa del Consejo Oleícola Internacional (COI). Las relaciones funcionales entre las variables respuesta y la temperatura definidas podrán incorporarse en modelos de simulación del rendimiento y calidad que se desarrollen en el futuro. Tomados en conjunto con la tendencia al aumento de la temperatura, producto del calentamiento global, los resultados permitirán estimar la magnitud de los impactos negativos que podrían afectar el rendimiento del cultivo y la calidad del aceite.

Inferencia bayesiana sobre los parámetros de dispersión genéticos y ambientales en modelos animales con efectos maternos

Los modelos 'modelos animales con efectos maternos' (MAM) son modelos lineales mixtos que se utilizan para ajustar registros de caracteres bajo la influencia de efectos maternos. Uno de los desafíos más importantes en el marco de los MAM es la estimación de los parámetros de dispersión o 'componentes de (co) varianza' (CVC). En esta tesis se introducen desde una perspectiva bayesiana contribuciones teóricas y metodológicas con relación a la estimación de CVC para MAM sujetos a estructuras de covarianza novedosas. En primer lugar, se describe una implementación del análisis bayesiano jerárquico vía el algoritmo del muestreo de Gibbs. Luego, se considera una especificación conjugada diferente para la distribución a priori de la matriz de covarianza genética, basada en la distribución Wishart invertida generalizada, y se presenta una estrategia para determinar los correspondientes hiperparámetros. Esta estrategia fue comparada contra otras especificaciones a priori mediante un estudio de simulación estocástica, y produjo estimaciones precisas de los parámetros genéticos, con menores errores estándares y mejor tasa de convergencia. En segundo lugar, se presenta una formulación alternativa del MAM que incluye un parámetro de correlación ambiental entre pares de observaciones madre-progenie, y se desarrolla un procedimiento de estimación basado en un algoritmo de muestreo por grilla. El procedimiento fue programado y ejecutado exitosamente, y se obtuvo la primera estimación del parámetro de correlación con datos de campo para peso al destete en bovinos de carne. Por último, se considera el problema de la estimación de CVC en una población multirracial, donde en general es necesario especificar una estructura de covarianza heterogénea para los valores de cría. En particular, se demuestra que el modelo basado en la descomposición de la matriz de covarianza genética es equivalente al que deriva de la teoría genética cuantitativa. Además, se extiende el modelo para incluir efectos maternos y se describe la implementación de un análisis bayesiano jerárquico con el objetivo de estimar los CVC. El procedimiento fue implementado con éxito en datos experimentales de peso al destete y se obtuvieron por primera vez estimaciones para el conjunto completo de CVC.

Efecto de variables ambientales sobre la productividad primaria neta aérea y la concentración de proteína bruta de gramíneas en sistemas silvopastoriles de ñire (Nothofagus antarctica) : creación de un modelo de simulación

Nothofagus antarctica (ñire) ocupa una importante superficie de los bosques de Patagonia Sur, siendo ampliamente usado como sistemas silvopastoriles con ganadería ovina y bovina. Sin embargo, en la mayoría de los casos el uso de estos bosques se realiza sin un manejo que permita optimizar la sustentabilidad del sistema. En esta tesis se estudió el efecto de distintas variables ambientales (temperatura, humedad del suelo, radiación) sobre la productividad y concentración de Proteína bruta (PB) de gramíneas en distintas situaciones de calidad forestal y coberturas de copas de bosques de ñire, así como también procesos asociados: descomposición y dinámica de nutrientes (mineralización-inmovilización) de residuos del estrato arbóreo y graminoso y mineralización de nitrógeno del suelo. La productividad de materia seca (MS) aérea de gramíneas (rango entre 135 a 1863 kg MS ha-1) estuvo positivamente asociada a la calidad forestal de los sitios, mientras que el efecto de las coberturas de copa dependió de la interacción con otras variables ambientales. Los valores de PB variaron principalmente según la cobertura de copa encontrándose los mayores valores en lugares más sombreados (transmisividades menor al 70 por ciento). La descomposición, tanto en hojas de ñire como en gramíneas, varió según el nivel de radiación, siendo mayores en los lugares con menor cobertura. Por otro lado, la dinámica de nutrientes no fue afectada por la CS o alguna variable estudiada. La mineralización de nitrógeno del suelo fue mayormente afectada por la humedad del suelo y el efecto de la cobertura de copa según la calidad del sitio forestal. Finalmente con los datos obtenidos e información de bibliografía se generó un modelo de simulación de productividad de MS y concentración de PB de gramíneas creciendo en bosques de ñire. El modelo generado después de ser sometido a diferentes pruebas se considera aceptable para ser usado como una herramienta de manejo sustentable en estos sistemas, aportando además nuevos conocimientos sobre las causales en la variación productiva y de calidad de los pastizales en bosques de ñire.

Un modelo de simulación de la productividad actual y potencial de un pastizal de la Depresión del Salado

p.95-101

Perspectivas de corto y mediano plazo de la oferta de papa con relación al comportamiento de la demanda : un modelo econométrico de simulación dinámica

p.169

Efecto de variables ambientales sobre la productividad primaria neta aérea y la concentración de proteína bruta de gramíneas en sistemas silvopastoriles de ñire (Nothofagus antarctica): creación de un modelo de simulación

Nothofagus antarctica (ñire)ocupa una importante superficie de los bosques de Patagonia Sur, siendo ampliamente usado como sistemas silvopastoriles con ganadería ovina y bovina. Sin embargo, en la mayoría de los casos el uso de estos bosques se realiza sin un manejo que permita optimizar la sustentabilidad del sistema. En esta tesis se estudió el efecto de distintas variables ambientales (temperatura, humedad del suelo, radiación)sobre la productividad y concentración de Proteína bruta (PB)de gramíneas en distintas situaciones de calidad forestal y coberturas de copas de bosques de ñire, así como también procesos asociados: descomposición y dinámica de nutrientes (mineralización-inmovilización)de residuos del estrato arbóreo y graminoso y mineralización de nitrógeno del suelo. La productividad de materia seca (MS)aérea de gramíneas (rango entre 135 a 1863 kg MS ha-1)estuvo positivamente asociada a la calidad forestal de los sitios, mientras que el efecto de las coberturas de copa dependió de la interacción con otras variables ambientales. Los valores de PB variaron principalmente según la cobertura de copa encontrándose los mayores valores en lugares más sombreados (transmisividades menor al 70 por ciento). La descomposición, tanto en hojas de ñire como en gramíneas, varió según el nivel de radiación, siendo mayores en los lugares con menor cobertura. Por otro lado, la dinámica de nutrientes no fue afectada por la CS o alguna variable estudiada. La mineralización de nitrógeno del suelo fue mayormente afectada por la humedad del suelo y el efecto de la cobertura de copa según la calidad del sitio forestal. Finalmente con los datos obtenidos e información de bibliografía se generó un modelo de simulación de productividad de MS y concentración de PB de gramíneas creciendo en bosques de ñire. El modelo generado después de ser sometido a diferentes pruebas se considera aceptable para ser usado como una herramienta de manejo sustentable en estos sistemas, aportando además nuevos conocimientos sobre las causales en la variación productiva y de calidad de los pastizales en bosques de ñire.

Inferencia bayesiana sobre los parámetros de dispersión genéticos y ambientales en modelos animales con efectos maternos

Los modelos 'modelos animales con efectos maternos' (MAM)son modelos lineales mixtos que se utilizan para ajustar registros de caracteres bajo la influencia de efectos maternos. Uno de los desafíos más importantes en el marco de los MAM es la estimación de los parámetros de dispersión o 'componentes de (co)varianza' (CVC). En esta tesis se introducen desde una perspectiva bayesiana contribuciones teóricas y metodológicas con relación a la estimación de CVC para MAM sujetos a estructuras de covarianza novedosas. En primer lugar, se describe una implementación del análisis bayesiano jerárquico vía el algoritmo del muestreo de Gibbs. Luego, se considera una especificación conjugada diferente para la distribución a priori de la matriz de covarianza genética, basada en la distribución Wishart invertida generalizada, y se presenta una estrategia para determinar los correspondientes hiperparámetros. Esta estrategia fue comparada contra otras especificaciones a priori mediante un estudio de simulación estocástica, y produjo estimaciones precisas de los parámetros genéticos, con menores errores estándares y mejor tasa de convergencia. En segundo lugar, se presenta una formulación alternativa del MAM que incluye un parámetro de correlación ambiental entre pares de observaciones madre-progenie, y se desarrolla un procedimiento de estimación basado en un algoritmo de muestreo por grilla. El procedimiento fue programado y ejecutado exitosamente, y se obtuvo la primera estimación del parámetro de correlación con datos de campo para peso al destete en bovinos de carne. Por último, se considera el problema de la estimación de CVC en una población multirracial, donde en general es necesario especificar una estructura de covarianza heterogénea para los valores de cría. En particular, se demuestra que el modelo basado en la descomposición de la matriz de covarianza genética es equivalente al que deriva de la teoría genética cuantitativa. Además, se extiende el modelo para incluir efectos maternos y se describe la implementación de un análisis bayesiano jerárquico con el objetivo de estimar los CVC. El procedimiento fue implementado con éxito en datos experimentales de peso al destete y se obtuvieron por primera vez estimaciones para el conjunto completo de CVC.

Promedio bayesiano de modelos para la predicción del consumo máximo de oxígeno en deportistas con datos de no-ejercicio

Históricamente, los modelos de no-ejercicio para predecir el consumo máximo de oxígeno (VO2max) han sido construidos mediante regresión lineal frecuentista, usando técnicas estándar de selección de modelos. Sin embargo, existe incertidumbre acerca de la estructura estadística en el proceso de selección del modelo. En este estudio se propuso construir un modelo de no-ejercicio para predecir el VO2max en deportistas orientados al rendimiento, considerando la incertidumbre de modelo a través del Promedio Bayesiano de Modelos (BMA). Un objetivo adicional fue comparar la performance predictiva del BMA con las de los modelos derivados de varias técnicas frecuentistas usuales de selección de variables. Con tal fin, se implementó un submuestreo aleatorio estratificado repetido. Los datos incluyeron observaciones de la variable respuesta (en L·min-1), así como registros de Género, Deporte, Edad, Peso, Talla e Índice de masa corporal (BMI) (Edad = 22.1 ± 4.9 años, media ± SD; n = 272). Se propuso una clasificación de deportes con el objetivo de incluirla dentro del proceso de construcción del modelo: Combate, Juego, Resistencia 1 y Resistencia 2. El enfoque BMA se implementó en base a dos métodos: Occam's window y Composición de Modelo mediante el método de Monte Carlo con Cadenas de Markov (MC²). Se observaron discrepancias en la selección de variables entre los procedimientos frecuentistas. Ambos métodos de BMA produjeron resultados muy similares. Los modelos que incluyeron Género y las variables dummies para Resistencia 1 y Resistencia 2 acumularon virtualmente toda la probabilidad de modelo a posteriori. El Peso fue el predictor continuo con la más alta probabilidad de inclusión a posteriori (menor a 0.8). Las combinaciones de variables que involucraron predictores con un alto nivel de multicolinealidad fueron desacreditadas. Los modelos con sustancial contribución para el BMA presentaron un ajuste apreciable (R² ajustado menor a 0.8). Entre los modelos seleccionados por estrategias frecuentistas, el obtenido mediante el método de regresión por pasos (Stepwise regression method) con alfa igual a 0.05 fue el más respaldado por los datos, en términos de probabilidad de modelo a posteriori. En concordancia con la literatura, el BMA tuvo mejor performance predictiva de los datos fuera de la muestra que los modelos seleccionados por técnicas frecuentistas, medida por la cobertura del intervalo de predicción de 90 por ciento. La clasificación de deportes reveló resultados consistentes.

Influencia de los incendios en la dinámica poblacional de dos arbustos dominantes en la sabana mesopotámica

Los incendios son disturbios frecuentes en muchos ecosistemas terrestres como pastizales, sabanas, arbustales y bosques. Las perturbaciones del régimen natural de incendios pueden tener profunda influencia sobre la estructura y la dinámica de las poblaciones y comunidades y provocar invasiones o extinciones locales de especies. Estudios realizados en las sabanas mesopotámicas indican que los incendios serían parte del régimen natural de disturbios de esos ecosistemas. Desde la creación del Parque Nacional El Palmar se produjeron cambios en la fisonomía de la vegetación notables, entre los cuales el aumento en la abundancia de Baccharis dracunculifolia y Eupatorium buniifolium fue uno de los más notables. En esta tesis evalué las influencias de los incendios en interacción con otros factores como variaciones en el suelo, en la conectividad entre poblaciones y en las interacciones entre individuos sobre la distribución y demografía estas dos especies de arbustos. Realicé censos en 80 sitios de muestreo distribuidos en el parque con los que caractericé la distribución, abundancia y estructura poblacional de las dos especies en relación con la fisonomía de la vegetación, la textura del suelo y la fecha del último incendio. Estimé el reclutamiento y la mortalidad de las 80 poblaciones locales y construí modelos matriciales de la dinámica poblacional de un gran número de ellas. Con estos datos describí los controles de las tasas y realicé experimentos de simulación con los que evalué los efectos de la frecuencia de incendios y del aislamiento de las poblaciones locales sobre las tendencias demográficas a mediano plazo de cada especie. En una porción del paisaje, estudié la disposición espacial de los individuos de estas especies antes y después de un incendio. Mediante análisis de patrones de puntos puse a prueba hipótesis sobre la presencia de interacciones entre individuos mediadas por el fuego. Baccharis dracunculifolia es una especie ampliamente distribuida en el parque, tanto en sitios arenosos como arcillosos. La estructura de las poblaciones está controlada principalmente por el tiempo desde el último incendio y por la fisonomía de la vegetación. Todos los individuos mueren luego de un incendio y el reclutamiento se produce principalmente luego de los incendios; estas dos características definen la estrategia semillantes de respuesta a los incendios. Mis experimentos de simulación sugieren que es necesaria que ocurran incendios cada 5-8 años así como una migración importante de semillas entre poblaciones para mantener un sitio ocupado en el mediano y largo plazo. Esto sugiere que la dinámica poblacional en el paisaje del parque sigue un modelo de fuente sumidero, con poblaciones con alta densidad de adultos que aportan semillas a poblaciones vecinas. Eupatorium buniifolium tiene una distribución más acotada en el paisaje del Parque y está restringida a los sitios relativamente arcillosos. Con respecto a su estrategia de respuesta a los incendios, esta especie se comporta como rebrotante facultativa, con una gran proporción de individuos que sobrevive los incendios y un pulso de reclutamiento posterior al mismo. En los experimentos de simulación, la permanencia de las poblaciones Eupatorium en el largo plazo estuvo asegurada por un rango amplio de frecuencia de incendios de entre 3 y 25 años. Además, cuando el tiempo entre incendios era menor a 25 años, la migración de propágulos no tuvo efecto sobre la dinámica de las poblaciones locales. Los análisis de patrones de puntos permitieron demostrar la existencia de interacciones mediadas por el fuego, donde la cercanía de individuos de Baccharis aumenta la mortalidad a causa del fuego de los individuos de Eupatorium. El espacio liberado por la muerte de los individuos de Eupatorium es parcialmente ocupado por plántulas de Baccharis.

Modelos animales para caracteres con efectos maternos, madres no identificadas e información incompleta

La evaluación genética para caracteres de crecimiento pre - destete requiere ajustar modelos animales con efectos maternos (MAM). Tanto la estimación paramétrica de la variabilidad como la evaluación genética mediante MAM son realizadas empleando datos de campo, muchos de los cuales no poseen información completa para todas las variables explicativas maternas. Es común no contar con la identificación de madres (biológicas y/o receptoras), de abuelas maternas y, consecuentemente, de la edad de la madre (EM). Este problema es bien marcado en razas compuestas como Brangus y Braford que tienen políticas para registrar animales de pedigrí "abierto". Además, no existe un consenso sobre cuál es el mejor modelo de predicción, y existen interrogantes sobre la magnitud de los componentes de (co) varianza genético-aditivos y ambientales del modelo de evaluación. La primera investigación de esta tesis consistió en la estimación, mediante métodos bayesianos de los parámetros de dispersión en MAMs con distintas estructuras de (co) varianza, para datos de peso al destete de animales Angus de pedigrí. El análisis se caracterizó por la originalidad en los muestreos de las distribuciones marginales posteriores de las covarianzas genéticas aditivas y de la correlación entre los efectos ambientales maternos permanentes de una vaca y sus hijas también madres. Con el objeto de especificar correctamente la fracción aditiva de las (co) varianzas cuando se desconocen las madres y/o abuelas maternas de los animales con datos, en otro capítulo se desarrollaron MAMs equivalentes que no requieren alargar los vectores de los valores de cría con madres o abuelas fantasmas. Finalmente, se desarrolló un modelo mixto que atenúa el sesgo por error de medición clásico en el efecto EM, e introduce splines penalizadas y una estructura de (co) variación autoregresiva de orden 1 para suavizar las covarianzas residuales Este modelo es apropiado para ajustar datos de animales nacidos por transplante embrionario con madres receptoras desconocidas

Calibración del modelo Cropsyst para simular el balance hídrico en cultivos invernales del centro oeste de Entre Ríos

El trabajo tiene por objetivo calibrar el modelo CropSyst para simular la producción de materia seca y el balance hídrico de los cultivos de trigo (Triticum aestivum L.) ciclo largo y ciclo corto, colza (Brassica napus), lino (Linum usitatissimum) y arveja proteica (Pisum sativum) en el centro oeste de Entre Ríos. Adicionalmente, se incluyó un análisis de sensibilidad de parámetros que podrían tener, a priori, un mayor impacto sobre el crecimiento y desarrollo y en especial sobre el balance hídrico. La optimización de Cropsyst se realizó utilizando datos obtenidos de un experimento a campo, realizado con un diseño en bloques completos al azar, con cuatro repeticiones. Los parámetros medidos fueron: fenología, humedad del suelo, radiación interceptada por el cultivo, biomasa, índice de área foliar, área foliar específica y rendimiento. Durante el proceso de calibración se utilizó un procedimiento iterativo para optimizar las salidas de crecimiento y desarrollo de los cultivos. Una vez optimizado el modelo agronómico se procedió a evaluar las predicciones relacionadas con el contenido de agua (cm-3 cm-3) por capas y la lámina de agua total (mm) hasta una profundidad de 50 cm. En relación con el contenido de agua se tomaron dos fechas aproximadas de observación y registro, seleccionadas en base al período crítico de cada especie y la lámina de agua total se evaluó durante todo el ciclo de los cultivos. A los efectos de simular la infiltración de agua en el suelo se usaron dos métodos incluidos en modelo: el método en cascada y el método de las diferencias finitas. En general, Cropsyst brindó estimaciones razonables en cuanto a la biomasa, intercepción de radiación y rendimiento. Las predicciones del contenido de agua en el suelo para cada profundidad resultaron poco satisfactorias. Sin embargo, cuando se consideró la lámina de agua a 50 cm. las predicciones fueron más aceptables. Con relación al análisis de sensibilidad, los parámetros evaluados no mostraron cambios significativos en la simulación del agua en el suelo ni en el resto de los indicadores analizados, sugiriendo que mejores ajustes en la predicción del contenido de agua deberían lograrse seleccionando otro tipo de parámetros. La revisión y adaptación de las funciones de pedo-transferencia para estimar los límites máximos y mínimos de almacenaje de agua en el suelo probablemente mejore los resultados obtenidos en este trabajo.