ALBERT

Description: Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10 000 BP to present) and into the last glacial period, extending back to 60 000 BP. Tephra records of past volcanism, compared to ice cores, are less complete but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 140 000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record, we include stochastically generated synthetic eruptions assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice-core-based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate and the mechanisms leading to abrupt transitions in the palaeoclimate record with low- to high-complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey and Schindlbeck-Belo, 2023).

Description: Much of our understanding of Cenozoic climate is based on the record of δ18O measured in benthic foraminifera. However, this measurement reflects a combined signal of global temperature and sea level, thus preventing a clear understanding of the interactions and feedbacks of the climate system in causing global temperature change. Our new reconstruction of temperature change over the past 4.5 million years includes two phases of long-term cooling, with the second phase of accelerated cooling during the Middle Pleistocene Transition (1.5 to 0.9 million years ago) being accompanied by a transition from dominant 41,000-year low-amplitude periodicity to dominant 100,000-year high-amplitude periodicity. Changes in the rates of long-term cooling and variability are consistent with changes in the carbon cycle driven initially by geologic processes, followed by additional changes in the Southern Ocean carbon cycle. 〈/jats:p〉

Description: Antarctic krill, crucial to the Southern Ocean ecosystem and a vital fisheries resource, is endangered by climate change. Identifying drivers of krill biomass is therefore essential for determining catch limits and designating protection zones. We present a modeling approach to pinpointing effects of sea surface temperature, ice cover, chlorophyll levels, climate indices, and intraspecific competition. Our study reveals that larval recruitment is driven by both competition among age classes and chlorophyll levels. In addition, while milder ice and temperature in spring and summer favor reproduction and early larval survival, both larvae and juveniles strongly benefit from heavier ice and colder temperatures in winter. We conclude that omitting top-down control of resources by krill is only acceptable for retrospective or single-year prognostic models that use field chlorophyll data but that incorporating intraspecific competition is essential for longer-term forecasts. Our findings can guide future krill modeling strategies, reinforcing the sustainability of this keystone species.

Description: 〈jats:p〉Arctic Ocean gateway fluxes play a crucial role in linking the Arctic with the global ocean and affecting climate and marine ecosystems. We reviewed past studies on Arctic–Subarctic ocean linkages and examined their changes and driving mechanisms. Our review highlights that radical changes occurred in the inflows and outflows of the Arctic Ocean during the 2010s. Specifically, the Pacific inflow temperature in the Bering Strait and Atlantic inflow temperature in the Fram Strait hit record highs, while the Pacific inflow salinity in the Bering Strait and Arctic outflow salinity in the Davis and Fram straits hit record lows. Both the ocean heat convergence from lower latitudes to the Arctic and the hydrological cycle connecting the Arctic with Subarctic seas were stronger in 2000–2020 than in 1980–2000. CMIP6 models project a continuing increase in poleward ocean heat convergence in the 21st century, mainly due to warming of inflow waters. They also predict an increase in freshwater input to the Arctic Ocean, with the largest increase in freshwater export expected to occur in the Fram Strait due to both increased ocean volume export and decreased salinity. Fram Strait sea ice volume export hit a record low in the 2010s and is projected to continue to decrease along with Arctic sea ice decline. We quantitatively attribute the variability of the volume, heat, and freshwater transports in the Arctic gateways to forcing within and outside the Arctic based on dedicated numerical simulations and emphasize the importance of both origins in driving the variability.〈/jats:p〉

Description: Abstract. Clouds are assumed to play an important role in the Arctic amplification process. This motivated a detailed investigation of cloud processes, including radiative and turbulent fluxes. Data from the aircraft campaign ACLOUD were analyzed with a focus on the mean and turbulent structure of the cloudy boundary layer over the Fram Strait marginal sea ice zone in late spring and early summer 2017. Vertical profiles of turbulence moments are presented from contrasting atmospheric boundary layers (ABLs) from 4 d. They differ by the magnitude of wind speed, boundary-layer height, stability, the strength of the cloud-top radiative cooling and the number of cloud layers. Turbulence statistics up to third-order moments are presented, which were obtained from horizontal-level flights and from slanted profiles. It is shown that both of these flight patterns complement each other and form a data set that resolves the vertical structure of the ABL turbulence well. The comparison of the 4 d shows that especially during weak wind, even in shallow Arctic ABLs with mixing ratios below 3 g kg-1, cloud-top cooling can serve as a main source of turbulent kinetic energy (TKE).Well-mixed ABLs are generated where TKE is increased and vertical velocity variance shows pronounced maxima in the cloud layer. Negative vertical velocity skewness points then to upside-down convection. Turbulent heat fluxes are directed upward in the cloud layer as a result of cold downdrafts. In two cases with single-layer stratocumulus, turbulent transport of heat flux and of temperature variance are both negative in the cloud layer, suggesting an important role of large eddies. In contrast, in a case with weak cloud-top cooling, these quantities are positive in the ABL due to the heating from the surface. Based on observations and results of a mixed-layer model it is shown that the maxima of turbulent fluxes are, however, smaller than the jump of the net terrestrial radiation flux across the upper part of a cloud due to the (i) shallowness of the mixed layer and (ii) the presence of a downward entrainment heat flux. The mixed-layer model also shows that the buoyancy production of TKE is substantially smaller in stratocumulus over the Arctic sea ice compared to subtropics due to a smaller surface moisture flux and smaller decrease in specific humidity (or even humidity inversions) right above the cloud top. In a case of strong wind, wind shear shapes the ABL turbulent structure, especially over rough sea ice, despite the presence of a strong cloud-top cooling. In the presence of mid-level clouds, cloud-top radiative cooling and thus also TKE in the lowermost cloud layer are strongly reduced, and the ABL turbulent structure becomes governed by stability, i.e., by the surface–air temperature difference and wind speed. A comparison of slightly unstable and weakly stable cases shows a strong reduction of TKE due to increased stability even though the absolute value of wind speed was similar. In summary, the presented study documents vertical profiles of the ABL turbulence with a high resolution in a wide range of conditions. It can serve as a basis for turbulence closure evaluation and process studies in Arctic clouds.

Description: Springtime Arctic mixed-phase convection over open water in the Fram Strait as observed during the recent ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) field campaign is simulated at turbulence-resolving resolutions. The first objective is to assess the skill of large-eddy simulation (LES) in reproducing the observed mixed-phase convection. The second goal is to then use the model to investigate how aerosol modulates the way in which turbulent mixing and clouds transform the low-level air mass. The focus lies on the low-level thermal structure and lapse rate, the heating efficiency of turbulent entrainment, and the low-level energy budget. A composite case is constructed based on data collected by two research aircraft on 18 June 2017. Simulations are evaluated against independent datasets, showing that the observed thermodynamic, cloudy, and turbulent states are well reproduced. Sensitivity tests on cloud condensation nuclei (CCN) concentration are then performed, covering a broad range between pristine polar and polluted continental values. We find a significant response in the resolved mixed-phase convection, which is in line with previous LES studies. An increased CCN substantially enhances the depth of convection and liquid cloud amount, accompanied by reduced surface precipitation. Initializing with the in situ CCN data yields the best agreement with the cloud and turbulence observations, a result that prioritizes its measurement during field campaigns for supporting high-resolution modeling efforts. A deeper analysis reveals that CCN significantly increases the efficiency of radiatively driven entrainment in warming the boundary layer. The marked strengthening of the thermal inversion plays a key role in this effect. The low-level heat budget shifts from surface driven to radiatively driven. This response is accompanied by a substantial reduction in the surface energy budget, featuring a weakened flow of solar radiation into the ocean. Results are interpreted in the context of air–sea interactions, air mass transformations, and climate feedbacks at high latitudes.

Description: The presence of reactive bromine in polar regions is a widespread phenomenon that plays an important role in the photochemistry of the Arctic and Antarctic lower troposphere, including the destruction of ozone, the disturbance of radical cycles, and the oxidation of gaseous elemental mercury. The chemical mechanisms leading to the heterogeneous release of gaseous bromine compounds from saline surfaces are in principle well understood. There are, however, substantial uncertainties about the contribution of different potential sources to the release of reactive bromine, such as sea ice, brine, aerosols, and the snow surface, as well as about the seasonal and diurnal variation and the vertical distribution of reactive bromine. Here we use continuous long-term measurements of the vertical distribution of bromine monoxide (BrO) and aerosols at the two Antarctic sites Neumayer (NM) and Arrival Heights (AH), covering the periods of 2003–2021 and 2012–2021, respectively, to investigate how chemical and physical parameters affect the abundance of BrO. We find the strongest correlation between BrO and aerosol extinction (R=0.56 for NM and R=0.28 for AH during spring), suggesting that the heterogeneous release of Br2 from saline airborne particles (blowing snow and aerosols) is a dominant source for reactive bromine. Positive correlations between BrO and contact time of air masses, both with sea ice and the Antarctic ice sheet, suggest that reactive bromine is not only emitted by the sea ice surface but by the snowpack on the ice shelf and in the coastal regions of Antarctica. In addition, the open ocean appears to represent a source for reactive bromine during late summer and autumn when the sea ice extent is at its minimum. A source–receptor analysis based on back trajectories and sea ice maps shows that main source regions for BrO at NM are the Weddell Sea and the Filchner–Ronne Ice Shelf, as well as coastal polynyas where sea ice is newly formed. A strong morning peak in BrO frequently occurring during summer and that is particularly strong during autumn suggests a night-time build-up of Br2 by heterogeneous reaction of ozone on the saline snowpack in the vicinity of the measurement sites. We furthermore show that BrO can be sustained for at least 3 d while travelling across the Antarctic continent in the absence of any saline surfaces that could serve as a source for reactive bromine.

Description: 〈jats:p〉Comprehensive sampling of natural genetic diversity with metagenomics enables highly resolved insights into the interplay between ecology and evolution. However, resolving adaptive, neutral, or purifying processes of evolution from intrapopulation genomic variation remains a challenge, partly due to the sole reliance on gene sequences to interpret variants. Here, we describe an approach to analyze genetic variation in the context of predicted protein structures and apply it to a marine microbial population within the SAR11 subclade 1a.3.V, which dominates low-latitude surface oceans. Our analyses reveal a tight association between genetic variation and protein structure. In a central gene in nitrogen metabolism, we observe decreased occurrence of nonsynonymous variants from ligand-binding sites as a function of nitrate concentrations, revealing genetic targets of distinct evolutionary pressures maintained by nutrient availability. Our work yields insights into the governing principles of evolution and enables structure-aware investigations of microbial population genetics.〈/jats:p〉

Description: 〈jats:p〉Rapidly shrinking Arctic sea ice has had substantial impacts on the Earth system. Therefore, reliably estimating the Arctic sea-ice thickness (SIT) using a combination of available observations and numerical modeling is urgently needed. Here, for the first time, we assimilate the latest CryoSat-2 summer SIT data into a coupled ice-ocean model. In particular, an incremental analysis update scheme is implemented to overcome the discontinuity resulting from the combined assimilation of biweekly SIT and daily sea-ice concentration (SIC) data. Along with improved estimates of sea-ice volume, our SIT estimates corrected the overestimation of SIT produced by the reanalysis that assimilates only SIC in summer in areas where the sea ice is roughest and experiences strong deformation, e.g., around the Fram Strait and Greenland. This study suggests that the newly developed CryoSat-2 SIT product, when assimilated properly using our approach, has great potential for Arctic sea-ice simulation and prediction.〈/jats:p〉

Description: The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.

