ALBERT

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: Previous studies have concluded that the wind-input vorticity in ocean gyres is balanced by bottom pressure torques (BPT), when integrated over latitude bands. However, the BPT must vanish when integrated over any area enclosed by an isobath. This constraint raises ambiguities regarding the regions over which BPT should close the vorticity budget, and implies that BPT generated to balance a local wind stress curl necessitates the generation of a compensating, non-local BPT and thus non-local circulation. This study aims to clarify the role of BPT in wind-driven gyres using an idealized isopycnal model. Experiments performed with a single-signed wind stress curl in an enclosed, sloped basin reveal that BPT balances the winds only when integrated over latitude bands. Integrating over other, dynamically-motivated definitions of the gyre, such as barotropic streamlines, yields a balance between wind stress curl and bottom frictional torques. This implies that bottom friction plays a non-negligible role in structuring the gyre circulation. Non-local bottom pressure torques manifest in the form of along-slope pressure gradients associated with a weak basin-scale circulation, and are associated with a transition to a balance between wind stress and bottom friction around the coasts. Finally, a suite of perturbation experiments is used to investigate the dynamics of BPT. To predict the BPT, the authors extend previous theory that describes propagation of surface pressure signals from the gyre interior toward the coast along planetary potential vorticity contours. This theory is shown to agree closely with the diagnosed contributions to the vorticity budget across the suite of model experiments.

Description: Sufficient and accurate tide data are essential for analyzing physical processes in the ocean. A method is developed to spatially fit the tidal amplitude and phase lag data along satellite altimeter tracks near Hawaii and construct reliable cotidal charts by using the Chebyshev polynomials. The method is completely dependent on satellite altimeter data. By using the cross-validation method, the optimal orders of Chebyshev polynomials are determined and the polynomial coefficients are calculated by the least squares method. The tidal amplitudes and phase lags obtained by the method are compared with those from the Finite Element Solutions 2014 (FES2014), National Astronomical Observatory 99b (NAO.99b) and TPXO9 models. Results indicate that the method yields accurate results as its fitting results are consistent with the harmonic constants of the three models. The feasibility of this method is also validated by the harmonic constants from tidal gauges near Hawaii.

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: 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: The temporal-spatial variations of the static stability of dry air and the relative importance of their influencing quantities are explored. Derivation shows that while it links to the vertical difference of temperature, static stability also relates to the temperature itself. The static stability is expressed as a nonlinear function of temperature and the vertical difference of temperature. The relative importance of the two influencing quantities is assessed with the linear regression. Tests show that the linear fitting method is robust. The results of the dominance rely on the data examined, which include an interannual variation, a seasonal variation, and a spatial variation that consists of the grid points over the globe. It is revealed that in lower troposphere, while the temporal variations of static stability are dominated by the vertical difference of temperature, the temperature itself may also have considerable influence, especially over the high latitudes of the two hemispheres. In stratosphere, temperature tends to have more contributions. Over Antarctic, temperature dominates the seasonal and interannual variations of the static stability. The spatial variation of the static stability of July is influenced by both temperature and its vertical difference before 1980, but after that it is dominated by temperature.

Description: In recent winters, there have been repeated observations of extreme warm and cold spells in the mid-latitude countries. This has evoked questions regarding how winter temperature extremes are induced. In this study, we demonstrate that abnormally warm winter weather in East Asia can drive the onset of extremely cold weather in North America approximately one week forward. These seesawing extremes across the basin are mediated by the North Pacific Oscillation (NPO), one of the recurrent atmospheric patterns over the North Pacific. Budget analysis of the quasi-geostrophic geopotential tendency equation shows that intense thermal advection over East Asia is able to trigger the growth of the NPO. Vorticity fluxes associated with the upper-level stationary trough then strengthen and maintain the NPO against thermal damping following the onset of the NPO. Differential diabatic heating accompanied by changes in circulation also positively contribute to the growth and maintenance of the NPO. These results imply that recurrent cold extremes, seemingly contrary to global warming, may be an inherent feature resulting from strengthening warm extremes.

