Copper oxide nanoparticles (CuO NPs), as an antimicrobial nanomaterial, have found many applications in agriculture. Ubiquitous and complex root exudates (RE) in the plant root zone motivates the determination of how specific components of RE interact with CuO NPs. This work aims to reveal the role of maize (Zea mays L.)-derived RE and their components on the aggregation and dissolution of CuO NPs in the rhizosphere. We observed that RE significantly inhibited the aggregation of CuO NPs regardless of ionic strength and electrolyte type. In the presence of RE, the CCC of CuO NPs in NaCl shifted from 30 to 125 mM and the value in CaCl2 shifted from 4 to 20 mM. Furthermore, this inhibition was correlated with molecular weight (MW) of RE fractions. Higher MW fraction (>10 kDa) reduced the aggregation most. We also discovered that RE significantly promoted the dissolution of CuO NPs and lower MW fraction (<3 kDa) RE mainly contributed to this process. Additionally, phytotoxicity of CuO NPs in the presence of RE and different fractions of RE was evaluated. The addition of 20 mg/L RE reduced the seedlings growth rate to 1.89% after 7 days exposure to 25 mg/L CuO NPs, which were significantly lower than the control group (4.82%). Notably, Cu accumulation in plant root tissues was significantly enhanced by 20 mg/L RE. This study provides useful insights into the interactions between RE and CuO NPs, which is of significance for the safe use of CuO NPs-based antimicrobial products in agricultural production.

Similar to soils, soilless media are composed of three phases: solid, aqueous, and gaseous. The substantial differences between the former and the latter are related to the solid phase composition and associated texture and structure. While the solid phase of soils is mostly minerals, it is a mixture of minerals and organic matter for soilless culture, resulting in substantial differences in their physical and hydraulic properties. The physical and hydraulic properties of soil have been intensively studied, both experimentally and theoretically, the application and implementation of these knowledge and methods to soilless substrates is limited. This chapter aims to close the gaps in physics and hydraulic characterization and application between mineral soils and soilless media by applying terminology, methods, and approaches that are commonly used for the former to the latter. The current version of this Chapter extents the topic of particle wettability and its effect on the physical properties of soilless media. Although research on sales-media wettability and its effect on water retention and flow in the media has only recently initiated, its relevance is essential owing to the high organic matter content of these media and the relationship between the two.

Secondary treated wastewater, a commonly used water resource in agriculture in (semi-)arid areas, often contains salts, sodium, and organic matter which may affect soil structure and hydraulic properties. The main objective of this study was to jointly analyse the effects of long-term irrigation with treated wastewater on physicochemical soil characteristics, soil structure, and soil water dynamics in undisturbed soils. X-ray microtomography was used to determine changes in macro-porosity (> 19 μm), pore size distribution, and pore connectivity of a sandy clay loam and a loamy sand. Differences in the pore network among soils irrigated with treated wastewater, fresh water that replaced treated wastewater, and non-irrigated control plots could be explained by changes in textural composition, soil physicochemical parameters, and hydraulic properties. In this study we showed that irrigation led to the development of a connected macro-pore network, independent of the studied water quality. The leaching of silt and clay particles in the sandy soil due to treated wastewater irrigation resulted in an increase of pores < 130 μm. While this change in texture reduced water retention, the unsaturated hydraulic conductivity was diminished by physicochemical alteration, i.e. induced water repellency and clay mineral swelling. Overall, the fine textured sandy clay loam was much more resistant to soil alteration by treated wastewater irrigation than the loamy sand.

