Large Spatial Variability Research Articles

Site-specific effects of fertilizer on hay and grain yields of oats: evidence from large-scale field experiments.

Oat (Avena sativa L.) is a valuable crop due to its strong adaptability to marginal environments, making it an important component of agricultural systems in regions where other cereals may not thrive. The application of chemical fertilizer can influence oat hay and grain yield significantly. However, large-scale meta-analytical studies of the size and variability of oat hay and grain yields in response to fertilizer addition are still lacking. Based on 83 studies worldwide, this meta-analysis quantifies the impact of the addition of fertilizer on oat hay and grain yields under varying environmental conditions (e.g., soil nutrient levels, texture, and climate). The results confirmed that the fertilizer application increased oat hay yield by 48.9% and grain yield by 36.2%. This study demonstrated that balanced fertilization with nitrogen, phosphorus, and potassium generally enhances oat hay and grain yield despite large temporal and spatial variations. Boosted regression tree (BRT) models suggest that changes in hay and grain yield were primarily dominated by soil pH and nitrogen fertilizer. The response ratio (the natural logarithm of the mean values of hay yield or grain yield with and without fertilization, respectively) of hay yield declined linearly with soil pH. Elevation was the second most important factor affecting the change in response ratio of hay yield and the third most important factor affecting the change in response ratio of grain yield but climatic conditions were not the dominant factors affecting changes in oat hay or grain yield. Overall, these results will benefit producers considering site-specific fertilization management of oat. They could increase yields and save investment in fertilizer, and help to facilitate the genetic breeding of oat varieties with high nutrient use efficiency. © 2024 The Author(s). Journal of the Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Bioaccumulation and tissue distribution of pharmaceuticals and their transformation products in fish along the pollution gradients of a wastewater-impacted river

A field study on the occurrence and distribution of forty-three pharmaceutically active compounds (PhACs) in water and fish samples from anthropogenically impacted section of the Sava River (Croatia) was performed to estimate the importance of bioaccumulation for the environmental risk assessment of PhACs. The study was performed using a highly specific LC-MS/MS method, tailored to include the most prominent PhACs from different therapeutic categories as well as their major metabolites and/or transformation products (TPs). The results revealed a widespread occurrence of PhAC residues both in water and fish samples with a large spatial variability reflecting the distance from the dominant wastewater discharges. The most prominent PhAC categories in less polluted upstream part of the river were common psychostimulants caffeine and cotinine, therapeutic opioids and cardiovascular drugs, while in the river section affected by the local municipal and industrial wastewater inputs, antibiotic drugs became clearly predominant, especially in fish tissue samples. The apparent bioconcentration factors (BCFs) of investigated PhACs varied over several orders of magnitude, from 0.02 ± 0.01 Lkg−1 for O-desmethyl tramadol in fish muscle to 784 ± 260 Lkg−1 for terbinafine in fish liver, indicating rather large differences in their bioconcentration potential and affinity to different tissues, with the tissue-specific BCFs increasing in the following order: muscle < gills < gonads < heart < liver < kidneys. The bioconcentration potential of most of the PhACs included in this study was only low to moderate however moderately high BCFs of certain PhACs (e.g. sertraline, terbinafine, loratadine, diazepam and azithromycin) in some tissues should be taken into consideration when assessing their potential environmental risks. Moreover, it was shown that BCFs could be strongly affected by biotransformation in fish. Risk prioritization based on risk quotient (RQ) and ToxPi index, revealed antibiotics, in particular azithromycin, and therapeutic psychoactive substances as the most hazardous pharmaceutical contaminants in the Sava River.

Characterization of soil organic carbon at profile scale in two forest soils under pine and holm oak

