Assessing salinity leaching efficiency in three soils by the HYDRUS-1D and -2D simulations

Salinity leaching is necessary to sustain agricultural production in irrigated croplands. Improving salinity leaching efficiency not only conserves water but also reduces groundwater contamination. Current leaching requirement (LR) calculations are based on steady-state and one-dimensional (1D) approaches, and consequently, this LR concept may not be applicable to drip irrigation (approximately 2D), which is becoming more common due to its higher water use efficiency. The aims of this study were to assess the salinity leaching fraction (LF) in clay, loam, and sand soils under 1D (to mimic sprinkler irrigation) and 2D (to mimic drip irrigation) transient conditions with a numerical model (HYDRUS). Water applications used the actual irrigation scheme in an almond orchard located in central California without considering precipitation. Model simulations showed that soil salinity at the lower boundary (depth of 150 cm) reached steady-state in 10 years in HYDRUS-1D simulations. The leaching fractions calculated from the ratio of drainage-water depth to irrigation-water depth (LFw = Ddw/Diw) and irrigation-water salinity to drainage-water salinity (LFEC = ECiw/ECdw) from HYDRUS-1D were similar among different textured soils. However, they were much higher under drip irrigation (2D) than under sprinkler irrigation (1D) when the same amount of water was applied, and LFEC values were much greater than the LFw values under 2D simulations. Salt balance (SB) and leaching efficiency (LE) indicated that sprinkler irrigation (1D) is more effective for salinity leaching than drip irrigation (2D). To improve salinity leaching efficiency, further evaluation of LRs under drip irrigation is needed

Optimizing drip irrigation for eggplant crops in semi-arid zones using evolving thresholds

AbstractField experiments were combined with a numerical model to optimize drip irrigation management based on soil matric potential (SMP) measurements. An experimental crop of eggplant was grown in Burkina Faso from December 2014 to March 2015 and plant response to water stress was investigated by applying four different irrigation treatments. Treatments consisted in using two different irrigation depths (low or high), combined with a water provision of 150%, 100% or 66% (150/100/66) of the maximum crop evapotranspiration (T150low, T66low, T100high, T66high). Soil matric potential measurements at 5, 10 and 15cm depth were taken using a wireless sensor network and were compared with measurements of plant and root biomass and crop yields. Field data were used to calibrate a numerical model to simulate triggered drip irrigation. Different simulations were built using the software HYDRUS 2D/3D to analyze the impact of the irrigation depth and frequency, the irrigation threshold and the soil texture on plant transpiration and water losses. Numerical results highlighted the great impact of the root distribution on the soil water dynamics and the importance of the sensor location to define thresholds. A fixed optimal sensor depth of 10 cm was found to manage irrigation from the vegetative state to the end of fruit development. Thresholds were defined to minimize water losses while allowing a sufficient soil water availability for optimal crop production. A threshold at 10cm depth of −15kPa is recommended for the early growth stage and −40kPa during the fruit formation and maturation phase. Simulations showed that those thresholds resulted in optimal transpiration regardless of the soil texture so that this management system can constitute the basis of an irrigation schedule for eggplant crops and possibly other vegetable crops in semi-arid regions

Spatial distribution of photoelectrons participating in formation of x-ray absorption spectra

Interpretation of x-ray absorption near-edge structure (XANES) experiments is often done via analyzing the role of particular atoms in the formation of specific peaks in the calculated spectrum. Typically, this is achieved by calculating the spectrum for a series of trial structures where various atoms are moved and/or removed. A more quantitative approach is presented here, based on comparing the probabilities that a XANES photoelectron of a given energy can be found near particular atoms. Such a photoelectron probability density can be consistently defined as a sum over squares of wave functions which describe participating photoelectron diffraction processes, weighted by their normalized cross sections. A fine structure in the energy dependence of these probabilities can be extracted and compared to XANES spectrum. As an illustration of this novel technique, we analyze the photoelectron probability density at the Ti K pre-edge of TiS2 and at the Ti K-edge of rutile TiO2.Comment: Journal abstract available on-line at http://link.aps.org/abstract/PRB/v65/e20511

