ALBERT

Description: Porous poly(N-isopropylacrylamide) (PNIPAAm) hydrogels with pore diameters in the nanometer and the micrometer range were synthesized using two variations of the surfactant-based template method. We showed that smaller pore diameters lead to faster swelling and deswelling. A graphic representation of a model describing the swelling kinetics explains the assumption that the release and the uptake of water is faster for a larger specific surface area, that is smaller pore diameters. Additionally, the open-porous channel structure benefits the water transport compared to a nonporous PNIPAAm hydrogel. Sensor measurements result in response times between 36 s and 4 min also showing the importance of mechanical stability of porous hydrogels.

Description: Fixed feed-in tariffs based on the Renewable Energy Act grant secure revenues from selling electricity for wind turbine operators in Germany. Anyhow, the level of federal financial support is being reduced consecutively. Plant operators must trade self-sufficiently in the future and hence generate revenue by selling electricity directly on electricity markets. Therefore, uncertain future market price developments will influence investment considerations and may lead to stagnation in the expansion of renewable energies. This study estimates future revenue potentials of non-subsidized wind turbines in Germany to reduce this risk. The paper introduces and analyses a forecasting model that generates electricity price time series suited for revenue estimation of wind turbines based on the electricity exchange market. Revenues from the capacity market are neglected. The model is based on openly accessible data and applies a merit-order approach in combination with a simple agent-based approach to forecast long-term day-ahead prices at an hourly resolution. The hourly generation profile of wind turbines can be mapped over several years in conjunction with fluctuations in the electricity price. Levelized revenue of energy is used to assess both dynamic variables (electricity supply and price). The merit-order effect from the expansion of renewables as well as the phasing out of nuclear energy and coal are assessed in a scenario analysis. Based on the assumptions made, the opposing effects could result in a constant average price level for Germany over the next 20 years. The influence of emission prices is considered in a sensitivity analysis and correlates with the share of fossil generation capacities in the generation mix. In a brief case study, it was observed that current average wind turbines are not able to yield financial profit over their lifetime without additional subsidies for the given case. This underlines a need for technical development and new business models like power purchase agreements. The model results can be used for setting and negotiating appropriate terms, such as energy price schedule or penalties for those agreements.

Description: Vegetation optical depth (VOD) retrieved from microwave radiometry correlates with the total amount of water in vegetation, based on theoretical and empirical evidence. Because the total amount of water in vegetation varies with relative water content (as well as with biomass), this correlation further suggests a possible relationship between VOD and plant water potential, a quantity that drives plant hydraulic behavior. Previous studies have found evidence for that relationship on the scale of satellite pixels tens of kilometers across, but these comparisons suffer from significant scaling error. Here we used small-scale remote sensing to test the link between remotely sensed VOD and plant water potential. We placed an L-band radiometer on a tower above the canopy looking down at red oak forest stand during the 2019 growing season in central Massachusetts, United States. We measured stem xylem and leaf water potentials of trees within the stand and retrieved VOD with a single-channel algorithm based on continuous radiometer measurements and measured soil moisture. VOD exhibited a diurnal cycle similar to that of leaf and stem water potential, with a peak at approximately 05:00 eastern daylight time (UTC−4). VOD was also positively correlated with both the measured dielectric constant and water potentials of stem xylem over the growing season. The presence of moisture on the leaves did not affect the observed relationship between VOD and stem water potential. We used our observed VOD–water-potential relationship to estimate stand-level values for a radiative transfer parameter and a plant hydraulic parameter, which compared well with the published literature. Our findings support the use of VOD for plant hydraulic studies in temperate forests.

Description: The case of the heavy precipitation event on 14 and 15 October 2018 which has led to severe flash flooding in the Aude watershed in south-western France is studied from a meteorological point of view using deterministic and probabilistic numerical weather prediction systems, as well as a unique combination of observations from both standard and personal weather stations. This case features typical characteristics of Mediterranean heavy precipitation events such as its classic synoptic situation and its quasi-stationary convective precipitation that regenerates continuously, as well as some peculiarities such as the presence of a former hurricane and a pre-existing cold air mass close to the ground. Mediterranean Sea surface temperature and soil moisture anomalies are briefly reviewed, as they are known to play a role in this type of hydrometeorological events. A study of rainfall forecasts shows that the event had limited predictability, in particular given the small size of the watersheds involved. It is shown that the stationarity of precipitation, whose estimation benefits from data from personal stations, is linked to the presence near the ground of a trough and a strong potential virtual temperature gradient, the stationarity of both of which is highlighted by a combination of observations from standard and personal stations. The forecast that comes closest to the rainfall observations contains the warmest, wettest, and fastest low-level jet and also simulates near the ground a trough and a marked boundary between cold air in the west and warm air in the east, both of which are stationary.