Context It is essential to promote soil carbon sequestration as a means to mitigate climate change. Thus, it is fundamental to know the distribution of C in the soil at profile scale, the characteristics of humic substances as indicators of soil organic matter (SOM) turnover and their relationship with other soil properties. Aims Two forest ecosystems under pine and holm oak developed under a Meso-Mediterranean climate in Spain were used to characterize SOM through the complete sequence of layers of the soil profile. Methods General soil analysis, infrared spectroscopic analysis and soil color measurements were conducted for the characterization. Key results Humus form under oak was found to be Mull mesotrophic–Mull acid while humus under pine Moder oligotrophic. The infrared spectrum determined that oxidation of the humic acids was more complete in the deeper horizons. Relationships between intensities of the main spectral bands in both soils followed similar maximum and minimum sequence values. The total humic extract (THE) color measured by reflection was found inversely related to the THE color measured by transmission. In the same way, the color spectrum between 350 and 800 nm in the THE showed an inverse relation between hue/brightness and absorbance values. Conclusions Infrared analysis and color measurements provided evidence of a different level of stabilization of humic substances from each soil, and between the different horizons. Large spatial variability in soil organic carbon quantity and quality was observed. Implications Better understanding of carbon sequestration behavior in different soil ecosystems in its crucial role within the global carbon cycle

Damage Analysis of High-Speed Railway Bridge-Track System with Unequal Pier Heights Under Near-Fault Earthquake Based on FK-FE Method

ABSTRACT Canyons have a significant influence on near-fault earthquakes, and near-fault earthquakes exhibit strong amplitude and large spatial variability. In turn, this causes the damage and deformation of the high-speed railway (HSR) bridge-track system with unequal pier heights, which affects the normal service of the unequal pier heights bridge and the safe running of the train. This study investigates the damage and deformation of HSR bridge-track systems with unequal pier height in canyon areas under near-fault seismic based on FK (frequency wavenumber)-FE (finite element) method. It aims to offer suggestions for seismic design in canyon regions of HSR unequal pier height bridges. First, a finite fault kinematic hybrid source model and a canyon site model were established, and the FK-FE hybrid method was employed to simulate the three-dimensional (3D) site motion considering the full process of the seismic source-propagation medium-complex site (SSPC). Then, the influence of near-fault earthquakes and the canyon topography on the spatial distribution of near-fault ground motions were analyzed, and the non-uniform near-fault seismic excitation at the bottom of each pier considering the canyon topography effect was obtained. Finally, a nonlinear dynamic FE model of the HSR bridge-track system was established, and its damage patterns and mechanism in canyon near-fault earthquakes were analyzed. The results indicate that near-fault earthquakes in the canyon exhibit significant vertical and sliding effects with significant spatial variability and the FK-FE method can provide more accurate ground motion input. The seismic damage to HSR bridge-track systems increases as the fault distance decreases, and topography effects further aggravates seismic damage to HSR bridges. Moreover, vulnerable areas in the track system are often found at girder ends. Therefore, for HSR bridges in canyon regions, seismic input needs to comprehensively consider the effects of near-fault earthquakes and topography, with post-earthquake focus on damage assessment at girder ends.

Multi-objective calibration and evaluation of the ORCHIDEE land surface model over France at high resolution

Abstract. Here we present a strategy to obtain a reliable hydrological simulation over France with the ORCHIDEE land surface model. The model is forced by the SAFRAN atmospheric reanalysis at 8 km resolution and hourly time steps from 1959 to 2020 and by a high-resolution DEM (around 1.3 km in France). Each SAFRAN grid cell is decomposed into a graph of hydrological transfer units (HTUs) based on the higher-resolution DEM to better describe lateral water movements. In particular, it is possible to accurately locate 3507 stations among the 4081 stations collected from the national hydrometric network HydroPortail (filtered to drain an upstream area larger than 64 km2). A simple trial-and-error calibration is conducted by modifying selected parameters of ORCHIDEE to reduce the biases of the simulated water budget compared to the evapotranspiration products (the GLEAM and FLUXCOM datasets) and the HydroPortail observations of river discharge. The simulation that is eventually preferred is extensively assessed with classic goodness-of-fit indicators complemented by trend analysis at 1785 stations (filtered to have records for at least 8 entire years) across France. For example, the median bias of evapotranspiration is −0.5 % against GLEAM (−4.3 % against FLUXCOM), the median bias of river discharge is 6.3 %, and the median Kling–Gupta efficiency (KGE) of square-rooted river discharge is 0.59. These indicators, however, exhibit a large spatial variability, with poor performance in the Alps and the Seine sedimentary basin. The spatial contrasts and temporal trends of river discharge across France are well represented with an accuracy of 76.4 % for the trend sign and an accuracy of 62.7 % for the trend significance. Although it does not yet integrate human impacts on river basins, the selected parameterization of ORCHIDEE offers a reliable historical overview of water resources and a robust configuration for climate change impact analysis at the nationwide scale of France.