Exploring the effects of topoisomerase II inhibitor XK469 on anthracycline cardiotoxicity and DNA damage

Anthracyclines, such as doxorubicin (adriamycin), daunorubicin, or epirubicin, rank among the most effective agents in classical anticancer chemotherapy. However, cardiotoxicity remains the main limitation of their clinical use. Topoisomerase IIβ has recently been identified as a plausible target of anthracyclines in cardiomyocytes. We examined the putative topoisomerase IIβ selective agent XK469 as a potential cardioprotective and designed several new analogues. In our experiments, XK469 inhibited both topoisomerase isoforms (α and β) and did not induce topoisomerase II covalent complexes in isolated cardiomyocytes and HL-60, but induced proteasomal degradation of topoisomerase II in these cell types. The cardioprotective potential of XK469 was studied on rat neonatal cardiomyocytes, where dexrazoxane (ICRF-187), the only clinically approved cardioprotective, was effective. Initially, XK469 prevented daunorubicin-induced toxicity and p53 phosphorylation in cardiomyocytes. However, it only partially prevented the phosphorylation of H2AX and did not affect DNA damage measured by Comet Assay. It also did not compromise the daunorubicin antiproliferative effect in HL-60 leukemic cells. When administered to rabbits to evaluate its cardioprotective potential in vivo, XK469 failed to prevent the daunorubicin induced cardiac toxicity in either acute or chronic settings. In the following in vitro analysis, we found that prolonged and continuous exposure of rat neonatal cardiomyocytes to XK469 led to significant toxicity. In conclusion, this study provides important evidence on the effects of XK469 and its combination with daunorubicin in clinically relevant doses in cardiomyocytes. Despite its promising characteristics, long-term treatments and in vivo experiments have not confirmed its cardioprotective potential

Surface softening in metal-ceramic sliding contacts: An experimental and numerical investigation

This study investigates the tribolayer properties at the interface of ceramic/metal (i.e., WC/W) sliding contacts using various experimental approaches and classical atomistic simulations. Experimentally, nanoindentation and micropillar compression tests, as well as adhesion mapping by means of atomic force microscopy, are used to evaluate the strength of tungsten?carbon tribolayers. To capture the influence of environmental conditions, a detailed chemical and structural analysis is performed on the worn surfaces by means of XPS mapping and depth profiling along with transmission electron microscopy of the debris particles. Experimentally, the results indicate a decrease in hardness and modulus of the worn surface compared to the unworn one. Atomistic simulations of nanoindentation on deformed and undeformed specimens are used to probe the strength of the WC tribolayer and despite the fact that the simulations do not include oxygen, the simulations correlate well with the experiments on deformed and undeformed surfaces, where the difference in behavior is attributed to the bonding and structural differences of amorphous and crystalline W-C. Adhesion mapping indicates a decrease in surface adhesion, which based on chemical analysis is attributed to surface passivation

Stochastic upscaling of hydrodynamic dispersion and retardation factor in a physically and chemically heterogeneous tropical soil

