Quantifying NO_x emissions in Egypt using TROPOMI observations

Urban areas and industrial facilities, which concentrate most human activity and industrial production, are major sources of air pollutants, with serious implications for human health and global climate. For most of these pollutants, emission inventories are often highly uncertain, especially in developing countries. Spaceborne observations from the TROPOMI instrument, onboard the Sentinel-5 Precursor satellite, are used to measure nitrogen dioxide (NO2) slant column densities with a high spatial resolution. Here, we use two years of TROPOMI retrievals to map nitrogen oxides (NOx = NO + NO2) emissions in Egypt with a top-down model based on the continuity equation in steady state. Emissions are expressed as the sum of a transport term and a sink term representing the three-body reaction comprising NO2 and OH. This sink term requires information on the lifetime of NO2, which is calculated with the use of CAMS near-real-time temperature and hydroxyl radical (OH) concentration fields. The applicability of the OH concentration field is evaluated by comparing the lifetime it provides with the lifetime inferred from the fitting of NO2 line density profiles with an exponentially modified Gaussian function. This comparison, which is conducted for 39 samples of NO2 patterns above the city of Riyadh, provides information on the reliability of the CAMS near-real-time OH concentration fields; It also provides the location of the most appropriate vertical level to represent typical pollution sources in industrial areas and megacities in the Middle East. In Egypt, total derived emissions of NOx are dominated by the sink term. However, they can be locally dominated by wind transport, especially along the Nile where human activities are concentrated. Megacities and industrial regions clearly appear as the largest sources of NOx emissions in the country. Our top-down model produces emissions whose annual variability is consistent with the national electricity consumption. It is also able to detect lower emissions on Fridays, which are inherent to the social norm of the country, and to quantify the drop in emissions due to the COVID-19 pandemic. Overall, our indications of NOx emissions for Egypt are found to be 25.0 % higher than the CAMS-GLOB-ANT_v4.2 inventory, but significantly differ in terms of seasonality.

Optimizing Water Consumption in Richards’ Equation Framework with Step-Wise Root Water Uptake: A Simplified Model

This work arises from the need of exploring new features for modeling and optimizing water consumption in irrigation processes. In particular, we focus on water flow model in unsaturated soils, accounting also for a root water uptake term, which is assumed to be discontinuos in the state variable. We investigate the possibility of accomplishing such optimization by computing the steady solutions of a $$\theta$$ θ -based Richards equation revised as equilibrium points of the ODEs system resulting from a numerical semi-dicretization in the space; after such semi-discretization, these equilibrium points are computed exactly as the solutions of a linear system of algebraic equations: the case in which the equilibrium lies on the threshold for the uptake term is of particular interest, since the system considerably simplifies. In this framework, the problem of minimizing the water waste below the root level is investigated. Numerical simulations are provided for representing the obtained results. Article Highlights Root water uptake is modelled in a Richards’ equation framework with a discontinuous sink term. After a proper semidiscretization in space, equilibrium points of the resulting nonlinear ODE system are computed exactly. The proposed approach simplifies a control problem for optimizing water consumption.

Momentum budget of low level westerly jet during monsoon onset

In this paper a diagnostic study is carried Out to the source and sink terms for the formation and acceleration/deceleration of low level westerly Jet during monsoon onset. For this purpose momentum budget technique is used. The budget equation is derived in the (x.y.p.t) system. The area is confined to a small box the boundaries over the central Arabian Sea the westerly flow prominent during the onset of monsoon. Each tem in the budget equation is calculated separately. FGG E III b. 1200 UTC data set is used for the analysis. The Coriolis force term is found to be sink term rather than a source tern Tram. Transient north-south term is prominent source terms when time averaged momentum budget is considered. When the time averaged momentum budget for each pressure slab is considered. it is concluded that, north-south force terms are prominent source terms of momentum for all slabs and large .cumulus type convection may contribute to frictional dissipation of momentum for the upper pressure slabs. Frictional force is the main sink term when one examines the momentum budget for each day and the source term are varying day-by-day. On the average, the jet is accelerated during the period. It is also found that the net momentum tendency is small and oscillatory in nature. it is also found that at Minicoy. rainfall is inversely related to momentum tendency and whenever westerly jet is strong (weak) the rainfall is less (more). Distribution of U-momentum is also found to be oscillatory ill nature.