Forest structure explains spatial heterogeneity of decadal carbon dynamics in a cool-temperate forest

Abstract Accurate evaluation of forest biomass distribution and its long-term change over wide areas is required for effective forest carbon management and prediction of landscape-scale forest dynamics. We evaluated a landscape-scale (225 km2) decadal forest carbon budget at a 1 ha spatial resolution in a cool-temperate forest, by repeating airborne laser observations 10 years apart and partitioning net forest biomass change (FBC) into growth and mortality. Using &gt;10 000 samples, we revealed that naturally regenerated forests have large spatial heterogeneity in net biomass change, and 3/4 of the photosynthetically acquired carbon stock moved to necromass even without anthropogenic disturbances. Actual carbon residence time as living tree biomass was estimated by dividing biomass by growth or mortality rates. The residence time was 107 and 106 years, respectively with large spatial variation among stands (48 and 42 years, respectively, as the difference between 25 and 75 percentile), although studied forest stands have small variation in the forest functional type in a landscape-scale. The best predictors of subsequent decadal biomass changes were two forest structural factors, mean canopy height and canopy height variation in addition to one environmental factor, elevation. Considering the long lifetime of trees, these structural factors may be an indicator of forest soundness rather than a cause of forest growth or mortality. However, in any cases, these structural factors can be powerful predictors of subsequent FBC.

A Field‐Based Estimation of the Variability of Particle Entrainment in Coarse‐Bed Rivers

AbstractThe determination of critical shear stresses is fundamental to bedload sediment transport prediction in gravel‐bed rivers. Due to the heterogeneous shape and arrangement of the individual clasts in a riverbed, critical shear stresses typically show a large spatial variability, which is not adequately captured by the reach‐averaged description followed in common studies. In this regard, there is a general paucity of field data on this spatial variability of the critical shear stress, largely due to the lack of a standardized measurement method. In an attempt to fill this gap, we propose a field‐based workflow to estimate the frequency distribution of dimensionless critical shear stress (also named critical Shields number), which is based on the measurement of a series of variables related to the position, orientation and resistance to motion of individual clasts in a gravel‐bed river, combined with a probabilistic approximation to drag and lift coefficients. Following this workflow, the patch‐scale variability of particle incipient‐motion conditions was determined in a gravel bar of the Upper Cinca River, Spain. The results are consistent with what is known about sediment entrainment in gravel‐bed rivers. We consider this method to have great potential to advance our understanding of particle initiation of motion in gravel‐bed rivers as it provides valuable systematic field information.

Brightening and Fading in the Youngest Galactic Supernova Remnant G1.9+0.3: 13 Years of Monitoring with the Chandra X-Ray Observatory

We report 13 years of Chandra monitoring of the youngest Galactic supernova remnant G1.9+0.3, the only remnant known to be increasing in brightness. We confirm the spatially integrated rate of flux increase of (1.2 ± 0.2)% yr−1 (1–7 keV), but find large spatial variations, from −3% yr−1 to +7% yr−1, over length scales as small as 10″ or smaller. We observe relatively little change in spectral slope, though one region shows significant hardening as it brightens by 1% yr−1. Such rates of change can be accommodated by any of several explanations, including steady blast-wave evolution, expansion or compression of discrete plasma blobs, magnetic turbulence, or variations in magnetic-field aspect angle. Our results do not constrain the mean magnetic-field strength, but a self-consistent picture can be produced in which the maximum particle energies are limited by the remnant age (applying both to electrons and to ions) to about 20 TeV, and the remnant-averaged magnetic-field strength is about 30 μG. The deceleration parameter m (average shock radius varying as t m ) is about 0.7, consistent with estimates from overall expansion dynamics and confirming an explosion date of about 1900 CE. Shock-efficiency factors ϵ e and ϵ B (fractions of shock energy in relativistic electrons and magnetic field) are 0.003 and 0.0002 in this picture. However, the large range of rates of brightness change indicates that such a global model is oversimplified. Temporal variations of photon index, expected to be small but measurable with longer time baselines, can discriminate among possible models.