Description: This study revisits the long-term variabilities of East Asian summer monsoon (EASM) in 1958-2017 through examining diurnal cycles. We group monsoon days into four dynamic quadrants, with emphasis on the strong daily southerlies coupled with a large (Q1) or small (Q4) diurnal amplitude over Southeast China. The occurrence day of Q1 increases in June-July with the seasonal progress of EASM. It is most pronounced in 1960s-1970s and declines to the lowest in 1980s-1990s, while the Q4 occurrence increases notably from 1970s to 1990s; both groups return to normal in recent years. The interdecadal decrease (increase) of Q1 (Q4) occurrence corresponds well to the known weakening of EASM in the 20th century, and it also coincides with the rainfall anomalies over China shifting from “North flooding and South drought” to “North drought and South flooding” modes. The rainfall under Q1 (Q4) can account for ∼60% of the interannual variance of summer rainfall in northern (southern) China. The contrasting effects of Q1 and Q4 on rainfall are due to their remarkably different regulation on water vapor transports and convergence. The interannual/interdecadal variations of Q1 (Q4) occurrence determine the anomalous water vapor transports to northern (southern) China, in association with the various expansion of the western Pacific subtropical high. In particular, Q1 condition can greatly intensify nighttime moisture convergence, which is responsible for the long-term variations of rainfall in northern China. The results highlight that the diurnal cycles in monsoon flow act as a key regional process working with large-scale circulations to regulate the spatial distributions and long-term variabilities of EASM rainfall.

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 lightning risk assessment for application to human safety was created and applied in 10 West Texas locations from 2 May 2016 to 30 September 2016. The method combined spatial lightning mapping data, probabilistic risk calculation adapted from the International Electrotechnical Commission Standard 62305-2, and weighted average interpolation to produce risk magnitudes that were compared to tolerability thresholds to issue lightning warnings. These warnings were compared to warnings created for the same dataset using a more standard lightning safety approach based on National Lightning Detection Network (NLDN) total lightning within 5 nautical miles of each location. Four variations of the calculation as well as different units of risk were tested to find the optimal configuration to calculate risk to an isolated human outdoors.The best performing risk configuration using risk 10min−1 or larger produced the most comparable results to the standard method, such as number of failures, average warning duration, and total time under warnings. This risk configuration produced fewer failures than the standard method, but longer total time under warnings and higher false alarm ratios. Median lead times associated with the risk configuration were longer than the standard method for all units considered, while median down times were shorter for risk 10min−1 and risk 15min−1. Overall, the risk method provides a baseline framework to quantify the changing lightning hazard on the storm-scale, and could be a useful tool to aid in lightning decision support scenarios.

Description: Baroclinic waves drive both regional variations in weather and large-scale variability in the extratropical general circulation. They generally do not exist in isolation, but rather often form into coherent wave packets that propagate to the east via a mechanism called downstream development. Downstream development has been widely documented and explored. Here we document a novel but also key aspect of baroclinic waves: the downstream suppression of baroclinic activity that occurs in the wake of eastward propagating disturbances. Downstream suppression is apparent not only in the Southern Hemisphere storm track as shown in previous work, but also in the North Pacific and North Atlantic storm tracks. It plays an essential role in driving subseasonal periodicity in extratropical eddy activity in both hemispheres, and gives rise to the observed quiescence of the North Atlantic storm track 1–2 weeks following pronounced eddy activity in the North Pacific sector. It is argued that downstream suppression results from the anomalously low baroclinicity that arises as eastward propagating wave packets convert potential to kinetic energy. In contrast to baroclinic wave packets, which propagate to the east at roughly the group velocity in the upper troposphere, the suppression of baroclinic activity propagates eastward at a slower rate that is comparable to that of the lower to midtropospheric flow. The results have implications for understanding subseasonal variability in the extratropical troposphere of both hemispheres.

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: 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: Persistent multiyear cold states of the tropical Pacific Ocean drive hydroclimate anomalies worldwide, including persistent droughts in the extratropical Americas. Here, the atmosphere and ocean dynamics and thermodynamics of multiyear cold states of the tropical Pacific Ocean are investigated using European Centre for Medium-Range Weather Forecasts reanalyses and simplified models of the ocean and atmosphere. The cold states are maintained by anomalous ocean heat flux divergence and damped by increased surface heat flux from the atmosphere to ocean. The anomalous ocean heat flux divergence is contributed to by both changes in the ocean circulation and thermal structure. The keys are an anomalously shallow thermocline that enhances cooling by upwelling and anomalous westward equatorial currents that enhance cold advection. The thermocline depth anomalies are shown to be a response to equatorial wind stress anomalies. The wind stress anomalies are shown to be a simple dynamical response to equatorial SST anomalies as mediated by precipitation anomalies. The cold states are fundamentally maintained by atmosphere-ocean coupling in the equatorial Pacific. The physical processes that maintain the cold states are well approximated by linear dynamics for ocean and atmosphere and simple thermodynamics.

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

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

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

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

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

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

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

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

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

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

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

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