[EN] Stochastic upscaling of flow and reactive solute transport in a tropical soil is performed using real data collected in the laboratory. Upscaling of hydraulic conductivity, longitudinal hydrodynamic dispersion, and retardation factor were done using three different approaches of varying complexity. How uncertainty propagates after upscaling was also studied. The results show that upscaling must be taken into account if a good reproduction of the flow and transport behavior of a given soil is to be attained when modeled at larger than laboratory scales. The results also show that arrival time uncertainty was well reproduced after solute transport upscaling. This work represents a first demonstration of flow and reactive transport upscaling in a soil based on laboratory data. It also shows how simple upscaling methods can be incorporated into daily modeling practice using commercial flow and transport codes.The authors thank the financial support by the Brazilian National Council for Scientific and Technological Development (CNPq) (Project 401441/2014-8). The doctoral fellowship award to the first author by the Coordination of Improvement of Higher Level Personnel (CAPES) is acknowledged. The first author also thanks the international mobility grant awarded by CNPq, through the Sciences Without Borders program (Grant Number: 200597/2015-9). The international mobility grant awarded by Santander Mobility in cooperation with the University of Sao Paulo is also acknowledged. DHI-WASI is gratefully thanked for providing a FEFLOW license.Almeida De-Godoy, V.; Zuquette, L.; Gómez-Hernández, JJ. (2019). Stochastic upscaling of hydrodynamic dispersion and retardation factor in a physically and chemically heterogeneous tropical soil. Stochastic Environmental Research and Risk Assessment. 33(1):201-216. https://doi.org/10.1007/s00477-018-1624-zS201216331Ahuja LR, Naney JW, Green RE, Nielsen DR (1984) Macroporosity to characterize spatial variability of hydraulic conductivity and effects of land management. Soil Sci Soc Am J 48:699. https://doi.org/10.2136/sssaj1984.03615995004800040001xBellin A, Lawrence AE, Rubin Y (2004) Models of sub-grid variability in numerical simulations of solute transport in heterogeneous porous formations: three-dimensional flow and effect of pore-scale dispersion. Stoch Environ Res Risk Assess 18:31–38. https://doi.org/10.1007/s00477-003-0164-2Brent RP (1973) Algorithms for minimization without derivatives. Prentice Hall, Englewood CliffsBrusseau ML (1998) Non-ideal transport of reactive solutes in heterogeneous porous media: 3. model testing and data analysis using calibration versus prediction. J Hydrol 209:147–165. https://doi.org/10.1016/S0022-1694(98)00121-8Brusseau ML, Srivastava R (1999) Nonideal transport of reactive solutes in heterogeneous porous media: 4. Analysis of the cape cod natural-gradient field experiment. Water Resour Res 35:1113–1125. https://doi.org/10.1029/1998WR900019Brutsaert W (1967) Some methods of calculating unsaturated permeability. Trans ASAE 10:400–404Cadini F, De Sanctis J, Bertoli I, Zio E (2013) Upscaling of a dual-permeability Monte Carlo simulation model for contaminant transport in fractured networks by genetic algorithm parameter identification. Stoch Environ Res Risk Assess 27:505–516. https://doi.org/10.1007/s00477-012-0595-8Cambardella CA, Moorman TB, Parkin TB, Karlen DL, Novak JM, Turco RF, Konopka AE (1994) Field-scale variability of soil properties in central iowa soils. Soil Sci Soc Am J 58:1501. https://doi.org/10.2136/sssaj1994.03615995005800050033xCapilla JE, Rodrigo J, Gómez-Hernández JJ (1999) Simulation of non-Gaussian transmissivity fields honoring piezometric data and integrating soft and secondary information. Math Geol 31:907–927. https://doi.org/10.1023/A:1007580902175Cassiraga EF, Fernàndez-Garcia D, Gómez-Hernández JJ (2005) Performance assessment of solute transport upscaling methods in the context of nuclear waste disposal. Int J Rock Mech Min Sci 42:756–764. https://doi.org/10.1016/j.ijrmms.2005.03.013Corey AT (1977) Mechanics of heterogeneous fluids in porous media. Water Resources Publications, Fort Collins, CO, p 259Dagan G (1989) Flow and transport in porous formations. Springer, Berlin. https://doi.org/10.1007/978-3-642-75015-1Dagan G (2004) On application of stochastic modeling of groundwater flow and transport. Stoch Environ Res Risk Assess. https://doi.org/10.1007/s00477-004-0191-7de Azevedo AAB, Pressinotti MMN, Massoli M (1981) Sedimentological studies of the Botucatu and Pirambóia formations in the region of Santa Rita do Passa Quatro (In portuguese). Rev do Inst Geológico 2:31–38. https://doi.org/10.5935/0100-929X.19810003Deng H, Dai Z, Wolfsberg AV, Ye M, Stauffer PH, Lu Z, Kwicklis E (2013) Upscaling retardation factor in hierarchical porous media with multimodal reactive mineral facies. Chemosphere 91:248–257. https://doi.org/10.1016/j.chemosphere.2012.10.105Diersch H-JG (2014) Finite element modeling of flow, mass and heat transport in porous and fractured media. Springer, Berlin. https://doi.org/10.1007/978-3-642-38739-5Dippenaar MA (2014) Porosity reviewed: quantitative multi-disciplinary understanding, recent advances and applications in vadose zone hydrology. Geotech Geol Eng 32:1–19. https://doi.org/10.1007/s10706-013-9704-9Fagundes JRT, Zuquette LV (2011) Sorption behavior of the sandy residual unconsolidated materials from the sandstones of the Botucatu Formation, the main aquifer of Brazil. Environ Earth Sci 62:831–845. https://doi.org/10.1007/s12665-010-0570-yFenton GA, Griffiths DV (2008) Risk assessment in geotechnical engineering. Wiley, p 463Fernàndez-Garcia D, Gómez-Hernández JJ (2007) Impact of upscaling on solute transport: Traveltimes, scale dependence of dispersivity, and propagation of uncertainty. Water Resour Res. https://doi.org/10.1029/2005WR004727Fernàndez-Garcia D, Llerar-Meza G, Gómez-Hernández JJ (2009) Upscaling transport with mass transfer models: mean behavior and propagation of uncertainty. Water Resour Res. https://doi.org/10.1029/2009WR007764Feyen L, Gómez-Hernández JJ, Ribeiro PJ, Beven KJ, De Smedt F (2003a) A Bayesian approach to stochastic capture zone delineation incorporating tracer arrival times, conductivity measurements, and hydraulic head observations. Water Resour Res. https://doi.org/10.1029/2002WR001544Feyen L, Ribeiro PJ, Gómez-Hernández JJ, Beven KJ, De Smedt F (2003b) Bayesian methodology for stochastic capture zone delineation incorporating transmissivity measurements and hydraulic head observations. J Hydrol 271:156–170. https://doi.org/10.1016/S0022-1694(02)00314-1Forsythe GE, Malcolm MA, Moler CB (1976) Computer methods for mathematical computations. Prentice-Hall, Englewood Cliffs, p 259Freeze R, Cherry J (1979) Groundwater. PrenticeHall Inc, Englewood cliffs, p 604Frippiat CC, Holeyman AE (2008) A comparative review of upscaling methods for solute transport in heterogeneous porous media. J Hydrol 362:150–176. https://doi.org/10.1016/j.jhydrol.2008.08.015Fu J, Gómez-Hernández JJ (2009) Uncertainty assessment and data worth in groundwater flow and mass transport modeling using a blocking Markov chain Monte Carlo method. J Hydrol 364:328–341. https://doi.org/10.1016/j.jhydrol.2008.11.014Gelhar LW, Axness CL (1983) Three-dimensional stochastic analysis of macrodispersion in aquifers. Water Resour Res 19:161–180. https://doi.org/10.1029/WR019i001p00161Gelhar LW, Welty C, Rehfeldt KR (1992) A critical review of data on field-scale dispersion in aquifers. Water Resour Res 28:1955–1974. https://doi.org/10.1029/92WR00607Giacheti HL, Rohm SA, Nogueira JB, Cintra JCA (1993) Geotechnical properties of the Cenozoic sediment (in protuguese). In: Albiero JH, Cintra JCA (eds) Soil from the interior of São Paulo. ABMS, Sao Paulo, pp 143–175Gómez-Hernandez JJ (1990) A stochastic approach to the simulation of block conductivity fields conditional upon data measured at a smaller scale. Stanford University, StanfordGómez-Hernández JJ, Gorelick SM (1989) Effective groundwater model parameter values: influence of spatial variabiity of hydraulic conductivity, leackance, and recharge. Water Resour Res 25:405–419Gómez-Hernández JJ, Journel A (1993) Joint sequential simulation of multigaussian fields. In: Geostatistics Tróia’92. pp 85–94. https://doi.org/10.1007/978-94-011-1739-5_8Gómez-Hernández JJ, Wen X-H (1994) Probabilistic assessment of travel times in groundwater modeling. Stoch Hydrol Hydraul 8:19–55. https://doi.org/10.1007/BF01581389Gómez-Hernández JJ, Fu J, Fernandez-Garcia D (2006) Upscaling retardation factors in 2-D porous media. In: Bierkens MFP, Gehrels JC, Kovar K (eds) Calibration and reliability in groundwater modelling: from uncertainty to decision making: proceedings of the ModelCARE 2005 conference held in The Hague, The Netherlands, 6–9 June, 2005. IAHS Publication, pp 130–136Goovaerts P (1999) Geostatistics in soil science: state-of-the-art and perspectives. Geoderma 89:1–45. https://doi.org/10.1016/S0016-7061(98)00078-0Jarvis NJ (2007) A review of non-equilibrium water fl ow and solute transport in soil macropores: principles, controlling factors and consequences for water quality. Eur J Soil Sci 58:523–546. https://doi.org/10.4141/cjss2011-050Jellali S, Diamantopoulos E, Kallali H, Bennaceur S, Anane M, Jedidi N (2010) Dynamic sorption of ammonium by sandy soil in fixed bed columns: evaluation of equilibrium and non-equilibrium transport processes. J Environ Manag 91:897–905. https://doi.org/10.1016/j.jenvman.2009.11.006Journel AG, Gomez-Hernandez JJ (1993) Stochastic imaging of the wilmington clastic sequence. SPE Form Eval 8:33–40. https://doi.org/10.2118/19857-PAJournel A, Deutsch C, Desbarats A (1986) Power averaging for block effective permeability. Proc SPE Calif Reg Meet. https://doi.org/10.2118/15128-MSKronberg BI, Fyfe WS, Leonardos OH, Santos AM (1979) The chemistry of some Brazilian soils: element mobility during intense weathering. Chem Geol 24:211–229. https://doi.org/10.1016/0009-2541(79)90124-4Lake LW (1988) The origins of anisotropy (includes associated papers 18394 and 18458). J Pet Technol 40:395–396. https://doi.org/10.2118/17652-PALawrence AE, Rubin Y (2007) Block-effective macrodispersion for numerical simulations of sorbing solute transport in heterogeneous porous formations. Adv Water Resour 30:1272–1285. https://doi.org/10.1016/j.advwatres.2006.11.005Lemke LD, Barrack WA II, Abriola LM, Goovaerts P (2004) Matching solute breakthrough with deterministic and stochastic aquifer models. Groundwater 42:920–934Li L, Zhou H, Gómez-Hernández JJ (2011a) A comparative study of three-dimensional hydraulic conductivity upscaling at the macro-dispersion experiment (MADE) site, Columbus Air Force Base, Mississippi (USA). J Hydrol 404:278–293. https://doi.org/10.1016/j.jhydrol.2011.05.001Li L, Zhou H, Gómez-Hernández JJ (2011b) Transport upscaling using multi-rate mass transfer in three-dimensional highly heterogeneous porous media. Adv Water Resour 34:478–489. https://doi.org/10.1016/j.advwatres.2011.01.001Logsdon Keller KE, Moorman TB (2002) Measured and predicted solute leaching from multiple undisturbed soil columns. Soil Sci Soc Am J 66:686–695. https://doi.org/10.2136/sssaj2002.6860Lourens A, van Geer FC (2016) Uncertainty propagation of arbitrary probability density functions applied to upscaling of transmissivities. Stoch Environ Res Risk Assess 30:237–249. https://doi.org/10.1007/s00477-015-1075-8Mahapatra IC, Singh KN, Pillai KG, Bapat SR (1985) Rice soils and their management. Indian J Agron 30:R1–R41Morakinyo JA, Mackay R (2006) Geostatistical modelling of ground conditions to support the assessment of site contamination. Stoch Environ Res Risk Assess 20:106–118. https://doi.org/10.1007/s00477-005-0015-4Moslehi M, de Barros FPJ, Ebrahimi F, Sahimi M (2016) Upscaling of solute transport in disordered porous media by wavelet transformations. Adv Water Resour 96:180–189. https://doi.org/10.1016/j.advwatres.2016.07.013Osinubi KJ, Nwaiwu CM (2005) Hydraulic conductivity of compacted lateritic soil. J Geotech Geoenviron Eng 131:1034–1041. https://doi.org/10.1061/(ASCE)1090-0241(2005)131:8(1034)Remy N (2004) SGeMS: stanford geostatistical modeling software. Softw Man. https://doi.org/10.1007/978-1-4020-3610-1_89Renard P, de Marsily G (1997) Calculating equivalent permeability: a review. Adv Water Resour 20:253–278. https://doi.org/10.1016/S0309-1708(96)00050-4Robin MJL, Sudicky EA, Gillham RW, Kachanoski RG (1991) Spatial variability of strontium distribution coefficients and their correlation with hydraulic conductivity in the Canadian forces base borden aquifer. Water Resour Res 27:2619–2632. https://doi.org/10.1029/91WR01107Salamon P, Fernàndez-Garcia D, Gómez-Hernández JJ (2007) Modeling tracer transport at the MADE site: the importance of heterogeneity. Water Resour Res. https://doi.org/10.1029/2006WR005522Sánchez-Vila X, Carrera J, Girardi JP (1996) Scale effects in transmissivity. J Hydrol 183:1–22. https://doi.org/10.1016/S0022-1694(96)80031-XScheibe T, Yabusaki S (1998) Scaling of flow and transport behavior in heterogeneous groundwater systems. Adv Water Resour 22:223–238. https://doi.org/10.1016/S0309-1708(98)00014-1Selvadurai PA, Selvadurai APS (2014) On the effective permeability of a heterogeneous porous medium: the role of the geometric mean. Philos Mag 94:2318–2338. https://doi.org/10.1080/14786435.2014.913111Shackelford CD (1994) Critical concepts for column testing. J Geotech Eng 120:1804–1828. https://doi.org/10.1016/0148-9062(95)96996-OŠimůnek J, van Genuchten MT, Šejna M, Toride N, Leij FJ (1999) The STANMOD computer software for evaluating solute transport in porous media using analytical solutions of convection-dispersion equation. Riverside, CaliforniaTaskinen A, Sirviö H, Bruen M (2008) Modelling effects of spatial variability of saturated hydraulic conductivity on autocorrelated overland flow data: linear mixed model approach. Stoch Environ Res Risk Assess 22:67–82. https://doi.org/10.1007/s00477-006-0099-5Tuli A, Hopmans JW, Rolston DE, Moldrup P (2005) Comparison of air and water permeability between disturbed and undisturbed soils. Soil Sci Soc Am J 69:1361. https://doi.org/10.2136/sssaj2004.0332Tyukhova AR, Willmann M (2016) Conservative transport upscaling based on information of connectivity. Water Resour Res 52:6867–6880. https://doi.org/10.1002/2015WR018331van Genuchten MTh (1980) Determining transport parameters from solute displacement experiments. Research Report 118. U.S. Salinity Lab., Riverside, CAVanderborght J, Timmerman A, Feyen J (2000) Solute transport for steady-state and transient flow in soils with and without macropores. Soil Sci Soc Am J 64:1305–1317. https://doi.org/10.2136/sssaj2000.6441305xVanmarcke E (2010) Random fields: analysis and synthesis. World Scientific. MIT Press, Cambridge, MA, p 364Vishal V, Leung JY (2017) Statistical scale-up of 3D particle-tracking simulation for non-Fickian dispersive solute transport modeling. Environ Res Risk Assess, Stoch. https://doi.org/10.1007/s00477-017-1501-1Wen X-H, Gómez-Hernández JJ (1996) Upscaling hydraulic conductivities in heterogeneous media: an overview. J Hydrol 183:ix–xxxii. https://doi.org/10.1016/S0022-1694(96)80030-8Wen XH, Gómez-Hernández JJ (1998) Numerical modeling of macrodispersion in heterogeneous media: a comparison of multi-Gaussian and non-multi-Gaussian models. J Contam Hydrol 30:129–156. https://doi.org/10.1016/S0169-7722(97)00035-1Wen XH, Capilla JE, Deutsch CV, Gómez-Hernández JJ, Cullick AS (1999) A program to create permeability fields that honor single-phase flow rate and pressure data. Comput Geosci 25:217–230. https://doi.org/10.1016/S0098-3004(98)00126-5Wilding LP, Drees LR (1983) Spatial variability and pedology. In: Wilding LP, Smeck NE, Hall GF (eds) Pedogenesis and soil taxonomy: the soil orders. Elsevier, Amsterdam, pp 83–116Willmann M, Carrera J, Guadagnini A (2006) Block-upscaling of transport in heterogeneous aquifers. h2ogeo.upc.edu 1–7Xu Z, Meakin P (2013) Upscaling of solute transport in heterogeneous media with non-uniform flow and dispersion fields. Appl Math Model 37:8533–8542. https://doi.org/10.1016/j.apm.2013.03.070Zech A, Attinger S, Cvetkovic V, Dagan G, Dietrich P, Fiori A, Rubin Y, Teutsch G (2015) Is unique scaling of aquifer macrodispersivity supported by field data? Water Resour Res 51:7662–7679. https://doi.org/10.1002/2015WR017220Zhou H, Li L, Gómez-Hernández JJ (2010) Three-dimensional hydraulic conductivity upscaling in groundwater modeling. Comput Geosci 36:1224–1235. https://doi.org/10.1016/j.cageo.2010.03.008Zhou H, Li L, Hendricks Franssen H-J, Gómez-Hernández JJ (2012) Pattern recognition in a bimodal aquifer using the normal-score ensemble Kalman filter. Math Geosci 44:169–185. https://doi.org/10.1007/s11004-011-9372-