Assessment of Essential Trace Elements Cu, Zn, Mn, and Ni in Mediterranean Mussels (Mytilus galloprovincialis) from the Algerian Coast.

The aim of this study was to determine the levels of essential trace elements copper, zinc, manganese, and nickel in two size groups of mussels, Mytilus galloprovincialis, sampled from three sites along the Algerian coast with different levels of pollution. The results of this work are very interesting for determining the effect of mussel size on the variation of contaminant levels in environmental monitoring studies. Thirty individuals from 2 size groups, less than 4-cm and more than 4-cm shell length, were collected over 4 seasons at each site studied. The samples were dissected, and copper, zinc, manganese, and nickel were measured in the whole soft tissues by inductively coupled plasma mass spectrometry (ICP-MS). The maximum concentrations of copper, zinc, manganese, and nickel expressed in mg/kg dry weight at the three sites were 54.01, 146.25, 20.69, and 8.85, respectively. A clear distribution pattern of the studied elements was observed in the two size classes. In fact, significantly higher concentrations (P < 0.05) were found in individuals with a shell length of less than 4 cm. Significantly, higher concentrations (P < 0.05) were found in individuals with a shell length of less than 4 cm due to their high filtration and low excretion activities compared to long individuals. In addition, metal pollution index (MPI) values were higher in small bivalve individuals. The trace element contents found in the mussel samples of the present study showed large spatial and temporal variations among the different sites studied, highlighting the typical exposure of marine sites to different marine pollutants. Notably, elevated concentrations of copper were observed in the port of Zemmouri, likely due to port activities. Conversely, the highest levels of manganese and nickel were recorded at Boudouaou El Bahri, potentially attributable to the use of fertilisers by agricultural activities. Some shellfish samples had exceeded the permissible concentration levels of copper and manganese, thereby posing a danger to bivalve consumers.

Projection of precipitation under RCP4.5 and RCP 8.5 in central and southern regions of Iraq

Climate change has significant impacts on natural systems, especially in terms of altering hydrological patterns due to changing precipitation and melting snow and ice. This study examines projected changes in precipitation for central and southern Iraq under the RCP4.5 and RCP8.5 scenarios for the periods 2046–2065 and 2081–2100 using CMIP5 climate model HadGEM2-AO output. The results revealed that under RCP4.5, there was large spatial variation in the projected precipitation over the region during 2064-2065, ranging from as low as 27 mm in western part to as high as 1541 mm in the eastern part. The spatial variability as well as precipitation amount decreased considerably during 2081-2100 periods. Under RCP8.5 the projected precipitation was only 14-164 mm in 2046-2065 and 36-532 mm in 2081-2100 periods. Thus, under RCP8.5 scenario anticipated precipitation will be quite low.

Impact of cloud properties and aerosols on the spatial distribution of Indian summer monsoon rainfall: A case study

Quantitative seasonal Indian summer monsoon (ISM) rainfall forecasting has been a challenge for the climate models as it is sensitive to cloud properties, which are modulated by large-scale dynamics. Proper understanding of the role of dynamics, aerosols, and cloud microphysical processes is important for the precipitation calculation in dynamical models. The all-India averaged seasonal ISM rainfall (ISMR) from June to September had been normal or above normal, however there were large spatial variations. Some meteorological subdivisions, the majority of which lie in the Gangetic Plains, received deficient seasonal rainfall during recent years 2022 and 2010. The quantitative seasonal prediction of the spatial distribution of ISMR from global coupled models (GCMs) has also shown a hard time for accurate spatial rainfall distribution. Four monsoon years (2022, 2021, 2011, and 2010) are considered, in which the country as a whole (all India land points) was ‘above normal’ or ‘normal’. However, in two specific years (2022 and 2010), some subdivisions received deficient rainfall (17% and 15% of the total area). We have used several datasets from the ground, and satellite observations, along with reanalysis products and computed climatology and anomalies for monsoon years of 2022 and 2010. The low-level wind and vertical velocity along with cloud microphysics played a pivotal role in the heterogeneous distribution of ISMR during 2022 and 2010 as compared to years 2021 and 2011. Here, we have demonstrated that the spatial heterogeneity of the positive and negative precipitation anomalies was consistent with cloud properties. Aerosols might have an impact on the cloud microphysical processes and modulate cloud properties and rainfall over different regions due to the availability of water vapor. Therefore, the inclusion of aerosol indirect effects in climate models is crucial for better interactions among dynamics, aerosol, and cloud microphysics. Proper representation of cloud droplet size distribution and microphysical conversion rates might be helpful for quantitative seasonal ISMR forecasting at a smaller spatial scale.

