The present study deals with groundwater pollution in multilayer aquifer. The model is based on decomposition of finite layers in semi-infinite groundwater reservoir. A constant pollutant source is injected at the input boundary of the uppermost layer (UML) of the landfill. At the intermediate inlet boundary, some average value for the longitudinal exchange of the input source concentration in each sub-layer is considered from the previous layer. Initially, the aquifer is not solute free in each sub layer that means some constant background contaminant concentration exists. In each sub layer, concentration gradient is assumed to be zero at the extreme boundary. The linear sorption and first orders decay terms are considered to model the groundwater pollution in multilayer aquifer. The Laplace transform technique is adopted to solve one-dimensional (1D) advection-dispersion equation (ADE). This approach is helpful to understand the solute migration in finite sub layers. The results are elucidated for the different time periods to examine the peak of pollutant concentration level in geological formations.

Vertical measurements of ozone (O₃) within the 3000-m lower troposphere were obtained using an O₃ lidar to investigate the contribution of the interactions between the transport and boundary layer processes to the surface O₃ levels in urban Shanghai, China during July 23–28, 2017. An extremely severe pollution episode with a maximum hourly O₃ mixing ratio of 160.4 ppb was observed. In addition to enhanced local photochemical production, both downward and advection transport in the lower troposphere may have played important roles in forming the pollution episode. The O₃-rich air masses in the lower free troposphere primarily originated from central China and the northern Yangtze River Delta (YRD) region. The downward transport of O₃ from the lower free troposphere may have an average contribution of up to 49.1% to the daytime (09:00–16:00 local time) surface O₃ in urban Shanghai during the pollution episode (July 23–26, 2017). As for the advection transport, large amounts of O₃ were transported outward from Shanghai in the planetary boundary layer under the influence of southeasterly winds during the field study. In this condition, the boundary-layer O₃ that was transported downward from the free troposphere in Shanghai could be transported back to the northern YRD region and accumulated therein, leading to the occurrence of severe O₃ pollution events over the whole YRD region. Our results indicate that effective regional emission control measures are urgently required to mitigate O₃ pollution in the YRD region.

In this study, we investigated the daily variability of PM₁₀ concentrations in January in Korea during the past 19 years (2001–2019), as well as the associated atmospheric circulation patterns. The daily PM₁₀ concentrations were classified into three cases: low (L; < 50 μg/m³), high (H; 50–100 μg/m³), and extremely high (EH; ≥ 100 μg/m³). We found that the strength of the East Asian winter monsoon influenced the PM₁₀ variability in the L and H cases. However, the EH cases were strongly influenced by the rapid growth of barotropic warming (anticyclonic anomaly) over the eastern North Atlantic and Northern Europe (ENE), and the stationary Rossby waves grew rapidly over Eurasia within only four days. Analysis of the quasi-geostrophic geopotential tendency budget revealed that the anticyclonic anomaly over the ENE was enhanced by vorticity advection. Linear baroclinic model experiments confirmed that vorticity forcing over the ENE induces favorable atmospheric conditions for the occurrence of EH PM₁₀ events in East Asia. As a result, the PM₁₀ concentration sharply increased sharply by approximately three times over four days. This study suggests that understanding atmospheric teleconnections between the ENE and East Asia can effectively predict the occurrence of EH PM₁₀ events in Korea, helping to reduce the human health risks from atmospheric pollution.

Oil spills occurring either in oceans or inland waterways may cause serious economic losses and ecological damage. Previous studies pertaining to oil spills and their consequences are primarily based on marine environments, whereas few have focused on oil spills occurring in inland waterways characterised by pronounced flow advection transport effects, which differ from the marine environment. A generalised flume experiment is performed to investigate the spread and transport of oil spills, and the relationships between the area and thickness of oil slick over time are analysed parametrically. An oil spill model combined with a depth–integrated two–dimensional non–uniform flow model, which is suitable for modelling inland waterways based on the Lagrangian method, is established; it is calibrated and verified using measured data from the flume experiment. The model is applied to three scenarios on the Luoqi reach of the Yangtze River, and spilled oil drifting trajectory maps are obtained and analysed considering the field wind parameters. The results show that the drift distance of the oil slick in the inland waterway is primarily controlled by the flow velocity with effects of advection transport; however, the oil spill trajectory spreads toward the wind direction when the flow velocity is relatively small compared with the wind speed. The results of this study serve as a reference for predicting the spread and transport of oil spills in inland waterways.

A Lagrangian model which describes radionuclide transport in the northern Indian Ocean is described. Water circulation is obtained from HYCOM ocean model for year 2017. The model includes advection by currents, turbulent mixing and radionuclide interactions between water and sediments, described in a dynamic way using kinetic transfer coefficients. Hypothetical releases from five coastal nuclear power plants operating in the northern Indian Ocean were simulated. Releases were supposed to start both during the winter and summer monsoons, to study reversing circulation effects. Age distributions of releases were calculated, which adds information about circulation and radionuclide pathways. It was found that, for some of the NPPs, radionuclide distributions resulting from releases starting in both seasons were not as different as could be expected from the opposed circulation schemes during each monsoon. Effective ¹³⁷Cs half-lives in the ocean surface were calculated and results were two orders of magnitude below previous estimations.

