Search

Acute exacerbation of chronic obstruction pulmonary disease (AECOPD) as a respiratory disease, is considered to be related to air pollution by more and more studies. However, the evidence on how air pollution affect the incidence of AECOPD and whether there are population differences is still insufficient. Therefore, we select PM10, PM2.5, SO2, NO2, CO, and O3 as representatives combined with daily AECOPD admission data from 1 January 2015 to 26 June 2016 in the rural areas of Qingyang, northwestern China to explore the associations of air pollution with AECOPD. Based on a time‐stratified case‐crossover design, we constructed a distributed lag nonlinear model to qualify the single and cumulative lagged effects of air pollution on AECOPD. Stratified related risks by sex and age were also reported. The cumulative exposure‐response curves were approximately linear for PM2.5, “V”‐shaped for PM10, “U”‐shaped for NO2 and inverted‐“V” for SO2, CO and O3. Exposure to high‐PM2.5 (42 μg/m3), high‐PM10 (91 μg/m3), high‐SO2 (58 μg/m3), low‐NO2 (12 μg/m3), and high‐CO (1.55 mg/m3) increased the risk of AECOPD. Females aged 15–64 were more susceptible under extreme concentrations of PM2.5, SO2, CO, and low‐PM10 than other subgroups. In addition, adults aged 15–64 were more sensitive to extreme concentrations of NO2 compared with the elderly ≥65 years old, while the latter were more sensitive to high‐PM10. High‐SO2, high‐NO2, and extreme concentrations of PM2.5 had the greatest effects on the day of exposure, while low‐SO2 and low‐CO had lagged effects on AECOPD. Precautionary measures should be taken with a focus on vulnerable subgroups, to control hospitalization for AECOPD associated with air pollutants., Key Points Exposure to high‐PM2.5, high‐PM10, high‐SO2, low‐NO2, and high‐CO increased the risk of acute exacerbation of chronic obstruction pulmonary disease (COPD)The cumulative curves were approximately linear for PM2.5, “V”‐shaped for PM10, “U”‐shaped for NO2 and inverted‐”V” for SO2, CO and O3 The nonlinear effects on acute exacerbation of COPD at different lags varied based on the air pollutants, involved gender and age

Recurrent and large forest fires negatively impact ecosystem, air quality, and human health. Moderate Resolution Imaging Spectroradiometer fire product is used to identify forest fires over central India domain, an extremely fire prone region. The study finds that from 2001 to 2020, ∼70% of yearly forest fires over the region occurred during March (1,857.5 counts/month) and April (922.8 counts/month). Some years such as 2009, 2012, and 2017 show anomalously high forest fires. The role of persistent warmer temperatures and multiple climate extremes in increasing forest fire activity over central India is comprehensively investigated. Warmer period from 2006 to 2020 showed doubling and tripling of forest fire activity during forest fire (February–June; FMAMJ) and non‐fire (July–January; JASONDJ) seasons, respectively. From 2015 JASONDJ to 2018 FMAMJ, central India experienced a severe heatwave, a rare drought and an extremely strong El Niño, the combined effect of which is linked to increased forest fires. Further, the study assesses quinquennial spatiotemporal changes in forest fire characteristics such as fire count density and average fire intensity. Deciduous forests of Jagdalpur‐Gadchiroli Range and Indravati National Park in Chhattisgarh state are particularly fire prone (>61 fire counts/grid) during FMAMJ and many forest fires are of high intensity (>45 MW). Statistical associations link high near surface air temperature and low precipitation during FMAMJ to significantly high soil temperature, low soil moisture content, low evapotranspiration and low normalized difference vegetation index. This creates a significantly drier environment, conducive for high forest fire activity in the region., Key Points Compared to 2001–2005, central India domain forest fire activity during 2006–2020 doubled in forest fire season and tripled in non‐fire seasonRole of persistent warmer temperatures and multiple climate extremes in increasing forest fire activity over central India is highlightedDeciduous forests of Jagdalpur‐Gadchiroli Range and Indravati National Park are extremely fire prone and fires are of high intensity