Description: Surface snowfall rate estimates from the Global Precipitation Measurement (GPM) mission’s Core Observatory sensors and the CloudSat radar are compared to those from the Multi-Radar Multi-Sensor (MRMS) radar composite product over the continental United States during the period from November 2014 to September 2020. The analysis includes: the Dual-Frequency Precipitation Radar (DPR) retrieval and its single frequency counterparts, the GPM Combined Radar Radiometer Algorithm (CORRA), the CloudSat Snow Profile product (2C-SNOW-PROFILE) and two passive microwave retrievals, i.e., the Goddard PROFiling algorithm (GPROF) and the Snow retrievaL ALgorithm fOr gMi (SLALOM). The 2C-SNOW retrieval has the highest Heidke Skill Score (HSS) for detecting snowfall among the products analysed. SLALOM ranks second; it outperforms GPROF and the other GPM algorithms, all detecting only 30% of the snow events. Since SLALOM is trained with 2C-SNOW, it suggests that the optimal use of the information content in the GMI observations critically depends on the precipitation training dataset. All the retrievals underestimate snowfall rates by a factor of two compared to MRMS. Large discrepancies (RMSE of 0.7 to 1.5 mm h-1) between space-borne and ground-based snowfall rate estimates are attributed to the complexity of the ice scattering properties and to the limitations of the remote sensing systems: the DPR instrument has low sensitivity, while the radiometric measurements are affected by the confounding effects of the background surface emissivity and of the emission of supercooled liquid droplet layers.

Description: Quantitative precipitation estimates (QPE) at high spatiotemporal resolution are essential for flash flood forecasting, especially in urban environments and headwater areas. An accurate quantification of precipitation is directly related to the temporal and spatial sampling of the precipitation system. The advent of phased array radar (PAR) technology, a potential next-generation weather radar, can provide updates that are at least 4-5 times faster than the conventional WSR-88D scanning rate. In this study, data collected by the KOUN WSR-88D radar with ~1 minute temporal resolution is used as an approximation of data that a future PAR system could provide to force the Ensemble Framework for Flash Flood Forecasting (EF5) hydrologic model. To assess the effect of errors resulting from temporal and spatial sampling of precipitation on flash flood warnings, KOUN precipitation data (1-km/1-min) is used to generate precipitation products at other spatial/temporal resolutions commonly used in hydrologic models, such as those provided by conventional WSR-88D radar (1-km/5-min), spaced-based observations (10-km/30-min), and hourly rainfall products (1-km/60-min). The effect of precipitation sampling errors on flash flood warnings are then examined and quantified by using discharge simulated from KOUN (1-km/1-min) as truth to assess simulations conducted using other generated coarser spatial/temporal resolutions of other precipitation products. Our results show that: 1) observations with coarse spatial and temporal sampling can cause large errors in quantification of the amount, intensity, and distribution of precipitation, 2) time series of precipitation products show that QPE peak values decrease as the temporal resolution gets coarser, and 3) the effect of precipitation sampling error on flash flood forecasting is large in headwater areas and decrease quickly as drainage area increases.