Low-permeability media such as clay appear in nearly all hydrogeological systems. To date, although significant efforts have been put forward by hydrologists, transport mechanism is still not well understood in such media, especially in an aquifer-aquitard system. In this study, two-dimensional experiments of groundwater flow and solute transport were conducted in a clay-sand two-layer system to investigate the characteristics of flow and transport in such a system. Sodium chloride (NaCl) (a conservative tracer) from a tank was injected after passing by the pre-inlet reservoir where the mixing effect and flow transiency were analyzed. A new numerical model considering the mixing effect and flow transiency was developed to interpret the experimental data based on the finite-element COMSOL Multiphysics platform. Transport parameters were assessed by best fitting the observed breakthrough curves (BTCs). Several important results were obtained. Firstly, aquitard advection was found to be non-negligible and should be considered in a proper mathematical model for describing the transport process. Secondly, advective velocities were temporally variable and showed decreasing trends in the sand and clay layers, mainly due to the impacts of physical and biological clogging. Thirdly, the mixing effect in the pre-inlet reservoir led to a lower tracer concentration in the sand layer at early times. Finally, the observed BTCs exhibited early arrivals in the clay layer, possibly resulting from preferential flow pathways. These findings can provide hints for contamination remediation works in aquifer-aquitard systems.

Scientific evidence suggests that Saharan dust intrusions in Southern Europe contribute to the worsening of multiple pathologies and increase the concentrations of particulate matter (PM) and other pollutants. However, few studies have examined whether Saharan dust intrusions influence the incidence and severity of COVID-19 cases. To address this question, in this study we carried out generalized linear models with Poisson link between incidence rates and daily hospital admissions and average daily concentrations of PM₁₀, NO₂, and O₃ in nine Spanish regions for the period from February 1, 2020 to December 31, 2020. The models were adjusted by maximum daily temperature and average daily absolute humidity. Furthermore, we controlled for trend, seasonality, and the autoregressive nature of the series. The variable relating to Saharan dust intrusions was introduced using a dichotomous variable, NAF, averaged across daily lags in ranges of 0–7 days, 8–14 days, 14–21 days, and 22–28 days. The results obtained in this study suggest that chemical air pollutants, and especially NO₂, are related to the incidence and severity of COVID-19 in Spain. Furthermore, Saharan dust intrusions have an additional effect beyond what is attributable to the variation in air pollution; they are related, in different lags, to both the incidence and hospital admissions rates for COVID-19. These results serve to support public health measures that minimize population exposure on days with particulate matter advection from the Sahara.

Meaningful assessment of pesticide fate in soils and plants is based on fate models that represent all relevant processes. With mechanistic models, these processes can be simulated based on soil, substance, and plant properties. We present a mechanistic model that simulates pesticide uptake from soil and investigate how it is influenced, depending on the governing uptake process, by root and substance properties and by distributions of the substance and water in the soil profile. A new root solute uptake model based on a lumped version of the Trapp model (Trapp, 2000) was implemented in a coupled version of R-SWMS-ParTrace models for 3-D water flow and solute transport in soil and root systems. Solute uptake was modeled as two individual processes: advection with the transpiration stream and diffusion through the root membrane. We set up the model for a FOCUS scenario used in the European Union (EU) for pesticide registration. Considering a single vertical root and advective uptake only, the root hydraulic properties could be defined so that water and substance uptake and substance fate in soil showed a good agreement with the results of the 1D PEARL model, one of the reference models used in the EU for pesticide registration. Simulations with a complex root system and using root hydraulic parameters reported in the literature predicted larger water uptake from the upper root zone, leading to larger pesticide uptake when pesticides are concentrated in the upper root zone. Dilution of root water concentrations at the top root zone with water with low pesticide concentration taken up from the bottom of the root zone leads to larger uptake of solute when uptake was simulated as a diffusive process. This illustrates the importance of modeling uptake mechanistically and considering root and solute physical and chemical properties, especially when root-zone pesticide concentrations are non-uniform.

To overcome the weaknesses of traditional landfills, a modified aerobic landfill concept with intermediate covers of coarse material between waste layers functioning as facilities of drainage and aeration has been proposed recently. In this study, a one-dimensional coupled model, including aerobic biodegradation, oxygen diffusion, and advection, is proposed to describe oxygen distribution in this new type of landfill. Homotopy analysis method and perturbation method are applied to solve this model at passive aeration and active aeration, respectively. The model has six input variables, that is, oxygen diffusion coefficient, gas permeability, maximum oxygen consumption rate, layer thickness of waste, and injection pressure and extraction pressure. A combination of their typical values gives rise to over 700,000 scenarios which can be calculated by the proposed solution. The coupled effect of the above variables on oxygen migration is quantitatively investigated, followed by an estimation formula of the minimum oxygen concentration in waste layer. The maximum waste layer thickness is defined as a function of other variables for a given aeration target of oxygen volume concentration larger than 5%. A generalized design method of waste layer thickness, injection pressure, and extraction pressure is then developed for the newly proposed modified layered aerobic landfill, which can promote its popularization and application.