Following extreme drought during the 2019–2020 bushfire summer, the eastern part of Australia suffered from a week‐long intense rainfall and extensive flooding in March 2021. Understanding how much water storage changes in response to these climate extremes is critical for developing timely water management strategies. To quantify prompt water storage changes associated with the 2021 March flooding, we processed the low‐latency (1–3 days), high‐precision intersatellite laser ranging measurements from GRACE Follow‐On spacecraft and determined instantaneous gravity changes along spacecraft orbital passes. Such new data processing detected an abrupt surge of water storage approaching 60–70 trillion liters (km3 of water) over a week in the region, which concurrently caused land subsidence of ∼5 mm measured by a network of ground GPS stations. This was the highest speed of ground water recharge ever recorded in the region over the last two decades. Compared to the condition in February 2020, the amount of recharged water was similar but the recharge speed was much faster in March 2021. While these two events together replenished the region up to ∼80% of the maximum storage over the last two decades, the wet antecedent condition of soils in 2021 was distinctly different from the dry conditions in 2020 and led to generating extensive runoff and flooding in 2021., Key Points New GRACE Follow‐On gravity data processing quantified prompt changes in water storage by the heavy rainfall and flooding in March 2021The eastern Australia experienced the highest speed of ground water recharge and wet antecedent soils were responsible for extensive flooding in 2021The study demonstrated the feasibility of rapid processing (1–3 days) for immediate assessment of mass changes from extreme events

Estimates of Amazon rainforest gross primary productivity (GPP) differ by a factor of 2 across a suite of three statistical and 18 process models. This wide spread contributes uncertainty to predictions of future climate. We compare the mean and variance of GPP from these models to that of GPP at six eddy covariance (EC) towers. Only one model's mean GPP across all sites falls within a 99% confidence interval for EC GPP, and only one model matches EC variance. The strength of model response to climate drivers is related to model ability to match the seasonal pattern of the EC GPP. Models with stronger seasonal swings in GPP have stronger responses to rain, light, and temperature than does EC GPP. The model to data comparison illustrates a trade‐off inherent to deterministic models between accurate simulation of a mean (average) and accurate responsiveness to drivers. The trade‐off exists because all deterministic models simplify processes and lack at least some consequential driver or interaction. If a model's sensitivities to included drivers and their interactions are accurate, then deterministically predicted outcomes have less variability than is realistic. If a GPP model has stronger responses to climate drivers than found in data, model predictions may match the observed variance and seasonal pattern but are likely to overpredict GPP response to climate change. High or realistic variability of model estimates relative to reference data indicate that the model is hypersensitive to one or more drivers., Key Points Regression logic is reason to doubt the accurate climate sensitivity of predictions whose variability is realistic or higherA suite of models poorly reproduces eddy covariance estimates of Amazon rainforest gross primary productivityHighly seasonal models predict stronger primary productivity responsiveness to meteorology than is likely to be true

GNSS Reflectometry (GNSS‐R) measurements over inland water bodies, such as lakes, rivers, and wetlands exhibit strong coherent signals. The strength of the coherent reflections is highly sensitive to small‐scale surface roughness. For inland waters, this roughness is primarily due to wind‐driven surface waves. The sensitivity of the coherent reflections to surface roughness can be leveraged to estimate wave height profiles across inland waters. Coupled with a wind wave model, an approach to retrieve a wind vector is described using a forward model, which is potentially able to predict scattered power profiles for different wind speeds and directions and choosing the minimum‐squared error solution. The ability for spaceborne or airborne GNSS‐R to measure an inland water wind vector and wave heights could contribute to scientific applications focused on understanding nearshore ecosystems, monitoring climate change effects on inland waters, sediment resuspension, biomass production, fish habitat, and others. This paper presents a novel approach to potentially retrieve wind vector and wave heights over inland waters using GNSS‐R and discusses the issues with performing such retrievals using simulation and very few available raw signals recorded from CYGNSS satellites., Key Points Global Navigation Satellite Systems (GNSS) signals reflected from inland waters exhibit coherent scattering properties that make them highly sensitive to water surface roughnessWind‐induced roughness will vary across a water body, with a strong dependence on wind speed, wind direction, and the water depthDiscuss potentials of a novel model‐based approach for retrieval a of wind vector and surface wave heights over lakes using CYGNSS data