Analytical Solution of 1D Advection-Dispersion Equation In Finite Groundwater Reservoir With Spatially Dependent Dispersivity And Two Inputs Sources With Source-Sink Term Impact

This study derives an analytical solution of a one-dimensional (1D) advection-dispersion equation (ADE) for solute transport with two contaminant sources that takes into account the source term. For a heterogeneous medium, groundwater velocity is considered as a linear function while the dispersion as a nth-power of linear function of space and analytical solutions are obtained for and . The solution in a heterogeneous finite domain with unsteady coefficients is obtained using the Generalized Integral Transform Technique (GITT) with a new regular Sturm-Liouville Problem (SLP). The solutions are validated with the numerical solutions obtained using MATLAB pedpe solver and the existing solution from the proposed solutions. We exanimated the influence of the source term, the heterogeneity parameters and the unsteady coefficient on the solute concentration distribution. The results show that the source term produces a solute build-up while the heterogeneity level decreases the concentration level in the medium. As an illustration, model predictions are used to estimate the time histories of the radiological doses of uranium at different distances from the sources boundary in order to understand the potential radiological impact on the general public.

Representation of the Autoconversion from Cloud to Rain Using a Weighted Ensemble Approach

Cloud and precipitation processes remain among the largest sources of uncertainties in weather and climate modeling, and considerable attention has been paid to improve the representation of the cloud and precipitation processes in numerical models in the last several decades. In this study, we develop a weighted ensemble (named as EN) scheme by employing several widely used autoconversion (ATC) schemes to represent the ATC from cloud water to rainwater. One unique feature of the EN approach is that ATC rate is a weighted mean value based on the calculations from several ATC schemes within a microphysics scheme with a negligible increase of computation cost. The EN scheme is compared with the several commonly used ATC schemes by performing a real case simulations. In terms of accumulated rainfall and extreme hourly rainfall rate, the EN scheme provides better simulations than that by using the single Berry-Reinhardt scheme which was originally used in the Thompson scheme. It is worth emphasizing, in the present study, we only pay our attention to the ATC process from cloud water into rainwater with the purpose to improve the modeling of the extreme rainfall events over southern China. Actually, any (source/sink) term in a cloud microphysics scheme can be dealt with the same approach. The ensemble method proposed herein appears to have important implications for developing cloud microphysics schemes in numerical models, especially for the models with variable grid resolution, which would be expected to improve of the representation of cloud microphysical processes in the weather and climate models.

Evaluation of NOAA National Water Model Parameter Calibration in Semi-Arid Environments Prone to Channel Infiltration

The NOAA National Water Model (NWM), maintained and executed by the NOAA National Weather Service (NWS) Office of Water Prediction, provides operational hydrological guidance throughout the Contiguous United States. Based on the WRF-Hydro model architecture developed by the National Center for Atmospheric Research (NCAR), the NWM was recently modified for semi-arid domains, by permitting it to explicitly resolve infiltration from ephemeral channels into the underlying channel bed as an added model sink term. To analyze the added value of channel infiltration in semi-arid environments, we calibrated NWM v2.1 (with the channel infiltration function) to 56 independent basins in the western CONUS, following identical calibration methods as the pre-operational NWM v2.1 (not including channel infiltration). Calibration of the model consists of two parts, including 1) calibration of channel infiltration only with other parameters set to the calibrated parameters used for pre-operational NWM v2.1 and 2) calibration of all parameters including channel infiltration with settings otherwise equivalent to the calibration of NWM v2.1. The calibrated channel-infiltration enhanced NWM improves predictive skill compared to the control NWM in 85% of evaluated basins, for the calibration period. The current NWM settings for physical processes and the biases of the calibration scheme limit model performance in semi-arid environments. To explore whether channel infiltration paired with an alternative calibration scheme could address these limitations, NWM v2.1 was calibrated with a new objective function in selected basins. We found that this updated objective function could ameliorate model biases in some semi-arid environments.

Investigation of several proxies to estimate sulfuric acid concentration under volcanic plume conditions