Separating the impact of climate change and human activities on the connection between meteorological and hydrological drought

AbstractUnderstanding the transition from meteorological to hydrological drought is essential for accurately predicting hydrological droughts. However, the factors driving this transition are intricate, and a comprehensive understanding of how direct human activities influence this shift in drought is lacking. In this study, we initially explored the spatiotemporal correlation between the occurrence of meteorological and hydrological droughts. Subsequently, we formulated multiple hydrological replenishment scenarios using the soil and water assessment tool (SWAT) model to assess the environmental impact of the transition from meteorological to hydrological droughts. The Xijiang River Basin (XRB), the primary tributary of the Pearl River basin, was selected as the study area. Our results identified 91 meteorological droughts in the XRB, and only 66 hydrological droughts from 1978 to 2018. The transition rate from meteorological to hydrological drought demonstrated large spatial variability, with a basin‐average rate of 56% and the lowest transition rate of 45% in the headstream The transition from meteorological to hydrological drought was mostly rapid between November and December (~2 months), but can be prolonged during spring (about 3–5 months) and winter (around 7–9 months). Additionally, analysis of regeneration scenarios indicated that human activity has mitigated drought severity over recent decades. The primary driver affecting drought duration and frequency during the transition from meteorological to hydrological droughts shows conspicuous spatial disparities in densely populated areas. Although human activity significantly contributes to drought duration and severity compared to climate change during the transition in the headstream, its effects are more pronounced downstream in terms of drought duration and frequency. One plausible explanation is that increased water consumption downstream has considerably prolonged drought progression, whereas water management has counteractive effects on drought progression due to climate change in the headstream. Our findings offer valuable insights into the transition process from meteorological to hydrological droughts in the presence of extensive human activities.

Phosphate Influx and Dust Deposition Create Zonal and Meridional Biogeochemical Gradients in Trichodesmium Abundance

AbstractTrichodesmium plays a key role in the biogeochemical cycling of carbon, nitrogen and phosphorus in the Tropical Atlantic Ocean. A complex interplay of physicochemical factors control the growth of Trichodesmium. However, owing to the large spatial and temporal variability, the relative influence of these factors in controlling Trichodesmium distribution and abundance remains unclear. In this study, we examined the basin‐scale distribution pattern of Trichodesmium in the upper 200 m water column of the Atlantic Ocean (25°N–30°S and 70°W–20°E) using a large data set (n = 33,235) and tried to constrain the distribution based on various physicochemical parameters. We suggest that the combined effect of warm temperatures and phosphate (PO43−) availability determines the zonal spatial extent and the abundance of Trichodesmium in the Tropical North Atlantic Ocean. However, the availability of dissolved iron, along with high sea surface temperatures and meteorological parameters such as the wind direction and precipitation, likely govern the meridional distribution of Trichodesmium across the Atlantic Ocean. Excess PO43− at the surface rules out the possibility of PO43− limitation in regulating the meridional distribution of the Trichodesmium. Depth‐integrated nitrogen fixation rates, based on a multiple linear regression, vary from 0.07 to 306 μmol N m−2 d−1. The presence of Trichodesmium colonies down to a depth of 200 m and the depth‐integrated nitrogen fixation rates reflect the pivotal role of Trichodesmium in the nitrogen budget of this region.

Higher risk of spring frost under future climate change across China's apple planting regions