Description: Marine stratocumulus cloud properties over the Southeast Atlantic Ocean are impacted by contact between above-cloud biomass burning aerosols and cloud tops. Different vertical separations (0 to 2000 m) between the aerosol layer and cloud tops were observed on six research flights in September 2016 during the NASA ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) field campaign. There were 30 contact profiles, where an aerosol layer with aerosol concentration (Na) 〉 500 cm−3 was within 100 m of cloud tops, and 41 separated profiles, where the aerosol layer with Na 〉 500 cm−3 was located more than 100 m above cloud tops. For contact profiles, the average cloud droplet concentration (Nc) in the cloud layer was up to 68 cm−3 higher, the effective radius (Re) up to 1.3 µm lower, and the liquid water content (LWC) within 0.01 g m−3 compared to separated profiles. Free-tropospheric humidity was higher in the presence of biomass burning aerosols, and contact profiles had a smaller decrease in humidity (and positive buoyancy) across cloud tops with higher median above-cloud Na (895 cm−3) compared to separated profiles (30 cm−3). Due to droplet evaporation from entrainment mixing of warm, dry free-tropospheric air into the clouds, the median Nc and LWC for contact profiles decreased with height by 21 and 9 % in the top 20 % of the cloud layer. The impact of droplet evaporation was stronger during separated profiles as a greater decrease in humidity (and negative buoyancy) across cloud tops led to greater decreases in median Nc (30 %) and LWC (16 %) near cloud tops. Below-cloud Na was sampled during 61 profiles, and most contact profiles (20 out of 28) were within high-Na (〉 350 cm−3) boundary layers, while most separated profiles (22 out of 33) were within low-Na (

Description: Mass loss from the Antarctic Ice Sheet is the main source of uncertainty in projections of future sea-level rise, with important implications for coastal regions worldwide. Central to ongoing and future changes is the marine ice sheet instability: once a critical threshold, or tipping point, is crossed, ice internal dynamics can drive a self-sustaining retreat committing a glacier to irreversible, rapid and substantial ice loss. This process might have already been triggered in the Amundsen Sea region, where Pine Island and Thwaites glaciers dominate the current mass loss from Antarctica, but modelling and observational techniques have not been able to establish this rigorously, leading to divergent views on the future mass loss of the West Antarctic Ice Sheet. Here, we aim at closing this knowledge gap by conducting a systematic investigation of the stability regime of Pine Island Glacier. To this end we show that early warning indicators in model simulations robustly detect the onset of the marine ice sheet instability. We are thereby able to identify three distinct tipping points in response to increases in ocean-induced melt. The third and final event, triggered by an ocean warming of approximately 1.2 ∘C from the steady-state model configuration, leads to a retreat of the entire glacier that could initiate a collapse of the West Antarctic Ice Sheet.

Description: Dimethyl sulfide (DMS), dimethylsulfoniopropionate (DMSP) and dimethyl sulfoxide (DMSO) were measured at the Boknis Eck Time Series Station (BE, Eckernförde Bay, SW Baltic Sea) during the period February 2009–December 2018. Our results show considerable interannual and seasonal variabilities in the mixed-layer concentrations of DMS, total DMSP (DMSPt) and total DMSO (DMSOt). Positive correlations were found between particulate DMSP (DMSPp) and particulate DMSO (DMSOp) as well as DMSPt and DMSOt in the mixed layer, suggesting a similar source for both compounds. The decreasing long-term trends, observed for DMSPt and DMS in the mixed layer, were linked to the concurrent trend of the sum of 19′-hexanoyloxyfucoxanthin and 19′-butanoyloxy-fucoxanthin, which are the marker pigments of prymnesiophytes and chrysophytes, respectively. Major Baltic inflow (MBI) events influenced the distribution of sulfur compounds due to phytoplankton community changes, and sediment might be a potential source for DMS in the bottom layer during seasonal hypoxia/anoxia at BE. A modified algorithm based on the phytoplankton pigments reproduces the DMSPp : Chl a ratios well during this study and could be used to estimate future surface (5 m) DMSPp concentrations at BE.

Description: Seismometers have detected the social response to lockdown measures implemented following the onset of the COVID-19 pandemic in cities around the world. This long-lasting pandemic has been a particular challenge in countries such as Mexico, where the informal economy constitutes most of the working population. This context motivated the monitoring of the mobility of populations throughout the various phases of lockdown measures independently of people's access to the internet and mobile technology. Here we use the variation of anthropogenic seismic noise in the city of Querétaro (central Mexico) recorded by a network of low-cost Raspberry Shake seismic stations to study the spatial and temporal variation of human activity in the city throughout the pandemic and during sporting events. The results emphasize the importance of densifying urban seismic networks and of tracking human activities without the privacy concerns associated with mobile technology.

Description: Virioplankton are a key component of the marine biosphere in maintaining diversity of microorganisms and stabilizing ecosystems. They also contribute greatly to nutrient cycles/cycling by releasing organic matter after lysis of hosts. In this study, we constructed the first global viral oceanography database (gVOD) by collecting 10 931 viral abundance (VA) data and 727 viral production (VP) data, along with host and relevant oceanographic parameters when available. Most VA data were obtained in the North Atlantic (32 %) and North Pacific (29 %) oceans, while the southeast Pacific and Indian oceans were quite undersampled. The VA in the global ocean was 1.17(±3.31)×107 particles mL−1. The lytic and lysogenic VP in the global ocean was 9.87(±24.16)×105 and 2.53(±8.64)×105 particlesmL-1h-1, respectively. Average VA in coastal oceans was higher than that in surface open oceans (3.61(±6.30)×107 versus 0.73(±1.24)×107 particles mL−1), while average VP in coastal and surface open oceans was close. Vertically, VA, lytic VP and lysogenic VP deceased from surface to deep oceans by about 1 order of magnitude. The total number of viruses in the global ocean estimated by bin-averaging and the random forest method was 1.56×1030 and 1.49×1030 particles, leading to an estimate of global ocean viral biomass at 35.9 and 34.4 Tg C, respectively. We expect that the gVOD will be a fundamental and very useful database for laboratory, field and modeling studies in marine ecology and biogeochemistry. The full gVOD database (Xie et al., 2020) is stored in PANGAEA (https://doi.org/10.1594/PANGAEA.915758).

Description: Oil seed crops, especially oil palm, are among the most rapidly expanding agricultural land uses, and their expansion is known to cause significant environmental damage. Accordingly, these crops often feature in public and policy debates which are hampered or biased by a lack of accurate information on environmental impacts. In particular, the lack of accurate global crop maps remains a concern. Recent advances in deep-learning and remotely sensed data access make it possible to address this gap. We present a map of closed-canopy oil palm (Elaeis guineensis) plantations by typology (industrial versus smallholder plantations) at the global scale and with unprecedented detail (10 m resolution) for the year 2019. The DeepLabv3+ model, a convolutional neural network (CNN) for semantic segmentation, was trained to classify Sentinel-1 and Sentinel-2 images onto an oil palm land cover map. The characteristic backscatter response of closed-canopy oil palm stands in Sentinel-1 and the ability of CNN to learn spatial patterns, such as the harvest road networks, allowed the distinction between industrial and smallholder plantations globally (overall accuracy =98.52±0.20 %), outperforming the accuracy of existing regional oil palm datasets that used conventional machine-learning algorithms. The user's accuracy, reflecting commission error, in industrial and smallholders was 88.22 ± 2.73 % and 76.56 ± 4.53 %, and the producer's accuracy, reflecting omission error, was 75.78 ± 3.55 % and 86.92 ± 5.12 %, respectively. The global oil palm layer reveals that closed-canopy oil palm plantations are found in 49 countries, covering a mapped area of 19.60 Mha; the area estimate was 21.00 ± 0.42 Mha (72.7 % industrial and 27.3 % smallholder plantations). Southeast Asia ranks as the main producing region with an oil palm area estimate of 18.69 ± 0.33 Mha or 89 % of global closed-canopy plantations. Our analysis confirms significant regional variation in the ratio of industrial versus smallholder growers, but it also confirms that, from a typical land development perspective, large areas of legally defined smallholder oil palm resemble industrial-scale plantings. Since our study identified only closed-canopy oil palm stands, our area estimate was lower than the harvested area reported by the Food and Agriculture Organization (FAO), particularly in West Africa, due to the omission of young and sparse oil palm stands, oil palm in nonhomogeneous settings, and semi-wild oil palm plantations. An accurate global map of planted oil palm can help to shape the ongoing debate about the environmental impacts of oil seed crop expansion, especially if other crops can be mapped to the same level of accuracy. As our model can be regularly rerun as new images become available, it can be used to monitor the expansion of the crop in monocultural settings. The global oil palm layer for the second half of 2019 at a spatial resolution of 10 m can be found at https://doi.org/10.5281/zenodo.4473715 (Descals et al., 2021).

Description: Recent advances in soil moisture remote sensing have produced satellite data sets with improved soil moisture mapping under vegetation and with higher spatial and temporal resolutions. In this study, we evaluate the potential of a new, experimental version of the Advanced Scatterometer (ASCAT) soil water index data set for multiple objective calibrations of a conceptual hydrologic model. The analysis is performed in 213 catchments in Austria for the period 2000–2014. An HBV (Hydrologiska Byråns Vattenbalansavdelning)-type hydrologic model is calibrated based on runoff data, ASCAT soil moisture data, and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover data for various calibration variants. Results show that the inclusion of soil moisture data in the calibration mainly improves the soil moisture simulations, the inclusion of snow data mainly improves the snow simulations, and the inclusion of both of them improves both soil moisture and snow simulations to almost the same extent. The snow data are more efficient at improving snow simulations than the soil moisture data are at improving soil moisture simulations. The improvements of both runoff and soil moisture model efficiencies are larger in low elevation and agricultural catchments than in others. The calibrated snow-related parameters are strongly affected by including snow data and, to a lesser extent, by soil moisture data. In contrast, the soil-related parameters are only affected by the inclusion of soil moisture data. The results indicate that the use of multiple remote sensing products in hydrological modeling can improve the representation of hydrological fluxes and prediction of runoff hydrographs at the catchment scale.

Description: Quantum gravimeters are a promising new development allowing for continuous absolute gravity monitoring while remaining user-friendly and transportable. In this study, we present experiments carried out to assess the capacity of the AQG#B01 in view of future deployment as a field gravimeter for hydrogeophysical applications. The AQG#B01 is the field version follow-up of the AQG#A01 portable absolute quantum gravimeter developed by the French quantum sensor company Muquans. We assess the instrument's performance in terms of stability (absence of instrumental drift) and sensitivity in relation to other gravimeters. No significant instrumental drift was observed over several weeks of measurement. We discuss the observations concerning the accuracy of the AQG#B01 in comparison with a state-of-the-art absolute gravimeter (Micro-g-LaCoste, FG5#228). We report the repeatability to be better than 50 nm s−2. This study furthermore investigates whether changes in instrument tilt and external temperature and a combination of both, which are likely to occur during field campaigns, influence the measurement of gravitational attraction. We repeatedly tested external temperatures between 20 and 30 ∘C and did not find any significant effect. As an example of a geophysical signal, a 100 nm s−2 gravity change is detected with the AQG#B01 after a rainfall event at the Larzac geodetic observatory (southern France). The data agreed with the gravity changes measured with a superconducting relative gravimeter (GWR, iGrav#002) and the expected gravity change simulated as an infinite Bouguer slab approximation. We report 2 weeks of stable operation under semi-terrain conditions in a garage without temperature-control. We close with operational recommendations for potential users and discuss specific possible future field applications. While not claiming completeness, we nevertheless present the first characterization of a quantum gravimeter carried out by future users. Selected criteria for the assessment of its suitability in field applications have been investigated and are complemented with a discussion of further necessary experiments.

Description: Air pollution is one of the major challenges in urban areas. It can have a major impact on human health and society and is currently a subject of several litigations in European courts. Information on the level of air pollution is based on near-surface measurements, which are often irregularly distributed along the main traffic roads and provide almost no information about the residential areas and office districts in the cities. To further enhance the process understanding and give scientific support to decision makers, we developed a prototype for an air quality forecasting system (AQFS) within the EU demonstration project “Open Forecast”. For AQFS, the Weather Research and Forecasting model together with its coupled chemistry component (WRF-Chem) is applied for the Stuttgart metropolitan area in Germany. Three model domains from 1.25 km down to a turbulence-permitting resolution of 50 m were used, and a single-layer urban canopy model was active in all domains. As a demonstration case study, 21 January 2019 was selected, which was a heavily polluted day with observed PM10 concentrations exceeding 50 µg m−3. Our results show that the model is able to reasonably simulate the diurnal cycle of surface fluxes and 2 m temperatures as well as evolution of the stable and shallow boundary layer typically occurring in wintertime in Stuttgart. The simulated fields of particulates with a diameter of less than 10 µm (PM10) and nitrogen dioxide (NO2) allow a clear statement about the most heavily polluted areas apart from the irregularly distributed measurement sites. Together with information about the vertical distribution of PM10 and NO2 from the model, AQFS will serve as a valuable tool for air quality forecasting and has the potential of being applied to other cities around the world.

Description: Multiple global change drivers affect plant productivity of grasslands and thus ecosystem services like forage production and the soil carbon sink. Subalpine grasslands seem particularly affected and may serve as a proxy for the cold, continental grasslands of the Northern Hemisphere. Here, we conducted a 4-year field experiment (AlpGrass) with 216 turf monoliths, subjected to three global change drivers: warming, moisture, and N deposition. Monoliths from six different subalpine pastures were transplanted to a common location with six climate scenario sites (CSs). CSs were located along an altitudinal gradient from 2360 to 1680 m a.s.l., representing an April–October mean temperature change of −1.4 to +3.0 ∘C, compared to CSreference with no temperature change and with climate conditions comparable to the sites of origin. To uncouple temperature effects along the altitudinal gradient from soil moisture and soil fertility effects, an irrigation treatment (+12 %–21 % of ambient precipitation) and an N-deposition treatment (+3 kg and +15 kg N ha−1 a−1) were applied in a factorial design, the latter simulating a fertilizing air pollution effect. Moderate warming led to increased productivity. Across the 4-year experimental period, the mean annual yield peaked at intermediate CSs (+43 % at +0.7 ∘C and +44 % at +1.8 ∘C), coinciding with ca. 50 % of days with less than 40 % soil moisture during the growing season. The yield increase was smaller at the lowest, warmest CS (+3.0 ∘C) but was still 12 % larger than at CSreference. These yield differences among CSs were well explained by differences in soil moisture and received thermal energy. Irrigation had a significant effect on yield (+16 %–19 %) in dry years, whereas atmospheric N deposition did not result in a significant yield response. We conclude that productivity of semi-natural, highly diverse subalpine grassland will increase in the near future. Despite increasingly limiting soil water content, plant growth will respond positively to up to +1.8 ∘C warming during the growing period, corresponding to +1.3 ∘C annual mean warming.

Description: How global gridded crop models (GGCMs) differ in their simulation of potential yield and reasons for those differences have never been assessed. The GGCM Intercomparison (GGCMI) offers a good framework for this assessment. Here, we built an emulator (called SMM for simple mechanistic model) of GGCMs based on generic and simplified formalism. The SMM equations describe crop phenology by a sum of growing degree days, canopy radiation absorption by the Beer–Lambert law, and its conversion into aboveground biomass by a radiation use efficiency (RUE). We fitted the parameters of this emulator against gridded aboveground maize biomass at the end of the growing season simulated by eight different GGCMs in a given year (2000). Our assumption is that the simple set of equations of SMM, after calibration, could reproduce the response of most GGCMs so that differences between GGCMs can be attributed to the parameters related to processes captured by the emulator. Despite huge differences between GGCMs, we show that if we fit both a parameter describing the thermal requirement for leaf emergence by adjusting its value to each grid-point in space, as done by GGCM modellers following the GGCMI protocol, and a GGCM-dependent globally uniform RUE, then the simple set of equations of the SMM emulator is sufficient to reproduce the spatial distribution of the original aboveground biomass simulated by most GGCMs. The grain filling is simulated in SMM by considering a fixed-in-time fraction of net primary productivity allocated to the grains (frac) once a threshold in leaves number (nthresh) is reached. Once calibrated, these two parameters allow for the capture of the relationship between potential yield and final aboveground biomass of each GGCM. It is particularly important as the divergence among GGCMs is larger for yield than for aboveground biomass. Thus, we showed that the divergence between GGCMs can be summarized by the differences in a few parameters. Our simple but mechanistic model could also be an interesting tool to test new developments in order to improve the simulation of potential yield at the global scale.

Description: Most precipitation from deep clouds over the continents and in the intertropical convergence zone is strongly influenced by the presence of ice crystals whose formation requires the presence of ice nucleating particles (INPs). Although there are a large number of INP sources, the ice nucleating abilities of aerosol particles originating from oceans, deserts, and wildfires are poorly described at tropical latitudes. To fill this gap in knowledge, the National Autonomous University of Mexico micro-orifice uniform deposit impactor droplet freezing technique (UNAM-MOUDI-DFT) was constructed to measure the ice nucleating activity of aerosol samples that were collected in Sisal and Mérida, Yucatán (Mexico) under the influence of cold fronts, biomass burning (BB), and African dust (AD) intrusions during five short-term field campaigns between January 2017 and July 2018. The three different aerosol types were distinguished by their physicochemical properties. Marine aerosol (MA), BB, and AD air masses were found to contain INPs; the highest concentrations were in AD (from 0.071 to 36.07 L−1 at temperatures between −18 and −27 ∘C), followed by MA (from 0.068 to 18.90 L−1 at temperatures between −15 and −28 ∘C) and BB (from 0.063 to 10.21 L−1 at temperatures between −20 and −27 ∘C). However, MA had the highest surface active site densities (ns) between −15 and −30 ∘C. Additionally, supermicron particles contributed more than 72 % of the total INP concentration independent of aerosol type.

Description: For the first time in the North China Plain (NCP) region, we investigated the seasonal variations in submicron particles (NR-PM1) and their chemical composition at a background mountainous site of Xinglong using an Aerodyne high-resolution time-of-flight aerosol mass spectrometer. The average concentration of NR-PM1 was highest in autumn (15.1 µg m−3) and lowest in summer (12.4 µg m−3), with a greater abundance of nitrate in spring (34 %), winter (31 %) and autumn (34 %) and elevated organics (40 %) and sulfate (38 %) in summer. PM1 in Xinglong showed higher acidity in summer and moderate acidity in spring, autumn and winter, with average pH values of 2.7±0.6, 4.2±0.7, 3.5±0.5 and 3.7±0.6, respectively, which is higher than those estimated in the United States and Europe. The size distribution of all PM1 species showed a consistent accumulation mode peaking at approximately 600–800 nm (dva), indicating a highly aged and internally mixed nature of the background aerosols, which was further supported by the source appointment results using positive matrix factorization and multilinear engine analysis. Significant contributions of aged secondary organic aerosol (SOA) in organic aerosol (OA) were resolved in all seasons (〉77 %), especially in summer. The oxidation state and the process of evolution of OAs in the four seasons were further investigated, and an enhanced carbon oxidation state (−0.45–0.10) and O/C (0.54–0.75) and OM/OC (1.86–2.13) ratios – compared with urban studies – were observed, with the highest oxidation state appearing in summer, likely because of the relatively stronger photochemical processing that dominated the formation processes of both less oxidized OA (LO-OOA) and more oxidized OA (MO-OOA). Aqueous-phase processing also contributed to the SOA formation and prevailed in winter, with the share to MO-OOA being more important than that to LO-OOA. In addition, regional transport also played an important role in the variations in SOA. Especially in summer, continuous increases in SOA concentration as a function of odd oxygen (Ox=O3+NO2) were found to be associated with the increases in wind speed. Furthermore, backward trajectory analysis showed that higher concentrations of submicron particles were associated with air masses transported short distances from the southern regions in all four seasons, while long-range transport from Inner Mongolia (western and northern regions) also contributed to summertime particulate pollution in the background areas of the NCP. Our results illustrate that the background particles in the NCP are influenced significantly by aging processes and regional transport, and the increased contribution of aerosol nitrate highlights how regional reductions in nitrogen oxide emissions are critical for remedying occurrence of nitrate-dominated haze events over the NCP.

Description: The sea surface microlayer (SML) represents the boundary layer at the air–sea interface. Microbial eukaryotes in the SML potentially influence air–sea gas exchange directly by taking up and producing gases and indirectly by excreting and degrading organic matter, which may modify the viscoelastic properties of the SML. However, little is known about the distribution of microbial eukaryotes in the SML. We studied the composition of the microbial community, transparent exopolymer particles and polysaccharides in the SML during the PEACETIME cruise along a west–east transect in the Mediterranean Sea, covering the western basin, Tyrrhenian Sea and Ionian Sea. At the stations located in the Ionian Sea, fungi – likely of continental origin and delivered by atmospheric deposition – were found in high relative abundances, making up a significant proportion of the sequences recovered. Concomitantly, bacterial and picophytoplankton counts decreased from west to east, while transparent exopolymer particle (TEP) abundance and total carbohydrate (TCHO) concentrations remained constant in all basins. Our results suggest that the presence of substrates for fungi, such as Cladosporium, known to take up phytoplankton-derived polysaccharides, in combination with decreased substrate competition by bacteria, might favor fungal dominance in the neuston of the Ionian Sea and other low-nutrient, low-chlorophyll (LNLC) regions.

Description: Ice-nucleating particles (INPs) influence the formation of ice crystals in clouds and many types of precipitation. This study reports unique properties of INPs collected from 42 precipitation samples in the Texas Panhandle region from June 2018 to July 2019. We used a cold stage instrument called the West Texas Cryogenic Refrigerator Applied to Freezing Test system to estimate INP concentrations per unit volume of air (nINP) through immersion freezing in our precipitation samples with our detection capability of 〉 0.006 INP L−1. A disdrometer was used for two purposes: (1) to characterize the ground-level precipitation type and (2) to measure the precipitation intensity as well as size of precipitating particles at the ground level during each precipitation event. While no clear seasonal variations of nINP values were apparent, the analysis of yearlong ground-level precipitation observation as well as INPs in the precipitation samples showed some INP variations, e.g., the highest and lowest nINP values at −25 ∘C both in the summer for hail-involved severe thunderstorm samples (3.0 to 1130 INP L−1), followed by the second lowest at the same temperature from one of our snow samples collected during the winter (3.2 INP L−1). Furthermore, we conducted bacteria community analyses using a subset of our precipitation samples to examine the presence of known biological INPs. In parallel, we also performed metagenomics characterization of the bacterial microbiome in suspended ambient dust samples collected at commercial open-lot livestock facilities (cattle feedyards hereafter) in the Texas Panhandle (i.e., the northernmost counties of Texas, also known as “West Texas”) to ascertain whether local cattle feedyards can act as a source of bioaerosol particles and/or INPs found in the precipitation samples. Some key bacterial phyla present in cattle feedyard samples appeared in precipitation samples. However, no known ice nucleation active species were detected in our samples. Overall, our results showed that cumulative nINP in our precipitation samples below −20 ∘C could be high in the samples collected while observing 〉 10 mm h−1 precipitation with notably large hydrometeor sizes and an implication of cattle feedyard bacteria inclusion.

Description: The comparative reactivity method (CRM) was developed more than a decade to measure OH reactivity (i.e., OH loss frequency) in both laboratory and field studies. However, accurate OH reactivity quantification remains challenging under real ambient conditions, especially for OH reactivity measurements in high-NOx (e.g., 〉 10 ppbv) environments, as ambient NO enhances the regeneration of OH radicals in the CRM reactor. To solve this problem, we design a new and improved CRM reactor (ICRM) and add NO into the system continuously so that the HO2 radical concentration is suppressed. We confirmed the appropriate level of NO by determining the maximum decrease in the pyrrole level caused by regenerated OH radicals from NO + HO2. RO2 radicals induced by volatile organic compounds (VOCs) in the ICRM reactor were also found to react with NO, which led to the regeneration of OH radicals and thus the underestimation of OH reactivity. This effect was quantified by the calibration of representative VOC species at different NO levels, and the correction coefficients obtained were used to correct the measured OH reactivity. All these efforts resulted in reducing the uncertainty of the NO-artifact correction by at least an order of magnitude compared to the original CRM system. Additionally, these technological improvements also considerably reduced the systematic errors from pyrrole photolysis in the original system. A new operation mode was proposed for the ICRM, which is able to avoid the interference resulting from OH radicals produced by photolysis of residual humidity and save time for ambient measurement. The ICRM system was employed in a field campaign to measure OH reactivity and performed well with ambient NO levels ranging from 0 to 50 ppbv, which are typically observed in the urban and suburban atmosphere.

Description: Precipitation is a crucial driver of hydrological processes. Ironically, a reliable characterization of its spatiotemporal variability is challenging. Ground-based rainfall measurement using rain gauges is more accurate. However, installing a dense gauging network to capture rainfall variability can be impractical. Satellite-based rainfall estimates (SREs) could be good alternatives, especially for data-scarce basins like in Ethiopia. However, SRE rainfall is plagued with uncertainties arising from many sources. The objective of this study was to evaluate the performance of the latest versions of several SRE products (i.e., CHIRPS2, IMERG6, TAMSAT3 and 3B42/3) for the Dhidhessa River Basin (DRB). Both statistical and hydrological modeling approaches were used for the performance evaluation. The Soil and Water Analysis Tool (SWAT) was used for hydrological simulations. The results showed that whereas all four SRE products are promising to estimate and detect rainfall for the DRB, the CHIRPS2 dataset performed the best at annual, seasonal and monthly timescales. The hydrological simulation-based evaluation showed that SWAT's calibration results are sensitive to the rainfall dataset. The hydrological response of the basin is found to be dominated by the subsurface processes, primarily by the groundwater flux. Overall, the study showed that both CHIRPS2 and IMERG6 products could be reliable rainfall data sources for the hydrological analysis of the DRB. Moreover, the climatic season in the DRB influences rainfall and streamflow estimation. Such information is important for rainfall estimation algorithm developers.

Description: Extreme events in the ocean severely impact marine organisms and ecosystems. Of particular concern are compound events, i.e., when conditions are extreme for multiple potential ocean ecosystem stressors such as temperature and chlorophyll. Yet, little is known about the occurrence, intensity, and duration of such compound high-temperature (a.k.a. marine heatwaves – MHWs) and low-chlorophyll (LChl) extreme events, whether their distributions have changed in the past decades, and what the potential drivers are. Here we use satellite-based sea surface temperature and chlorophyll concentration estimates to provide a first assessment of such compound extreme events. We reveal hotspots of compound MHW and LChl events in the equatorial Pacific, along the boundaries of the subtropical gyres, in the northern Indian Ocean, and around Antarctica. In these regions, compound events that typically last 1 week occur 3 to 7 times more often than expected under the assumption of independence between MHWs and LChl events. The occurrence of compound MHW and LChl events varies on seasonal to interannual timescales. At the seasonal timescale, most compound events occur in summer in both hemispheres. At the interannual timescale, the frequency of compound MHW and LChl events is strongly modulated by large-scale modes of natural climate variability such as the El Niño–Southern Oscillation, whose positive phase is associated with increased compound event occurrence in the eastern equatorial Pacific and in the Indian Ocean by a factor of up to 4. Our results provide a first understanding of where, when, and why compound MHW and LChl events occur. Further studies are needed to identify the exact physical and biological drivers of these potentially harmful events in the ocean and their evolution under global warming.

Description: Crop yield is reduced by heat and water stress and even more when these conditions co-occur. Yet, compound effects of air temperature and water availability on crop heat stress are poorly quantified. Existing crop models, by relying at least partially on empirical functions, cannot account for the feedbacks of plant traits and response to heat and water stress on canopy temperature. We developed a fully mechanistic model, coupling crop energy and water balances, to determine canopy temperature as a function of plant traits, stochastic environmental conditions, and irrigation applications. While general, the model was parameterized for wheat. Canopy temperature largely followed air temperature under well-watered conditions. But, when soil water potential was more negative than −0.14 MPa, further reductions in soil water availability led to a rapid rise in canopy temperature – up to 10 ∘C warmer than air at soil water potential of −0.62 MPa. More intermittent precipitation led to higher canopy temperatures and longer periods of potentially damaging crop canopy temperatures. Irrigation applications aimed at keeping crops under well-watered conditions could reduce canopy temperature but in most cases were unable to maintain it below the threshold temperature for potential heat damage; the benefits of irrigation in terms of reduction of canopy temperature decreased as average air temperature increased. Hence, irrigation is only a partial solution to adapt to warmer and drier climates.

Description: Ammonia (NH3) plays an important role in particulate matter formation; hence, its atmospheric level is relevant to human health and climate change. Due to different relative distributions of NH3 sources, concentrations of atmospheric NH3 may behave differently in urban and rural areas. However, few parallel long-term observations of NH3 exist to reveal the different behaviors of NH3 concentrations at urban and rural sites in a same region. In this study, online ammonia analyzers were used to continuously observe atmospheric NH3 concentrations at an urban site and a suburban site in Beijing from 13 January 2018 to 13 January 2019. The observed mixing ratio of NH3 averaged 21±14 ppb (range of 1.6–133 ppb) at the urban site and 22±15 ppb (range of 0.8–199 ppb) at the suburban site. The NH3 mixing ratios at the urban and suburban sites exhibited similar seasonal variations, with high values in summer and spring and low values in autumn and winter. The hourly mean NH3 mixing ratios at the urban site were highly correlated (R=0.849, P

Description: 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.

Description: Increasing evidence from experimental studies suggests that the losses of semi-volatile vapors to chamber walls could be responsible for the underestimation of organic aerosol (OA) in air quality models that use parameters obtained from chamber experiments. In this study, a box model with a volatility basis set (VBS) scheme was developed, and the secondary organic aerosol (SOA) yields with vapor wall loss correction were optimized by a genetic algorithm based on advanced chamber experimental data for biomass burning. The vapor wall loss correction increases the SOA yields by a factor of 1.9–4.9 and leads to better agreement with measured OA for 14 chamber experiments under different temperatures and emission loads. To investigate the influence of vapor wall loss correction on regional OA simulations, the optimized parameterizations (SOA yields, emissions of intermediate-volatility organic compounds from biomass burning, and enthalpy of vaporization) were implemented in the regional air quality model CAMx (Comprehensive Air Quality Model with extensions). The model results from the VBS schemes with standard (VBS_BASE) and vapor-wall-loss-corrected parameters (VBS_WLS), as well as the traditional two-product approach, were compared and evaluated by OA measurements from five Aerodyne aerosol chemical speciation monitor (ACSM) or aerosol mass spectrometer (AMS) stations in the winter of 2011. An additional reference scenario, VBS_noWLS, was also developed using the same parameterization as VBS_WLS except for the SOA yields, which were optimized by assuming there is no vapor wall loss. The VBS_WLS generally shows the best performance for predicting OA among all OA schemes and reduces the mean fractional bias from −72.9 % (VBS_BASE) to −1.6 % for the winter OA. In Europe, the VBS_WLS produces the highest domain average OA in winter (2.3 µg m−3), which is 106.6 % and 26.2 % higher than VBS_BASE and VBS_noWLS, respectively. Compared to VBS_noWLS, VBS_WLS leads to an increase in SOA by up to ∼80 % (in the Balkans). VBS_WLS also leads to better agreement between the modeled SOA fraction in OA (fSOA) and the estimated values in the literature. The substantial influence of vapor wall loss correction on modeled OA in Europe highlights the importance of further improvements in parameterizations based on laboratory studies for a wider range of chamber conditions and field observations with higher spatial and temporal coverage.

Description: The ongoing phenomenon of climate change is leading to an upsurge in the number of extreme events. Territories must adapt to these modifications in order to protect their populations and the properties present in coastal areas. The adaptation of coastal areas also aims to make them more resilient to future events. In this article, we examine two strategies for adapting to coastal risks: holding the coastal line through hard constructions such as seawalls or ripraps and the managed retreat of activities and populations to a part of the territory not exposed to hazards. In France, these approaches are financed by a solidarity insurance system at the national level as well as local taxes. These solidarity systems aim to compensate the affected populations and finance implementation of the strategies chosen by local authorities. However, the French mainland coast generally attracts affluent residents, the price of land being higher than inland. This situation induces the presence of inequalities in these territories, inequalities which can be maintained or reinforced in the short and medium term when a defense strategy based on hard constructions is implemented. In such a trajectory, it appears that these territories would be less resilient in the long term because of the maintenance costs of the structures and the uncertainties relating to the hazards (submersion, rising sea levels, erosion). Conversely, with a managed-retreat strategy, inequalities would instead be done away with since property and populations would no longer be exposed to hazards, which would cost society less and would lead these territories towards greater resilience in the long term. Only one social group would be strongly impacted by this strategy in the short term when they are subjected to a managed retreat to another part of the territory.

Description: G2DC-PL+, a gridded 2 km daily climate dataset for the union of the Polish territory and the Vistula and Odra basins, is an update and extension of the CHASE-PL Forcing Data – Gridded Daily Precipitation and Temperature Dataset – 5 km (CPLFD-GDPT5). The latter was the first publicly available, high-resolution climate forcing dataset in Poland, used for a range of purposes including hydrological modelling and bias correction of climate projections. While the spatial coverage of the new dataset remained the same, it has undergone several major changes: (1) the time coverage was increased from 1951–2013 to 1951–2019; (2) its spatial resolution increased from 5 to 2 km; (3) the number of stations used for interpolation of temperature and precipitation approximately doubled; and (4) in addition to precipitation and temperature, the dataset consists of relative humidity and wind speed data. The main purpose for developing this product was the need for long-term areal climate data for earth-system modelling, and particularly hydrological modelling. Geostatistical methods (kriging) were used for interpolation of the studied climate variables. The kriging cross-validation revealed improved performance for precipitation compared to the original dataset expressed by the median of the root mean squared errors standardized by standard deviation of observations (0.59 vs. 0.79). Kriging errors were negatively correlated with station density only for the period 1951–1970. Values of the root mean squared error normalized to the standard deviation (RMSEsd) were equal to 0.52 and 0.4 for minimum and maximum temperature, respectively, suggesting a small to moderate improvement over the original dataset. Relative humidity and wind speed exhibited lower performance, with median RMSEsd equal to 0.82 and 0.87, respectively. The dataset is openly available from the 4TU Centre for Research Data at https://doi.org/10.4121/uuid:a3bed3b8-e22a-4b68-8d75-7b87109c9feb (Piniewski et al., 2020).

Description: Continental shelf sediments are places of both rapid organic carbon turnover and accumulation, while at the same time increasingly subjected to human-induced disturbances. Recent research suggests that shelf sediments might have a role to play as a natural climate solution, e.g. by storing organic carbon if left undisturbed from anthropogenic activity. However, we have an incomplete understanding about the centres of organic carbon accumulation and storage on continental shelves. To better constrain the rate of accumulation and the mass of organic carbon that is stored in sediments, we developed and applied a spatial modelling framework that allows us to estimate those quantities from sparse observations and predictor variables known or suspected to influence the spatial patterns of these parameters. This paper presents spatial distribution patterns of organic carbon densities and accumulation rates in the North Sea and Skagerrak. We found that organic carbon stocks and accumulation rates are highest in the Norwegian Trough, while large parts of the North Sea are characterised by low stocks and zero net accumulation. The total stock of organic carbon that is stored in the upper 0.1 m of sediments amounted to 230.5 ± 134.5 Tg C, of which approximately 26 % is stored in the Norwegian Trough. Rates of organic carbon accumulation in the Norwegian Trough are comparable with those reported from nearby fjords. We provide baseline datasets that could be used in marine management, e.g. for the establishment of “carbon protection zones”. Additionally, we highlight the complex nature of continental shelves with zones of rapid carbon cycling and accumulation juxtaposed, which will require further detailed and spatially explicit analyses to constrain sedimentary organic carbon stocks and accumulation rates globally.

Description: The MARINE (Model of Anticipation of Runoff and INundations for Extreme events) hydrological model is a distributed model dedicated to flash flood simulation. Recent developments of the MARINE model are explored in this work. On one hand, transfers of water through the subsurface, formerly relying on water height, now take place in a homogeneous soil column based on the soil saturation degree (SSF model). On the other hand, the soil column is divided into two layers, which represent, respectively, the upper soil layer and the deep weathered rocks (SSF–DWF model). The aim of the present work is to assess the accuracy of these new representations for the simulation of soil moisture during flash flood events. An exploration of the various products available in the literature for soil moisture estimation is performed. The efficiency of the models for soil saturation degree simulation is estimated with respect to several products either at the local scale or spatially distributed: (i) the gridded soil moisture product provided by the operational modeling chain SAFRAN-ISBA-MODCOU; (ii) the gridded soil moisture product provided by the LDAS-Monde assimilation chain, which is based on the ISBA-A-gs land surface model and assimilating satellite derived data; (iii) the upper soil water content hourly measurements taken from the SMOSMANIA observation network; and (iv) the Soil Water Index provided by the Copernicus Global Land Service (CGLS), which is derived from Sentinel-1 C-SAR and ASCAT satellite data. The case study is performed over two French Mediterranean catchments impacted by flash flood events over the 2017–2019 period. The local comparison of the MARINE outputs with the SMOSMANIA measurements, as well as the comparison at the basin scale of the MARINE outputs with the gridded LDAS-Monde and CGLS data, lead to the following conclusion: both the dynamics and the amplitudes of the soil saturation degree simulated with the SSF and SSF–DWF models are better correlated with both the SMOSMANIA measurements and the LDAS-Monde data than the outputs of the base model. Finally, the soil saturation degree simulated by the two-layers model for the deep layer is compared to the soil saturation degree provided by the LDAS-Monde product at corresponding depths. In conclusion, the developments presented for the representation of subsurface flow in the MARINE model enhance the soil saturation degree simulation during flash floods with respect to both gridded data and local soil moisture measurements.

Description: Aerosol sulfate is a major component of submicron particulate matter (PM1). Sulfate can be present as inorganic (mainly ammonium sulfate, AS) or organosulfate (OS). Although OS is thought to be a smaller fraction of total sulfate in most cases, recent literature argues that this may not be the case in more polluted environments. Aerodyne aerosol mass spectrometers (AMSs) measure total submicron sulfate, but it has been difficult to apportion AS vs. OS as the detected ion fragments are similar. Recently, two new methods have been proposed to quantify OS separately from AS with AMS data. We use observations collected during several airborne field campaigns covering a wide range of sources and air mass ages (spanning the continental US, marine remote troposphere, and Korea) and targeted laboratory experiments to investigate the performance and validity of the proposed OS methods. Four chemical regimes are defined to categorize the factors impacting sulfate fragmentation. In polluted areas with high ammonium nitrate concentrations and in remote areas with high aerosol acidity, the decomposition and fragmentation of sulfate in the AMS is influenced by multiple complex effects, and estimation of OS does not seem possible with current methods. In regions with lower acidity (pH 〉 0) and ammonium nitrate (fraction of total mass

Description: The early Permian Meisenheim Formation of the Saar–Nahe Basin (Germany) is famous for its richness in vertebrate fossils, among which the temnospondyls were present with microvores and fish-eating apex predators. The latter trophic guild was occupied exclusively by the genus Sclerocephalus in that basin within a long time interval up to M8, whereas in M9, a new taxon, Glanochthon lellbachae, appeared. This taxon is defined by (1) a preorbital region 1.8–2.0 times as long as the postorbital skull table, (2) dermal ornament with tall radial ridges, (3) a prefrontal anteriorly wider with straight lateral margin, (4) a squamosal posteriorly only half as wide as the quadratojugal, (5) phalanges of manus and pes long and gracile, (6) carpals unossified in adults, and (7) tail substantially longer than skull and trunk combined. Phylogenetic analysis finds that G. lellbachae forms the basal sister taxon of the stratigraphically younger G. angusta and G. latirostre and that this clade nests within the paraphyletic taxon Sclerocephalus, with S. nobilis forming the sister taxon of the genus Glanochthon (urn:lsid:zoobank.org:act:3038F794-17B9-4FCA-B241-CCC3F4423651; registration date: 15 March 2021).

Description: The Brewer ozone spectrophotometer (the Brewer) was designed at Environment and Climate Change Canada (ECCC) in the 1970s to make accurate automated total ozone column measurements. Since the 1980s, the Brewer instrument has become a World Meteorological Organization (WMO) Global Atmosphere Watch (GAW) standard ozone monitoring instrument. Now, more than 230 Brewers have been produced. To assure the quality of the Brewer measurements, a calibration chain is maintained, i.e., first, the reference instruments are independently absolutely calibrated, and then the calibration is transferred from the reference instrument to the travelling standard, and subsequently from the travelling standard to field instruments. ECCC has maintained the world Brewer reference instruments since the 1980s to provide transferable calibration to field instruments at monitoring sites. Three single-monochromator (Mark II) type instruments (serial numbers 008, 014, and 015) formed this world Brewer reference triad (BrT) and started their service in Toronto, Canada, in 1984. In the 1990s, the Mark III type Brewer (known as the double Brewer) was developed, which has two monochromators to reduce the internal instrumental stray light. The double-Brewer world reference triad (BrT-D) was formed in 2013 (serial numbers 145, 187 and 191), co-located with the BrT. The first assessment of the BrT's performance was made in 2005, covering the period between 1984 and 2004 (Fioletov et al., 2005). The current work provides an updated assessment of the BrT's performance (from 1999 to 2019) and the first comprehensive assessment of the BrT-D. The random uncertainties of individual reference instruments are within the WMO/GAW requirement of 1 % (WMO, 2001): 0.49 % and 0.42 % for BrT and BrT-D, respectively, as estimated in this study. The long-term stability of the reference instruments is also evaluated in terms of uncertainties of the key instrument characteristics: the extraterrestrial calibration constant (ETC) and effective ozone absorption coefficients (both having an effect of less than 2 % on total column ozone). Measurements from a ground-based instrument (Pandora spectrometer), satellites (11 datasets, including the most recent high-resolution satellite, TROPOspheric Monitoring Instrument), and reanalysis model (the second Modern-Era Retrospective analysis for Research and Applications, MERRA-2) are used to further assess the performance of world Brewer reference instruments and to provide a context for the requirements of stratospheric ozone observations during the last two decades.

Description: Atmospheric submicrometer aerosols have a great effect on air quality and human health, while their formation and evolution processes are still not fully understood. Herein, the crucial role of atmospheric oxidation capacity, as characterized by OH exposure dose in the formation and evolution of secondary submicrometer aerosols, was systematically investigated based on a highly time-resolved chemical characterization of PM1 in a southern suburb of Beijing in summertime from 25 July to 21 August 2019. The averaged concentration of PM1 was 19.3 ± 11.3 µg m−3, and nearly half (48.3 %) of the mass was organic aerosols (OAs) during the observation period. The equivalent photochemical age (ta) estimated from the ratios of toluene to benzene was applied to characterize the OH exposure dose of the air mass, in which an observation period with the similar sources and minimal influence of fresh emission was adopted. The relationships of non-refractory PM1 species, OA factors (i.e., one hydrocarbon-like and three oxygenated organic aerosol factors) and elemental compositions (e.g., H∕C, O∕C, N∕C, S∕C, OM∕OC, and OSc) to ta were analyzed in detail. It was found that higher PM1 concentration accompanied longer ta, with an average increase rate of 0.8 µgm-3h-1. Meanwhile, the formation of sulfate and more oxidized oxygenated OA were most sensitive to the increase in ta, and their contributions to PM1 were enhanced from 22 % to 28 % and from 29 % to 48 %, respectively, as ta increased. In addition, OSc and the ratios of O∕C and OM∕OC increased with the increase in ta. These results indicated that photochemical aging is a key factor leading to the evolution of OA and the increase in PM1 in summertime.

Description: Terrestrial and airborne laser scanning and structure from motion techniques have emerged as viable methods to map snow depths. While these systems have advanced snow hydrology, these techniques have noted limitations in either horizontal or vertical resolution. Lidar on an unpiloted aerial vehicle (UAV) is another potential method to observe field- and slope-scale variations at the vertical resolutions needed to resolve local variations in snowpack depth and to quantify snow depth when snowpacks are shallow. This paper provides some of the earliest snow depth mapping results on the landscape scale that were measured using lidar on a UAV. The system, which uses modest-cost, commercially available components, was assessed in a mixed deciduous and coniferous forest and open field for a thin snowpack (

Description: The Gridpoint Statistical Interpolation data assimilation (DA) system was developed for the four size bin sectional Model for Simulating Aerosol Interactions and Chemistry (MOSAIC) aerosol mechanism in the Weather Research and Forecasting-Chemistry (WRF-Chem) model. The forward and tangent linear operators for the aerosol optical depth (AOD) analysis were derived from WRF-Chem aerosol optical code. We applied three-dimensional variational DA to assimilate the multi-wavelength AOD, ambient aerosol scattering coefficient, and aerosol absorption coefficient, measured by the sun–sky photometer, nephelometer, and aethalometer, respectively. These measurements were undertaken during a dust observation field campaign at Kashi in northwestern China in April 2019. The results showed that the DA analyses decreased the model aerosols' low biases; however, it had some deficiencies. Assimilating the surface particle concentration increased the coarse particles in the dust episodes, but AOD and the coefficients for aerosol scattering and absorption were still lower than those observed. Assimilating aerosol scattering coefficient separately from AOD improved the two optical quantities. However, it caused an overestimation of the particle concentrations at the surface. Assimilating the aerosol absorption coefficient yielded the highest positive bias in the surface particle concentration, aerosol scattering coefficient, and AOD. The positive biases in the DA analysis were caused by the forward operator underestimating aerosol mass scattering and absorption efficiency. As compensation, the DA system increased particle concentrations excessively to fit the observed optical values. The best overall improvements were obtained from the simultaneous assimilation of the surface particle concentration and AOD. The assimilation did not substantially change the aerosol chemical fractions. After DA, the clear-sky aerosol radiative forcing at Kashi was −10.4 W m−2 at the top of the atmosphere, which was 55 % higher than the radiative forcing value before DA.

Description: The Southern Ocean has a profound impact on the Earth's climate system. Its strong winds, intense currents, and fierce waves are critical components of the air–sea interface and contribute to absorbing, storing, and releasing heat, moisture, gases, and momentum. Owing to its remoteness and harsh environment, this region is significantly undersampled, hampering the validation of prediction models and large-scale observations from satellite sensors. Here, an unprecedented data set of simultaneous observations of winds, surface currents, and ocean waves is presented, to address the scarcity of in situ observations in the region – https://doi.org/10.26179/5ed0a30aaf764 (Alberello et al., 2020c) and https://doi.org/10.26179/5e9d038c396f2 (Derkani et al., 2020). Records were acquired underway during the Antarctic Circumnavigation Expedition (ACE), which went around the Southern Ocean from December 2016 to March 2017 (Austral summer). Observations were obtained with the wave and surface current monitoring system WaMoS-II, which scanned the ocean surface around the vessel using marine radars. Measurements were assessed for quality control and compared against available satellite observations. The data set is the most extensive and comprehensive collection of observations of surface processes for the Southern Ocean and is intended to underpin improvements of wave prediction models around Antarctica and research of air–sea interaction processes, including gas exchange and dynamics of sea spray aerosol particles. The data set has further potentials to support theoretical and numerical research on lower atmosphere, air–sea interface, and upper-ocean processes.

Description: Assessing the risk of a historical-level flood is essential for regional flood protection and resilience establishment. However, due to the limited spatiotemporal coverage of observations, the impact assessment relies on model simulations and is thus subject to uncertainties from cascade physical processes. This study assesses the flood hazard map with uncertainties subject to different combinations of runoff inputs, variables for flood frequency analysis and fitting distributions based on estimations by the CaMa-Flood global hydrodynamic model. Our results show that deviation in the runoff inputs is the most influential source of uncertainties in the estimated flooded water depth and inundation area, contributing more than 80 % of the total uncertainties investigated in this study. Global and regional inundation maps for floods with 1-in-100 year return periods show large uncertainty values but small uncertainty ratios for river channels and lakes, while the opposite results are found for dry zones and mountainous regions. This uncertainty is a result of increasing variation at tails among various fitting distributions. In addition, the uncertainty between selected variables is limited but increases from the regular period to the rarer floods, both for the water depth at points and for inundation area over regions. The uncertainties in inundation area also lead to uncertainties in estimating the population and economy exposure to the floods. In total, inundation accounts for 9.1 % [8.1 %–10.3 %] of the land area for a 1-in-100 year flood, leading to 13.4 % [12.1 %–15 %] of population exposure and 13.1 % [11.8 %–14.7 %] of economic exposure for the globe. The flood exposure and uncertainties vary by continent and the results in Africa have the largest uncertainty, probably due to the limited observations to constrain runoff simulations, indicating a necessity to improve the performance of different hydrological models especially for data-limited regions.

Description: The evaluation of models in general is a nontrivial task and can, due to epistemological and practical reasons, never be considered complete. Due to this incompleteness, a model may yield correct results for the wrong reasons, i.e., via a different chain of processes than found in observations. While guidelines and strategies exist in the atmospheric sciences to maximize the chances that models are correct for the right reasons, these are mostly applicable to full physics models, such as numerical weather prediction models. The Intermediate Complexity Atmospheric Research (ICAR) model is an atmospheric model employing linear mountain wave theory to represent the wind field. In this wind field, atmospheric quantities such as temperature and moisture are advected and a microphysics scheme is applied to represent the formation of clouds and precipitation. This study conducts an in-depth process-based evaluation of ICAR, employing idealized simulations to increase the understanding of the model and develop recommendations to maximize the probability that its results are correct for the right reasons. To contrast the obtained results from the linear-theory-based ICAR model to a full physics model, idealized simulations with the Weather Research and Forecasting (WRF) model are conducted. The impact of the developed recommendations is then demonstrated with a case study for the South Island of New Zealand. The results of this investigation suggest three modifications to improve different aspects of ICAR simulations. The representation of the wind field within the domain improves when the dry and the moist Brunt–Väisälä frequencies are calculated in accordance with linear mountain wave theory from the unperturbed base state rather than from the time-dependent perturbed atmosphere. Imposing boundary conditions at the upper boundary that are different to the standard zero-gradient boundary condition is shown to reduce errors in the potential temperature and water vapor fields. Furthermore, the results show that there is a lowest possible model top elevation that should not be undercut to avoid influences of the model top on cloud and precipitation processes within the domain. The method to determine the lowest model top elevation is applied to both the idealized simulations and the real terrain case study. Notable differences between the ICAR and WRF simulations are observed across all investigated quantities such as the wind field, water vapor and hydrometeor distributions, and the distribution of precipitation. The case study indicates that the precipitation maximum calculated by the ICAR simulation employing the developed recommendations is spatially shifted upwind in comparison to an unmodified version of ICAR. The cause for the shift is found in influences of the model top on cloud formation and precipitation processes in the ICAR simulations. Furthermore, the results show that when model skill is evaluated from statistical metrics based on comparisons to surface observations only, such an analysis may not reflect the skill of the model in capturing atmospheric processes like gravity waves and cloud formation.

Description: Complex microbial communities facilitate iron and methane transformations in anoxic methanic sediments of freshwater lakes, such as Lake Kinneret (the Sea of Galilee, Israel). The phylogenetic and functional diversity of these consortia are not fully understood, and it is not clear which lineages perform iron reduction and anaerobic oxidation of methane (AOM). Here, we investigated microbial communities from both natural Lake Kinneret iron-rich methanic sediments (〉20 cm depth) and iron-amended slurry incubations from this zone using metagenomics, focusing on functions associated with iron reduction and methane cycling. Analyses of the phylogenetic and functional diversity indicate that consortia of archaea (mainly Bathyarchaeia, Methanomicrobia, Thermoplasmata, and Thermococci) and bacteria (mainly Chloroflexi (Chloroflexota), Nitrospirae (Nitrospirota), and Proteobacteria) perform key metabolic reactions such as amino acid uptake and dissimilation, organic matter fermentation, and methanogenesis. The Deltaproteobacteria, especially Desulfuromondales (Desulfuromonadota), have the potential to transfer electrons extracellularly either to iron mineral particles or to microbial syntrophs, including methanogens. This is likely via transmembrane cytochromes, outer-membrane hexaheme c-type cytochrome (OmcS) in particular, or pilin monomers (PilA), all of which were attributed to this lineage. Bona fide anaerobic oxidizers of methane (ANME) and denitrifying methanotrophs Methylomirabilia (NC10) may mediate AOM in these methanogenic sediments; however we also consider the role of methanogens in active AOM or back flux of methanogenesis. Putative aerobes, such as methane-oxidizing bacteria Methylomonas and their methylotrophic syntrophs Methylotenera, are found among the anaerobic lineages in Lake Kinneret iron-amended slurries and are also involved in the oxidation of methane or its intermediates, as suggested previously. We propose a reaction model for the metabolic interactions in these sediments, linking the potential players that interact via intricate metabolic tradeoffs and direct electron transfer between species. Our results highlight the metabolic complexity of microbial communities in an energy-limited environment, where aerobe and anaerobe communities may co-exist and facilitate AOM as one strategy for survival.

Description: In recent years, ammonia emissions have been continuously increasing, being almost 4 times higher than in the 20th century. Although an important species, as its use as a fertilizer sustains human living, ammonia has major consequences for both humans and the environment because of its reactive gas-phase chemistry that makes it easily convertible to particles. Despite its pronounced importance, ammonia emissions are highly uncertain in most emission inventories. However, the great development of satellite remote sensing nowadays provides the opportunity for more targeted research on constraining ammonia emissions. Here, we used satellite measurements to calculate global ammonia emissions over the period 2008–2017. Then, the calculated ammonia emissions were fed to a chemistry transport model, and ammonia concentrations were simulated for the period 2008–2017. The simulated concentrations of ammonia were compared with ground measurements from Europe, North America and Southeastern Asia, as well as with satellite measurements. The satellite-constrained ammonia emissions represent global concentrations more accurately than state-of-the-art emissions. Calculated fluxes in the North China Plain were seen to be more increased after 2015, which is not due to emission changes but due to changes in sulfate emissions that resulted in less ammonia neutralization and hence in larger atmospheric loads. Emissions over Europe were also twice as much as those in traditional datasets with dominant sources being industrial and agricultural applications. Four hot-spot regions of high ammonia emissions were seen in North America, which are characterized by high agricultural activity, such as animal breeding, animal farms and agricultural practices. South America is dominated by ammonia emissions from biomass burning, which causes a strong seasonality. In Southeastern Asia, ammonia emissions from fertilizer plants in China, Pakistan, India and Indonesia are the most important, while a strong seasonality was observed with a spring and late summer peak due to rice and wheat cultivation. Measurements of ammonia surface concentrations were better reproduced with satellite-constrained emissions, such as measurements from CrIS (Cross-track Infrared Sounder).

Description: In this study, we present an estimate of the gravity signal of the slabs beneath the Alpine mountain belt. Estimates of the gravity effect of the subducting slabs are often omitted or simplified in crustal-scale models. The related signal is calculated here for alternative slab configurations at near-surface height and at a satellite altitude of 225 km. We apply three different modelling approaches in order to estimate the gravity signal from the subducting slab segments: (i) direct conversion of upper mantle seismic velocities to density distribution, which are then forward calculated to obtain the gravity signal; (ii) definition of slab geometries based on seismic crustal thickness and high-resolution upper mantle tomography for two competing slab configurations – the geometries are then forward calculated by assigning a constant density contrast and slab thickness; (iii) accounting for compositional and thermal variations with depth within the predefined slab geometry. Forward calculations predict a gravity signal of up to 40 mGal for the Alpine slab configuration. Significant differences in the gravity anomaly patterns are visible for different slab geometries in the near-surface gravity field. However, different contributing slab segments are not easily separated, especially at satellite altitude. Our results demonstrate that future studies addressing the lithospheric structure of the Alps should have to account for the subducting slabs in order to provide a meaningful representation of the geodynamic complex Alpine area.

Description: Boreal forest soils are globally an important sink for methane (CH4), while these soils are also capable of emitting CH4 under favourable conditions. Soil wetness is a well-known driver of CH4 flux, and the wetness can be estimated with several terrain indices developed for the purpose. The aim of this study was to quantify the spatial variability of the forest floor CH4 flux with a topography-based upscaling method connecting the flux with its driving factors. We conducted spatially extensive forest floor CH4 flux and soil moisture measurements, complemented by ground vegetation classification, in a boreal pine forest. We then modelled the soil moisture with a random forest model using digital-elevation-model-derived topographic indices, based on which we upscaled the forest floor CH4 flux. The modelling was performed for two seasons: May–July and August–October. Additionally, we evaluated the number of flux measurement points needed to get an accurate estimate of the flux at the whole study site merely by averaging. Our results demonstrate high spatial heterogeneity in the forest floor CH4 flux resulting from the soil moisture variability as well as from the related ground vegetation. The mean measured CH4 flux at the sample points was −5.07 µmol m−2 h−1 in May–July and −8.67 µmol m−2 h−1 in August–October, while the modelled flux for the whole area was −7.42 and −9.91 µmol m−2 h−1 for the two seasons, respectively. The spatial variability in the soil moisture and consequently in the CH4 flux was higher in the early summer (modelled range from −12.3 to 6.19 µmol m−2 h−1) compared to the autumn period (range from −14.6 to −2.12 µmol m−2 h−1), and overall the CH4 uptake rate was higher in autumn compared to early summer. In the early summer there were patches emitting high amounts of CH4; however, these wet patches got drier and smaller in size towards the autumn, changing their dynamics to CH4 uptake. The mean values of the measured and modelled CH4 fluxes for the sample point locations were similar, indicating that the model was able to reproduce the results. For the whole site, upscaling predicted stronger CH4 uptake compared to simply averaging over the sample points. The results highlight the small-scale spatial variability of the boreal forest floor CH4 flux and the importance of soil chamber placement in order to obtain spatially representative CH4 flux results. To predict the CH4 fluxes over large areas more reliably, the locations of the sample points should be selected based on the spatial variability of the driving parameters, in addition to linking the measured fluxes with the parameters.

Description: Enhancing resilience is critical for coastal urban systems to cope with and minimize flood disaster risks. This issue is certainly more important in Africa, where the increase in flood frequency is a significant concern for many areas. In this context, urban planners need accurate approaches to set up a standard for measuring the resilience to floods. In Morocco, this issue is still not fully covered by the scientific community despite the obvious need for a new approach adapted to local conditions. This study applied a composite index and geographic-information-system approach to measure and map resilience to floods in three northern coastal municipalities. The approach is also based on a linear ranking of resilience parameters, offering a more optimal classification of spatial resilience variation. The results allowed us to identify specific areas with different resilience levels and revealed the relationship between urban dimensions and the flood resilience degree. This approach provides an efficient decision-support tool to facilitate flood risk management, especially in terms of prioritizing protective actions.

Description: The power generation and loading dynamic responses of a 2.2 m diameter horizontal axis wind turbine (HAWT) under some of the IEC 61400-1 transient extreme operational conditions, more specifically extreme wind shears (EWSs) and extreme operational gust (EOG), that were reproduced at the WindEEE Dome at Western University were investigated. The global forces were measured by a multi-axis force balance at the HAWT tower base. The unsteady horizontal shear induced a significant yaw moment on the rotor with a dynamic similar to that of the extreme event without affecting the power generation. The EOG severely affected all the performance parameters of the turbine.

Description: Compound dry and hot conditions frequently cause large impacts on ecosystems and societies worldwide. A suite of indices is available for the assessment of droughts and heatwaves, yet there is no index available for incorporating the joint variability of dry and hot conditions at the sub-monthly scale. Here we introduce a daily-scale index, called the standardized compound drought and heat index (SCDHI), to assess compound dry-hot conditions. The SCDHI is based on a daily drought index (the standardized antecedent precipitation evapotranspiration index – SAPEI), the daily-scale standardized temperature index (STI), and a joint probability distribution method. The new index is verified against real-world compound dry and hot events and associated observed vegetation impacts in China. The SCDHI can not only capture compound dry and hot events at both monthly and sub-monthly scales, but is also a good indicator for associated vegetation impacts. Using the SCDHI, we quantify the frequency, severity, duration, and intensity of compound dry-hot events during the historical period and assess the ability of climate models to reproduce these characteristics in China. We find that compound events whose severity is at least light and which last longer than 2 weeks generally persisted for 20–35 d in China. Southern China suffered from compound events most frequently, and the most severe compound events were mainly detected in this region. Climate models generally overestimate the frequency, duration, severity, and intensity of compound events in China, especially for western regions, which can be attributed to a too strong dependence between the SAPEI and STI in those models. The SCDHI provides a new tool to quantify sub-monthly characteristics of compound dry and hot events and to monitor their initiation, development, and decay. This is important information for decision-makers and stakeholders to release early and timely warnings.

Description: Sugar compounds (SCs) are major water-soluble constituents in atmospheric aerosols. In this study, we investigated their molecular compositions and abundances in the northern receptor site (Mangshan) of Beijing, China, to better understand the contributions from biogenic and anthropogenic sources using a gas chromatography–mass spectrometry technique. The sampling site receives anthropogenic air mass transported from Beijing by southerly winds, while northerly winds transport relatively clean air mass from the forest areas. Day- and nighttime variations were analyzed for anhydrosugars, primary sugars, and sugar alcohols in autumn 2007. We found that biomass burning (BB) tracers were more abundant at nighttime than daytime, while other SCs showed different diurnal variations. Levoglucosan was found to be dominant sugar among the SCs observed, indicating an intense influence of local BB for cooking and space heating at the surroundings of the Mangshan site. The high levels of arabitol and mannitol in daytime suggest a significant contribution of locally emitted fungal spores and long-range-transported bioaerosols from the Beijing area. The plant emissions from Mangshan forest park significantly control the diurnal variations of glucose, fructose, and mannitol. The meteorological parameters (relative humidity, temperature, and rainfall) significantly affect the concentrations and diurnal variations of SCs. Sucrose (pollen tracer) showed a clear diurnal variation, peaking in the daytime due to higher ambient temperature and wind speed, which influences the pollen release from the forest plants. We found the contribution of trehalose from soil dust in daytime, while microbial and fungal spores were responsible for nighttime. Anhydrosugar and primary sugars are prime carbon sources of the Mangshan aerosols. The high ratios of levoglucosan in organic carbon and water-soluble organic carbon at nighttime suggest a significant contribution of BB to organic aerosols at night. Levoglucosan / mannosan ratios demonstrate that low-temperature burning of hardwood is dominant in Mangshan. The positive matrix factorization analysis concluded that forest vegetation, fungal species, and local BB are the significant sources of SCs.

Description: As critical precursors of ozone (O3) and secondary organic aerosols, volatile organic compounds (VOCs) play a vital role in air quality, human health, and climate change. In this study, a campaign of comprehensive field observations and VOC grid sampling was conducted in Xi'an, China, from 20 June to 20 July 2019 to identify the spatiotemporal concentration levels, sources, and secondary transformation potential of VOCs. During the observation period, the average VOC concentrations at the Chanba (CB), Di Huan Suo (DHS), Qinling (QL), and gridded sampling sites were 27.8 ± 8.9, 33.8 ± 10.5, 15.5 ± 5.8, and 29.1 ± 8.4 ppb, respectively. Vehicle exhaust was the primary source of VOC emissions in Xi'an, and the contributions of vehicle exhaust to VOCs at the CB, DHS, and QL sites were 41.3 %, 30.6 %, and 23.6 %–41.4 %, respectively. While industrial emissions were the second-largest source of VOCs in urban areas, contributions from aging sources were high in rural areas. High potential source contribution function values primarily appeared in eastern and southern Xi'an near the sampling site, which indicates that Xi'an exhibits a strong local VOC source. Moreover, alkenes, aromatics, and oxygenated VOCs played a dominant role in secondary transformation, which is a major concern in reducing O3 pollution in Xi'an.

Description: This paper discusses a systematic vehicle design process in which light weight is taken as the vehicle design objective, and the designed frame is analyzed in detail. The load condition of a vehicle under different circumstances is calculated according to the distances from the front and rear wheels to the centroid position. The stress on the components in the condition is analyzed by finite element analysis, the steering geometry of the vehicle is analyzed, and the vehicle's turning angle and radius are designed. The displacement of the vehicle under a load is calculated by rigidity analysis to determine the stability of the vehicle in motion. The experimental modal analysis of the real frame and the finite element method are verified mutually for the electric vehicle body-in-white (BIW) manufacturing process to determine the consistency of model formation and the real frame. In terms of the circuit design, we used no-fuse switches and fuses to provide overcurrent protection for the main power supply, and the chip is combined with an optically coupled circuit and current sensor, which is driven by a restriction controller for protection. Moreover, a solid-state relay (SSR) is used for current protection and for controlling the forward/reverse rotation of the motor.

Description: The integration of surface and subsurface geoscience data is critical for efficient and effective mineral exploration and mining. Publicly accessible data sets to evaluate the various geoscience analytical tools and their effectiveness for characterisation of mineral assemblages and lithologies or discrimination of ore from waste are however scarce. The open-access Rocklea Dome 3D Mineral Mapping Test Data Set (Laukamp, 2020; https://doi.org/10.25919/5ed83bf55be6a) provides an opportunity for evaluating proximal and remote sensing data, validated and calibrated by independent geochemical and mineralogical analyses, for exploration of channel iron deposits (CIDs) through cover. We present hyperspectral airborne, surface, and drill core reflectance spectra collected in the visible–near-infrared and shortwave infrared wavelength ranges (VNIR–SWIR; 350 to 2500 nm), as well as whole-rock geochemistry obtained by means of X-ray fluorescence analysis and loss-on-ignition measurements of drill core samples. The integration of surface with subsurface hyperspectral data collected in the frame of previously published Rocklea Dome 3D Mineral Mapping case studies demonstrated that about 30 % of exploration drill holes were sunk into barren ground and could have been of better use, located elsewhere, if airborne hyperspectral imagery had been consulted for drill hole planning. The remote mapping of transported Tertiary detritals (i.e. potential hosts of channel iron ore resources) versus weathered in situ Archaean bedrock (i.e. barren ground) has significant implications for other areas where “cover” (i.e. regolith and/or sediments covering bedrock hosting mineral deposits) hinders mineral exploration. Hyperspectral remote sensing represents a cost-effective method for regolith landform mapping required for planning drilling programmes. In the Rocklea Dome area, vegetation unmixing methods applied to airborne hyperspectral data, integrated with subsurface data, resulted in seamless mapping of ore zones from the weathered surface to the base of the CID – a concept that can be applied to other mineral exploration and mineral deposit studies. Furthermore, the associated, independent calibration data allowed the quantification of iron oxide phases and associated mineralogy from hyperspectral data. Using the Rocklea Dome data set, novel geostatistical clustering methods were applied to the drill core data sets for ore body domaining that introduced scientific rigour to a traditionally subjective procedure, resulting in reproducible objective domains that are critical for the mining process. Beyond the previously published case studies, the Rocklea Dome 3D Mineral Mapping Test Data Set has the potential to develop new methods for advanced resource characterisation and develop new applications that aid exploration for mineral deposits through cover. The white mica and chlorite abundance maps derived from airborne hyperspectral, presented here for the first time, highlight the additional applications of remote sensing for geological mapping and could help to evaluate newly launched hyper- and multispectral spaceborne systems for geoscience and mineral exploration.

Description: Land-use models and integrated assessment models provide scenarios of land-use and land-cover (LULC) changes following pathways or storylines related to different socioeconomic and environmental developments. The large diversity of available scenario projections leads to a recognizable variability in impacts on land ecosystems and the levels of services provided. We evaluated 16 projections of future LULC until 2040 that reflected different assumptions regarding socioeconomic demands and modeling protocols. By using these LULC projections in a state-of-the-art dynamic global vegetation model, we simulated their effect on selected ecosystem service indicators related to ecosystem productivity and carbon sequestration potential, agricultural production and the water cycle. We found that although a common trend for agricultural expansion exists across the scenarios, where and how particular LULC changes are realized differs widely across models and scenarios. They are linked to model-specific considerations of some demands over others and their respective translation into LULC changes and also reflect the simplified or missing representation of processes related to land dynamics or other influencing factors (e.g., trade, climate change). As a result, some scenarios show questionable and possibly unrealistic features in their LULC allocations, including highly regionalized LULC changes with rates of conversion that are contrary to or exceed rates observed in the past. Across the diverging LULC projections, we identified positive global trends of net primary productivity (+10.2 % ± 1.4 %), vegetation carbon (+9.2 % ± 4.1 %), crop production (+31.2 % ± 12.2 %) and water runoff (+9.3 % ± 1.7 %), and a negative trend of soil and litter carbon stocks (−0.5 % ± 0.4 %). The variability in ecosystem service indicators across scenarios was especially high for vegetation carbon stocks and crop production. Regionally, variability was highest in tropical forest regions, especially at current forest boundaries, because of intense and strongly diverging LULC change projections in combination with high vegetation productivity dampening or amplifying the effects of climatic change. Our results emphasize that information on future changes in ecosystem functioning and the related ecosystem service indicators should be seen in light of the variability originating from diverging projections of LULC. This is necessary to allow for adequate policy support towards sustainable transformations.

Description: Advanced aeroelastically optimized tip extensions are among rotor innovation concepts which could contribute to the higher performance and lower cost of wind turbines. A novel design optimization framework for wind turbine blade tip extensions based on surrogate aeroelastic modeling is presented. An academic wind turbine is modeled in an aeroelastic code equipped with a near-wake aerodynamic module, and tip extensions with complex shapes are parametrized using 11 design variables. The design space is explored via full aeroelastic simulations in extreme turbulence, and a surrogate model is fitted to the data. Direct optimization is performed based on the surrogate model seeking to maximize the power of the retrofitted turbine within the ultimate load constraints. The presented optimized design achieves a load-neutral gain of up to 6 % in annual energy production. Its performance is further evaluated in detail by means of the near-wake model used for the generation of the surrogate model and compared with a higher-fidelity aerodynamic module comprising a hybrid filament-particle-mesh vortex method with a lifting-line implementation. A good agreement between the solvers is obtained at low turbulence levels, while differences in predicted power and flapwise blade root bending moment grow with increasing turbulence intensity.

Description: Changes in the vegetation of Brazilian Cerrado may occur over time. However, long-term dynamics are not fully understood yet, especially woody plant encroachment (WPE). The objective of this study was to examine changes in vegetation structure in a preserved area in Triângulo Mineiro region, within the southern Brazilian Cerrado domain, over 32 years (1987, 2005, and 2019). We based the study on field and literature surveys, as well as satellite imagery, and hypothesized that, due to the absence of periodic fires or grazing, Cerrado open formations (i.e., grassland or savanna) tend to become denser due to WPE. Shrubby grassland cover assessed in 1987 disappeared in the following periods (from 30.0 % to 0.0 % in 2019) while forest formations increased (from 7.0 % in 1987 to 11.0 % in 2019). Changes between 2005 and 2019 occurred within the stricto sensu cerrado subdivisions, with reduction of sparse cerrado (from 34.2 % to 7.7 %) and an increase in dense cerrado (from 6.9 % to 39.8 %). Normalized difference vegetation index (NDVI) applied for similar periods indicates a progressive increase of values over time (from 1986 (0.61±0.10) to 2004 (0.65±0.06) and 2018 (0.78±0.05)) and corroborates the WPE process. These patterns imply the loss of biodiversity in open plant formation. Another major consequence was the reduction of wetlands and possible impact on water supply. Such patterns are important to support plant management plans for the threatened Cerrado open plant formations.

Description: Nitrogen (N) availability can be highly variable in tropical forests on regional and local scales. While environmental gradients influence N cycling on a regional scale, topography is known to affect N availability on a local scale. We compared natural abundance of 15N isotopes of soil profiles in tropical lowland forest, tropical montane forest, and subtropical Miombo woodland within the Congo Basin as a proxy to assess ecosystem-level differences in N cycling. Soil δ15N profiles indicated that N cycling in the montane forest is relatively more closed and dominated by organic N turnover, whereas the lowland forest and Miombo woodland experienced a more open N cycle dominated by inorganic N. Furthermore, we examined the effect of slope gradient on soil δ15N within forest types to quantify local differences induced by topography. Our results show that slope gradient only affects the soil δ15N in the Miombo forest, which is prone to erosion due to a lower vegetation cover and intense rainfall at the onset of the wet season. Lowland forest, on the other hand, with a flat topography and protective vegetation cover, showed no influence of topography on soil δ15N in our study site. Despite the steep topography, slope angles do not affect soil δ15N in the montane forest, although stable isotope signatures exhibited higher variability within this ecosystem. A pan-tropical analysis of soil δ15N values (i.e., from our study and literature) reveals that soil δ15N in tropical forests is best explained by factors controlling erosion, namely mean annual precipitation, leaf area index, and slope gradient. Erosive forces vary immensely between different tropical forest ecosystems, and our results highlight the need for more spatial coverage of N cycling studies in tropical forests, to further elucidate the local impact of topography on N cycling in this biome.

Description: We present a set of observations on meltwater meandering rivulets on ice and compare them (qualitatively and quantitatively) to morphologies commonly found in meandering channels in different media. The observations include data from planned centimeter-scale experiments and from incidental self-formed millimeter-scale rivulets. Our data show pulsed lateral migration features, undercut banks and overhangs, meander bend skewness, and meander bend cutoffs. The data also compare well with planform characteristics of alluvial meandering rivers (sinuosity, wavelength-to-width ratios, and meander bend fatness and skewness). We discuss the (ir)relevance of scale in our experiments, which, in spite of being in the laminar flow regime and likely affected by surface tension effects, are capable of shedding light into the processes driving formation and evolution of supraglacial meltwater meandering channels. Our observations suggest that sinuosity growth in meltwater meandering channels on ice is a function of flow velocity and the interplay between vertical and lateral incision driven by temperature differences between flow and ice. In the absence of recrystallization (depositional analog to alluvial rivers), bends are more likely to be downstream-skewed and channels show lower sinuosities.

Description: The El Niño‐-Southern Oscillation (ENSO) influences the global climate and the variability in the terrestrial carbon cycle on interannual timescales. Two different expressions of El Niño have recently been identified: (i) central Pacific (CP) and (ii) eastern Pacific (EP). Both types of El Niño are characterised by above-average sea surface temperature anomalies at the respective locations. Studies exploring the impact of these expressions of El Niño on the carbon cycle have identified changes in the amplitude of the concentration of interannual atmospheric carbon dioxide (CO2) variability following increased tropical near-surface air temperature and decreased precipitation. We employ the dynamic global vegetation model LPJ-GUESS (Lund–Potsdam–Jena General Ecosystem Simulator) within a synthetic experimental framework to examine the sensitivity and potential long-term impacts of these two expressions of El Niño on the terrestrial carbon cycle. We manipulated the occurrence of CP and EP events in two climate reanalysis datasets during the latter half of the 20th and early 21st century by replacing all EP with CP and separately all CP with EP El Niño events. We found that the different expressions of El Niño affect interannual variability in the terrestrial carbon cycle. However, the effect on longer timescales was small for both climate reanalysis datasets. We conclude that capturing any future trends in the relative frequency of CP and EP El Niño events may not be critical for robust simulations of the terrestrial carbon cycle.

Description: A primary mercury gas standard was developed at Van Swinden Laboratory (VSL) to establish an International System of Units (SI)-traceable reference point for mercury concentrations at emission and background levels in the atmosphere. The majority of mercury concentration measurements are currently made traceable to the empirically determined vapour pressure of mercury. The primary mercury gas standard can be used for the accurate and precise calibration of analytical systems used for measuring mercury concentrations in air. It has been especially developed to support measurements related to ambient air monitoring (1–2 ng m−3), indoor and workplace-related mercury concentration levels according to health standards (from 50 ng m−3 upwards) as well as stationary source emissions (from 1 µg m−3 upwards). The primary mercury gas standard is based on diffusion according to ISO 6154-8. Calibration gas mixtures are obtained by combining calibrated mass flows of nitrogen and air through a generator holding diffusion cells containing elemental mercury. In this paper, we present the results of comparisons between the primary gas standard and mercury calibration methods maintained by NPL (National Physical Laboratory in the United Kingdom), a National Metrology Institute (NMI), and the Jozef Stefan Institute (JSI), a Designated Institute (DI). The calibration methods currently used at NPL and JSI are based on the bell-jar calibration apparatus in combination with the Dumarey equation or a NIST (National Institute of Standards and Technology in the United States) reference material. For the comparisons, mercury was sampled on sorbent traps to obtain transfer standards with levels between 2 and 1000 ng with an expanded uncertainty not exceeding 3 % (k=2). The comparisons performed show that the results for the primary gas standard and the NIST reference material are comparable, whereas a difference of −8 % exists between results traceable to the primary gas standard and the Dumarey equation.

Description: Analysis of short sediment cores collected in Grand Lake, Labrador, revealed that this lake is an excellent candidate for the preservation of a laminated sediment record. The great depth of Grand Lake, the availability of fine sediments along its tributaries and its important seasonal river inflow have favoured the formation of a 160-year-long clastic varved sequence. Each varve represents 1 hydrological year. Varve formation is mainly related to spring discharge conditions with contributions from summer and autumn rainfall events. The statistically significant relation between varve parameters and the Naskaupi River discharge observations provided the opportunity to develop local hydrological reconstructions beyond the instrumental period. The combined detrital layer thickness and the particle size (99th percentile) series extracted from each varve yield the strongest correlations with instrumental data (r=0.68 and 0.75 respectively) and have been used to reconstruct the respective Naskaupi River mean and maximum annual discharges over the 1856–2016 period. The reconstructed Q-mean series suggest that high Q-mean years occurred during the 1920–1960 period, and a slight decrease in Q-mean took place during the second half of the 20th century. Independent reconstructions based on rainfall–runoff modelling of the watershed from historical reanalysis of global geopotential height fields display a significant correlation with the reconstructed Naskaupi River discharge based on varve physical parameters. The Grand Lake varved sequence contains a regional hydrological signal, as suggested by the statistically significant relation between the combined detrital layer thickness series and the observed Labrador region Q-mean series extracted from five watersheds of different sizes.

Description: We investigate the contributions of local and regional emission sources to air pollution in Beijing to inform the design of short-term emission control strategies for mitigating major pollution episodes. We use a well-evaluated version of the WRF-Chem model at 3 km horizontal resolution to determine the daily accumulation of pollution over Beijing from local and regional sources in October 2014 under a range of meteorological conditions. Considering feasible emission reductions across residential, transport, power, and industrial sectors, we find that 1 d controls on local emissions have an immediate effect on PM2.5 (particulate matter with diameter less than 2.5 µm) concentrations on the same day but can have lingering effects as much as 5 d later under stagnant conditions. One-day controls in surrounding provinces have the greatest effect in Beijing on the day following the controls but may have negligible effects under northwesterly winds when local emissions dominate. To explore the contribution of different emission sectors and regions, we perform simulations with each source removed in turn. We find that residential and industrial sectors from neighbouring provinces dominate PM2.5 levels in Beijing during major pollution episodes but that local residential emissions and industrial or residential emissions from more distant provinces can also contribute significantly during some episodes. We then perform a structured set of perturbed emission simulations to allow us to build statistical emulators that represent the relationships between emission sources and air pollution in Beijing over the period. We use these computationally fast emulators to determine the sensitivity of PM2.5 concentrations to different emission sources and the interactions between them, including for secondary PM, and to create pollutant response surfaces for daily average PM2.5 concentrations in Beijing. We use these surfaces to identify the short-term emission controls needed to meet the national air quality target of daily average PM2.5 less than 75 µg m−3 for pollution episodes of different intensities. We find that for heavily polluted days with daily mean PM2.5 higher than 225 µg m−3, even emission reductions of 90 % across all sectors over Beijing and surrounding provinces may be insufficient to meet the national air quality standards. These results highlight the regional nature of PM pollution and the challenges of tackling it during major pollution episodes.

Description: Climate is critically affected by aerosols, which alter cloud lifecycles and precipitation distribution through radiative and microphysical effects. In this study, aerosol and cloud property datasets from MODIS (Moderate Resolution Imaging Spectroradiometer), onboard the Aqua satellite, and surface observations, including aerosol concentrations, raindrop size distribution, and meteorological parameters, were used to statistically quantify the effects of aerosols on low-level warm-cloud microphysics and drizzle over northern Taiwan during multiple fall seasons (from 15 October to 30 November of 2005–2017). Our results indicated that northwestern Taiwan, which has several densely populated cities, is dominated by low-level clouds (e.g., warm, thin, and broken clouds) during the fall season. The observed effects of aerosols on warm clouds indicated aerosol indirect effects (i.e., increased aerosol loading caused a decrease in cloud effective radius (CER)), an increase in cloud optical thickness, an increase in cloud fraction, and a decrease in cloud-top temperature under a fixed cloud water path. Quantitatively, aerosol–cloud interactions (ACI=-∂ln⁡CER∂ln⁡α|CWP, changes in CER relative to changes in aerosol amounts) were 0.07 for our research domain and varied between 0.09 and 0.06 in the surrounding remote (i.e., ocean) and polluted (i.e., land) areas, respectively, indicating aerosol indirect effects were stronger in the remote area. From the raindrop size distribution analysis, high aerosol loading resulted in a decreased frequency of drizzle events, redistribution of cloud water to more numerous and smaller droplets, and reduced collision–coalescence rates. However, during light rain (≤1 mm h−1), high aerosol concentrations drove raindrops towards smaller droplet sizes and increased the appearance of drizzle drops. This study used long-term surface and satellite data to determine aerosol variations in northern Taiwan, effects on clouds and precipitation, and observational strategies for future research on aerosol–cloud–precipitation interactions.

Description: In this study, the detection and characteristics of dry/wet spells (defined as episodes when precipitation is abnormally low or high compared to usual climatology) drawn from several datasets are compared for Senegal. Here, four datasets are based on satellite data (TRMM-3B42 V7, CMORPH V1.0, TAMSAT V3, and CHIRPS V2. 0), two on reanalysis products (NCEP-CFSR and ERA5), and three on rain gauge observations (CPC Unified V1.0/RT and a 65-rain-gauge network regridded by using two kriging methods, namely ordinary kriging, OK, and block kriging, BK). All datasets were converted to the same spatio-temporal resolution: daily cumulative rainfall on a regular 0.25∘ grid. The BK dataset was used as a reference. Despite strong agreement between the datasets on the spatial variability in cumulative seasonal rainfall (correlations ranging from 0.94 to 0.99), there were significant disparities in dry/wet spells. The occurrence of dry spells is less in products using infrared measurement techniques than in products coupling infrared and microwave, pointing to more frequent dry spell events. All datasets show that dry spells appear to be more frequent at the start and end of rainy seasons. Thus, dry spell occurrences have a major influence on the duration of the rainy season, in particular through the “false onset” or “early cessation” of seasons. The amplitude of wet spells shows the greatest variation between datasets. Indeed, these major wet spells appear more intense in the OK and Tropical Rainfall Measuring Mission (TRMM) datasets than in the others. Lastly, the products indicate a similar wet spell frequency occurring at the height of the West African monsoon. Our findings provide guidance in choosing the most suitable datasets for implementing early warning systems (EWSs) using a multi-risk approach and integrating effective dry/wet spell indicators for monitoring and detecting extreme events.

Description: The Mg/Ca and Sr/Ca ratios of marine shells have been widely used in environmental paleoreconstructions to understand past marine conditions. Temperature calibrations to ostracod Mg/Ca ratios are known to be species-specific but only available for a few species, despite the large number of known ostracod species. Here, we develop temperature calibrations for two shallow marine ostracods, Sinocytheridea impressa and Neomonoceratina delicata, using modern sediment samples. Our results show that adult specimens of these two species might be useful as a paleothermometer. We observed significant correlations using the Mg/Ca ratios of both species to the annual (Mg/CaS. impressa=3.7 ⋅ T−62.7; Mg/CaN. delicata=1.6 ⋅ T−16.8) and April (Mg/CaS. impressa=2.8 ⋅ T−39.2; Mg/CaN. delicata=1.6 ⋅ T−15.7) temperatures. The correlation of temperature to the Mg/Ca ratio of S. impressa is more significant and therefore should be preferred for paleoreconstructions. Re-analysis from satellite data allows us to validate our temperature calibration to an extended area around the Pearl River estuary. Our results show that Mg/Ca of S. impressa and N. delicata ostracods can be used to reconstruct water temperature at a regional scale, which provides information on the oceanic circulation in coastal areas of the South China Sea. Sr/Ca ratios of both species do not correlate with any of the 24 water parameters recorded by the Environmental Protection Department of Hong Kong, including temperature (21.7–24.1 ∘C), salinity (23.8–33.7 PSU), dissolved oxygen (4.3–7.1 mg L−1), suspended solids (1.9–35.4 mg L−1) and pH (7.7–8.2).

Description: The continental divide along the spine of the Canadian Rockies in southwestern Canada is a critical headwater region for hydrological drainages to the Pacific, Arctic, and Atlantic oceans. Major flooding events are typically attributed to heavy precipitation on its eastern side due to upslope (easterly) flows. Precipitation can also occur on the western side of the divide when moisture originating from the Pacific Ocean encounters the west-facing slopes of the Canadian Rockies. Often, storms propagating across the divide result in significant precipitation on both sides. Meteorological data over this critical region are sparse, with few stations located at high elevations. Given the importance of all these types of events, the Storms and Precipitation Across the continental Divide Experiment (SPADE) was initiated to enhance our knowledge of the atmospheric processes leading to storms and precipitation on either side of the continental divide. This was accomplished by installing specialized meteorological instrumentation on both sides of the continental divide and carrying out manual observations during an intensive field campaign from 24 April–26 June 2019. On the eastern side, there were two field sites: (i) at Fortress Mountain Powerline (2076 m a.s.l.) and (ii) at Fortress Junction Service, located in a high-elevation valley (1580 m a.s.l.). On the western side, Nipika Mountain Resort, also located in a valley (1087 m a.s.l.), was chosen as a field site. Various meteorological instruments were deployed including two Doppler light detection and ranging instruments (lidars), three vertically pointing micro rain radars, and three optical disdrometers. The three main sites were nearly identically instrumented, and observers were on site at Fortress Mountain Powerline and Nipika Mountain Resort during precipitation events to take manual observations of precipitation type and microphotographs of solid particles. The objective of the field campaign was to gather high-temporal-frequency meteorological data and to compare the different conditions on either side of the divide to study the precipitation processes that can lead to catastrophic flooding in the region. Details on field sites, instrumentation used, and collection methods are discussed. Data from the study are publicly accessible from the Federated Research Data Repository at https://doi.org/10.20383/101.0221 (Thériault et al., 2020). This dataset will be used to study atmospheric conditions associated with precipitation events documented simultaneously on either side of a continental divide. This paper also provides a sample of the data gathered during a precipitation event.

Description: A leading hypothesis for the mechanism of fast supraglacial lake drainages is that transient extensional stresses briefly allow crevassing in otherwise compressional ice flow regimes. Lake water can then hydrofracture a crevasse to the base of the ice sheet, and river inputs can maintain this connection as a moulin. If future ice sheet models are to accurately represent moulins, we must understand their formation processes, timescales, and locations. Here, we use remote-sensing velocity products to constrain the relationship between strain rates and lake drainages across ∼ 1600 km2 in Pâkitsoq, western Greenland, between 2002–2019. We find significantly more extensional background strain rates at moulins associated with fast-draining lakes than at slow-draining or non-draining lake moulins. We test whether moulins in more extensional background settings drain their lakes earlier, but we find insignificant correlation. To investigate the frequency at which strain-rate transients are associated with fast lake drainage, we examined Landsat-derived strain rates over 16 and 32 d periods at moulins associated with 240 fast-lake-drainage events over 18 years. A low signal-to-noise ratio, the presence of water, and the multi-week repeat cycle obscured any resolution of the hypothesized transient strain rates. Our results support the hypothesis that transient strain rates drive fast lake drainages. However, the current generation of ice sheet velocity products, even when stacked across hundreds of fast lake drainages, cannot resolve these transients. Thus, observational progress in understanding lake drainage initiation will rely on field-based tools such as GPS networks and photogrammetry.

Description: To prevent the spread of the COVID-19 epidemic, restrictions such as “lockdowns” were conducted globally, which led to a significant reduction in fossil fuel emissions, especially in urban areas. However, CO2 concentrations in urban areas are affected by many factors, such as weather, biological sinks and background CO2 fluctuations. Thus, it is difficult to directly observe the CO2 reductions from sparse ground observations. Here, we focus on urban ground transportation emissions, which were dramatically affected by the restrictions, to determine the reduction signals. We conducted six series of on-road CO2 observations in Beijing using mobile platforms before (BC), during (DC) and after (AC) the implementation of COVID-19 restrictions. To reduce the impacts of weather conditions and background fluctuations, we analyze vehicle trips with the most similar weather conditions possible and calculated the enhancement metric, which is the difference between the on-road CO2 concentration and the “urban background” CO2 concentration measured at the tower of the Institute of Atmospheric Physics (IAP), Chinese Academy of Sciences. The results showed that the DC CO2 enhancement was decreased by 41 (±1.3) parts per million (ppm) and 26 (±6.2) ppm compared to those for the BC and AC trips, respectively. Detailed analysis showed that, during COVID-19 restrictions, there was no difference between weekdays and weekends during working hours (09:00–17:00 local standard time; LST). The enhancements during rush hours (07:00–09:00 and 17:00–20:00 LST) were almost twice those during working hours, indicating that emissions during rush hours were much higher. For DC and BC, the enhancement reductions during rush hours were much larger than those during working hours. Our findings showed a clear CO2 concentration decrease during COVID-19 restrictions, which is consistent with the CO2 emissions reductions due to the pandemic. The enhancement method used in this study is an effective method to reduce the impacts of weather and background fluctuations. Low-cost sensors, which are inexpensive and convenient, could play an important role in further on-road and other urban observations.

Description: Dated to approximately 13 000 years ago, the Laacher See (East Eifel volcanic zone) eruption was one of the largest midlatitude Northern Hemisphere volcanic events of the Late Pleistocene. This eruptive event not only impacted local environments and human communities but probably also affected Northern Hemispheric climate. To better understand the impact of a Laacher See-type eruption on NH circulation and climate, we have simulated the evolution of its fine ash and sulfur cloud with an interactive stratospheric aerosol model. Our experiments are based around a central estimate for the Laacher See aerosol cloud of 15 Tg of sulfur dioxide (SO2) and 150 Tg of fine ash, across the main eruptive phases in May and a smaller one in June with 5 Tg SO2 and 50 Tg of fine ash. Additional sensitivity experiments reflect the estimated range of uncertainty of the injection rate and altitude and assess how the solar-absorptive heating from the fine ash emitted in the first eruptive phase changed the volcanic clouds' dispersion. The chosen eruption dates were determined by the stratospheric wind fields to reflect the empirically observed ash lobes as derived from geological, paleoecological and archeological evidence linked directly to the prehistoric Laacher See eruption. Whilst our simulations are based on present-day conditions, and we do not seek to replicate the climate conditions that prevailed 13 000 years ago, we consider our experimental design to be a reasonable approximation of the transport pathways in the midlatitude stratosphere at this time of year. Our simulations suggest that the heating of the ash plays an important role for the transport of ash and sulfate. Depending on the altitude of the injection, the simulated volcanic cloud begins to rotate 1 to 3 d after the eruption. This mesocyclone, as well as the additional radiative heating of the fine ash, then changes the dispersion of the cloud itself to be more southward compared to dispersal estimated without fine ash heating. This ash-cloud-generated southerly migration process may at least partially explain why, as yet, no Laacher See tephra has been found in Greenland ice cores. Sulfate transport is similarly impacted by the heating of the ash, resulting in stronger transport to low latitudes, later arrival of the volcanic cloud in the Arctic regions and a longer lifetime compared to cases without injection of fine ash. Our study offers new insights into the dispersion of volcanic clouds in midlatitudes and addresses a likely behavior of the ash cloud of the Laacher See eruption that darkened European skies at the end of the Pleistocene. In turn, this study can also serve as significant input for scenarios that consider the risks associated with re-awakened volcanism in the Eifel.

Description: In this paper, we compare column-averaged dry-air mole fractions of water vapor (XH2O) retrievals from the COllaborative Carbon Column Observing Network (COCCON) with retrievals from two co-located high-resolution Fourier transform infrared (FTIR) spectrometers as references at two boreal sites, Kiruna, Sweden, and Sodankylä, Finland, from 6 March 2017 to 20 September 2019. In the framework of the Network for the Detection of Atmospheric Composition Change (NDACC), an FTIR spectrometer is operated at Kiruna. The H2O product derived from these observations has been generated with the MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) processor. In Sodankylä, a Total Carbon Column Observing Network (TCCON) spectrometer is operated, and the official XH2O data as provided by TCCON are used for this study. The datasets are in good overall agreement, with COCCON data showing a wet bias of (49.20±58.61) ppm ((3.33±3.37) %, R2=0.9992) compared with MUSICA NDACC and (56.32±45.63) ppm ((3.44±1.77) %, R2=0.9997) compared with TCCON. Furthermore, the a priori H2O volume mixing ratio (VMR) profiles (MAP) used as a priori information in the TCCON retrievals (also adopted for COCCON retrievals) are evaluated with respect to radiosonde (Vaisala RS41) profiles at Sodankylä. The MAP and radiosonde profiles show similar shapes and a good linear correlation of integrated XH2O, indicating that MAP is a reasonable approximation of the true atmospheric state and an appropriate choice for the scaling retrieval methods as applied by COCCON and TCCON. COCCON shows a reduced dry bias (−14.96 %) in comparison with TCCON (−19.08 %) with respect to radiosonde XH2O. Finally, we investigate the quality of satellite data at high latitudes. For this purpose, the COCCON XH2O is compared with retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) generated with the MUSICA processor (MUSICA IASI) and with retrievals from the TROPOspheric Monitoring Instrument (TROPOMI). Both paired datasets generally show good agreement and similar correlations at the two sites. COCCON measures 4.64 % less XH2O at Kiruna and 3.36 % less at Sodankylä with respect to MUSICA IASI, whereas COCCON measures 9.71 % more XH2O at Kiruna and 7.75 % more at Sodankylä compared with TROPOMI. Our study supports the assumption that COCCON also delivers a well-characterized XH2O data product. This emphasizes that this approach might complement the TCCON network with respect to satellite validation efforts. This is the first published study where COCCON XH2O has been compared with MUSICA NDACC and TCCON retrievals and has been used for MUSICA IASI and TROPOMI validation.

Description: We derive recent surface mass balance (SMB) estimates from airborne radar observations along the iSTAR traverse (2013, 2014) at Pine Island Glacier (PIG), West Antarctica. Ground-based neutron probe measurements provide information of snow and firn density with depth at 22 locations and were used to date internal annual reflection layers. The 2005 layer was traced for a total distance of 2367 km to determine annual mean SMB for the period 2005–2014. Using complementary SMB estimates from two regional climate models, RACMO2.3p2 and MAR, and a geostatistical kriging scheme, we determine a regional-scale SMB distribution with similar main characteristics to that determined for the period 1985–2009 in previous studies. Local departures exist for the northern PIG slopes, where the orographic precipitation shadow effect appears to be more pronounced in our observations, and the southward interior, where the SMB gradient is more pronounced in previous studies. We derive total mass inputs of 79.9±19.2 and 82.1±19.2 Gt yr−1 to the PIG basin based on complementary ASIRAS–RACMO and ASIRAS–MAR SMB estimates, respectively. These are not significantly different to the value of 78.3±6.8 Gt yr−1 for the period 1985–2009. Thus, there is no evidence of a secular trend at decadal scales in total mass input to the PIG basin. We note, however, that our estimated uncertainty is more than twice the uncertainty for the 1985–2009 estimate on total mass input. Our error analysis indicates that uncertainty estimates on total mass input are highly sensitive to the selected krige methodology and assumptions made on the interpolation error, which we identify as the main cause for the increased uncertainty range compared to the 1985–2009 estimates.

Description: Uncertain or inaccurate parameters in sea ice models influence seasonal predictions and climate change projections in terms of both mean and trend. We explore the feasibility and benefits of applying an ensemble Kalman filter (EnKF) to estimate parameters in the Los Alamos sea ice model (CICE). Parameter estimation (PE) is applied to the highly influential dry snow grain radius and combined with state estimation in a series of perfect model observing system simulation experiments (OSSEs). Allowing the parameter to vary in space improves performance along the sea ice edge but degrades in the central Arctic compared to requiring the parameter to be uniform everywhere, suggesting that spatially varying parameters will likely improve PE performance at local scales and should be considered with caution. We compare experiments with both PE and state estimation to experiments with only the latter and have found that the benefits of PE mostly occur after the data assimilation period, when no observations are available to assimilate (i.e., the forecast period), which suggests PE's relevance for improving seasonal predictions of Arctic sea ice.

Description: Evolution of organic carbon content in soils has the potential to be a major driver of atmospheric greenhouse gas concentrations over the next century. Understanding soil carbon dynamics is a challenge due to a wide range of residence times of soil organic matter and limited constraints on the mechanisms influencing its persistence. In particular, large uncertainties exist regarding the persistence of pyrogenic organic carbon in soils. In order to characterize organic matter with varying degrees of persistence and to distinguish pyrogenic organic carbon, we combined Rock-Eval analysis, a thermo-chemical method, with the benzene polycarboxylic acid molecular marker method and Raman spectroscopy to characterize samples from long-term bare-fallow experiments, progressively depleted in the most labile organic carbon over time. Considering the heterogeneity of soil samples, size fractions have been separated to distinguish pools of organic carbon with distinct properties. We observe that organic carbon dynamics is dependent on granulometry. A pool of organic carbon with intermediate residence times, from years to a few decades, representing ca. 65 % of the bulk soil organic carbon stock, is mainly associated with fine fractions ( 20 µm) are rich in centennially persistent organic carbon, representing ca. 20 % of the initial organic carbon stock, due to the chemical recalcitrance of organic matter in these fractions, dominated by pyrogenic organic carbon. A second pool of persistent organic carbon, representing ca. 15 % of the initial organic carbon stock, is associated with the clay fraction, indicating mechanisms of protection occurring at the submicron scale (

Description: Geoengineering has been discussed as a potential option to offset the global impacts of anthropogenic climate change and at the same time reach the global temperature targets of the Paris Agreement. Before any implementation of geoengineering, however, the complex natural responses and consequences of such methods should be fully understood to avoid any unexpected and potentially degrading impacts. Here we assess the changes in ecosystem carbon exchange and storage among different terrestrial biomes under three aerosol-based radiation management methods with the baseline of RCP8.5 using an Earth system model (NorESM1-ME). All three methods used in this study (stratospheric aerosol injection, marine sky brightening, cirrus cloud thinning) target the global mean radiation balance at the top of the atmosphere to reach that of the RCP4.5 scenario. The three radiation management (RM) methods investigated in this study show vastly different precipitation patterns, especially in the tropical forest biome. Precipitation differences from the three RM methods result in large variability in global vegetation carbon uptake and storage. Our findings show that there are unforeseen regional consequences under geoengineering, and these consequences should be taken into account in future climate policies as they have a substantial impact on terrestrial ecosystems. Although changes in temperature and precipitation play a large role in vegetation carbon uptake and storage, our results show that CO2 fertilization also plays a considerable role. We find that the effects of geoengineering on vegetation carbon storage are much smaller than the effects of mitigation under the RCP4.5 scenario (e.g., afforestation in the tropics). Our results emphasize the importance of considering multiple combined effects and responses of land biomes while achieving the global temperature targets of the Paris Agreement.

Description: We conduct a global inverse analysis of 2010–2018 GOSAT observations to better understand the factors controlling atmospheric methane and its accelerating increase over the 2010–2018 period. The inversion optimizes anthropogenic methane emissions and their 2010–2018 trends on a 4∘×5∘ grid, monthly regional wetland emissions, and annual hemispheric concentrations of tropospheric OH (the main sink of methane). We use an analytical solution to the Bayesian optimization problem that provides closed-form estimates of error covariances and information content for the solution. We verify our inversion results with independent methane observations from the TCCON and NOAA networks. Our inversion successfully reproduces the interannual variability of the methane growth rate inferred from NOAA background sites. We find that prior estimates of fuel-related emissions reported by individual countries to the United Nations are too high for China (coal) and Russia (oil and gas) and too low for Venezuela (oil and gas) and the US (oil and gas). We show large 2010–2018 increases in anthropogenic methane emissions over South Asia, tropical Africa, and Brazil, coincident with rapidly growing livestock populations in these regions. We do not find a significant trend in anthropogenic emissions over regions with high rates of production or use of fossil methane, including the US, Russia, and Europe. Our results indicate that the peak methane growth rates in 2014–2015 are driven by low OH concentrations (2014) and high fire emissions (2015), while strong emissions from tropical (Amazon and tropical Africa) and boreal (Eurasia) wetlands combined with increasing anthropogenic emissions drive high growth rates in 2016–2018. Our best estimate is that OH did not contribute significantly to the 2010–2018 methane trend other than the 2014 spike, though error correlation with global anthropogenic emissions limits confidence in this result.

Description: Social work in prisons not only works with and for people who are confined; it also constitutes a form of social work that is carried out under conditions of confinement. This article draws on carceral geography to understand the corporeal and spatial aspects of social work in prison settings. Based on insights from two prisons in Switzerland, we argue that understanding carceral social work as a spatial and materially situated practice helps to gain deep insight into the intricate layers of meaning and powerful modes of functioning of prisons and of the people involved. In particular, it shows how the way social work is carried out in prison is supported and strongly structured by the spatiality of the prison itself, allowing for counselling, desk-type social work, rather than for social work that actively initiates and creates spaces for encounters or activities.

Description: Relative to their surface area, estuaries make a disproportionately large contribution of dissolved organic carbon (DOC) to the global carbon cycle, but it is unknown how this will change under a future climate. As such, the response of DOC fluxes from microbially dominated unvegetated sediments to individual and combined future climate stressors of temperature change (from Δ−3 to Δ+5 ∘C compared to ambient mean temperatures) and ocean acidification (OA, ∼ 2× current CO2 partial pressure, pCO2) was investigated ex situ. Warming alone increased sediment heterotrophy, resulting in a proportional increase in sediment DOC uptake; sediments became net sinks of DOC (3.5 to 8.8 mmol C m−2 d−1) at warmer temperatures (Δ+3 and Δ+5 ∘C, respectively). This temperature response changed under OA conditions, with sediments becoming more autotrophic and a greater sink of DOC (up to 4× greater than under current pCO2 conditions). This response was attributed to the stimulation of heterotrophic bacteria with the autochthonous production of labile organic matter by microphytobenthos. Extrapolating these results to the global area of unvegetated subtidal estuarine sediments, we find that the future climate of warming (Δ+3 ∘C) and OA may decrease estuarine export of DOC by ∼ 80 % (∼ 150 Tg C yr−1) and have a disproportionately large impact on the global DOC budget.

Description: A new integrated mass-flux adjustment filter is introduced, which uses the analyzed integrated mass-flux divergence field to correct the analyzed wind field. The filter has been examined in twin experiments with rapid update cycling using an idealized setup for convective-scale radar data assimilation. It is found that the new filter slightly reduces the accuracy of background and analysis states; however, it preserves the main structure of cold pools and primary mesocyclone properties of supercells. More importantly, it considerably diminishes spurious mass-flux divergence and the high surface pressure tendency, and it thus results in more dynamically balanced analysis states. For the ensuing 3 h forecasts, the experiment that employs the filter becomes more skillful after 1 h. These preliminary results show that the filter is a promising tool to alleviate the imbalance problem caused by data assimilation, especially for convective-scale applications.

Description: The origin of intraplate volcanism is not explained by plate tectonic theory, and several models have been put forward for explanation. One of these models involves edge-driven convection (EDC), in which cold and thick continental lithosphere is juxtaposed with warm and thin oceanic lithosphere to trigger convective instability. To test whether EDC can produce long-lived high-volume magmatism, we run numerical models of EDC for a wide range of mantle properties and edge (i.e., the oceanic–continental transition) geometries. We find that the most important parameters that govern EDC are the rheological parameters mantle viscosity η0 and activation energy Ea. However, even the maximum melting volumes predicted by our most extreme cases are insufficient to account for island-building volcanism on old seafloor, such as at the Canary Islands and Cabo Verde. Also, beneath old seafloor, localized EDC-related melting commonly transitions into widespread melting due to small-scale sublithospheric convection, inconsistent with the distribution of volcanism at these volcano chains. In turn, EDC is a good candidate to sustain the formation of small seamounts on young seafloor, as it is a highly transient phenomenon that occurs in all our models soon after initiation. In a companion paper, we investigate the implications of interaction of EDC with mantle plume activity (Manjón-Cabeza Córdoba and Ballmer, 2021).

Description: Infrared spectroscopic observations have shown that crystalline ammonium nitrate (AN) particles are an abundant constituent of the upper tropospheric aerosol layer which is formed during the Asian summer monsoon period, the so-called Asian Tropopause Aerosol Layer (ATAL). At upper tropospheric temperatures, the thermodynamically stable phase of AN is different from that at 298 K, meaning that presently available room-temperature optical constants of AN, that is, the real and imaginary parts of the complex refractive index, cannot be applied for the quantitative analysis of these infrared measurements. In this work, we have retrieved the first low-temperature data set of optical constants for crystalline AN in the 800–6000 cm−1 wavenumber range with a spectral resolution of 0.5 cm−1. The optical constants were iteratively derived from an infrared extinction spectrum of 1 µm sized AN particles suspended in a cloud chamber at 223 K. The uncertainties of the new data set were carefully assessed in a comprehensive sensitivity analysis. We show that our data accurately fit aircraft-borne infrared measurements of ammonium nitrate particles in the ATAL.

Description: Air pollution, in particular high concentrations of particulate matter smaller than 1 µm in diameter (PM1), continues to be a major health problem, and meteorology is known to substantially influence atmospheric PM concentrations. However, the scientific understanding of the ways in which complex interactions of meteorological factors lead to high-pollution episodes is inconclusive. In this study, a novel, data-driven approach based on empirical relationships is used to characterize and better understand the meteorology-driven component of PM1 variability. A tree-based machine learning model is set up to reproduce concentrations of speciated PM1 at a suburban site southwest of Paris, France, using meteorological variables as input features. The model is able to capture the majority of occurring variance of mean afternoon total PM1 concentrations (coefficient of determination (R2) of 0.58), with model performance depending on the individual PM1 species predicted. Based on the models, an isolation and quantification of individual, season-specific meteorological influences for process understanding at the measurement site is achieved using SHapley Additive exPlanation (SHAP) regression values. Model results suggest that winter pollution episodes are often driven by a combination of shallow mixed layer heights (MLHs), low temperatures, low wind speeds, or inflow from northeastern wind directions. Contributions of MLHs to the winter pollution episodes are quantified to be on average ∼5 µg/m3 for MLHs below

Description: The Organization of Tropical East Pacific Convection (OTREC) field campaign investigated the dynamical structure of convection in the tropical east Pacific and Caribbean. One of the central data sets for this field campaign is the thermodynamic structure of the atmosphere measured by dropsondes released from the NSF/NCAR G-V research aircraft. Between 7 August and 2 October 2019, 648 dropsondes were successfully released from 22 research flights. Soundings were launched in a grid pattern with a typical spacing of 1∘ longitude and 1.2∘ latitude and provided profiles of pressure, temperature, humidity, and winds between the surface and on average 13.3 km. Of these soundings, 636 provided complete vertical profiles of all parameters with a nominal vertical resolution between 6 to 12 m from the surface to almost flight altitude. OTREC deployed the new NRD41 dropsonde, which is the most advanced model that has been developed at NCAR. Here, we describe the data set, the processing of the measurements, and general statistics of all dropsonde observations. The data set is available at https://doi.org/10.26023/EHRT-TN96-9W04 (UCAR/NCAR and Vömel, 2019).

Description: The Tempest, the last work entirely attributed to William Shakespeare, has been subject to many studies and interpretations, ranging from adventure and Shakespeare's biography to colonialism and the cultural revolution, and is studied in this paper in the context of naturally occurring hazards. The play tells the story of a magician, Prospero, and his daughter who are shipwrecked on an unknown island where they encounter strange creatures and beings. But is it a fantastic island or was the author inspired by real places? Literary scholars proposed several hypotheses through the years, based on historical sources. Here, we analyse the play in the light of geosciences and mythology supporting the hypothesis that the playwright was inspired by the Mediterranean. Our goal is not to identify the island but rather to examine the various geographical and philosophical–political factors that may have influenced Shakespeare's literary creation. Nevertheless, some verses in the play suggest volcanism, placing the island in the Sicilian sea. This underlines once again how deep the playwright's knowledge of Italy was. It also suggests that this part of the Mediterranean was known, at the time of Shakespeare, as the theatre of phenomena originated in the volcanism of the area. One implication is that he could have used historical sources, still unknown and precious, to reconstruct geological events that occurred off the Sicilian coast.

Description: Air pollution is a pressing issue that is associated with adverse effects on human health, ecosystems, and climate. Despite many years of effort to improve air quality, nitrogen dioxide (NO2) limit values are still regularly exceeded in Europe, particularly in cities and along streets. This study explores how concentrations of nitrogen oxides (NOx = NO + NO2) in European urban areas have changed over the last decades and how this relates to changes in emissions. To do so, the incremental approach was used, comparing urban increments (i.e. urban background minus rural concentrations) to total emissions, and roadside increments (i.e. urban roadside concentrations minus urban background concentrations) to traffic emissions. In total, nine European cities were assessed. The study revealed that potentially confounding factors like the impact of urban pollution at rural monitoring sites through atmospheric transport are generally negligible for NOx. The approach proves therefore particularly useful for this pollutant. The estimated urban increments all showed downward trends, and for the majority of the cities the trends aligned well with the total emissions. However, it was found that factors like a very densely populated surrounding or local emission sources in the rural area such as shipping traffic on inland waterways restrict the application of the approach for some cities. The roadside increments showed an overall very diverse picture in their absolute values and trends and also in their relation to traffic emissions. This variability and the discrepancies between roadside increments and emissions could be attributed to a combination of local influencing factors at the street level and different aspects introducing inaccuracies to the trends of the emission inventories used, including deficient emission factors. Applying the incremental approach was evaluated as useful for long-term pan-European studies, but at the same time it was found to be restricted to certain regions and cities due to data availability issues. The results also highlight that using emission inventories for the prediction of future health impacts and compliance with limit values needs to consider the distinct variability in the concentrations not only across but also within cities.

Description: Chronological uncertainty is a hallmark of the paleoenvironmental sciences and geosciences. While many tools have been made available to researchers to quantify age uncertainties suitable for various settings and assumptions, disparate tools and output formats often discourage integrative approaches. In addition, associated tasks like propagating age-model uncertainties to subsequent analyses, and visualizing the results, have received comparatively little attention in the literature and available software. Here, we describe geoChronR, an open-source R package to facilitate these tasks. geoChronR is built around an emerging data standard (Linked PaleoData, or LiPD) and offers access to four popular age-modeling techniques (Bacon, BChron, OxCal, BAM). The output of these models is used to conduct ensemble data analysis, quantifying the impact of chronological uncertainties on common analyses like correlation, regression, principal component, and spectral analyses by repeating the analysis across a large collection of plausible age models. We present five real-world use cases to illustrate how geoChronR may be used to facilitate these tasks, visualize the results in intuitive ways, and store the results for further analysis, promoting transparency and reusability.

Description: Gender affects all aspects of life, and the working and learning environments of science, technology, engineering and geosciences present no exception. Gender issues concerning the access, permanence and ascension of women in exact sciences and Earth sciences careers in general are related to a variety of causes. The underrepresentation of women in science communications, sexual or moral harassment caused by professors and colleagues during undergraduate and graduate ages or the overloading of girls, when compared to boys, with housework during early school ages are some examples mentioned in the literature. In other words, the gender imbalance in science and technology careers may be seen as the result of a series of structured oppression suffered by women of all ages. In this context, we propose the development of an education package that is designed to understand these processes at different levels. One of the tools of this package is known as the “Theatre of the Oppressed”. Elaborated on by Augusto Boal in the 1970s, the Theatre of the Oppressed uses theatre techniques as a means of promoting social and political changes. Usually, a scene takes place that reveals a situation of oppression. The audience become what is called “spect-actors”, where they become active by exploring, showing and transforming the reality in which they are living. In the context of gender issues in exact sciences careers, the students can stage situations that reveal the subtle actions of power relations that usually put women in subservient positions. Our experience showed that, even though the acting is based on fiction, the spectators learn a great deal from the enactment because the simulation of real-life situations, problems and solutions stimulates the practice of resisting oppression in reality from within a setting that offers a safe space to practise making a change.

Description: Unconventional sources of data that enhance our understanding of internal interactions between socio-economic and hydrological processes are central to modeling human–water systems. Participatory modeling (PM) departs from conventional modeling tools by informing and conceptualizing human–water systems through stakeholder engagement. However, the implementation of many PM processes remains biased, particularly in regions where marginalized communities are present. Many PM processes are not cognizant of differentiation and diversity within a society and tend to treat communities as homogeneous units with similar capabilities, needs, and interests. This undifferentiation leads to the exclusion of key actors, many of whom are associated with marginalized communities. In this study, a participatory model-building framework (PMBF), aiming to ensure the inclusiveness of marginalized stakeholders – who (1) have low literacy, (2) are comparatively powerless, and/or (3) are associated with a marginalized language – in participatory modeling, is proposed. The adopted approach employs interdisciplinary storylines to inform and conceptualize human–water systems. The suggested method is underpinned by the multi-level perspective (MLP) framework, which was developed by Geels et al. (2002) to conceptualize socio-technical transitions and modified in this study to accommodate the development of interdisciplinary storylines. A case study was conducted in Atitlán Basin, Guatemala, to understand the relationships that govern the lake's cultural eutrophication problem. This research integrated key stakeholders from the Indigenous Mayan community, associated with diverse literacy ranges, and emerging from three different marginalized linguistic backgrounds (Kaqchikel, Tz'utujil, and K'iche'), in the PM activity. The proposed approach facilitated the participation of marginalized stakeholders. Moreover, it (1) helped develop an understanding of mechanisms governing the eutrophication of the lake, (2) initiated a dialogue between Indigenous Peoples and non-Indigenous stakeholders, and (3) extracted potential solutions targeting the system's leverage points. The participatory model-building activity generated three submodules: (1) agriculture, (2) tourism, and (3) environmental awareness. Each submodule contained socioculturally specific mechanisms associated with nutrient discharge to Lake Atitlán. The delineation of such nuanced relationships helps develop well-targeted policies and best management practices (BMPs). Additionally, the suggested process helped decrease the impact of power imbalances in water resources management and empowered community-based decision-making.

Description: The spatiotemporal distributions of aerosol optical properties and major aerosol types, along with the vertical distribution of major aerosol types over Australia, are investigated based on multi-year Aerosol Robotic Network (AERONET) observations at nine sites, the Moderate Resolution Imaging Spectroradiometer (MODIS), Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), and back-trajectory analysis from the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT). During the observation period from 2001–2020, the annual aerosol optical depth (AOD) at most sites showed increasing trends (0.002–0.029 yr−1), except for that at three sites, Canberra, Jabiru, and Lake Argyle, which showed decreasing trends (−0.004 to −0.014 yr−1). In contrast, the annual Ångström exponent (AE) showed decreasing tendencies at most sites (−0.045 to −0.005 yr−1). The results showed strong seasonal variations in AOD, with high values in the austral spring and summer and relatively low values in the austral fall and winter, and weak seasonal variations in AE, with the highest mean values in the austral spring at most sites. Monthly average AOD increases from August to December or the following January and decreases during March–July. Spatially, the MODIS AOD showed obvious spatial heterogeneity, with high values appearing over the Australian tropical savanna regions, Lake Eyre Basin, and southeastern regions of Australia, while low values appeared over the arid regions in western Australia. MERRA-2 showed that carbonaceous aerosol over northern Australia, dust over central Australia, sulfate over densely populated northwestern and southeastern Australia, and sea salt over Australian coastal regions are the major types of atmospheric aerosols. The nine ground-based AERONET sites over Australia showed that the mixed type of aerosols (biomass burning and dust) is dominant in all seasons. Moreover, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) showed that polluted dust is the dominant aerosol type detected at heights 0.5–5 km over the Australian continent during all seasons. The results suggested that Australian aerosol has similar source characteristics due to the regional transport over Australia, especially for biomass burning and dust aerosols. However, the dust-prone characteristic of aerosol is more prominent over central Australia, while the biomass-burning-prone characteristic of aerosol is more prominent in northern Australia.

Description: Rock salt has remarkable mechanical properties and high economic importance; however, the strength of salt compared to other rocks makes it a rather vulnerable material. Human activities could lead to acceleration of the dissolution of soluble rock salt and collapse of subsurface caverns. Although sinkhole development can be considered a local geological disaster regarding the characteristic size of surface depressions, the deformations can result in catastrophic events. In this study we report the spatiotemporal evolution of surface deformation in the Solotvyno salt mine area in Ukraine based on Sentinel-1 interferometric synthetic aperture radar measurements. Although the mining operations were finished in 2010, several sinkholes have been opened up since then. Our results show that despite the enormous risk management efforts, the sinkholes continue to expand with a maximum line-of-sight deformation rate of 5 cm/yr. The deformation time series show a rather linear feature, and unfortunately no slowdown of the processes can be recognized based on the investigated 4.5-year-long data set. We utilized both ascending and descending satellite passes to discriminate the horizontal and vertical deformations, and our results revealed that vertical deformation is much more pronounced in the area. Analytical source modeling confirmed that the complex deformation pattern observed by Sentinel-1 radar interferometry has a direct connection to the former mining activity and is confined to the mining territory. With the 6 d repetition time of Sentinel-1 observations, the evolution of surface changes can be detected in quasi real time, which can facilitate disaster response and recovery.

Description: Tropical ecosystems contribute significantly to global emissions of methane (CH4), and landscape topography influences the rate of CH4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O2) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O2 concentrations increased in the valley. We simulated the dynamics of CH4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O2, acetate, and associated changes in soil pH with topography regulated simulated CH4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH4 production (acetate, CO2, and H2) and oxidation (O2 and CH4) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH4. A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH4 emissions. This study enhanced the predictive capability for CH4 emissions associated with complex topography and drought in wet tropical forest soils.