Sulfuric acid (H2SO4) is commonly accepted as a key precursor for atmospheric new particle formation (NPF). However, direct measurements of [H2SO4] remain challenging, thereby preventing the determination of this important quantity, and, consequently, a complete understanding of its contribution to the NPF process. Several proxies have been developed to bridge the gaps, but their ability to predict [H2SO4] under very specific conditions, such as those encountered in volcanic plumes (including, in particular, high sulfur dioxide mixing ratios), has not been evaluated so far. In this context, the main objective of the present study was to develop new proxies for daytime [H2SO4] under volcanic plume conditions and compare their performance to that of the proxies available in the literature. Specifically, the data collected at Maïdo during the OCTAVE (Oxygenated organic Compounds in the Tropical Atmosphere: variability and atmosphere–biosphere Exchanges) 2018 campaign, in the volcanic eruption plume of the Piton de la Fournaise, were first used to derive seven proxies based on knowledge of the sulfur dioxide (SO2) mixing ratio, global radiation, condensation sink (CS) and relative humidity (RH). A specific combination of some or all of these variables was tested in each of the seven proxies. In three of them (F1–F3), all considered variables were given equal weight in the prediction of [H2SO4], whereas adjusted powers were allowed (and determined during the fitting procedure) for the different variables in the other four proxies (A1–A4). Overall, proxies A1–A4 were found to perform better than proxies F1–F3, with, in particular, improved predictive ability for [H2SO4] > 2 × 108 cm−3. The CS was observed to play an important role in regulating [H2SO4], whereas the inclusion of RH did not improve the predictions. A last expression accounting for an additional sink term related to cluster formation, S1, was also tested and showed a very good predictive ability over the whole range of measured [H2SO4]. In a second step, the newly developed proxies were further evaluated using airborne measurements performed in the passive degassing plume of Etna during the STRAP (Synergie Transdisciplinaire pour Répondre aux Aléas liés aux Panaches volcaniques) 2016 campaign. Increased correlations between observed and predicted [H2SO4] were obtained when the dependence of predicted [H2SO4] on the CS was the lowest and when the dependence on [SO2] was concurrently the highest. The best predictions were finally retrieved by the simple formulation of F2 (in which [SO2] and radiation alone were assumed to explain the variations in [H2SO4] with equal contributions), with a pre-factor adapted to the STRAP data. All in all, our results illustrate the fairly good capacity of the proxies available in the literature to describe [H2SO4] under volcanic plume conditions, but they concurrently highlight the benefit of the newly developed proxies for the prediction of the highest concentrations ([H2SO4] > 2–3 × 108 cm−3). Moreover, the contrasting behaviours of the new proxies in the two investigated datasets indicate that in volcanic plumes, like in other environments, the relevance of a proxy can be affected by changes in environmental conditions and that location-specific coefficients do logically improve the predictions.

The Cumulative Amount of Exuded Citrate Controls Its Efficiency to Mobilize Soil Phosphorus

Root exudation of citrate is discussed as mechanism to mobilize P from the soils' solid phase. Microbial processes can mitigate the mobilization efficiency of citrate. Due to higher microbial activity in topsoils compared to subsoils, we hypothesized a lower mobilization efficiency of exuded citrate in topsoils than in the subsoils. As a model system we used microdialysis (MD) probes and we followed diffusive fluxes of citrate from the perfusate into the soil and of phosphate from the soil into the dialysate in three soil horizons (Oa, Ah, Bw) of a Fagus sylvatica L. stand Cambisol. Three different MD perfusates with a KCl background concentration have been used: control, 1, and 3 mmol L−1 citric acid. Fluxes have been measured after 24, 48, and 144 h. The high-citrate perfusate increased the cumulative 144 h P-influx by a factor of 8, 13, and 113 in the Oa, Ah, and Bw horizon, respectively. With the high-citrate treatment, P mobilization efficiency decreased over time, whereas for the low citrate, P mobilization efficiency had a maximum at day 2. Minimum P mobilization efficiency of citrate was 1:25,000 mol phosphate per mol citrate in the Oa horizon between days 2 and 6, and maximum was 1:286 in the Bw-horizon during day 2. An increasing citrate efflux over time indicated an increasing sink term for citrate in the soil due to microbial decay or immobilization processes. Cumulative phosphate influx could be fitted to cumulative citrate efflux and soil horizon in a logarithmic model explaining 87% of the variability. For the first time, we could follow the localized P-uptake with citrate exudation over several days. Cumulative citrate efflux as the main control of P-mobilization has been barely discussed yet, however, it could explain some gaps in the role of carboxylates in the rhizosphere. Batch experiments are not capable to elucidate microscale dynamic competition for phosphate and carboxylates. MD is a promising tool for spatially explicit investigation of phosphate–citrate exchange, since such detailed insights in are not possible with batch experiments. In combination with the analysis of microbial properties, this technique has a huge potential to identify mobilization processes in soils as induced by citrate.