Description: Measuring rainfall is complex, due to the high temporal and spatial variability of precipitation, especially in a changing climate, but it is of great importance for all the scientific and operational disciplines dealing with rainfall effects on the environment, human activities, and economy.Microwave (MW) telecommunication links carry information on rainfall rates along their path, through signal attenuation caused by raindrops, and can become measurements of opportunity, offering inexpensive chances to augment information without deploying additional infrastructures, at the cost of some smart processing. Processing satellite telecom signals bring some specific complexities related to the effects of rainfall boundaries, melting layer, and non-weather attenuations, but with the potential to provide worldwide precipitation data with high temporal and spatial samplings. These measurements have to be processed according to the probabilistic nature of the information they carry. An EnKF-based (Ensemble Kalman Filter) method has been developed to dynamically retrieve rainfall fields in gridded domains, which manages such probabilistic information and exploits the high sampling rate of measurements. The paper presents the EnKF method with some representative tests from synthetic 3D experiments. Ancillary data are assumed as from worldwide-available operational meteorological satellites and models, for advection, initial and boundary conditions, rain height. The method reproduces rainfall structures and quantities in a correct way, and also manages possible link outages. It results computationally viable also for operational implementation and applicable to different link observation geometries and characteristics.

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

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

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

Description: Tropical cyclones (TCs) propagating into baroclinic midlatitude environments can transform into extratropical cyclones, in some cases resulting in high-impact weather conditions far from the tropics. This study extends analysis of extratropical transition (ET) changes in multi-seasonal global simulations using the Model for Prediction Across Scales-Atmosphere (MPAS-A) under present-day and projected future conditions. High-resolution (15 km) covers the Northern Hemisphere; TCs and ET events are tracked based on sea-level pressure minima accompanied by a warm core and use of a cyclone phase space method. Previous analysis of these simulations showed large changes in ET over the North Atlantic (NATL) basin, with ET events exhibiting a 4–5° northward latitudinal shift and a ~6 hPa strengthening of the post-transition extratropical cyclone. Storm-relative composites, primarily representing post-transformation cold-core events, indicate that this increase in post-transition storm intensity is associated with an intensification of the neighboring upper-level trough and downstream ridge, and a poleward shift in the storm center, conducive to enhanced trough-TC interactions after ET completion. Additionally, the future composite ET event is located in the right-jet entrance of an outflow jet that is strengthened relative to its present-day counterpart. Localized impacts associated with ET events, such as heavy precipitation and strong near-surface winds, are significantly enhanced in the future-climate simulations; 6-hourly precipitation for NATL events increases at a super-Clausius-Clapeyron rate with area-average precipitation increasing over 30%. Furthermore, intensified precipitation contributes to enhanced lower-tropospheric potential vorticity and stronger upper-tropospheric outflow, implying the potential for more extreme downstream impacts under the future climate scenario.

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

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

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

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

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

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

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

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

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

Description: Assimilation of dual-polarization (dual-pol) observations provides more accurate storm-scale analyses to initialize forecasts of severe convective thunderstorms. This study investigates the impact assimilating experimental sector-scan dual-pol observations has on storm-scale ensemble forecasts and how this impact changes over different data assimilation (DA) windows using the ensemble Kalman filter (EnKF). Ensemble forecasts are initialized after 30, 45, and 60 minutes of DA for two sets of experiments that assimilate either reflectivity and radial velocity only (EXPZ) or reflectivity and radial velocity plus differential reflectivity (EXPZZDR). This study uses the 31 May 2013 Oklahoma event which included multiple storms that produced tornadoes and severe hail, with focus placed on two storms that impacted El Reno and Stillwater during the event.The earliest initialized forecast of EXPZZDR better predicts the evolution of the El Reno storm than EXPZ, but the two sets of experiments become similar at subsequent forecast times. However, the later EXPZZDR forecasts of the Stillwater storm, which organized towards the end of the DA window, produce improved results compared to EXPZ, in which the storm is less intense and weakens. Evaluation of forecast products for supercell mesocyclones (updraft helicity [UH]) and hail show similar results with earlier EXPZZDR forecasts better predicting the UH swaths of the El Reno storm and later forecasts producing improved UH and hail swaths for the Stillwater storm. The results indicate that the assimilation of ZDR over fewer DA cycles can produce improved forecasts when DA windows sufficiently cover storms during their initial development and organization.