By injecting a mixture of gas and pyroclasts into the atmosphere, explosive volcanic eruptions frequently generate vortex rings, which are toroidal vortices formed by the jet's initial momentum. Here, we report high‐speed imaging and acoustic measurements of vortex rings sourcing from gas‐rich eruptive jets at Stromboli volcano (Italy). Volcanic vortex rings (VVRs) form at the vent together with an initial compression acoustic wave, VVRs maximum rise velocity being directly proportional to the amplitude and inversely proportional to the duration of the compression wave. The axial rise and acoustic signature of VVRs match well those predicted by recent fluid‐dynamic experiments. This good match allows using the high‐frequency (80–1,000 Hz) component of the jet sound and the time‐dependent rise of VVRs to retrieve two key eruption parameters: the Mach number of the eruptive jets (, Key Points Volcanic vortex rings are formed by gas‐rich, jet‐forming Strombolian‐style explosionsThe explosion acoustic signals bear the signature of the vortex rings and reveal supersonic eruptive jets with Mach number up to 1.5Vent diameter can be estimated from the time‐dependent rise of vortex rings at 0.7 m, in agreement with direct observation from Uncrewed Aerial Vehicle

The Mesozoic era (∼252 to 66 million years ago) was a key interval in Earth's evolution toward its modern state, witnessing the breakup of the supercontinent Pangaea and significant biotic innovations like the early evolution of mammals. Plate tectonic dynamics drove a fundamental climatic transition from the early Mesozoic supercontinent toward the Late Cretaceous fragmented continental configuration. Here, key aspects of Mesozoic long‐term environmental changes are assessed in a climate model ensemble framework. We analyze so far the most extended ensemble of equilibrium climate states simulated for evolving Mesozoic boundary conditions covering the period from 255 to 60 Ma in 5 Myr timesteps. Global mean temperatures are generally found to be elevated above the present and exhibit a baseline warming trend driven by rising sea levels and increasing solar luminosity. Warm (Triassic and mid‐Cretaceous) and cool (Jurassic and end‐Cretaceous) anomalies result from pCO2 changes indicated by different reconstructions. Seasonal and zonal temperature contrasts as well as continental aridity show an overall decrease from the Late Triassic‐Early Jurassic to the Late Cretaceous. Meridional temperature gradients are reduced at higher global temperatures and less land area in the high latitudes. With systematic sensitivity experiments, the influence of paleogeography, sea level, vegetation patterns, pCO2, solar luminosity, and orbital configuration on these trends is investigated. For example, long‐term seasonality trends are driven by paleogeography, but orbital cycles could have had similar‐scale effects on shorter timescales. Global mean temperatures, continental humidity, and meridional temperature gradients are, however, also strongly affected by pCO2., Key Points We assess global long‐term climate trends through the Mesozoic era with an ensemble of climate model simulationsVarying carbon dioxide levels cause anomalies around an overall warming trend due to changing paleogeography and increasing insolationSeasonal and zonal temperature contrasts as well as aridity decrease with time, while meridional gradients vary with paleogeography

Understanding the complex interrelationships between wildfire and its environmental and anthropogenic controls is crucial for wildfire modeling and management. Although machine learning (ML) models have yielded significant improvements in wildfire predictions, their limited interpretability has been an obstacle for their use in advancing understanding of wildfires. This study builds an ML model incorporating predictors of local meteorology, land‐surface characteristics, and socioeconomic variables to predict monthly burned area at grid cells of 0.25° × 0.25° resolution over the contiguous United States. Besides these predictors, we construct and include predictors representing the large‐scale circulation patterns conducive to wildfires, which largely improves the temporal correlations in several regions by 14%–44%. The Shapley additive explanation is introduced to quantify the contributions of the predictors to burned area. Results show a key role of longitude and latitude in delineating fire regimes with different temporal patterns of burned area. The model captures the physical relationship between burned area and vapor pressure deficit, relative humidity (RH), and energy release component (ERC), in agreement with the prior findings. Aggregating the contribution of predictor variables of all the grids by region, analyses show that ERC is the major contributor accounting for 14%–27% to large burned areas in the western US. In contrast, there is no leading factor contributing to large burned areas in the eastern US, although large‐scale circulation patterns featuring less active upper‐level ridge‐trough and low RH two months earlier in winter contribute relatively more to large burned areas in spring in the southeastern US., Key Points US fire burned areas are well predicted by a machine learning model and SHAP improves interpretation of the contributing factorsIncluding large‐scale circulation patterns conducive to wildfires as predictors improve prediction of burned areas in several regionsFire‐season fuel dryness and dry winters are important contributors to large burned areas in western and southeastern US, respectively