Description: Can experimental studies on the behavioural impacts of ocean acidification be trusted? That question was raised in early 2020 when a high-profile paper failed to corroborate previously observed responses of coral reef fish to high CO2. New information on the methodologies used in the “replicated” studies now provides a plausible explanation: the experimental conditions were substantially different. High sensitivity to test conditions is characteristic of ocean acidification research; such response variability shows that effects are complex, interacting with many other factors. Open-minded assessment of all research results, both negative and positive, remains the best way to develop process-based understanding. As in other fields, replication studies in ocean acidification are most likely to contribute to scientific advancement when carried out in a spirit of collaboration rather than confrontation.

Description: The present work analyzes the quality and reliability of an important class of general-purpose, second-order accurate finite-volume (FV) solvers for the large-eddy simulation of a neutrally stratified atmospheric boundary layer (ABL) flow. The analysis is carried out within the OpenFOAM® framework, which is based on a colocated grid arrangement. A series of open-channel flow simulations are carried out using a static Smagorinsky model for subgrid scale momentum fluxes in combination with an algebraic equilibrium wall-layer model. The sensitivity of the solution to variations in numerical parameters such as grid resolution (up to 1603 control volumes), numerical solvers, and interpolation schemes for the discretization of nonlinear terms is evaluated and results are contrasted against those from a well-established mixed pseudospectral–finite-difference code. Considered flow statistics include mean streamwise velocity, resolved Reynolds stresses, velocity skewness and kurtosis, velocity spectra, and two-point autocorrelations. A quadrant analysis along with the examination of the conditionally averaged flow field are performed to investigate the mechanisms responsible for momentum transfer in the flow. It is found that at the selected grid resolutions, the considered class of FV-based solvers yields a poorly correlated flow field and is not able to accurately capture the dominant mechanisms responsible for momentum transport in the ABL. Specifically, the predicted flow field lacks the well-known sweep and ejection pairs organized side by side along the cross-stream direction, which are representative of a streamwise roll mode. This is especially true when using linear interpolation schemes for the discretization of nonlinear terms. This shortcoming leads to a misprediction of flow statistics that are relevant for ABL flow applications and to an enhanced sensitivity of the solution to variations in grid resolution, thus calling for future research aimed at reducing the impact of modeling and discretization errors.