The recognition of treated wastewater (TWW) as an alternative water resource is expanding in areas with a shortage of freshwater (FW). While many studies have been devoted to the effects of long-term irrigation with TWW on soil wettability and spatial flow variations in the soil profile, much less attention has been given to the spatial distribution of soil water repellency in the soil surface layer. This is the objective of the current study. Undisturbed soil samples (5 cm thick) were taken at 15-cm intervals parallel to a drip lateral in two adjacent plots of a commercial citrus orchard in central Israel. Each soil sample was sectioned into five consecutive 1-cm layers for which soil water repellency was determined by water drop penetration time method, and soil organic matter by loss-on-ignition method. Geostatistics and multivariate empirical mode decomposition were used to investigate the overall and scale-specific spatial variation of soil water repellency and its dependence on dripper intervals along the lateral. A high degree of soil water repellency with strong spatial variation was found in the surface soil after 4–6 years of TWW irrigation. Weak to moderate spatial dependence of soil water repellency with maximum autocorrelation distance of around 30 cm was discovered by geostatistical analysis. The spatial distribution of soil water repellency was considered to be greatly affected by the location of the drippers, being higher between adjacent drippers and lower underneath them. This soil water repellency distribution is presumed to result from ongoing lateral displacement of the amphiphilic substances in the TWW toward the outer edge of the wetted plume periphery. Multivariate empirical mode decomposition of the overall spatial variation of soil water repellency yielded three scale-specific variations with corresponding characteristic scales of 30 cm, 110 cm and 200 cm. Most of the soil water repellency variation was separated into the 30 cm and 110 cm spatial scales, which were correlated to processes related to the drippers and trees. Replacing TWW with FW for the reclamation of water-repellent soils partially alleviated the intensity of TWW irrigation-induced soil water repellency. Moreover, an inconsistency between the hot spots of water-repellency development between adjacent drippers and the areas that are effectively ameliorated by FW irrigation below the drippers could be developed and affect the spatial distribution of flow pattern in an a priori unpredictable way.

The improvement of crop productivity under abiotic stress is one of the biggest challenges faced by the agricultural scientific community. Despite extensive research, the research-to-commercial transfer rate of abiotic stress-resistant crops remains very low. This is mainly due to the complexity of genotype◻×◻environment interactions and in particular, the ability to quantify the dynamic plant physiological response profile to a dynamic environment.Most existing phenotyping facilities collect information using robotics and automated image acquisition and analysis. However, their ability to directly measure the physiological properties of the whole plant is limited. We demonstrate a high-throughput functional phenotyping system (HFPS) that enables comparing plants’ dynamic responses to different ambient conditions in dynamic environments due to its direct and simultaneous measurement of yield-related physiological traits of plants under several treatments. The system is designed as one-to-one (1:1) plant–[sensors+controller] units, i.e., each individual plant has its own personalized sensor, controller and irrigation valves that enable (i) monitoring water-relation kinetics of each plant–environment response throughout the plant’s life cycle with high spatiotemporal resolution, (ii) a truly randomized experimental design due to multiple independent treatment scenarios for every plant, and (iii) reduction of artificial ambient perturbations due to the immobility of the plants or other objects. In addition, we propose two new resilience-quantifying-related traits that can also be phenotyped using the HFPS: transpiration recovery rate and night water reabsorption.We use the HFPS to screen the effects of two commercial biostimulants (a seaweed extract—ICL-SW, and a metabolite formula—ICL-NewFo1) on Capsicum annuum under different irrigation regimes. Biostimulants are considered an alternative approach to improving crop productivity. However, their complex mode of action necessitates cost-effective pre-field phenotyping. The combination of two types of treatment (biostimulants and drought) enabled us to evaluate the precision and resolution of the system in investigating the effect of biostimulants on drought tolerance. We analyze and discuss plant behavior at different stages, and assess the penalty and trade-off between productivity and survivability. In this test case, we suggest a protocol for the screening of biostimulants’ physiological mechanisms of action.