Mechanistic framework to link root growth models with weather and soil physical properties, including example applications to soybean growth in Brazil

Background and aimsRoot elongation is generally limited by a combination of mechanical impedance and water stress in most arable soils. However, dynamic changes of soil penetration resistance with soil water content are rarely included in models for predicting root growth. Better modelling frameworks are needed to understand root growth interactions between plant genotype, soil management, and climate. Aim of paper is to describe a new model of root elongation in relation to soil physical characteristics like penetration resistance, matric potential, and hypoxia.MethodsA new diagrammatic framework is proposed to illustrate the interaction between root elongation, soil management, and climatic conditions. The new model was written in Matlab®, using the root architecture model RootBox and a model that solves the 1D Richards equations for water flux in soil. Inputs: root architectural parameters for Soybean; soil hydraulic properties; root water uptake function in relation to matric flux potential; root elongation rate as a function of soil physical characteristics. Simulation scenarios: (a) compact soil layer at 16 to 20 cm; (b) test against a field experiment in Brazil during contrasting drought and normal rainfall seasons.Results(a) Soil compaction substantially slowed root growth into and below the compact layer. (b) Simulated root length density was very similar to field measurements, which was influenced greatly by drought. The main factor slowing root elongation in the simulations was evaluated using a stress reduction function.ConclusionThe proposed framework offers a way to explore the interaction between soil physical properties, weather and root growth. It may be applied to most root elongation models, and offers the potential to evaluate likely factors limiting root growth in different soils and tillage regimes