Integrodifference master equation describing actively growing blood vessels in angiogenesis

We study a system of particles in a two-dimensional geometry that move according to a reinforced random walk with transition probabilities dependent on the solutions of reaction-diffusion equations (RDEs) for the underlying fields. A birth process and a history-dependent killing process are also considered. This system models tumor-induced angiogenesis, the process of formation of blood vessels induced by a growth factor (GF) released by a tumor. Particles represent vessel tip cells, whose trajectories constitute the growing vessel network. New vessels appear and may fuse with existing ones during their evolution. Thus, the system is described by tracking the density of active tips, calculated as an ensemble average over many realizations of the stochastic process. Such density satisfies a novel discrete master equation with source and sink terms. The sink term is proportional to a space-dependent and suitably fitted killing coefficient. Results are illustrated studying two influential angiogenesis models.

Investigation of several proxies to estimate sulfuric acid concentration in volcanic plume conditions

Sulfuric acid (H2SO4) is commonly accepted as a key precursor for atmospheric new particle formation (NPF). However, direct measurements of [H2SO4] remain challenging, thus preventing the determination of this important quantity, and, consequently, a complete understanding of its contribution to the NPF process. Several proxies have been developed to bridge the gaps, but their ability to predict [H2SO4] in very specific conditions such as those encountered in volcanic plumes (including in particular high sulphur dioxide mixing ratios) has not been evaluated so far. In this context, the main objective of the present study was to develop new proxies for daytime [H2SO4] in volcanic plume conditions and compare their performance to that of the proxies available in the literature. In specific, the data collected at Maïdo during the OCTAVE 2018 campaign, in the volcanic eruption plume of the Piton de la Fournaise, were first used to derive seven proxies based on the knowledge of sulphur dioxide (SO2) mixing ratio, global radiation, condensation sink (CS) and relative humidity (RH). In three of the seven proxies (F1–F3), all variables were given equal weight in the prediction of [H2SO4], while adjusted powers were allowed for the different variables in the other four proxies (A1–A4). Proxies A1–A4 were overall found to perform better compared to F1–F3, with, in specific, improved predictive ability for [H2SO4] > 2 × 108 cm−3. The CS was observed to play an important role in regulating [H2SO4], while, in contrast, the inclusion of RH did not improve the predictions. A last expression accounting for an additional sink term related to cluster formation, S1, was also tested and showed a very good predictive ability over the whole range of measured [H2SO4]. The newly developed proxies were in a second step further evaluated using airborne measurements performed in the passive degassing plume of Etna during the STRAP 2016 campaign. Increased correlations between observed and predicted [H2SO4] were obtained when the dependence of predicted [H2SO4] over CS was the lowest, and when the dependence over [SO2] was concurrently the highest. The best predictions were finally retrieved by the simple formulation of F2 (in which [SO2] and radiation alone were assumed to explain the variations of [H2SO4] with equal contributions), with a pre factor adapted to the STRAP data. All in all, our results illustrate the fairly good capacity of the proxy available in the literature to describe [H2SO4] in volcanic plume conditions, but highlight at the same time the benefit of the newly developed proxies for the prediction of the highest concentrations ([H2SO4] > 2–3 × 108 cm−3). Also, the contrasting behaviours of the new proxies in the two investigated datasets indicate that in volcanic plumes like in other environments, the relevance of a proxy can be affected by changes in environmental conditions, and that location specific coefficients do logically improve the predictions.