The identification of preferential flow generation and the tracking of their pathways in agricultural soils are challenging tasks, in particular while keeping the soil profile undisturbed. The noninvasive electrical resistivity tomography (ERT) method seems, a priori, to be an ideal tool for these purposes, but the simultaneous dependence of soil electrical resistivity (ER) on soil moisture content and salt concentration poses a problem when high-salinity water is involved. The dependence of these variables is expressed, for example, by Archie’s law. Irrigation with brackish or treated wastewater are two such complicated cases. The latter has been found to render soils water-repellent with inherent flow in preferential pathways. A method that resolves this complication is developed and applied for in-situ-measured ERT data. This method combines a model that simulates the preferential ER(t) as a series of interconnected well-mixed units (WMUs) with frequent ERT scans during soil wetting and subsequent drying. Following validation by comparison of its output to simulations using the Richards equation for flow and convection-dispersion equation for transport, the WMU model was used to simulate scenarios of simultaneous variations in moisture content, salinity, and ER(t) in a well-mixed soil volume (unit). The characteristic patterns acquired by these simulations were compared with in-situ measured ER(z,t) to determine whether the observed preferential ER pathways coincide with preferential flow or salinity pathways, or both. From the eight sections in the profile where this analysis was made for, five exhibit preferential flow, characterized by fast wetting front propagation that includes by-passing of a certain soil volume, while the wetting front velocity in the other three was moderate. We, therefore, concluded, based on the water content/salinity effect on ER, that the measured preferential ER pathways coincide with preferential flow pathways. Salinity has, for the current case, a minor, if any, effect on the ER spatial variation. Therefore, ERT is an advantageous mean of mapping of preferential flow pathways in the soil profile, and the method developed herein can be used for studying the water distribution in the soil profile. Particular irrigation design and management can be then used to ease the negative effects of preferential flow pathways on soil water availability to plant roots and the enhanced leaching of agrochemical below the active root zone.

Alternative transportation fuels (ATFs) can reduce air pollution. However, the influence of conventional fuels—diesel and gasoline, and particularly ATFs on photochemical air pollution is not well-characterized, limiting assessments of ATFs and air quality. This is mainly due to frequent use of lumped chemical mechanisms by related atmospheric modeling. Here we hypothesized that applying a chemical mechanism that is specifically developed according to both emission fractions and photochemical ozone creation potential of volatile organic compounds (VOCs) is key to gaining reliable insights into the impact of transportation fuels on photochemistry. We used a heterogeneous chemical mechanism with 927 reactions and relatively detailed emission inventories to specifically meet the requirements for reliable simulation of the effect of exhaust emissions from vehicles fueled by selected model fuels—diesel, gasoline, and mixtures of 15% gasoline with 85% ethanol (E85) or 85% methanol (M85)—on photochemistry. These dispersion-box model simulations revealed a strong influence of atmospheric background balance between VOCs and nitrogen oxides (NOX = [NO] + [NO2]) on the impact of exhaust emissions on photochemistry, with higher tendency toward ozone (O3) formation or destruction for more VOC-limited or NOX-limited conditions, respectively. Accordingly, higher [NOX]/[VOC] exhaust emission, such as from diesel and M85, resulted in lower O3, not only locally but also downwind of the emission. This offers a new perspective and measure for transportation fuel assessment. Rapid conversion of O3 to hydroxyl and hydroperoxyl radicals downwind of the exhaust emission indicates the importance of simulating the impact of road transportation on photochemistry at high spatial and temporal resolution. Peroxyacetyl nitrate formation was more sensitive to VOC emission under VOC-limited conditions than to NOX emission under NOX-limited conditions. Secondary formaldehyde dominated over primary emitted formaldehyde several minutes after emission. These findings should be verified using a 3D modeling study under varying meteorological conditions.

