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: Mantle-derived noble gases in volcanic gases are powerful tracers of terrestrial volatile evolution, as they contain mixtures of both primordial (from Earth's accretion) and secondary (e.g., radiogenic) isotope signals that characterize the composition of deep Earth. However, volcanic gases emitted through subaerial hydrothermal systems also contain contributions from shallow reservoirs (groundwater, crust, atmosphere). Deconvolving deep and shallow source signals is critical for robust interpretations of mantle-derived signals. Here, we use a novel dynamic mass spectrometry technique to measure argon, krypton, and xenon isotopes in volcanic gas with ultrahigh precision. Data from Iceland, Germany, United States (Yellowstone, Salton Sea), Costa Rica, and Chile show that subsurface isotope fractionation within hydrothermal systems is a globally pervasive and previously unrecognized process causing substantial nonradiogenic Ar-Kr-Xe isotope variations. Quantitatively accounting for this process is vital for accurately interpreting mantle-derived volatile (e.g., noble gas and nitrogen) signals, with profound implications for our understanding of terrestrial volatile evolution.

Description: Published

Description: eadg2566

Description: OSV2: Complessità dei processi vulcanici: approcci multidisciplinari e multiparametrici

Description: JCR Journal

Description: TR17-08, a marine sedimentary core (14.6 m), was collected during 2017 from the Edisto Inlet (Ross Sea, Antarctica), a small fjord near Cape Hallett. The core is characterized by expanded laminated sedimentary sequences making it suitable for studying submillennial processes during the Early Holocene. By studying different well-known foraminifera species (Globocassidulina biora, G. subglobosa, Trifarina angulosa, Nonionella iridea, Epistominella exigua, Stainforthia feylingi, Miliammina arenacea, Paratrochammina bartrami and Portatrochammina antarctica), we were able to identify five different foraminiferal assemblages over the last ∼ 2000 years BP. Comparison with diatom assemblages and other geochemical proxies retrieved from nearby sediment cores in the Edisto Inlet (BAY05-20 and HLF17-1) made it possible to distinguish three different phases characterized by different environmental settings: (1) a seasonal phase (from 2012 to 1486 years BP) characterized by the dominance of calcareous species, indicating a seasonal opening of the inlet by more frequent events of melting of the sea-ice cover during the austral summer and, in general, a higher-productivity, more open and energetic environment; (2) a transitional phase (from 1486 to 696 years BP) during which the fjord experienced less extensive sea-ice melting, enhanced oxygen-poor conditions and carbonate dissolution conditions, indicated by the shifts from calcareous-dominated association to agglutinated-dominated association probably due to a freshwater input from the retreat of three local glaciers at the start of this period; and (3) a cooler phase (from 696 years BP to present) during which the sedimentation rate decreased and few to no foraminiferal specimens were present, indicating ephemeral openings or a more prolonged cover of the sea ice during the austral summer, affecting the nutrient supply and the sedimentation regime.

Description: Published

Description: 95–115

Description: OSA2: Evoluzione climatica: effetti e loro mitigazione

Description: JCR Journal

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: 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: 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: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Liu, C.-Z., Dick, H. J. B., Mitchell, R. N., Wei, W., Zhang, Z.-Y., Hofmann, A. W., Yang, J.-F., & Li, Y. Archean cratonic mantle recycled at a mid-ocean ridge. Science Advances, 8(22), (2022): eabn6749, https://doi.org/10.1126/sciadv.abn6749.

Description: Basalts and mantle peridotites of mid-ocean ridges are thought to sample Earth’s upper mantle. Osmium isotopes of abyssal peridotites uniquely preserve melt extraction events throughout Earth history, but existing records only indicate ages up to ~2 billion years (Ga) ago. Thus, the memory of the suspected large volumes of mantle lithosphere that existed in Archean time (〉2.5 Ga) has apparently been lost somehow. We report abyssal peridotites with melt-depletion ages up to 2.8 Ga, documented by extremely unradiogenic 187Os/188Os ratios (to as low as 0.1095) and refractory major elements that compositionally resemble the deep keels of Archean cratons. These oceanic rocks were thus derived from the once-extensive Archean continental keels that have been dislodged and recycled back into the mantle, the feasibility of which we confirm with numerical modeling. This unexpected connection between young oceanic and ancient continental lithosphere indicates an underappreciated degree of compositional recycling over time.

Description: This study was financially supported by the National Science Fund for Distinguished Young Scholars 42025201 (to C.-Z.L.), the National Key Research and Development Project of China 2020YFA0714801 (to C.-Z.L.), the Strategic Priority Research Program of the Chinese Academy of Sciences XDA13010106 (to C.-Z.L.), the Strategic Priority Research Program of the Chinese Academy of Sciences XDB42020301 (to C.-Z.L.), and NSF grants 2114652 and 1657983 (to H.J.B.D.).

Description: © The Author(s), 2021. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Zhang, Y., Gazel, E., Gaetani, G. A., & Klein, F. Serpentinite-derived slab fluids control the oxidation state of the subarc mantle. Science Advances, 7(48), (2021): eabj2515, https://doi.org/10.1126/sciadv.abj2515.

Description: Recent geochemical evidence confirms the oxidized nature of arc magmas, but the underlying processes that regulate the redox state of the subarc mantle remain yet to be determined. We established a link between deep subduction-related fluids derived from dehydration of serpentinite ± altered oceanic crust (AOC) using B isotopes and B/Nb as fluid proxies, and the oxidized nature of arc magmas as indicated by Cu enrichment during magma evolution and V/Yb. Our results suggest that arc magmas derived from source regions influenced by a greater serpentinite (±AOC) fluid component record higher oxygen fugacity. The incorporation of this component into the subarc mantle is controlled by the subduction system’s thermodynamic conditions and geometry. Our results suggest that the redox state of the subarc mantle is not homogeneous globally: Primitive arc magmas associated with flat, warm subduction are less oxidized overall than those generated in steep, cold subduction zones.

Description: Y.Z. acknowledges funding from the National Science Foundation of China (91958213), the Chinese Academy of Sciences (XDB42020402), and the Shandong Provincial Natural Science Foundation, China (ZR2020QD068). This study was supported in part by the U.S. National Science Foundation NSF EAR 1826673 to E.G. and G.A.G. and OCE 1756349 to E.G.

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: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Freeman, D. H., & Ward, C. P. Sunlight-driven dissolution is a major fate of oil at sea. Science Advances, 8(7), (2022): eabl7605, https://doi.org/10.1126/sciadv.abl7605.

Description: Oxygenation reactions initiated by sunlight can transform insoluble components of crude oil at sea into water-soluble products, a process called photo-dissolution. First reported a half century ago, photo-dissolution has never been included in spill models because key parameters required for rate modeling were unknown, including the wavelength and photon dose dependence. Here, we experimentally quantified photo-dissolution as a function of wavelength and photon dose, making possible a sensitivity analysis of environmental variables in hypothetical spill scenarios and a mass balance assessment for the 2010 Deepwater Horizon (DwH) spill. The sensitivity analysis revealed that rates were most sensitive to oil slick thickness, season/latitude, and wavelength and less sensitive to photon dose. We estimate that 3 to 17% (best estimate 8%) of DwH surface oil was subject to photo-dissolution, comparable in magnitude to other widely recognized fate processes. Our findings invite a critical reevaluation of surface oil budgets for both DwH and future spills at sea.

Description: This work was supported by the Fisheries and Oceans Canada Multi-Partner Research Initiative award to C.P.W. (project #1.06), the NSF Graduate Research Fellowship awarded to D.H.F. (award #174530), and NSF-OCE grant #1841092 to C.P.W.

Description: © The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Peng, Q., Xie, S.-P., Wang, D., Huang, R. X., Chen, G., Shu, Y., Shi, J.-R., & Liu, W. Surface warming-induced global acceleration of upper ocean currents. Science Advances, 8(16), (2022): eabj8394, https://doi.org/10.1126/sciadv.abj8394.

Description: How the ocean circulation changes in a warming climate is an important but poorly understood problem. Using a global ocean model, we decompose the problem into distinct responses to changes in sea surface temperature, salinity, and wind. Our results show that the surface warming effect, a robust feature of anthropogenic climate change, dominates and accelerates the upper ocean currents in 77% of the global ocean. Specifically, the increased vertical stratification intensifies the upper subtropical gyres and equatorial currents by shoaling these systems, while the differential warming between the Southern Ocean upwelling zone and the region to the north accelerates surface zonal currents in the Southern Ocean. In comparison, the wind stress and surface salinity changes affect regional current systems. Our study points a way forward for investigating ocean circulation change and evaluating the uncertainty.

Description: Q.P. is supported by the National Natural Science Foundation of China (42005035), the Science and Technology Planning Project of Guangzhou (202102020935), and the Independent Research Project Program of State Key Laboratory of Tropical Oceanography (LTOZZ2102). D.W. is supported by the National Natural Science Foundation of China (92158204), and the Innovation Group Project of Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (311020004). S.-P.X. is supported by the National Science Foundation (AGS-1934392). Y.S. is supported by the National Key Research and Development Program of China (2016YFC1401702). G.C. is supported by National Natural Science Foundation of China (41822602). The numerical simulation is supported by the High-Performance Computing Division and HPC managers of W. Zhou and D. Sui in the South China Sea Institute of Oceanology.

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.