Description: A large number of water masses are presented in the Atlantic Ocean, and knowledge of their distributions and properties is important for understanding and monitoring of a range of oceanographic phenomena. The characteristics and distributions of water masses in biogeochemical space are useful for, in particular, chemical and biological oceanography to understand the origin and mixing history of water samples. Here, we define the characteristics of the major water masses in the Atlantic Ocean as source water types (SWTs) from their formation areas, and map out their distributions. The SWTs are described by six properties taken from the biased-adjusted Global Ocean Data Analysis Project version 2 (GLODAPv2) data product, including both conservative (conservative temperature and absolute salinity) and non-conservative (oxygen, silicate, phosphate and nitrate) properties. The distributions of these water masses are investigated with the use of the optimum multi-parameter (OMP) method and mapped out. The Atlantic Ocean is divided into four vertical layers by distinct neutral densities and four zonal layers to guide the identification and characterization. The water masses in the upper layer originate from wintertime subduction and are defined as central waters. Below the upper layer, the intermediate layer consists of three main water masses: Antarctic Intermediate Water (AAIW), Subarctic Intermediate Water (SAIW) and Mediterranean Water (MW). The North Atlantic Deep Water (NADW, divided into its upper and lower components) is the dominating water mass in the deep and overflow layer. The origin of both the upper and lower NADW is the Labrador Sea Water (LSW), the Iceland–Scotland Overflow Water (ISOW) and the Denmark Strait Overflow Water (DSOW). The Antarctic Bottom Water (AABW) is the only natural water mass in the bottom layer, and this water mass is redefined as Northeast Atlantic Bottom Water (NEABW) in the north of the Equator due to the change of key properties, especially silicate. Similar with NADW, two additional water masses, Circumpolar Deep Water (CDW) and Weddell Sea Bottom Water (WSBW), are defined in the Weddell Sea region in order to understand the origin of AABW.