AbstractA reliable forecast of potential evapotranspiration (ET0) is key to precise irrigation scheduling toward reducing water and agrochemical use while optimizing crop yield. In this study, we examine the benefits of using the Weather Research and Forecasting (WRF) Model for ET0 and precipitation forecasts with simulations at a 3-km grid spatial resolution and an hourly temporal resolution output over Israel. The simulated parameters needed to calculate ET0 using the Penman?Monteith (PM) approach, as well as calculated ET0 and precipitation, were compared to observations from a network of meteorological stations. WRF forecasts of all PM meteorological parameters, except wind speed Ws, were significantly sensitive to seasonality and synoptic conditions, whereas forecasts of Ws consistently showed high bias associated with strong local effects, leading to high bias in the evaluated PM?ET0. Local Ws bias correction using observations on days preceding the forecast and interpolation of the resulting PM?ET0 to other locations led to significant improvement in ET0 forecasts over the investigated area. By using this hybrid forecast approach (WRFBC) that combines WRF numerical simulations with statistical bias corrections, daily ET0 forecast bias was reduced from an annual mean of 13% with WRF to 3% with WRFBC, while maintaining a high model?observation correlation. WRF was successful in predicting precipitation events on a daily event basis for all four forecast lead days. Considering the benefit of the hybrid approach for forecasting ET0, the WRF Model was found to be a high-potential tool for improving crop irrigation management.AbstractA reliable forecast of potential evapotranspiration (ET0) is key to precise irrigation scheduling toward reducing water and agrochemical use while optimizing crop yield. In this study, we examine the benefits of using the Weather Research and Forecasting (WRF) Model for ET0 and precipitation forecasts with simulations at a 3-km grid spatial resolution and an hourly temporal resolution output over Israel. The simulated parameters needed to calculate ET0 using the Penman?Monteith (PM) approach, as well as calculated ET0 and precipitation, were compared to observations from a network of meteorological stations. WRF forecasts of all PM meteorological parameters, except wind speed Ws, were significantly sensitive to seasonality and synoptic conditions, whereas forecasts of Ws consistently showed high bias associated with strong local effects, leading to high bias in the evaluated PM?ET0. Local Ws bias correction using observations on days preceding the forecast and interpolation of the resulting PM?ET0 to other locations led to significant improvement in ET0 forecasts over the investigated area. By using this hybrid forecast approach (WRFBC) that combines WRF numerical simulations with statistical bias corrections, daily ET0 forecast bias was reduced from an annual mean of 13% with WRF to 3% with WRFBC, while maintaining a high model?observation correlation. WRF was successful in predicting precipitation events on a daily event basis for all four forecast lead days. Considering the benefit of the hybrid approach for forecasting ET0, the WRF Model was found to be a high-potential tool for improving crop irrigation management.

Dry deposition of ozone (O3) to vegetation is an important removal pathway for tropospheric O3, while O3 uptake through plant stomata negatively affects vegetation and leads to climate change. Both processes are controlled by vegetation characteristics and ambient conditions via complex mechanisms. Recent studies have revealed that these processes can be fundamentally impacted by coastal effects, and by dry and warm conditions in ways that have not been fully characterized, largely due to lack of measurements under such conditions. Hence, we hypothesized that measuring dry deposition of O3 to vegetation along a sharp spatial climate gradient, and at different distances from the coast, can offer new insights into the characterization of these effects on O3 deposition to vegetation and stomatal uptake, providing important information for afforestation management and for climate and air-quality model improvement. To address these hypotheses, several measurement campaigns were performed at different sites, including pine, oak, and mixed Mediterranean forests, at distances of 20–59 km from the Eastern Mediterranean coast, under semiarid, Mediterranean and humid Mediterranean climate conditions. The eddy covariance technique was used to quantify vertical O3 flux (Ftot) and its partitioning to stomatal flux (Fst) and non-stomatal flux (Fns). Whereas Fst tended to peak around noon under humid Mediterranean and Mediterranean conditions in summer, it was strongly limited by drought under semiarid conditions from spring to early winter, with minimum average Fst/Ftot of 8–11% during the summer. Fns in the area was predominantly controlled by relative humidity (RH), whereas increasing Fns with RH for RH < 70% indicated enhancement of Fns by aerosols, via surface wetness stimulation. At night, efficient turbulence due to sea and land breezes, together with increased RH, resulted in strong enhancement of Ftot. Extreme dry surface events, some induced by dry intrusion from the upper troposphere, resulted in positive Fns events.