Coastal boulder deposits (CBD) provide what are sometimes the only remaining signatures of wave inundation on rocky coastlines; in recent decades, CBD combined with initiation of motion (IoM) analyses have repeatedly been used as primary evidence to infer the existence of ancient tsunamis. However, IoM storm wave heights inferred by these studies have been shown to be highly inaccurate, bringing some inferences into question. This work develops a dimensionless framework to relate CBD properties with storm‐wave hindcasts and measurements, producing data‐driven relations between wave climate and boulder properties. We present an elevation‐density‐size‐inland distance‐wave height analysis for individual storm‐transported boulders which delineates the dynamic space where storm‐wave CBD occur. Testing these new relations against presumed tsunami CBD demonstrates that some fall well within the capabilities of storm events, suggesting that some previous studies might be fruitfully reexamined within the context of this new framework., Key Points Size, elevation, and distance inland in storm‐wave‐generated coastal boulder deposits are all functions of climatological wave heightsClear data ranges exist where storm wave coastal boulder deposits do/do not existSome coastal boulder deposits previously identified as having tsunami origin may have been generated by storm waves

Six conceptually different transport models were applied to the macrodispersion experiment (MADE)-1 field tracer experiment as a first major attempt for model comparison. The objective was to show that complex mass distributions in heterogeneous aquifers can be predicted without calibration of transport parameters, solely making use of structural and flow data. The models differ in their conceptualization of the heterogeneous aquifer structure, computational complexity, and use of conductivity data obtained from various observation methods (direct push injection logging, DPIL, grain size analysis, pumping tests and flowmeter). They share the same underlying physical transport process of advection by the velocity field solely. Predictive capability is assessed by comparing results to observed longitudinal mass distributions of the MADE-1 experiment. The decreasing mass recovery of the observed plume is attributed to sampling and no physical process like mass transfer is invoked by the models. Measures like peak location and strength are used in comparing the modeled and measured plume mass distribution. Comparison of models reveals that the predictions of the solute plume agree reasonably well with observations, if the models are underlain by a few parameters of close values: mean velocity, a parameter reflecting log-conductivity variability, and a horizontal length scale related to conductivity spatial correlation. The robustness of the results implies that conservative transport models with appropriate conductivity upscaling strategies of various observation data provide reasonable predictions of plumes longitudinal mass distribution, as long as key features are taken into account., The authors thank Marco Bianchi and Boris Baeumer for their support during the development of the study. A. Fiori and A. Bellin acknowledge funding from the Italian Ministry of Education, University and Research (MIUR) in the frame of the Departments of Excellence Initiative 2018–2022 granted to Dept. of Engineering of Roma Tre University, and to the Dept. of Civil, Environmental and Mechanical Engineering of the University of Trento, respectively. M. Dentz acknowledges funding of the European Research Council (ERC) through the project MHetScale (contract number 617511), and the Spanish Research Agency (AEI) through the project HydroPore (contract number PID2019-106887GB-C31).

Intensive observations and WRF‐Chem simulations are applied in this study to investigate the adverse impacts of regional transport on the PM2.5 (fine particulate matter; diameter ≤2.5 μm) changes in Shanghai during the Coronavirus Disease 2019 lockdown. As the local atmospheric oxidation capacity was observed to be generally weakened, strong regional transport carried by the frequent westerly winds is suggested to be the main driver of the unexpected pollution episodes, increasing the input of both primary and secondary aerosols. Contributing 40%–80% to the PM2.5, the transport contributed aerosols are simulated to exhibit less decreases (13.2%–21.8%) than the local particles (37.1%–64.8%) in urban Shanghai due to the lockdown, which largely results from the less decreased industrial and residential emissions in surrounding provinces. To reduce the influence of the transport, synergetic emission control, especially synergetic ammonia control, measures are proved to be effective strategies, which need to be considered in future regulations., Key Points PM2.5 pollution in Shanghai were mostly attributed to the regional transport carried by frequent westerly winds during the COVID‐lockdownThe transport contributed aerosols in Shanghai are less affected by the lockdown measures due to the less decreased surrounding emissionsSynergetic emission control measures, especially those on ammonia emissions, are effective in mitigating the adverse transport impacts