Description: Dissolved iron (dFe) availability limits the uptake of atmospheric CO2 by the Southern Ocean (SO) biological pump. Hence, any change in bioavailable dFe in this region can directly influence climate. On the basis of Fe uptake experiments with Phaeocystis antarctica, we show that the range of dFe bioavailability in natural samples is wider (〈1 to ~200% compared to free inorganic Fe′) than previously thought, with higher bioavailability found near glacial sources. The degree of bioavailability varied regardless of in situ dFe concentration and depth, challenging the consensus that sole dFe concentrations can be used to predict Fe uptake in modeling studies. Further, our data suggest a disproportionately major role of biologically mediated ligands and encourage revisiting the role of humic substances in influencing marine Fe biogeochemical cycling in the SO. Last, we describe a linkage between in situ dFe bioavailability and isotopic signatures that, we anticipate, will stimulate future research.

Description: Stable water isotopologues of snow, firn and ice cores provide valuable information on past climate variations. Yet single profiles are generally not suitable for robust climate reconstructions. Stratigraphic noise, introduced by the irregular deposition, wind-driven erosion and redistribution of snow, impacts the utility of high-resolution isotope records, especially in low-Accumulation areas. However, it is currently unknown how stratigraphic noise differs across the East Antarctic Plateau and how it is affected by local environmental conditions. Here, we assess the amount and structure of stratigraphic noise at seven sites along a 120 km transect on the plateau of Dronning Maud Land, East Antarctica. Replicated oxygen isotope records of 1 m length were used to estimate signal-To-noise ratios as a measure of stratigraphic noise at sites characterised by different accumulation rates (43-64 mm w.e. a-1), snow surface roughnesses and slope inclinations. While we found a high level of stratigraphic noise at all sites, there was also considerable variation between sites. At sastrugi-dominated sites, greater stratigraphic noise coincided with stronger surface roughnesses, steeper slopes and lower accumulation rates, probably related to increased wind speeds. These results provide a first step to modelling stratigraphic noise and might guide site selection and sampling strategies for future expeditions to improve high-resolution climate reconstructions from low-Accumulation regions.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2022 is an update of the previous version, GLODAPv2.2021 (Lauvset et al., 2021). The major changes are as follows: data from 96 new cruises were added, data coverage was extended until 2021, and for the first time we performed secondary quality control on all sulphur hexafluoride (SF6) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF6 data, which are now subjected to secondary quality control, and carbon data measured onboard the RV Knorr in the Indian Ocean in 1994–1995 which are now adjusted using CRM measurements made at the time. GLODAPv2.2022 includes measurements from almost 1.4 million water samples from the global oceans collected on 1085 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, CCl4, and SF6) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 96 new cruises were derived by comparing those data with the data from the 989 quality controlled cruises in the GLODAPv2.2021 data product using crossover analysis. SF6 data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO2) chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 μmol kg-1 in dissolved inorganic carbon, 4 μmol kg-1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO2 fugacity (fCO2), were not subjected to bias comparison or adjustments.

Description: 〈jats:p〉 Next week, the Convention on the Conservation of Antarctic Marine Living Resources (CCAMLR) convenes in Hobart, Tasmania, to examine the state of marine life in the Southern Ocean. As part of the Antarctic Treaty System, this convention entered into force in 1982, and its focus on the region’s environmental integrity has never been more important, given the increasing effects of climate change and commercial fishing. An important focus over the past 40 years has been Antarctic krill, 〈jats:italic〉Euphausia superba〈/jats:italic〉 (hereafter krill), a keystone species that helps to hold this marine ecosystem together. Climate and fishing stresses should prompt the CCAMLR to address whether management of krill fishing is at a level that protects the Southern Ocean from losing its overall balance of marine life and the oceanic processes that regulate global climate. 〈/jats:p〉

Description: Changes in Southern Ocean export production have broad biogeochemical and climatic implications. Specifically, iron fertilization likely increased subantarctic nutrient utilization and enhanced the efficiency of the biological pump during glacials. However, past export production in the subantarctic Southeast Pacific is poorly documented, and its connection to Fe fertilization, potentially related to Patagonian Ice Sheet dynamics is unknown. We report on biological productivity changes over the past 400 ka, based on a combination of 230Thxs-normalized and stratigraphy-based mass accumulation rates of biogenic barium, organic carbon, biogenic opal, and calcium carbonate as indicators of paleo-export production in a sediment core upstream of the Drake Passage. In addition, we use fluxes of iron and lithogenic material as proxies for terrigenous matter, and thus potential micronutrient supply. Stratigraphy-based mass accumulation rates are strongly influenced by bottom-current dynamics, which result in variable sediment focussing or winnowing at our site. Carbonate is virtually absent in the core, except during peak interglacial intervals of the Holocene, and Marine Isotope Stages (MIS) 5 and 11, likely caused by transient decreases in carbonate dissolution. All other proxies suggest that export production increased during most glacial periods, coinciding with high iron fluxes. Such augmented glacial iron fluxes at the core site were most likely derived from glaciogenic input from the Patagonian Ice Sheet promoting the growth of phytoplankton. Additionally, glacial export production peaks are also consistent with northward shifts of the Subantarctic and Polar Fronts, which positioned our site south of the Subantarctic Front and closer to silicic acid-rich waters of the Polar Frontal Zone, as well as a with a decrease in the diatom utilization of Si relative to nitrate under Fe-replete conditions. However, glacial export production near the Drake Passage was lower than in the Atlantic and Indian sectors of the Southern Ocean, which may relate to complete consumption of silicic acid in the study area. Our results underline the importance of micro-nutrient fertilization through lateral terrigenous input from South America rather than aeolian transport, and exemplify the role of frontal shifts and nutrient limitation for past productivity changes in the Pacific entrance to the Drake Passage.

Description: We developed a new version of the Alfred Wegener Institute Climate Model (AWI-CM3), which has higher skills in representing the observed climatology and better computational efficiency than its predecessors. Its ocean component FESOM2 has the multi-resolution functionality typical for unstructured-mesh models while still featuring a scalability and efficiency similar to regular-grid models. The atmospheric component OpenIFS (CY43R3) enables the use of latest developments in the numerical weather prediction community in climate sciences. In this paper we describe the coupling of the model components and evaluate the model performance on a variable resolution (25–125 km) ocean mesh and a 61 km atmosphere grid, which serves as a reference and starting point for other on-going research activities with AWI-CM3. This includes the exploration of high and variable resolution, the development of a full Earth System Model as well as the creation of a new sea ice prediction system. At this early development stage and with the given coarse to medium resolutions, the model already features above CMIP6-average skills in representing the climatology and competitive model throughput. Finally we identify remaining biases and suggest further improvements to be made to the model.