Spring frost is one of major agro-meteorological disasters impacting apple production. Previous studies focused on the impacts of spring frost during the flowering period on apples, however, the frost events occurring at both pre- and post-flowering phases have significant impacts on apples. Our study employed four phenology models (parallel model, M1 model, sequential model, and growing degree hour model), along with their ensemble model driven by five global climate models (GCMs) to predict the key spring phenphases of apple across China’s apple planting regions under two emission scenarios (SSP245 and SSP585) during the 2050 s (2040–2069) and 2080 s (2070–2099). A new spring frost assessment index (FI) combining frequency and intensity of spring frost was developed to evaluate the impacts of climate change on spring frost risk of apple. The results showed a large spatial variation in projected spring phenophases across China’s apple planting regions. Projected budbreak (BBCH 09), initial flowering (BBCH 60), and final flowering (BBCH 69) would advance 5.6–15.9 days, 5.7–14.8 days, and 5.6–14.5 days, respectively, in regions I, II, and V of northern China apple planting regions. However, in regions III and IV of southern China apple planting regions, projected budbreak, initial flowering, and final flowering would delay –9.2–21.5 days, –7.5–20.6 days, and –6.8–20.8 days in the future climate scenarios, compared with the baseline (1980–2009). Projected last spring frost date considering the entire spring phenphases would advance 3.7–11.8 d in the future climate scenarios, compared with the baseline. Projected spring frost risk would increase in all the three growth stages including budbreak to initial flowering, flowering period, and final flowering to the last spring frost date across all the China's apple plating regions in future climate scenarios except southward China's apple plating region under the SSP585 in the 2080 s. Our study provides an important reference for early warning, prevention and mitigation of spring frost disasters, and the layout of the apple industry.

Preventing Child Maltreatment: Beneficial Side Effects of Public Childcare

Abstract We investigate the impact of public childcare provision on the incidence of severe child maltreatment. For identification, we exploit a reform that expanded early childcare in Germany, generating large temporal and spatial variation in childcare coverage at the county level. Using high-quality administrative data covering all reported cases of child maltreatment in Germany by county and year, we find that an increase in childcare slots by one percentage point in a county reduced child maltreatment cases leading to out-of-home placement by about 1%. Our results suggest that the provision of universal public childcare may be more cost effective than previously thought.

Modeling pan-Arctic seasonal and interannual landfast sea ice thickness and snow depth between 1979 and 2021

ABSTRACT Landfast sea ice (LFSI) is sensitive to local climate change, making it an important component of the cryosphere system. In this study, the LFSI around the pan-Arctic domain was simulated from 1979 to 2021 using a well-validated snow and ice thermodynamic model (HIGHTSI) under the framework of the Fast Ice Prediction System (FIPS), forced by the ERA5 reanalysis. The simulation results agree well with the in-situ observations in the Canadian Arctic, with a mean error of −0.06 ± 0.29 m for ice thickness and −0.04 ± 0.12 m for snow depth. A decrease of −2.8 ± 0.4 cm/10a in thickness and −16.2 ± 1.5 km3/a in volume for the Arctic LFSI was modeled during this period. There was significant spatial variability among the different domains, with the fastest decline found in the Vilkitsky Strait. The modeled snow depth shows large interannual and spatial variations, which was confirmed by other modeling results. The spatiotemporal variations in both air temperature and precipitation are the driving factors for the multi-decadal variations in LFSI thickness. The decreasing air temperature during the 2010s aligned with a slower thickness decrease and a slight volume increase for LFSI, which agreed with the pan-Arctic sea ice pattern.

Exploring how the heterogeneous urban landscape influences CO2 concentrations: The case study of the Metropolitan Area of Barcelona

Monitoring CO2 concentrations in urban areas is crucial for determining the efficacy of climate change mitigation policies. However, highly heterogeneous land use, local geography, and local convection patterns, which vary throughout the urban landscape, complicate this task. To establish continuous monitoring programs, it is important to first determine the heterogeneity of urban landscapes on the ground. To understand the role these factors play in the distribution of CO2 over an urban area, we conducted a CO2 measurement campaign over the Metropolitan Area of Barcelona (AMB) over four urban land uses: impervious, green, forest, and agricultural. There is a clear tendency for CO2 mixing ratios to decrease as the degree of urban vegetation increases, even in the midst of a developed boundary layer. For example, CO2 concentrations were 429 and 427 ppm at forest and agricultural sites, respectively, while 485 ppm was reported at urban sites. A decrease in atmospheric CO2 was observed from 458 to 428 ppm in the gradient from urban to suburban areas, in which the biosphere component increased. The biosphere component of the CO2 signal was significant and was observed in the gradient from urban to suburban areas, which averaged a reduction from 458 to 428 ppm. Our findings show that the large spatial variability in CO2 concentrations (ranging from 410 to 495 ppm) is best explained by anthropogenic activity. We propose increasing the spatiotemporal resolution of CO2 monitoring in the AMB to determine these trends more precisely over longer periods of time.