Abstract Wetting and drying affects soil structure and pesticide migration, which both may lead to land degradation. The effect on soil structure has been mainly addressed by classical methods at the macroaggregate scale, and the effect on herbicide leaching has not been thoroughly addressed. We aimed to characterize the effects of wetting and drying on soil microstructure, aggregate packing and stability, and subsequent effect on pesticide mobility in three agricultural soils. We developed advanced methods to quantitatively describe soil microstructure changes, induced by wetting and drying. Changes in soil packing, observed by micro-CT, indicate that large aggregates in a Clay soil disintegrate, whereas particles from a Sandy Loam form larger aggregates. To reflect these changes in terms of soil stability, we developed an aggregate durability index, based on changes in soil particle-size distribution measured by laser granulometry. For a Clay soil, the index decreased from 17.0% to 1.7%, indicating soil disaggregation upon wetting and drying. Whereas for a Sandy Clay Loam soil, the index increased from 0.4% to 2.6%, indicating formation of more durable aggregates, explained in terms of cementation by CaCO3. As expected, wetting and drying did not alter a Loamy Sand soil structure. The adverse effects of wetting and drying on soil structure correspond with the trends of atrazine mobility. As atrazine is trapped within the Clay soil aggregates, disaggregation leads to a 35% enhancement in pesticide mobility, whereas stabilization of a Sandy Clay Loam aggregates reduced atrazine leaching by 23%. Finally, wetting and drying directly affects the soil microstructure, which has an immense indirect effect on pollutant mobility, with both potentially leading to land degradation.

Polycation conformation upon adsorption to a negatively charged surface can be modulated by its charge density. At high charge density monomers directly interact with the surface in a ‘trains’ conformation and as charge density decreases a high degree of monomers dangle into solution in a ‘loops and tails’ conformation. In this study, the conformations of polycations upon adsorption to montmorillonite, as a function of polycation charge (20, 50 and 100% of the monomers, denoted as P-Q20, P-Q50 and P-Q100) were characterized and in accordance with their conformation, the adsorption of non-ionic and anionic molecules by the composite was tested. As expected, the adsorption of the nonionic pharmaceuticals increased to a composite with a ‘loops and tails’ conformation, due to both conformation and chemical properties. On the other hand, the anionic molecules, gemfibrozil and diclofenac, preferably adsorbed to composites with higher charge density (Q50 or Q100 composites). However, they showed different tendency toward the composites, i.e. higher adsorption of diclofenac by Q100 composite vs. higher adsorption of gemfibrozil by Q50 composite. To elucidate the differences in adsorption between these two pharmaceuticals, density functional theory calculations were employed. Both gemfibrozil and diclofenac were found to be better stabilized by methyl pyridinium sites (prevail in Q100 composite). The preferable adsorption of gemfibrozil by Q50 composite was explained by the presence of ‘loops and tails’ conformation enabling additional adsorption sites and diverse monomer-target molecule orientations.

Use of reclaimed wastewater for agricultural irrigation is seen as an attractive option to meet agricultural water demands of a growing number of countries suffering from water scarcity. However, reclaimed wastewater contains pollutants which are introduced to the agro-environment during the irrigation process. While water reuse guidelines do consider selected classes of pollutants, they do not account for the presence of pollutants of emerging concern such as pharmaceuticals and the potential risks these may pose. Here we use source–pathway–receptor analysis (S–P–R) to develop a holistic framework for evaluating the impacts of pharmaceuticals, present in wastewater used for agricultural irrigation, on human and ecosystem health and evaluate the data availability for the framework components. The developed framework comprised of 34 processes and compartments but a good level of knowledge was available for only five of these suggesting that currently it is not possible to fully establish the impacts of pharmaceuticals in wastewater irrigation systems. To address this, work is urgently needed to understand the fate and transport of pharmaceuticals in arable soil systems and the effects of chronic low-level exposure to these substances on microbes, invertebrates, plants, wildlife and humans. In addition, research pertaining to the fate, uptake and effects of pharmaceutical mixtures and metabolites is lacking as well as data on bio-accessibility of pharmaceuticals after ingestion. Scientific advancements in the five areas prioritised in terms of future research are needed before we are able to fully quantify the agricultural and human health risks associated with reclaimed wastewater use.