Description: Arctic river deltas and deltaic near-shore zones represent important land–ocean transition zones influencing sediment dynamics and nutrient fluxes from permafrost-affected terrestrial ecosystems into the coastal Arctic Ocean. To accurately model fluvial carbon and freshwater export from rapidly changing river catchments as well as assess impacts of future change on the Arctic shelf and coastal ecosystems, we need to understand the sea floor characteristics and topographic variety of the coastal zones. To date, digital bathymetrical data from the poorly accessible, shallow, and large areas of the eastern Siberian Arctic shelves are sparse. We have digitized bathymetrical information for nearly 75 000 locations from large-scale (1:25 000–1:500 000) current and historical nautical maps of the Lena Delta and the Kolyma Gulf region in northeastern Siberia. We present the first detailed and seamless digital models of coastal zone bathymetry for both delta and gulf regions in 50 and 200 m spatial resolution. We validated the resulting bathymetry layers using a combination of our own water depth measurements and a collection of available depth measurements, which showed a strong correlation (r〉0.9). Our bathymetrical models will serve as an input for a high-resolution coupled hydrodynamic–ecosystem model to better quantify fluvial and coastal carbon fluxes to the Arctic Ocean, but they may be useful for a range of other studies related to Arctic delta and near-shore dynamics such as modeling of submarine permafrost, near-shore sea ice, or shelf sediment transport. The new digital high-resolution bathymetry products are available on the PANGAEA data set repository for the Lena Delta (https://doi.org/10.1594/PANGAEA.934045; Fuchs et al., 2021a) and Kolyma Gulf region (https://doi.org/10.1594/PANGAEA.934049; Fuchs et al., 2021b), respectively. Likewise, the depth validation data are available on PANGAEA as well (https://doi.org/10.1594/PANGAEA.933187; Fuchs et al., 2021c).

Description: Hyrrokkin sarcophaga is a parasitic foraminifera that is commonly found in cold-water coral reefs where it infests the file clam Acesta excavata and the scleractinian coral Desmophyllum pertusum (formerly known as Lophelia pertusa). Here, we present measurements of the trace element and isotopic composition of these parasitic foraminifera, analyzed by inductively coupled optical emission spectrometry (ICP-OES), electron probe microanalysis (EPMA) and mass spectrometry (gas-source MS and inductively-coupled-plasma MS). Our results reveal that the geochemical signature of H. sarcophaga depends on the host organism it infests. Sr / Ca ratios are 1.1 mmol mol−1 higher in H. sarcophaga that infest D. pertusum, which could be an indication that dissolved host carbonate material is utilized in shell calcification, given that the aragonite of D. pertusum has a naturally higher Sr concentration compared to the calcite of A. excavata. Similarly, we measure 3.1 ‰ lower δ13C and 0.25 ‰ lower δ18O values in H. sarcophaga that lived on D. pertusum, which might be caused by the direct uptake of the host's carbonate material with a more negative isotopic composition or different pH regimes in these foraminifera (pH can exert a control on the extent of CO2 hydration/hydroxylation) due to the uptake of body fluids of the host. We also observe higher Mn / Ca ratios in foraminifera that lived on A. excavata but did not penetrate the host shell compared to specimen that penetrated the shell, which could be interpreted as a change in food source, changes in the calcification rate, Rayleigh fractionation or changing oxygen conditions. While our measurements provide an interesting insight into the calcification process of this unusual foraminifera, these data also indicate that the geochemistry of this parasitic foraminifera is unlikely to be a reliable indicator of paleoenvironmental conditions using Sr / Ca, Mn / Ca, δ18O or δ13C unless the host organism is known and its geochemical composition can be accounted for.

Description: In this work we have performed a detailed study of vectors to ore to a representative volcanic-rock-hosted replacive volcanogenic massive sulfide (VMS) deposit located in the northern Iberian Pyrite Belt (Spain), the Aguas Teñidas deposit. The investigated vectors include the following: (1) mineralogical zoning, (2) host sequence characterization and mineralized unit identification based on whole rock geochemistry discrimination diagrams, (3) study of the characteristics and behaviour of whole rock geochemical anomalies around the ore (e.g. alteration-related compositional changes, characteristics and extent of geochemical halos of indicative elements such as Cu, Zn, Pb, Sb, Tl, and Ba around the deposit), and (4) application of portable X-ray fluorescence (p-XRF) analysis to the detection of the previous vectors. In the footwall, a concentric cone-shaped hydrothermal alteration zone bearing the stockwork passes laterally, from core to edge, from quartz (only local) to chlorite–quartz, sericite–chlorite–quartz, and sericite–quartz alteration zones. The hydrothermal alteration is also found in the hanging wall despite being tectonically allochthonous to the orebody: a proximal sericite alteration zone is followed by a more distal albite-rich one. Whole rock major elements show an increase in alteration indexes (e.g. AI, CCPI) towards the mineralization, a general SiO2 enrichment, and FeO enrichment as well as K2O and Na2O depletion towards the centre of the hydrothermal system, with MgO showing a less systematic behaviour. K2O and Na2O leached from the centre of the system are transported and deposited in more external areas. Copper, Pb, and Zn produce proximal anomalies around mineralized areas, with the more mobile Sb, Tl, and Ba generating wider halos. Whereas Sb and Tl halos form around all mineralized areas, Ba anomalies are restricted to areas around the massive sulfide body. Our results show that proposed vectors, or adaptations designed to overcome p-XRF limitations, can be confidently used by analysing unprepared hand specimens, including the external rough curved surface of drill cores. The data presented in this work are not only applicable to VMS exploration in the Iberian Pyrite Belt, but on a broader scale they will also contribute to improving our general understanding of vectors to ore in replacive-type VMS deposits.

Description: We combine satellite data products to provide a first and general overview of the physical sea ice conditions along the drift of the international Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition and a comparison with previous years (2005–2006 to 2018–2019). We find that the MOSAiC drift was around 20 % faster than the climatological mean drift, as a consequence of large-scale low-pressure anomalies prevailing around the Barents–Kara–Laptev sea region between January and March. In winter (October–April), satellite observations show that the sea ice in the vicinity of the Central Observatory (CO; 50 km radius) was rather thin compared to the previous years along the same trajectory. Unlike ice thickness, satellite-derived sea ice concentration, lead frequency and snow thickness during winter months were close to the long-term mean with little variability. With the onset of spring and decreasing distance to the Fram Strait, variability in ice concentration and lead activity increased. In addition, the frequency and strength of deformation events (divergence, convergence and shear) were higher during summer than during winter. Overall, we find that sea ice conditions observed within 5 km distance of the CO are representative for the wider (50 and 100 km) surroundings. An exception is the ice thickness; here we find that sea ice within 50 km radius of the CO was thinner than sea ice within a 100 km radius by a small but consistent factor (4 %) for successive monthly averages. Moreover, satellite acquisitions indicate that the formation of large melt ponds began earlier on the MOSAiC floe than on neighbouring floes.

Description: Across the Qinghai–Tibet Plateau (QTP) there is a narrow engineering corridor with widely distributed slopes called the Qinghai–Tibet Engineering Corridor (QTEC), where a variety of important infrastructures are concentrated. These facilities are transportation routes for people, materials, energy, etc. from inland China to the Tibet Autonomous Region. From Golmud to Lhasa, the engineering corridor covers 632 km of permafrost containing the densely developed Qinghai–Tibet Railway and Qinghai–Tibet Highway, as well as power and communication towers. Slope failure in permafrost regions, caused by permafrost degradation, ground ice melting, etc., affects the engineering construction and permafrost environments in the QTEC. We implement a variety of sensors to monitor the hydrological and thermal deformation between permafrost slopes and permafrost engineering projects in the corridor. In addition to soil temperature and moisture sensors, the global navigation satellite system (GNSS), terrestrial laser scanning (TLS), and unmanned aerial vehicles (UAVs) were adopted to monitor the spatial distribution and changes in thermal deformation. An integrated dataset of hydrological and thermal deformation in permafrost engineering and slopes in the QTEC from the 1950s to 2020, including meteorological and ground observations, TLS point cloud data, and RGB and thermal infrared (TIR) images, can be of great value for estimating the hydrological and thermal impact and stability between engineering and slopes under the influence of climate change and engineering disturbance. The dataset and code were uploaded to the Zenodo repository and can be accessed through https://zenodo.org/communities/qtec (last access: 23 June 2021), including meteorological and ground observations at https://doi.org/10.5281/zenodo.5009871 (Luo et al., 2020d), TLS measurements at https://doi.org/10.5281/zenodo.5009558 (Luo et al., 2020a), UAV RGB and TIR images at https://doi.org/10.5281/zenodo.5016192 (Luo et al., 2020b), and R code for permafrost indices and visualisation at https://doi.org/10.5281/zenodo.5002981 (Luo et al., 2020c).

Description: We present the first incorporation of the Common Representative Intermediates version 2.2 tropospheric chemistry mechanism, CRI v2.2, combined with stratospheric chemistry, into the global chemistry–climate United Kingdom Chemistry and Aerosols (UKCA) model to give the CRI-Strat 2 mechanism. A rigorous comparison of CRI-Strat 2 with the earlier version, CRI-Strat, is performed in UKCA in addition to an evaluation of three mechanisms, CRI-Strat 2, CRI-Strat and the standard UKCA chemical mechanism, StratTrop v1.0, against a wide array of surface and airborne chemical data. CRI-Strat 2 comprises a state-of-the-art isoprene scheme, optimized against the Master Chemical Mechanism v3.3.1, which includes isoprene peroxy radical isomerization, HOx recycling through the addition of photolabile hydroperoxy aldehydes (HPALDs), and isoprene epoxy diol (IEPOX) formation. CRI-Strat 2 also features updates to several rate constants for the inorganic chemistry, including the reactions of inorganic nitrogen and O(1D). The update to the isoprene chemistry in CRI-Strat 2 increases OH over the lowest 500 m in tropical forested regions by 30 %–50 % relative to CRI-Strat, leading to an improvement in model–observation comparisons for surface OH and isoprene relative to CRI-Strat and StratTrop. Enhanced oxidants also cause a 25 % reduction in isoprene burden and an increase in oxidation fluxes of isoprene and other biogenic volatile organic compounds (BVOCs) at low altitudes with likely impacts on subsequent aerosol formation, atmospheric lifetime, and climate. By contrast, updates to the rate constants of O(1D) with its main reactants relative to CRI-Strat reduces OH in much of the free troposphere, producing a 2 % increase in the methane lifetime, and increases the tropospheric ozone burden by 8 %, primarily from reduced loss via O(1D)+H2O. The changes to inorganic nitrogen reaction rate constants increase the NOx burden by 4 % and shift the distribution of nitrated species closer to that simulated by StratTrop. CRI-Strat 2 is suitable for multi-decadal model integrations and the improved representation of isoprene chemistry provides an opportunity to explore the consequences of HOx recycling in the United Kingdom Earth System Model (UKESM1). This new mechanism will enable a re-evaluation of the impact of BVOCs on the chemical composition of the atmosphere and further probe the feedback between the biosphere and the climate.

Description: The use of mass flow simulations in volcanic hazard zonation and mapping is often limited by model complexity (i.e. uncertainty in correct values of model parameters), a lack of model uncertainty quantification, and limited approaches to incorporate this uncertainty into hazard maps. When quantified, mass flow simulation errors are typically evaluated on a pixel-pair basis, using the difference between simulated and observed (“actual”) map-cell values to evaluate the performance of a model. However, these comparisons conflate location and quantification errors, neglecting possible spatial autocorrelation of evaluated errors. As a result, model performance assessments typically yield moderate accuracy values. In this paper, similarly moderate accuracy values were found in a performance assessment of three depth-averaged numerical models using the 2012 debris avalanche from the Upper Te Maari crater, Tongariro Volcano, as a benchmark. To provide a fairer assessment of performance and evaluate spatial covariance of errors, we use a fuzzy set approach to indicate the proximity of similarly valued map cells. This “fuzzification” of simulated results yields improvements in targeted performance metrics relative to a length scale parameter at the expense of decreases in opposing metrics (e.g. fewer false negatives result in more false positives) and a reduction in resolution. The use of this approach to generate hazard zones incorporating the identified uncertainty and associated trade-offs is demonstrated and indicates a potential use for informed stakeholders by reducing the complexity of uncertainty estimation and supporting decision-making from simulated data.

Description: In this study, we analyse the climatology of ionosphere over Nepal based on GPS-derived vertical total electron content (VTEC) observed from four stations as defined in Table 1: KKN4 (27.80∘ N, 85.27∘ E), GRHI (27.95∘ N, 82.49∘ E), JMSM (28.80∘ N, 83.74∘ E) and DLPA (28.98∘ N, 82.81∘ E) during the years 2008 to 2018. The study illustrates the diurnal, monthly, annual, seasonal and solar cycle variations in VTEC during all times of solar cycle 24. The results clearly reveal the presence of equinoctial asymmetry in TEC, which is more pronounced in maximum phases of solar cycle in the year 2014 at KKN4 station, followed by descending, ascending and minimum phases. Diurnal variations in VTEC showed the short-lived day minimum which occurs between 05:00 to 06:00 LT (local time) at all the stations considered, with diurnal peaks between 12:00 and 15:00 LT. The maximum value of TEC is observed more often during the spring equinox than the autumn equinox, with a few asymmetries. Seasonal variation in TEC is observed to be a manifestation of variations in solar flux, particularly regarding the level of solar flux in consecutive solstices.

Description: We report on applications of the ultraviolet-light-emitting-diode-based incoherent broadband cavity-enhanced absorption spectroscopy (UV-LED-IBBCEAS) technique for optical monitoring of HONO, NO2 and CH2O in a simulation chamber. Performance intercomparison of UV-LED-IBBCEAS with a wet chemistry-based NitroMAC sensor and a Fourier transform infrared (FTIR) spectrometer has been carried out on real-time simultaneous measurement of HONO, NO2 and CH2O concentrations during the reaction of NO2 with H2O vapour in CESAM (French acronym for Experimental Multiphasic Atmospheric Simulation Chamber). The 1σ (signal-to-noise ratio (SNR) = 1) detection limits of 112 pptv for NO2, 56 pptv for HONO and 41 ppbv for CH2O over 120 s were found for the UV-LED-IBBCEAS measurement. On the contrary to many set-ups where cavities are installed outside the simulation chamber, we describe here an original in situ permanent installation. The intercomparison results demonstrate that IBBCEAS is a very well suitable technique for in situ simultaneous measurements of multiple chemically reactive species with high sensitivity and high precision even if the absorption bands of these species are overlapped. It offers excellent capacity for non-invasive optical monitoring of chemical reactions without any perturbation. For the application to simulation chambers, it has the advantage to provide a spatially integrated measurement across the reactor and hence to avoid point-sampling-related artefacts.

Description: Assimilation of weather radar measurements including radar reflectivity and radial wind data has been operational at the Deutscher Wetterdienst, with a diagonal observation error (OE) covariance matrix. For an implementation of a full OE covariance matrix, the statistics of the OE have to be a priori estimated, for which the Desroziers method has been often used. However, the resulted statistics consists of contributions from different error sources and are difficult to interpret. In this work, we use an approach that is based on samples for truncation error in radar observation space to approximate the representation error due to unresolved scales and processes (RE) and compare its statistics with the OE statistics estimated by the Desroziers method. It is found that the statistics of the RE help the understanding of several important features in the variances and correlation length scales of the OE for both reflectivity and radial wind data and the other error sources from the microphysical scheme, radar observation operator and the superobbing technique may also contribute, for instance, to differences among different elevations and observation types. The statistics presented here can serve as a guideline for selecting which observations are assimilated and for assignment of the OE covariance matrix that can be diagonal or full and correlated.

Description: Determining controls on the temperature sensitivity of heterotrophic soil respiration remains critical to incorporating soil–climate feedbacks into climate models. Most information on soil respiratory responses to temperature comes from laboratory incubations of isolated soils and typically subsamples of individual horizons. Inconsistencies between field and laboratory results may be explained by microbial priming supported by cross-horizon exchange of labile C or N. Such exchange is feasible in intact soil profiles but is absent when soils are isolated from surrounding depths. Here we assess the role of soil horizon connectivity, by which we mean the degree to which horizons remain layered and associated with each other as they are in situ, on microbial C and N substrate use and its relationship to the temperature sensitivity of respiration. We accomplished this by exploring changes in C : N, soil organic matter composition (via C : N, amino acid composition and concentration, and nuclear magnetic resonance spectroscopy), and the δ13C of respiratory CO2 during incubations of organic horizons collected across boreal forests in different climate regions where soil C and N compositions differ. The experiments consisted of two treatments: soil incubated (1) with each organic horizon separately and (2) as a whole organic profile, permitting cross-horizon exchange of substrates during the incubation. The soils were incubated at 5 and 15 ∘C for over 430 d. Enhanced microbial use of labile C-rich, but not N-rich, substrates were responsible for enhanced, whole-horizon respiratory responses to temperature relative to individual soil horizons. This impact of a labile C priming mechanism was most emergent in soils from the warmer region, consistent with these soils' lower C bioreactivity relative to soils from the colder region. Specifically, cross-horizon exchange within whole soil profiles prompted increases in mineralization of carbohydrates and more 13C-enriched substrates and increased soil respiratory responses to warming relative to soil horizons incubated in isolation. These findings highlight that soil horizon connectivity can impact microbial substrate use in ways that affect how soil effluxes of CO2 are controlled by temperature. The degree to which this mechanism exerts itself in other soils remains unknown, but these results highlight the importance of understanding mechanisms that operate in intact soil profiles – only rarely studied – in regulating a key soil–climate feedback.

Description: Plumes from the boreal spring biomass burning (BB) in northern peninsular Southeast Asia (nPSEA) are lifted into the subtropical jet stream and transported and deposited across nPSEA, South China, Taiwan and even the western North Pacific Ocean. This paper as part of the Seven SouthEast Asian Studies (7-SEAS) project effort attempts to improve the chemical weather prediction capability of the Weather Research and Forecasting coupled with the Community Multiscale for Air Quality (WRF–CMAQ) model over a vast region, from the mountainous near-source burning sites at nPSEA to its downwind region. Several sensitivity analyses of plume rise are compared in the paper, and it is discovered that the initial vertical allocation profile of BB plumes and the plume rise module (PLMRIM) are the main reasons causing the inaccuracies of the WRF–CMAQ simulations. The smoldering emission from the Western Regional Air Partnership (WRAP) empirical algorithm included has improved the accuracies of PM10, O3 and CO at the source. The best performance at the downwind sites is achieved with the inline PLMRIM, which accounts for the atmospheric stratification at the mountainous source region with the FINN burning emission dataset. Such a setup greatly improves not only the BB aerosol concentration prediction over near-source and receptor ground-based measurement sites but also the aerosol vertical distribution and column aerosol optical depth of the BB aerosol along the transport route. The BB aerosols from nPSEA are carried by the subtropical westerlies in the free troposphere to the western North Pacific, while BB aerosol has been found to interact with the local pollutants in the Taiwan region through three conditions: (a) overpassing western Taiwan and entering the central mountain area, (b) mixing down to western Taiwan, (c) transport of local pollutants upwards and mixing with a BB plume on higher ground. The second condition, which involves the prevailing high-pressure system from Asian cold surge, is able to impact most of the population in Taiwan.

Description: The German Bight located within the central North Sea is a hydro- and morphodynamically highly complex system of estuaries, barrier islands, and part of the world's largest coherent tidal flats, the Wadden Sea. To identify and understand challenges faced by coastal stakeholders, such as harbor operators or governmental agencies, to maintain waterways and employ numerical models for further analyses, it is imperative to have a consistent database for both bathymetry and surface sedimentology. Current commercial and public data products are insufficient in spatial and temporal resolution and coverage for recent analysis methods. Thus, this first part of a two-part publication series of the German joint project EasyGSH-DB describes annual bathymetric digital terrain models at a 10 m gridded resolution for the German North Sea coast and German Bight from 1996 to 2016 (Sievers et al., 2020a, https://doi.org/10.48437/02.2020.K2.7000.0001), as well as surface sedimentological models of discretized cumulative grain size distribution functions for 1996, 2006, and 2016 on 100 m grids (Sievers et al., 2020b, https://doi.org/10.48437/02.2020.K2.7000.0005). Furthermore, basic morphodynamic and sedimentological processing analyses, such as the estimation of, for example, bathymetric stability or surface maps of sedimentological parameters, are provided (Sievers et al., 2020a, b, see respective download links).

Description: The Phanerozoic tectonothermal evolution of the SW slope of the East European Platform (EEP) in Poland is reconstructed by means of thermal maturity, low-temperature thermochronometry, and thermal modelling. We provide a set of new thermochronometric data and integrate stratigraphic and thermal maturity information to constrain the burial and thermal history of sediments. Apatite fission track (AFT) analysis and zircon (U-Th)/He (ZHe) thermochronology have been carried out on samples of sandstones, bentonites, diabase, and crystalline basement rocks collected from 17 boreholes located in central and NE Poland. They penetrated sedimentary cover of the EEP subdivided from the north to south into the Baltic, Podlasie, and Lublin basins. The average ZHe ages from Proterozoic basement rocks as well as Ordovician to Silurian bentonites and Cambrian to lower Carboniferous sandstones range from 848 ± 81 to 255 ± 22 Ma with a single early Permian age of 288 Ma, corresponding to cooling after a thermal event. The remaining ZHe ages represent partial reset or source ages. The AFT ages of samples are dispersed in the range of 235.8 ± 17.3 Ma (Middle Triassic) to 42.1 ± 11.1 Ma (Paleogene) providing a record of Mesozoic and Cenozoic cooling. The highest frequency of the AFT ages is in the Jurassic and Early Cretaceous prior to Alpine basin inversion. Thermal maturity results are consistent with the SW-ward increase of the Paleozoic and Mesozoic sediments thickness. An important break in a thermal maturity profile exists across the base Permian–Mesozoic unconformity. Thermal modelling showed that significant heating of Ediacaran to Carboniferous sedimentary successions occurred before the Permian with maximum paleotemperatures in the earliest and latest Carboniferous for Baltic–Podlasie and Lublin basins, respectively. The results obtained suggest an important role of early Carboniferous uplift and exhumation at the SW margin of the EEP. The SW slope of the latter was afterward overridden in the Lublin Basin by the Variscan orogenic wedge. Its tectonic loading interrupted Carboniferous uplift and caused resumption of sedimentation in the late Viséan. Consequently, a thermal history of the Lublin Basin is different from that in the Podlasie and Baltic basins but similar to other sections of the Variscan foreland, characterized by maximum burial at the end of Carboniferous. The Mesozoic thermal history was characterized by gradual cooling from peak temperatures at the transition from Triassic to Jurassic due to decreasing heat flow. Burial caused maximum paleotemperatures in the SW part of the study area, where the EEP was covered by an extensive sedimentary pile. However, further NE, due to low temperatures caused by shallow burial, the impact of fluids can be detected by vitrinite reflectance, illite/smectite, and thermochronological data. Our new results emphasize the importance of using multiple low-temperature thermochronometers and thermal modelling in connection with thermal maturity analysis to elucidate the near-surface evolution of platform margins.

Description: Snow stands out from materials at the Earth’s surface owing to its unique optical properties. Snow optical properties are sensitive to the snow microstructure, triggering potent climate feedbacks. The impacts of snow microstructure on its optical properties such as reflectance are, to date, only partially understood. However, precise modelling of snow reflectance, particularly bidirectional reflectance, are required in many problems, e.g. to correctly process satellite data over snow-covered areas. This study presents a dataset that combines bidirectional reflectance measurements over 500–2500 nm and the X-ray tomography of the snow microstructure for three snow samples of two different morphological types. The dataset is used to evaluate the stereological approach from Malinka (2014) that relates snow optical properties to the chord length distribution in the snow microstructure. The mean chord length and specific surface area (SSA) retrieved with this approach from the albedo spectrum and those measured by the X-ray tomography are in excellent agreement. The analysis of the 3D images has shown that the random chords of the ice phase obey the gamma distribution with the shape parameter m taking the value approximately equal to or a little greater than 2. For weak and intermediate absorption (high and medium albedo), the simulated bidirectional reflectances reproduce the measured ones accurately but tend to slightly overestimate the anisotropy of the radiation. For such absorptions the use of the exponential law for the ice chord length distribution instead of the one measured with the X-ray tomography does not affect the simulated reflectance. In contrast, under high absorption (albedo of a few percent), snow microstructure and especially facet orientation at the surface play a significant role in the reflectance, particularly at oblique viewing and incidence.

Description: Attribution in its general definition aims to quantify drivers of change in a system. According to IPCC Working Group II (WGII) a change in a natural, human or managed system is attributed to climate change by quantifying the difference between the observed state of the system and a counterfactual baseline that characterizes the system's behavior in the absence of climate change, where “climate change refers to any long-term trend in climate, irrespective of its cause” (IPCC, 2014). Impact attribution following this definition remains a challenge because the counterfactual baseline, which characterizes the system behavior in the hypothetical absence of climate change, cannot be observed. Process-based and empirical impact models can fill this gap as they allow us to simulate the counterfactual climate impact baseline. In those simulations, the models are forced by observed direct (human) drivers such as land use changes, changes in water or agricultural management but a counterfactual climate without long-term changes. We here present ATTRICI (ATTRIbuting Climate Impacts), an approach to construct the required counterfactual stationary climate data from observational (factual) climate data. Our method identifies the long-term shifts in the considered daily climate variables that are correlated to global mean temperature change assuming a smooth annual cycle of the associated scaling coefficients for each day of the year. The produced counterfactual climate datasets are used as forcing data within the impact attribution setup of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a). Our method preserves the internal variability of the observed data in the sense that factual and counterfactual data for a given day have the same rank in their respective statistical distributions. The associated impact model simulations allow for quantifying the contribution of climate change to observed long-term changes in impact indicators and for quantifying the contribution of the observed trend in climate to the magnitude of individual impact events. Attribution of climate impacts to anthropogenic forcing would need an additional step separating anthropogenic climate forcing from other sources of climate trends, which is not covered by our method.

Description: Biocides used in film protection products leaching from facades are known to be a potential threat to the environment. This study identifies individual sources and entry pathways in a small-scale urban area. We investigate emissions of commonly used biocides (terbutryn, diuron, and octylisothiazolinone – OIT) and some of their transformation products (TPs; diuron-desmethyl, terbumeton, terbuthylazine-2-hydroxy, and terbutryn-desethyl) from a 2 ha residential area 13 years after construction has ended. Sampling utilizes existing urban water infrastructure representative for decentralized storm water management in central and northern Europe and applies a two-step approach to (a) determine the occurrence of biocides above water quality limits (i.e., predicted no-effect concentration, PNEC) and (b) identify source areas and characterize entry pathways into surface and groundwater. Monitoring focuses on the analysis of selected biocides and TPs by liquid chromatography–mass spectrometry/mass spectrometry (LC-MS/MS) in water samples taken from facades, rainwater pipes, drainage, and storm water infiltration systems. In standing water in a swale, we found high concentrations of diuron (174 ng L−1) and terbutryn (40 ng L−1) above PNEC for surface water. We confirmed expected sources, i.e., facades. Sampling of rain downpipes from flat roofs identified additional sources of all biocides and two TPs of terbutryn and one TP of diuron. Diuron and terbutryn were found in three drainage pipes representing different entry pathways of biocides. In one drainage pipe collecting road runoff, only diuron-desmethyl and terbutryn-desethyl were detected. In two other drainage pipes collecting infiltrated water through soil, terbuthylazine-2-hydroxy was additionally detected. One of the pipes collecting infiltrated water through soil concentration showed the highest concentrations of terbutryn and two of its TPs (terbutryn-desethyl and terbuthylazine-2-hydroxy). This suggests a high leaching potential of terbutryn. The applied two-step approach determined sources and pathways of biocide and their TPs. This study contributes to expanding knowledge on their entry and distribution and, thus, eventually towards reducing emissions.

Description: Regional anomalies of steric sea level are either due to redistribution of heat and freshwater anomalies or due to ocean–atmosphere buoyancy fluxes. Interannual to decadal variability in sea level across the tropical Pacific is mainly due to steric variations driven by wind stress anomalies. The importance of air–sea buoyancy fluxes is less clear. We use a global, eddy-permitting ocean model and a series of sensitivity experiments with quasi-climatological momentum and buoyancy fluxes to identify the contribution of buoyancy fluxes for interannual to decadal sea level variability in the tropical Pacific. We find their contribution on interannual timescales to be strongest in the central tropical Pacific at around a 10∘ latitude in both hemispheres and also relevant in the very east of the tropical domain. Buoyancy-flux-forced anomalies are correlated with variations driven by wind stress changes, but their effect on the prevailing anomalies and the importance of heat and freshwater fluxes vary locally. In the eastern tropical basin, interannual sea level variability is amplified by anomalous heat fluxes, while the importance of freshwater fluxes is small, and neither has any impact on decadal timescales. In the western tropical Pacific, the variability on interannual and decadal timescales is dampened by both heat and freshwater fluxes. The mechanism involves westward-propagating Rossby waves that are triggered during El Niño–Southern Oscillation (ENSO) events by anomalous buoyancy fluxes in the central tropical Pacific and counteract the prevailing sea level anomalies once they reach the western part of the basin.

Description: The Nioghalvfjerdsfjorden glacier (also known as the 79∘ North Glacier) drains approximately 8 % of the Greenland Ice Sheet. Supraglacial lakes (SGLs), or surface melt ponds, are a persistent summertime feature and are thought to drain rapidly to the base of the glacier and influence seasonal ice velocity. However, seasonal development and spatial distribution of SGLs in the north-east of Greenland are poorly understood, leaving a substantial error in the estimate of meltwater and its impacts on ice velocity. Using results from an automated detection of melt ponds, atmospheric and surface mass balance modelling, and reanalysis products, we investigate the role of specific climatic conditions in melt onset, extent, and duration from 2016 to 2019. The summers of 2016 and 2019 were characterised by above-average air temperatures, particularly in June, as well as a number of rainfall events, which led to extensive melt ponds to elevations up to 1600 m. Conversely, 2018 was particularly cold, with a large accumulated snowpack, which limited the development of lakes to altitudes less than 800 m. There is evidence of inland expansion and increases in the total area of lakes compared to the early 2000s, as projected by future global warming scenarios.

Description: Coincident wide-angle and multi-channel seismic data acquired within the scope of the PAMELA Moz3-5 project allow us to reconsider the formation mechanism of East African margins offshore of southern Mozambique. This study specifically focuses on the sedimentary and deep-crustal architecture of the Limpopo margin (LM) that fringes the eastern edge of the Mozambique’s Coastal Plain (MCP) and its offshore southern prolongation the North Natal Valley (NNV). It relies primarily on the MZ3 profile that runs obliquely from the northeastern NNV towards the Mozambique basin (MB) with additional inputs from a tectonostratigraphy analysis of industrial onshore–offshore seismic lines and nearby or crossing velocity models from companion studies. Over its entire N–S extension the LM appears segmented into (1) a western domain that shows the progressive eastward crustal thinning and termination of the MCP/NNV continental crust and its overlying pre-Neocomian volcano-sedimentary basement and (2) a central corridor of anomalous crust bounded to the east by the Mozambique fracture zone (MFZ) and the oceanic crust of the MB. A prominent basement high marks the boundary between these two domains. Its development was most probably controlled by a steep and deeply rooted fault, i.e., the Limpopo fault. We infer that strike-slip or slightly transtensional rifting occurred along the LM and was accommodated along this Limpopo fault. At depth we propose that ductile shearing was responsible for the thinning of the continental crust and an oceanward flow of lower crustal material. This process was accompanied by intense magmatism that extruded to form the volcanic basement and gave the corridor its peculiar structure and mixed nature. The whole region remained at a relative high level during the rifting period and a shallow marine environment dominated the pre-Neocomian period during the early phase of continent–ocean interaction. It is only some time after break-up in the MB and the initiation of the MFZ that decoupling occurred between the MCP/NNV and the corridor, allowing for the latter to subside and become covered by deep marine sediments. A scenario for the early evolution and formation of the LM is proposed taking into account both recent kinematic and geological constraints. It implies that no or little change in extensional direction occurred between the intra-continental rifting and subsequent phase of continent–ocean interaction.

Description: Models are an important tool to predict Earth system dynamics. An accurate prediction of future states of ecosystems depends on not only model structures but also parameterizations. Model parameters can be constrained by data assimilation. However, applications of data assimilation to ecology are restricted by highly technical requirements such as model-dependent coding. To alleviate this technical burden, we developed a model-independent data assimilation (MIDA) module. MIDA works in three steps including data preparation, execution of data assimilation, and visualization. The first step prepares prior ranges of parameter values, a defined number of iterations, and directory paths to access files of observations and models. The execution step calibrates parameter values to best fit the observations and estimates the parameter posterior distributions. The final step automatically visualizes the calibration performance and posterior distributions. MIDA is model independent, and modelers can use MIDA for an accurate and efficient data assimilation in a simple and interactive way without modification of their original models. We applied MIDA to four types of ecological models: the data assimilation linked ecosystem carbon (DALEC) model, a surrogate-based energy exascale earth system model: the land component (ELM), nine phenological models and a stand-alone biome ecological strategy simulator (BiomeE). The applications indicate that MIDA can effectively solve data assimilation problems for different ecological models. Additionally, the easy implementation and model-independent feature of MIDA breaks the technical barrier of applications of data–model fusion in ecology. MIDA facilitates the assimilation of various observations into models for uncertainty reduction in ecological modeling and forecasting.

Description: The northern Antarctic Peninsula (NAP) is a highly dynamic transitional zone between the subpolar-polar and oceanic-coastal environments, and it is located in an area affected by intense climate change, including intensification and spatial shifts of the westerlies as well as atmospheric and oceanic warming. In the NAP area, the water masses originate mainly from the Bellingshausen and Weddell seas, which create a marked regional dichotomy thermohaline characteristic. Although the NAP area has relatively easy access when compared to other Southern Ocean environments, our understanding of the water masses' distribution and the dynamical processes affecting the variability of the region is still limited. That limitation is closely linked to the sparse data coverage, as is commonly the case in most Southern Ocean environments. This work provides a novel seasonal three-dimensional high-resolution hydrographic gridded data set for the NAP (version 1), namely the NAPv1.0. Hydrographic measurements from 1990 to 2019 comprising data collected by conductivity, temperature, depth (CTD) casts; sensors from the Marine Mammals Exploring the Oceans Pole to Pole (MEOP) consortium; and Argo floats have been optimally interpolated to produce maps of in situ temperature, practical salinity, and dissolved oxygen at ∼ 10 km spatial resolution and 90 depth levels. The water masses and oceanographic features in this regional gridded product are more accurate than other climatologies and state estimate products currently available. The data sets are available in netCDF format at https://doi.org/10.5281/zenodo.4420006 (Dotto et al., 2021). The novel and comprehensive data sets presented here for the NAPv1.0 product are a valuable tool to be used in studies addressing climatological changes in the unique NAP region since they provide accurate initial conditions for ocean models and improve the end of the 20th- and early 21st-century ocean mean-state representation for that area.

Description: Time shapes every single human–environment relationship and is inherent in 21st-century global challenges such as climate change and the urgent move towards global sustainability. Nonetheless, the concept of time is still insufficiently addressed in climate justice debates. This paper aims to help fill this gap by presenting empirical results about experiences of climate change in farming communities in Tamil Nadu, South India, and fishing communities on the Lofoten Islands, Norway. With the help of the five dimensions of affectedness, rhythms and rituals, slow motion, care, and health and well-being, it exemplifies how time matters to issues of climate injustices faced by the communities. The paper promotes a qualitative understanding of time and climate change. Thereby, it may stimulate greater relatability to climate change, as well as discussion likely to lead to conceptual advances.

Description: A method is proposed for determining the height of the turbulent mixing layer on the basis of the vertical profiles of the dissipation rate of turbulent energy, which is estimated from lidar measurements of the radial wind velocity using conical scanning by a probe beam around the vertical axis. The accuracy of the proposed method is discussed in detail. It is shown that for the estimation of the mixing layer height (MLH) with the acceptable relative error not exceeding 20 %, the signal-to-noise ratio should be no less than −16 dB, when the relative error of lidar estimation of the dissipation rate does not exceed 30 %. The method was tested in a 6 d experiment in which the wind velocity turbulence was estimated in smog conditions due to forest fires in Siberia in summer 2019. The results of the experiment reveal that the relative error of determination of the MLH time series obtained by this method does not exceed 10 % in the period of turbulence development. The estimates of the turbulent mixing layer height by the proposed method are in a qualitative agreement with the MLH estimated from the distributions of the Richardson number in height and time obtained during the comparison experiment in spring 2020.

Description: In the 2019/2020 austral summer, the surface melt duration and extent on the northern George VI Ice Shelf (GVIIS) was exceptional compared to the 31 previous summers of distinctly lower melt. This finding is based on analysis of near-continuous 41-year satellite microwave radiometer and scatterometer data, which are sensitive to meltwater on the ice shelf surface and in the near-surface snow. Using optical satellite imagery from Landsat 8 (2013 to 2020) and Sentinel-2 (2017 to 2020), record volumes of surface meltwater ponding were also observed on the northern GVIIS in 2019/2020, with 23 % of the surface area covered by 0.62 km3 of ponded meltwater on 19 January. These exceptional melt and surface ponding conditions in 2019/2020 were driven by sustained air temperatures ≥0 ∘C for anomalously long periods (55 to 90 h) from late November onwards, which limited meltwater refreezing. The sustained warm periods were likely driven by warm, low-speed (≤7.5 m s−1) northwesterly and northeasterly winds and not by foehn wind conditions, which were only present for 9 h total in the 2019/2020 melt season. Increased surface ponding on ice shelves may threaten their stability through increased potential for hydrofracture initiation; a risk that may increase due to firn air content depletion in response to near-surface melting.

Description: In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ∼ 8–12 ∘C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Description: A ground-based network of more than 1200 Global Navigation Satellite System (GNSS) Continuously Operating Reference Stations (CORS) was analysed using GIPSY-OASIS II software package for the documentation of time and space variations of water vapor in atmosphere during the North Atlantic Waveguide and Downstream impact EXperiment (NAWDEX) during fall 2016. The network extends throughout the North Atlantic, from the Caribbeans to Morocco through Greenland. This paper presents the methodology used for GNSS data processing, screening, and conversion of Zenith Tropospheric Delay (ZTD) estimates to Integrated Water Vapor content (IWV) using surface parameters from reanalysis. The retrieved IWV are used to evaluate the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalyses ERAI and ERA5. ERA5 shows an overall improvement over ERAI in representing the spatial and temporal variability of IWV over the study area. The mean bias is decreased from 0.31±0.63 to 0.19±0.56 kg m−2 (mean ±1σ over all stations) and the standard deviation reduced from 2.17±0.67 to 1.64±0.53 kg m−2 combined with a slight improvement in correlation coefficient from 0.95 to 0.97. At regional scale, both reanalyses show a general wet bias at mid and northern latitudes but a dry bias in the Caribbeans. We hypothesize this results from the different nature of data being assimilated over the tropical oceans. This GNSS IWV data set is intended to be used for a better description of the high impact weather events that occurred during the NAWDEX experiment.

Description: The deep chlorophyll maximum (DCM) is a well-known feature of the global ocean. However, its description and the study of its formation are a challenge, especially in the peculiar environment that is the Black Sea. The retrieval of chlorophyll a (chl a) from fluorescence (Fluo) profiles recorded by Biogeochemical Argo (BGC-Argo) floats is not trivial in the Black Sea, due to the very high content of coloured dissolved organic matter (CDOM) which contributes to the fluorescence signal and produces an apparent increase in the chl a concentration with depth. Here, we revised Fluo correction protocols for the Black Sea context using co-located in situ high-performance liquid chromatography (HPLC) and BGC-Argo measurements. The processed set of chl a data (2014–2019) is then used to provide a systematic description of the seasonal DCM dynamics in the Black Sea and to explore different hypotheses concerning the mechanisms underlying its development. Our results show that the corrections applied to the chl a profiles are consistent with HPLC data. In the Black Sea, the DCM begins to form in March, throughout the basin, at a density level set by the previous winter mixed layer. During a first phase (April–May), the DCM remains attached to this particular layer. The spatial homogeneity of this feature suggests a hysteresis mechanism, i.e. that the DCM structure locally influences environmental conditions rather than adapting instantaneously to external factors. In a second phase (July–September), the DCM migrates upward, where there is higher irradiance, which suggests the interplay of biotic factors. Overall, the DCM concentrates around 45 % to 65 % of the total chlorophyll content within a 10 m layer centred around a depth of 30 to 40 m, which stresses the importance of considering DCM dynamics when evaluating phytoplankton productivity at basin scale.

Description: A normal contact stiffness model considering 3D topography and elastic–plastic contact of rough surfaces is presented in this paper. The asperities are generated from the measured surfaces using the watershed segmentation and a modified nine-point rectangle. The topography parameters, including the asperity locations, heights, and radii of the summit, are obtained. Asperity shoulder–shoulder contact is considered. The relationship of the contact parameters, such as the contact force, the deformation, and the mean separation of two surfaces, is modelled in the three different contact regimes, namely elastic, elastic–plastic and fully plastic. The asperity contact state is determined, and if the contact occurs, the stiffness of the single asperity pair is calculated and summed as the total normal stiffness of two contact surfaces. The developed model is validated using experimental tests conducted on two types of specimens and is compared with published theoretical models.

Description: In four study areas within different lithological settings and rockfall activity, lidar data were applied for a morphometric analysis of block sizes, block shapes and talus cone characteristics. This information was used to investigate the dependencies between block size, block shape and lithology on the one hand and runout distances on the other hand. In our study, we were able to show that lithology seems to have an influence on block size and shape and that gravitational sorting did not occur on all of the studied debris cones but that other parameters apparently control the runout length of boulders. Such a parameter seems to be the block shape, as it plays the role of a moderating parameter in two of the four study sites, while we could not confirm this for our other study sites. We also investigated the influence of terrain parameters such as slope inclination, profile curvature and roughness. The derived roughness values show a clear difference between the four study sites and seem to be a good proxy for block size distribution on the talus cones and thus could be used in further studies to analyse a larger sample of block size distribution on talus cones with different lithologies.

Description: Understanding the vertical distribution of aerosol helps to reduce the uncertainty in the aerosol life cycle and therefore in the estimation of the direct and indirect aerosol forcing. To improve our understanding, we use measurements from four deployments of the Atmospheric Tomography (ATom) field campaign (ATom1–4) which systematically sampled aerosol and trace gases over the Pacific and Atlantic oceans with near pole-to-pole coverage. We evaluate the UK Earth System Model (UKESM1) against ATom observations in terms of joint biases in the vertical profile of three variables related to new particle formation: total particle number concentration (NTotal), sulfur dioxide (SO2) mixing ratio and the condensation sink. The NTotal, SO2 and condensation sink are interdependent quantities and have a controlling influence on the vertical profile of each other; therefore, analysing them simultaneously helps to avoid getting the right answer for the wrong reasons. The simulated condensation sink in the baseline model is within a factor of 2 of observations, but the NTotal and SO2 show much larger biases mainly in the tropics and high latitudes. We performed a series of model sensitivity tests to identify atmospheric processes that have the strongest influence on overall model performance. The perturbations take the form of global scaling factors or improvements to the representation of atmospheric processes in the model, for example by adding a new boundary layer nucleation scheme. In the boundary layer (below 1 km altitude) and lower troposphere (1–4 km), inclusion of a boundary layer nucleation scheme (Metzger et al., 2010) is critical to obtaining better agreement with observations. However, in the mid (4–8 km) and upper troposphere (〉 8 km), sub-3 nm particle growth, pH of cloud droplets, dimethyl sulfide (DMS) emissions, upper-tropospheric nucleation rate, SO2 gas-scavenging rate and cloud erosion rate play a more dominant role. We find that perturbations to boundary layer nucleation, sub-3 nm growth, cloud droplet pH and DMS emissions reduce the boundary layer and upper tropospheric model bias simultaneously. In a combined simulation with all four perturbations, the SO2 and condensation sink profiles are in much better agreement with observations, but the NTotal profile still shows large deviations, which suggests a possible structural issue with how nucleation or gas/particle transport or aerosol scavenging is handled in the model. These perturbations are well-motivated in that they improve the physical basis of the model and are suitable for implementation in future versions of UKESM.

Description: Stratospheric ozone and water vapour are key components of the Earth system, and past and future changes to both have important impacts on global and regional climate. Here, we evaluate long-term changes in these species from the pre-industrial period (1850) to the end of the 21st century in Coupled Model Intercomparison Project phase 6 (CMIP6) models under a range of future emissions scenarios. There is good agreement between the CMIP multi-model mean and observations for total column ozone (TCO), although there is substantial variation between the individual CMIP6 models. For the CMIP6 multi-model mean, global mean TCO has increased from ∼ 300 DU in 1850 to ∼ 305 DU in 1960, before rapidly declining in the 1970s and 1980s following the use and emission of halogenated ozone-depleting substances (ODSs). TCO is projected to return to 1960s values by the middle of the 21st century under the SSP2-4.5, SSP3-7.0, SSP4-3.4, SSP4-6.0, and SSP5-8.5 scenarios, and under the SSP3-7.0 and SSP5-8.5 scenarios TCO values are projected to be ∼ 10 DU higher than the 1960s values by 2100. However, under the SSP1-1.9 and SSP1-1.6 scenarios, TCO is not projected to return to the 1960s values despite reductions in halogenated ODSs due to decreases in tropospheric ozone mixing ratios. This global pattern is similar to regional patterns, except in the tropics where TCO under most scenarios is not projected to return to 1960s values, either through reductions in tropospheric ozone under SSP1-1.9 and SSP1-2.6, or through reductions in lower stratospheric ozone resulting from an acceleration of the Brewer–Dobson circulation under other Shared Socioeconomic Pathways (SSPs). In contrast to TCO, there is poorer agreement between the CMIP6 multi-model mean and observed lower stratospheric water vapour mixing ratios, with the CMIP6 multi-model mean underestimating observed water vapour mixing ratios by ∼ 0.5 ppmv at 70 hPa. CMIP6 multi-model mean stratospheric water vapour mixing ratios in the tropical lower stratosphere have increased by ∼ 0.5 ppmv from the pre-industrial to the present-day period and are projected to increase further by the end of the 21st century. The largest increases (∼ 2 ppmv) are simulated under the future scenarios with the highest assumed forcing pathway (e.g. SSP5-8.5). Tropical lower stratospheric water vapour, and to a lesser extent TCO, shows large variations following explosive volcanic eruptions.

Description: High-quality and long-term soil moisture products are significant for hydrologic monitoring and agricultural management. However, the acquired daily Advanced Microwave Scanning Radiometer 2 (AMSR2) soil moisture products are incomplete in global land (just about 30 %–80 % coverage ratio), due to the satellite orbit coverage and the limitations of soil moisture retrieval algorithms. To solve this inevitable problem, we develop a novel spatio-temporal partial convolutional neural network (CNN) for AMSR2 soil moisture product gap-filling. Through the proposed framework, we generate the seamless daily global (SGD) AMSR2 long-term soil moisture products from 2013 to 2019. To further validate the effectiveness of these products, three verification methods are used as follows: (1) in situ validation, (2) time-series validation, and (3) simulated missing-region validation. Results show that the seamless global daily soil moisture products have reliable cooperativity with the selected in situ values. The evaluation indexes of the reconstructed (original) dataset are a correlation coefficient (R) of 0.685 (0.689), root-mean-squared error (RMSE) of 0.097 (0.093), and mean absolute error (MAE) of 0.079 (0.077). The temporal consistency of the reconstructed daily soil moisture products is ensured with the original time-series distribution of valid values. The spatial continuity of the reconstructed regions is in accordance with the spatial information (R: 0.963–0.974, RMSE: 0.065–0.073, and MAE: 0.044–0.052). This dataset can be downloaded at https://doi.org/10.5281/zenodo.4417458 (Zhang et al., 2021).

Description: The intertropical convergence zone (ITCZ) is an important component of the tropical rain belt. Climate models continue to struggle to adequately represent the ITCZ and differ substantially in its simulated response to climate change. Here we employ complex network approaches, which extract spatiotemporal variability patterns from climate data, to better understand differences in the dynamics of the ITCZ in state-of-the-art global circulation models (GCMs). For this purpose, we study simulations with 14 GCMs in an idealized slab-ocean aquaplanet setup from TRACMIP – the Tropical Rain belts with an Annual cycle and a Continent Model Intercomparison Project. We construct network representations based on the spatial correlation patterns of monthly surface temperature anomalies and study the zonal-mean patterns of different topological and spatial network characteristics. Specifically, we cluster the GCMs by means of the distributions of their zonal network measures utilizing hierarchical clustering. We find that in the control simulation, the distributions of the zonal network measures are able to pick up model differences in the tropical sea surface temperature (SST) contrast, the ITCZ position, and the strength of the Southern Hemisphere Hadley cell. Although we do not find evidence for consistent modifications in the network structure tracing the response of the ITCZ to global warming in the considered model ensemble, our analysis demonstrates that coherent variations of the global SST field are linked to ITCZ dynamics. This suggests that climate networks can provide a new perspective on ITCZ dynamics and model differences therein.

Description: Recent seismicity in Alberta and north-east British Columbia has been attributed to ongoing oil and gas development in the area, due to its temporal and spatial correlation. Prior to such development, the area was seismically quiescent. Here, we show evidence that latent seismicity may occur in areas where previous operations have occurred, even during a shutdown in operations. The global COVID-19 pandemic furnished the unique opportunity to study seismicity during a long period of anthropogenic quiescence. Within the Kiskatinaw area of British Columbia, 389 events were detected from April to August 2020, which encompasses a period with very little hydraulic fracturing operations. This reduction in operations was the result of a government-imposed lockdown severely restricting the movement of people as well as a downturn in the economic market causing industry stock prices to collapse. Except for a reduction in the seismicity rate and a lack of temporal clustering that is often characteristic of hydraulic fracturing induced sequences, the general characteristics of the observed seismicity were similar to the preceding time period of active operations. During the period of relative quiescence, event magnitudes were observed between ML −0.7 and ML 1.2, which is consistent with previous event magnitudes in the area. Hypocentres occurred in a corridor orientated NW–SE, just as seismicity had done in previous years, and were located at depths associated with the target Montney formation or shallower (

Description: Satellite observations (International Satellite Cloud Climatology Project (ISCCP), 1983–2009) of linear trends in cloud cover are compared to those in global precipitation (Global Precipitation Climatology Project (GPCP) pentad V2.2, 1983–2009), to investigate possible cause(s) of the linear trends in both cloud cover and precipitation. The spatial distributions of the linear trends in total cloud cover and precipitation are both characterized primarily by a broadening of the major ascending zone of Hadley circulation. Our correlation studies suggest that global warming, Atlantic Multidecadal Oscillation (AMO) and Pacific Decadal Oscillation (PDO) can explain 67 %, 49 % and 38 %, respectively, of the spatial variabilities in the linear trends in cloud cover, but causality is harder to establish. Further analysis of the broadening of the major ascending zone of Hadley circulation shows that the trend in global temperature, rather than those in AMO and PDO, is the primary contributor to the observed linear trends in total cloud cover and precipitation in 1983–2009. The underlying mechanism driving this broadening is proposed to be the moisture–convection–latent-heat feedback cycle under global warming conditions. The global analysis is extended by investigating connections between clouds and precipitation in China, based on a large number of long-running, high-quality surface weather stations in 1957–2005. This reveals a quantitative matching relationship between the reduction in light precipitation and the reduction in total cloud cover. Furthermore, our study suggests that the reduction in cloud cover in China is primarily driven by global temperature; PDO plays a secondary role, while the contribution from AMO and Niño3.4 is insignificant, consistent with the global analysis.

Description: This study was performed in the framework of HyMeX (Hydrological cycle in the Mediterranean Experiment), which aimed to study the heavy precipitation that regularly affects the Mediterranean area. A reanalysis with a convective-scale model AROME-WMED (Application of Research to Operations at MEsoscale western Mediterranean) was performed, which assimilated most of the available data for a 2-month period corresponding to the first special observation period of the field campaign (Fourrié et al., 2019). Among them, observations related to the low-level humidity flow were assimilated. Such observations are important for the description of the feeding of the convective mesoscale systems with humidity (Duffourg and Ducrocq, 2011; Bresson et al., 2012; Ricard et al., 2012). Among them there were a dense reprocessed network of high-quality Global Navigation Satellite System (GNSS) zenithal total delay (ZTD) observations, reprocessed data from wind profilers, lidar-derived vertical profiles of humidity (ground and airborne) and Spanish radar data. The aim of the paper is to assess the impact of the assimilation of these four observation types on the analyses and the forecasts from the 3 h forecast range (first guess) up to the 48 h forecast range. In order to assess this impact, several observing system experiments (OSEs) or so-called denial experiments, were carried out by removing one single data set from the observation data set assimilated in the reanalysis. Among the evaluated observations, it is found that the ground-based GNSS ZTD data set provides the largest impact on the analyses and the forecasts, as it represents an evenly spread and frequent data set providing information at each analysis time over the AROME-WMED domain. The impact of the reprocessing of GNSS ZTD data also improves the forecast quality, but this impact is not statistically significant. The assimilation of the Spanish radar data improves the 3 h precipitation forecast quality as well as the short-term (30 h) precipitation forecasts, but this impact remains located over Spain. Moreover, marginal impact from wind profilers was observed on wind background quality. No impacts have been found regarding lidar data, as they represent a very small data set, mainly located over the sea.

Description: Systematic measurements of the atmospheric Ar∕N2 ratio have been made at ground-based stations in Japan and Antarctica since 2012. Clear seasonal cycles of the Ar∕N2 ratio with summertime maxima were found at middle- to high-latitude stations, with seasonal amplitudes increasing with increasing latitude. Eight years of the observed Ar∕N2 ratio at Tsukuba (TKB) and Hateruma (HAT), Japan, showed interannual variations in phase with the observed variations in the global ocean heat content (OHC). We calculated secularly increasing trends of 0.75 ± 0.30 and 0.89 ± 0.60 per meg per year from the Ar∕N2 ratio observed at TKB and HAT, respectively, although these trend values are influenced by large interannual variations. In order to examine the possibility of the secular trend in the surface Ar∕N2 ratio being modified significantly by the gravitational separation in the stratosphere, two-dimensional model simulations were carried out by arbitrarily modifying the mass stream function in the model to simulate either a weakening or an enhancement of the Brewer–Dobson circulation (BDC). The secular trend of the Ar∕N2 ratio at TKB, corrected for gravitational separation under the assumption of weakening (enhancement) of BDC simulated by the 2-D model, was 0.60 ± 0.30 (0.88 ± 0.30) per meg per year. By using a conversion factor of 3.5 × 10−23 per meg per joule by assuming a one-box ocean with a temperature of 3.5 ∘C, average OHC increase rates of 17.1 ± 8.6 ZJ yr−1 and 25.1 ± 8.6 ZJ yr−1 for the period 2012–2019 were estimated from the corrected secular trends of the Ar∕N2 ratio for the weakened- and enhanced-BDC conditions, respectively. Both OHC increase rates from the uncorrected- and weakened-BDC secular trends of the Ar∕N2 ratio are consistent with 12.2 ± 1.2 ZJ yr−1 reported by ocean temperature measurements, while that from the enhanced-BDC is outside of the range of the uncertainties. Although the effect of the actual atmospheric circulation on the Ar∕N2 ratio is still unclear and longer-term observations are needed to reduce uncertainty of the secular trend of the surface Ar∕N2 ratio, the analytical results obtained in the present study imply that the surface Ar∕N2 ratio is an important tracer for detecting spatiotemporally integrated changes in OHC and BDC.

Description: Understanding the evolution of, and the interactions between, ice sheets and the global climate over geological timescales is important for being able to project their future evolution. However, direct observational evidence of past CO2 concentrations, and the implied radiative forcing, only exists for the past 800 000 years. Records of benthic δ18O date back millions of years but contain signals from both land ice volume and ocean temperature. In recent years, inverse forward modelling has been developed as a method to disentangle these two signals, resulting in mutually consistent reconstructions of ice volume, temperature, and CO2. We use this approach to force a hybrid ice-sheet–climate model with a benthic δ18O stack, reconstructing the evolution of the ice sheets, global mean sea level, and atmospheric CO2 during the late Pliocene and the Pleistocene, from 3.6 million years (Myr) ago to the present day. During the warmer-than-present climates of the late Pliocene, reconstructed CO2 varies widely, from 320–440 ppmv for warm periods to 235–250 ppmv for the early glacial excursion ∼3.3 million years ago. Sea level is relatively stable during this period, with maxima of 6–14 m and minima of 12–26 m during glacial episodes. Both CO2 and sea level are within the wide ranges of values covered by available proxy data for this period. Our results for the Pleistocene agree well with the ice-core CO2 record, as well as with different available sea-level proxy data. For the Early Pleistocene, 2.6–1.2 Myr ago, we simulate 40 kyr glacial cycles, with interglacial CO2 decreasing from 280–300 ppmv at the beginning of the Pleistocene to 250–280 ppmv just before the Mid-Pleistocene Transition (MPT). Peak glacial CO2 decreases from 220–250 to 205–225 ppmv during this period. After the MPT, when the glacial cycles change from 40 to 80 120 kyr cyclicity, the glacial–interglacial contrast increases, with interglacial CO2 varying between 250–320 ppmv and peak glacial values decreasing to 170–210 ppmv.

Description: In urban areas, road traffic is a dominant source of nitrogen oxides (NOx=NO+NO2). Although the emissions from individual vehicles are regulated by the European emission standards, real driving emissions often exceed these limits. In this study, two multi-axis differential optical absorption spectroscopy (MAX-DOAS) instruments on opposite sides of the motorway were used to measure the NO2 absorption caused by road traffic at the A60 motorway close to Mainz, Germany. In combination with wind data, the total NOx emissions for the occurring traffic volume can be estimated. Hereto, the ozone-dependent photochemical equilibrium between NO and NO2 is considered. We show that for 10 May 2019 the measured emissions exceed the maximum expected emissions calculated from the European emission standards for standardised test cycles by a factor of 11±7. One major advantage of the method used here is that MAX-DOAS measurements are very sensitive to the integrated NO2 concentration close to the surface. Thus, all emitted NO2 molecules are detected independently from their altitude, and therefore the whole emission plume originating from the nearby motorway is captured, which is a key advantage compared to other approaches such as in situ measurements.

Description: To provide long-term air pollutant exposure estimates for epidemiological studies, it is essential to test the feasibility of developing land-use regression (LUR) models using only routine air quality measurement data and to evaluate the transferability of LUR models between nearby cities. In this study, we developed and evaluated the intercity transferability of annual-average LUR models for ambient respirable suspended particulates (PM10), fine suspended particulates (PM2.5), nitrogen dioxide (NO2) and ozone (O3) in the Taipei–Keelung metropolitan area of northern Taiwan in 2019. Ambient PM10, PM2.5, NO2 and O3 measurements at 30 fixed-site stations were used as the dependent variables, and a total of 156 potential predictor variables in six categories (i.e., population density, road network, land-use type, normalized difference vegetation index, meteorology and elevation) were extracted using buffer spatial analysis. The LUR models were developed using the supervised forward linear regression approach. The LUR models for ambient PM10, PM2.5, NO2 and O3 achieved relatively high prediction performance, with R2 values of 〉 0.72 and leave-one-out cross-validation (LOOCV) R2 values of 〉 0.53. The intercity transferability of LUR models varied among the air pollutants, with transfer-predictive R2 values of 〉 0.62 for NO2 and

Description: The species Selenogonus narinoensis was described by Stirton (1947) based on a single specimen which comes from sediments cropping out in the Cocha Verde locality, Nariño Department (Colombia), tentatively referred to the late Pliocene–Pleistocene (MGN 931; IGM p002118, Museo Geológico Nacional, Servicio Geológico Colombiano, Bogotá). However, morphological studies and comparative morphometric observations of the specimen suggest that (1) no diagnostic character supports the validity of the species Selenogonus narinoensis (here considered species inquirenda); (2) a combination of features (e.g., the mandibular condyle located behind the posterior edge of the vertical mandibular ramus, the angular process which projects laterally outwards, a bunolophodont crown morphology, a mesodont crown height, and a simple crown morphology of the third lobe of m3) indicates it belongs to the genus Platygonus; (3) this specimen corresponds to one of the largest South American peccaries; (4) taking into account certain anatomical characters as well as its morphometric range, this specimen is assigned to Platygonus cf. marplatensis. Even though the stratigraphic provenance of the specimen is still doubtful, it can be proposed that (1) it could be one of the most ancient records of tayassuids in South America, as would be expected given its geographical position, and (2) considering the new taxonomic proposal, this specimen represents the first record of Platygonus cf. marplatensis in Colombia and represents one of the northernmost South American records of the genus. This new interpretation would be of great relevance in the Great American Biotic Interchange due to its strategic geographical proximity to the Isthmus of Panama.

Description: We experimentally quantified the impact of cloud fraction and cloud type on the heating rate (HR) of black and brown carbon (HRBC and HRBrC). In particular, we examined in more detail the cloud effect on the HR detected in a previous study (Ferrero et al., 2018). High-time-resolution measurements of the aerosol absorption coefficient at multiple wavelengths were coupled with spectral measurements of the direct, diffuse and surface reflected irradiance and with lidar–ceilometer data during a field campaign in Milan, Po Valley (Italy). The experimental set-up allowed for a direct determination of the total HR (and its speciation: HRBC and HRBrC) in all-sky conditions (from clear-sky conditions to cloudy). The highest total HR values were found in the middle of winter (1.43 ± 0.05 K d−1), and the lowest were in spring (0.54 ± 0.02 K d−1). Overall, the HRBrC accounted for 13.7 ± 0.2 % of the total HR, with the BrC being characterized by an absorption Ångström exponent (AAE) of 3.49 ± 0.01. To investigate the role of clouds, sky conditions were classified in terms of cloudiness (fraction of the sky covered by clouds: oktas) and cloud type (stratus, St; cumulus, Cu; stratocumulus, Sc; altostratus, As; altocumulus, Ac; cirrus, Ci; and cirrocumulus–cirrostratus, Cc–Cs). During the campaign, clear-sky conditions were present 23 % of the time, with the remaining time (77 %) being characterized by cloudy conditions. The average cloudiness was 3.58 ± 0.04 oktas (highest in February at 4.56 ± 0.07 oktas and lowest in November at 2.91 ± 0.06 oktas). St clouds were mostly responsible for overcast conditions (7–8 oktas, frequency of 87 % and 96 %); Sc clouds dominated the intermediate cloudiness conditions (5–6 oktas, frequency of 47 % and 66 %); and the transition from Cc–Cs to Sc determined moderate cloudiness (3–4 oktas); finally, low cloudiness (1–2 oktas) was mostly dominated by Ci and Cu (frequency of 59 % and 40 %, respectively). HR measurements showed a constant decrease with increasing cloudiness of the atmosphere, enabling us to quantify for the first time the bias (in %) of the aerosol HR introduced by the simplified assumption of clear-sky conditions in radiative-transfer model calculations. Our results showed that the HR of light-absorbing aerosol was ∼ 20 %–30 % lower in low cloudiness (1–2 oktas) and up to 80 % lower in completely overcast conditions (i.e. 7–8 oktas) compared to clear-sky ones. This means that, in the simplified assumption of clear-sky conditions, the HR of light-absorbing aerosol can be largely overestimated (by 50 % in low cloudiness, 1–2 oktas, and up to 500 % in completely overcast conditions, 7–8 oktas). The impact of different cloud types on the HR was also investigated. Cirrus clouds were found to have a modest impact, decreasing the HRBC and HRBrC by −5 % at most. Cumulus clouds decreased the HRBC and HRBrC by −31 ± 12 % and −26 ± 7 %, respectively; cirrocumulus–cirrostratus clouds decreased the HRBC and HRBrC by −60 ± 8 % and −54 ± 4 %, which was comparable to the impact of altocumulus (−60 ± 6 % and −46 ± 4 %). A higher impact on the HRBC and HRBrC suppression was found for stratocumulus (−63 ± 6 % and −58 ± 4 %, respectively) and altostratus (−78 ± 5 % and −73 ± 4 %, respectively). The highest impact was associated with stratus, suppressing the HRBC and HRBrC by −85 ± 5 % and −83 ± 3 %, respectively. The presence of clouds caused a decrease of both the HRBC and HRBrC (normalized to the absorption coefficient of the respective species) of −11.8 ± 1.2 % and −12.6 ± 1.4 % per okta. This study highlights the need to take into account the role of both cloudiness and different cloud types when estimating the HR caused by both BC and BrC and in turn decrease the uncertainties associated with the quantification of their impact on the climate.

Description: Extensive and rigorous model intercomparison is of great importance before model application due to the uncertainties in current land surface models (LSMs). Without considering the uncertainties in forcing data and model parameters, this study designed an ensemble of 55 296 experiments to evaluate the Noah LSM with multi-parameterization (Noah-MP) for snow cover events (SCEs), soil temperature (ST) and soil liquid water (SLW) simulation, and investigated the sensitivity of parameterization schemes at a typical permafrost site on the Qinghai–Tibet Plateau (QTP). The results showed that Noah-MP systematically overestimates snow cover, which could be greatly resolved when adopting the sublimation from wind and a semi-implicit snow/soil temperature time scheme. As a result of the overestimated snow, Noah-MP generally underestimates ST, which is mostly influenced by the snow process. A systematic cold bias and large uncertainties in soil temperature remain after eliminating the effects of snow, particularly in the deep layers and during the cold season. The combination of roughness length for heat and under-canopy (below-canopy) aerodynamic resistance contributes to resolving the cold bias in soil temperature. In addition, Noah-MP generally underestimates top SLW. The runoff and groundwater (RUN) process dominates the SLW simulation in comparison to the very limited impacts of all other physical processes. The analysis of the model structural uncertainties and characteristics of each scheme would be constructive to a better understanding of the land surface processes in the permafrost regions of the QTP as well as to further model improvements towards soil hydrothermal regime modeling using LSMs.

Description: Organic nitrates are secondary species in the atmosphere. Their fate is related to the chemical transport of pollutants from polluted areas to more distant zones. While their gas-phase chemistry has been studied, their reactivity in condensed phases is far from being understood. However, these compounds represent an important fraction of organic matter in condensed phases. In particular, their partition to the aqueous phase may be especially important for oxidized organic nitrates for which water solubility increases with functionalization. This work has studied for the first time the aqueous-phase ⚫OH-oxidation kinetics of four alkyl nitrates (isopropyl nitrate, isobutyl nitrate, 1-pentyl nitrate, and isopentyl nitrate) and three functionalized organic nitrates (α-nitrooxyacetone, 1-nitrooxy-2-propanol, and isosorbide 5-mononitrate) by developing a novel and accurate competition kinetic method. Low reactivity was observed, with kOH ranging from 8×107 to 3.1×109 L mol−1 s−1 at 296±2 K. Using these results, a previously developed aqueous-phase structure–activity relationship (SAR) was extended, and the resulting parameters confirmed the extreme deactivating effect of the nitrate group, up to two adjacent carbon atoms. The achieved extended SAR was then used to determine the ⚫OH-oxidation rate constants of 49 organic nitrates, including hydroxy nitrates, ketonitrates, aldehyde nitrates, nitrooxy carboxylic acids, and more functionalized organic nitrates such as isoprene and terpene nitrates. Their multiphase atmospheric lifetimes towards ⚫OH oxidation were calculated using these rate constants, and they were compared to their gas-phase lifetimes. Large differences were observed, especially for polyfunctional organic nitrates: for 50 % of the proposed organic nitrates for which the ⚫OH reaction occurs mainly in the aqueous phase (more than 50 % of the overall removal), their ⚫OH-oxidation lifetimes increased by 20 % to 155 % under cloud/fog conditions (liquid water content LWC = 0.35 g m−3). In particular, for 83 % of the proposed terpene nitrates, the reactivity towards ⚫OH occurred mostly (〉98 %) in the aqueous phase, while for 60 % of these terpene nitrates, their lifetimes increased by 25 % to 140 % compared to their gas-phase reactivity. We demonstrate that these effects are of importance under cloud/fog conditions but also under wet aerosol conditions, especially for the terpene nitrates. These results suggest that considering aqueous-phase ⚫OH-oxidation reactivity of biogenic nitrates is necessary to improve the predictions of their atmospheric fate.

Description: Accurate knowledge of cloud properties is essential to the measurement of atmospheric composition from space. In this work we assess the quality of the cloud data from three Copernicus Sentinel-5 Precursor (S5P) TROPOMI cloud products: (i) S5P OCRA/ROCINN_CAL (Optical Cloud Recognition Algorithm/Retrieval of Cloud Information using Neural Networks;Clouds-As-Layers), (ii) S5P OCRA/ROCINN_CRB (Clouds-as-Reflecting Boundaries), and (iii) S5P FRESCO-S (Fast Retrieval Scheme for Clouds from Oxygen absorption bands – Sentinel). Target properties of this work are cloud-top height and cloud optical thickness (OCRA/ROCINN_CAL), cloud height (OCRA/ROCINN_CRB and FRESCO-S), and radiometric cloud fraction (all three algorithms). The analysis combines (i) the examination of cloud maps for artificial geographical patterns, (ii) the comparison to other satellite cloud data (MODIS, NPP-VIIRS, and OMI O2–O2), and (iii) ground-based validation with respect to correlative observations (30 April 2018 to 27 February 2020) from the Cloudnet network of ceilometers, lidars, and radars. Zonal mean latitudinal variation of S5P cloud properties is similar to that of other satellite data. S5P OCRA/ROCINN_CAL agrees well with NPP VIIRS cloud-top height and cloud optical thickness and with Cloudnet cloud-top height, especially for the low (mostly liquid) clouds. For the high clouds, S5P OCRA/ROCINN_CAL cloud-top height is below the cloud-top height of VIIRS and of Cloudnet, while its cloud optical thickness is higher than that of VIIRS. S5P OCRA/ROCINN_CRB and S5P FRESCO cloud height are well below the Cloudnet cloud mean height for the low clouds but match on average better with the Cloudnet cloud mean height for the higher clouds. As opposed to S5P OCRA/ROCINN_CRB and S5P FRESCO, S5P OCRA/ROCINN_CAL is well able to match the lowest CTH mode of the Cloudnet observations. Peculiar geographical patterns are identified in the cloud products and will be mitigated in future releases of the cloud data products.

Description: Application of the environmentally friendly scaling inhibitor NC47.1 B in geothermal systems was studied in laboratory and field-scale experiments. Biodegradation was investigated under anaerobic, in situ-like conditions and a mass balance confirmed the almost complete conversion of the polycarboxylate to e.g. acetate, formate, methane and CO2. Much higher concentrations of inhibitor were chosen than applied in situ and rapid degradation was observed in biofilm-inoculated setups: A concentration of 100 mg/L of the inhibitor was degraded below detection limit within 8 d of incubation. Furthermore, the inhibitor was applied at the geothermal plant in Unterhaching, Germany. Monitoring of the microbial community in situ showed an increase in the abundance of Bacteria. Particularly, relatives of the fermenting Caldicellulosiruptor dominated the biocenosis after about six months of continuous inhibitor dosage (5–10 mg/L). However, in long-term laboratory experiments representatives of Caldicellulosiruptor were only detected in traces and the microbial community comprised a broader spectrum of fermentative bacteria. The different composition of the biocenosis in situ and in laboratory experiments is probably caused by the different inhibitor concentrations, temperatures as well as nutrient availability in situ compared to the closed system of the batch experiments.

Description: This work describes the latest design, calibration and application of a near-infrared laser diode-based photoacoustic (PA) hygrometer developed for total water content measurement in simulated atmospheric freezing precipitation and high ice water content conditions with relevance in fundamental icing research, aviation testing, and certification. The single-wavelength and single-pass PA absorption cell is calibrated for molar water vapor fractions with a two-pressure humidity generator integrated into the instrument. Laboratory calibration showed an estimated measurement accuracy better than 3.3 % in the water vapor mole fraction range of 510–12 360 ppm (5 % from 250–21 200 ppm) with a theoretical limit of detection (3σ) of 3.2 ppm. The hygrometer is examined in combination with a basic isokinetic evaporator probe (IKP) and sampling system designed for icing wind tunnel applications, for which a general description of total condensed water content (CWC) measurements and uncertainties are presented. Despite the current limitation of the IKP to a hydrometeor mass flux below 90 gm-2s-1, a CWC measurement accuracy better than 20 % is achieved by the instrument above a CWC of 0.14 g m−3 in cold air (−30 ∘C) with suitable background humidity measurement. Results of a comparison to the Cranfield University IKP instrument in freezing drizzle and rain show a CWC agreement of the two instruments within 20 %, which demonstrates the potential of PA hygrometers for water content measurement in atmospheric icing conditions.

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.