Improved Estimates of Sub‐Regional Groundwater Storage Anomaly Using Coordinated Forward Modeling

AbstractGroundwater storage anomaly (GWSA) can be estimated either at the large scale from the Gravity Recovery and Climate Experiment (GRACE) or at the local scale based on in situ observed groundwater level (GWL) and aquifer storage parameters. Yet, the accuracy of GRACE‐based estimate is affected by leakage errors, while that of local GWL‐based estimate requires the reliable specific yield (Sy) data that are usually not available. Here, we developed a novel approach, the coordinated forward modeling (CoFM), based on the iterative forward modeling to improve GWSA estimation at the sub‐regional scale smaller than the typical GRACE footprint. It is achieved by solving Sy through iterative comparisons between GRACE‐based and observation‐based GWSA at 0.5° grid scale, and then re‐calculating GWSA using the updated Sy and observed GWL. The utility of CoFM is explored by using the hypothetical experiments and a real case study in the Piedmont Plain (PP, ∼54,000 km2) and East‐central Plain (ECP, ∼86,000 km2) of North China Plain. Results show that CoFM can detect GWSA at 0.5° grid scale in the hypothetical experiments given the large spatial variability of GWL. While in the real case study, the CoFM distinguishes between the divergent unconfined GWSA trends (2005–2016) in PP (−41.80 ± 0.55 mm/yr) and ECP (−7.57 ± 0.60 mm/yr) caused by the differences in hydrogeological conditions and groundwater use. The improvement made by CoFM can be attributed to the use of the distributed GWL information to constrain GRACE leakage errors. This study highlights a practical important solution for improving sub‐regional GWSA estimation through the joint use of large‐scale GRACE data and local‐scale well observations.

Large-eddy simulation of aerosol concentrations in a realistic urban environment: Model validation and transport mechanism

Air pollution in urban environments exhibits large spatial and temporal variations due to high heterogeneous air flow and emissions. To address the complexity of local air pollutant dynamics, a comprehensive large-eddy simulation using the PALM model system v6.0 was conducted. The distribution of flow and vehicle emitted aerosol particles in a realistic urban environment in Malmö, Sweden, was studied and evaluated against on-site measurements made using portable instrumentation on a spring morning in 2021. The canyon transport mechanisms were investigated, and the convective and turbulent mass-transport rates compared to clarify their role in aerosol transport. The horizontal distribution of aerosols showed acceptable evaluation metrics for both mass and number. Flow and pollutant concentrations were more complex than those in idealized street canyon networks. Vertical turbulent mass-transport rate was found to dominate the mass transport process compared with the convective transport rate, contributing more than 70% of the pollutant transport process. Our findings highlight the necessity of examining various aerosol metric due their distinct dispersion behaviour. This study introduces a comprehensive high-resolution modelling framework that accounts for dynamic meteorological and aerosol background boundary conditions, real-time traffic emission, and detailed building features, offering a robust toll for local urban air quality assessment.

Recent trends and variability of temperature and atmospheric water vapor over South Asia

This study examines the recent trends and variability of temperature and atmospheric water vapor over South Asia during the summer monsoon season from 1970 to 2022. We employ gridded reanalysis and observational datasets to assess the trends and variability in surface temperature and further investigate the associated physical processes. Our results show that the surface temperature over South Asia has increased by 0.5 °C per decade. Notably, we find a statistically significant major transition in the surface temperature series in 1995, such that 1970–1994 (1995–2022) is associated with anomalous cooling (anomalous warming). Since 1995, not only have the temporal variations in regional surface temperature, atmospheric water vapor, and clear-sky downward longwave radiation been synchronous at interannual scales, but their largest spatial variations have been located mainly over the same region. A diagnosis of the atmospheric moisture budget reveals that the observed changes in water vapor convergence over South Asia are dominated by changes in mean circulation dynamics, while the change in thermodynamics is relatively small. The enhancement of the Somali low-level jet may contribute to increased water vapor in South Asia by modulating low-level circulations.