The intensification of heat stress reduces the labor capacity and hence poses a threat to socio‐economic development. The reliable projection of the changing climate and the development of sound adaptation strategies are thus desired for adapting to the decreasing labor productivity under climate change. In this study, an optimization modeling approach coupled with dynamical downscaling is proposed to design the optimal adaptation strategies for improving labor productivity under heat stress in China. The future changes in heat stress represented by the wet‐bulb globe temperature (WBGT) are projected with a spatial resolution of 25 × 25 km by a regional climate model (RCM) through the dynamical downscaling of its driving global climate model (GCM). Uncertain information such as system costs, environmental costs, and subsidies are also incorporated into the optimization process to provide reliable decision alternatives for improving labor productivity. Results indicate that the intensification of WBGT is overestimated by the GCM compared to the RCM. Such an overestimation can lead to more losses in working hours derived from the GCM than those from the RCM regardless of climate scenarios. Nevertheless, the overestimated heat stress does not alter the regional measures taken to adapt to decreasing labor productivity. Compared to inland regions, the monsoon‐affected regions tend to improve labor productivity by applying air conditioning rather than working overtime due to the cost differences. Consequently, decision‐makers need to optimally make a balance between working overtime and air conditioning measures to meet sustainable development goals., Key Points An optimization model coupled with dynamical downscaling is developed to improve labor productivity under intensified heat stress in ChinaCivilian workers in inland regions will be more vulnerable to the intensified heat stress in a changing climateCompared to the regional climate model, the global climate model overestimates heat stress and thus exaggerates related system costs to recover from the working hour loss

Supraglacial debris affects glacier mass balance as a thin layer enhances surface melting, while a thick layer reduces it. While many glaciers are debris‐covered, global glacier models do not account for debris because its thickness is unknown. We provide the first globally distributed debris thickness estimates using a novel approach combining sub‐debris melt and surface temperature inversion methods. Results are evaluated against observations from 22 glaciers. We find the median global debris thickness is ∼0.15 ± 0.06 m. In all regions, the net effect of accounting for debris is a reduction in sub‐debris melt, on average, by 37%, which can impact regional mass balance by up to 0.40 m water equivalent (w.e.) yr‐1. We also find recent observations of similar thinning rates over debris‐covered and clean ice glacier tongues is primarily due to differences in ice dynamics. Our results demonstrate the importance of accounting for debris in glacier modeling efforts., Key Points We produce the first distributed global debris thickness estimatesAccounting for debris significantly reduces regional glacier mass lossThe similar thinning rates of debris‐covered and clean ice glaciers in High Mountain Asia is primarily caused by differences in ice dynamics

We document an exceptional large‐spatial scale case of changes in tidal range in the North Sea, featuring pronounced trends between −2.3 mm/yr at tide gauges in the United Kingdom and up to 7 mm/yr in the German Bight between 1958 and 2014. These changes are spatially heterogeneous and driven by a superposition of local and large‐scale processes within the basin. We use principal component analysis to separate large‐scale signals appearing coherently over multiple stations from rather localized changes. We identify two leading principal components (PCs) that explain about 69% of tidal range changes in the entire North Sea including the divergent trend pattern along United Kingdom and German coastlines that reflects movement of the region’s semidiurnal amphidromic areas. By applying numerical and statistical analyses, we can assign a baroclinic (PC1) and a barotropic large‐scale signal (PC2), explaining a large part of the overall variance. A comparison between PC2 and tide gauge records along the European Atlantic coast, Iceland, and Canada shows significant correlations on time scales of less than 2 years, which points to an external and basin‐wide forcing mechanism. By contrast, PC1 dominates in the southern North Sea and originates, at least in part, from stratification changes in nearby shallow waters. In particular, from an analysis of observed density profiles, we suggest that an increased strength and duration of the summer pycnocline has stabilized the water column against turbulent dissipation and allowed for higher tidal elevations at the coast., Key Points 70 North Sea tide gauges evince contrasting trends in tidal range between the United Kingdom (−1.0 mm/yr) and the German Bight (3.3 mm/yr) since 1958We use principal component analysis (PCA) to separate local (e.g., building measures) from large‐scale effects of oceanographic originThe first PC explains 77% of variance in the German Bight and is linked to stability changes in shallow, seasonally stratified waters