Description: The characteristics of El-Niño-Southern Oscillation (ENSO) phase-locking in observations and CMIP5 and CMIP6 models are examined in this study. Two metrics based on the peaking month histogram for all El Niño and La Niña events are adopted to delineate the basic features of ENSO phase-locking in terms of the preferred calendar month and strength of this preference. It turns out that most models are poor at simulating the ENSO phase-locking, either showing little peak strengths or peaking at the wrong seasons. By deriving ENSO’s linear dynamics based on the conceptual recharge oscillator (RO) framework through the seasonal linear inverse model (sLIM) approach, various simulated phase-locking behaviors of CMIP models are systematically investigated in comparison with observations. In observations, phase-locking is mainly attributed to the seasonal modulation of ENSO’s SST growth rate. In contrast, in a significant portion of CMIP models, phase-locking is co-determined by the seasonal modulations of both SST growth and phase-transition rates. Further study of the joint effects of SST growth and phase-transition rates suggests that for simulating realistic winter peak ENSO phase-locking with the right dynamics, climate models need to have four key factors in the right combination: (1) correct phase of SST growth rate modulation peaking at the fall; (2) large enough amplitude for the annual cycle in growth rate; (3) amplitude of semi-annual cycle in growth rate needs to be small; and (4) amplitude of seasonal modulation in SST phase-transition rate needs to be small.

Description: Undisturbed trade-wind conditions comprise the most prevalent synoptic weather pattern in Hawai’i and produces a distinct pattern of orographic rainfall. Significant total rainfall contributions and extreme events are linked to four types of atmospheric disturbances: cold fronts, Kona lows, upper-tropospheric disturbances, and tropical cyclones. In this study, a 20- year (1990-2010) categorical disturbance time series is compiled and analyzed in relation to daily rainfall over the same period. The primary objective of this research is to determine how disturbances contribute to total wet season rainfall on the Island of O’ahu, Hawai’i. On average, 41% of wet seasonal rainfall occurs on disturbance days. Seventeen percent of seasonal rainfall can be directly attributed to disturbances (after a background signal is removed) and as much as 48% in a single season. The intensity of disturbance rainfall (mm/day) is a stronger predictor (r2 = 0.49; p 〈 0.001) of the total seasonal rainfall than the frequency of occurrence (r2 = 0.11; p = 0.153). Cold fronts are the most common disturbance type; however, the rainfall associated with fronts that cross the island is significantly higher than rainfall produced from non-crossing fronts. In fact, non-crossing fronts produce significantly less rainfall than under mean non-disturbance conditions 76% of the time. While the combined influence of atmospheric disturbances can account for almost half of the rainfall received during the wet season, the primary factor in determining a relatively wet or dry season/year on O’ahu are the frequency and rainfall intensity of Kona Low events.

Description: We investigate how sea ice decline in summer and warmer ocean and surface temperatures in winter affect sea ice growth in the Arctic. Sea ice volume changes are estimated from satellite observations during winter from 2002 to 2019 and partitioned into thermodynamic growth and dynamic volume change. Both components are compared to validated sea ice-ocean models forced by reanalysis data to extend observations back to 1980 and to understand the mechanisms that cause the observed trends and variability. We find that a negative feedback driven by the increasing sea ice retreat in summer yields increasing thermodynamic ice growth during winter in the Arctic marginal seas eastward from the Laptev Sea to the Beaufort Sea. However, in the Barents and Kara Seas, this feedback seems to be overpowered by the impact of increasing oceanic heat flux and air temperatures, resulting in negative trends in thermodynamic ice growth of -2 km3month-1yr-1 on average over 2002-2019 derived from satellite observations.

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

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