550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (902) (remove)
Has Fulltext
- yes (902)
Keywords
- Climate change (7)
- climate change (7)
- COSMO-CLM (6)
- Palaeoclimate (6)
- Atmospheric chemistry (5)
- Biogeochemistry (5)
- Palaeoceanography (5)
- precipitation (5)
- loess (4)
- Clausius–Clapeyron scaling (3)
Institute
- Geowissenschaften (546)
- Geowissenschaften / Geographie (106)
- Biodiversität und Klima Forschungszentrum (BiK-F) (59)
- Senckenbergische Naturforschende Gesellschaft (51)
- Extern (22)
- Biowissenschaften (21)
- Geographie (15)
- Präsidium (12)
- Institut für Ökologie, Evolution und Diversität (9)
- Physik (6)
- Institut für sozial-ökologische Forschung (ISOE) (5)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- Medizin (4)
- SFB 268 (4)
- Starker Start ins Studium: Qualitätspakt Lehre (3)
- Geschichtswissenschaften (2)
- Informatik und Mathematik (2)
- Zentrum für Interdisziplinäre Afrikaforschung (ZIAF) (2)
- Biochemie und Chemie (1)
- Center for Scientific Computing (CSC) (1)
- Gesellschaftswissenschaften (1)
- keine Angabe Institut (1)
Ice nucleating particle concentrations of the past: insights from a
600 year old Greenland ice core
(2020)
Ice nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 ◦C to −35 ◦C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to 1 sample per year. Additionally, a number of special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity and the wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L -1 at −25 ◦C. The INP variability was found to be about 1 – 2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to chemical tracers derived from continuous flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (dp > 1.2 µm). The correlation is higher for the seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) of single particles retrieved from selected samples found particles of soil origin to be the dominant fraction, verifying the significance of mineral dust particles as INPs. Overall, the concentrations compare reasonably well to present day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23 ◦C, −24 ◦C, −25 ◦C, −26 ◦C and −28 ◦C were significantly higher (and more variable) in the modern day period, which could indicate a potential anthropogenic impact or some post-coring contamination of the topmost, very porous firn.
Recently significant advances have been made in the collection, detection, and characterization of ice nucleating particles (INP). Ice nuclei are particles that facilitate the heterogeneous formation of ice within the atmospheric aerosol by lowering the free energy barrier to spontaneous nucleation and growth of ice from atmospheric water and/or vapor. The Frankfurt isostatic diiffusion chamber (FRIDGE) is an INP collection and offline detection system that has become widely deployed and shows additional potential for ambient measurements. Since its initial development FRIDGE has gone through several iterations and improvements. Here we describe improvements that have been made in the collection and analysis techniques. We detail the uncertainties inherent in the measurement method, and suggest a systematic method of error analysis for FRIDGE measurements. Thus what is presented herein should serve as a foundation for the dissemination of all current and future measurements using FRIDGE instrumentation.
Recently significant advances have been made in the collection, detection and characterization of ice nucleating particles (INPs). Ice nuclei are particles that facilitate the heterogeneous formation of ice within the atmospheric aerosol by lowering the free energy barrier to spontaneous nucleation and growth of ice from atmospheric water and/or vapor. The Frankfurt isostatic diffusion chamber (FRankfurt Ice nucleation Deposition freezinG Experiment: FRIDGE) is an INP collection and offline detection system that has become widely deployed and shows additional potential for ambient measurements. Since its initial development FRIDGE has gone through several iterations and improvements. Here we describe improvements that have been made in the collection and analysis techniques. We detail the uncertainties inherent in the measurement method and suggest a systematic method of error analysis for FRIDGE measurements. Thus what is presented herein should serve as a foundation for the dissemination of all current and future measurements using FRIDGE instrumentation.
The formation of particles from precursor vapors is an important source of atmospheric aerosol. Research at the Cosmics Leaving OUtdoor Droplets (CLOUD) facility at CERN tries to elucidate which vapors are responsible for this new-particle formation, and how in detail it proceeds. Initial measurement campaigns at the CLOUD stainless-steel aerosol chamber focused on investigating particle formation from ammonia (NH3) and sulfuric acid (H2SO4). Experiments were conducted in the presence of water, ozone and sulfur dioxide. Contaminant trace gases were suppressed at the technological limit. For this study, we mapped out the compositions of small NH3–H2SO4 clusters over a wide range of atmospherically relevant environmental conditions. We covered [NH3] in the range from < 2 to 1400 pptv, [H2SO4] from 3.3 × 106 to 1.4 × 109 cm−3 (0.1 to 56 pptv), and a temperature range from −25 to +20 °C. Negatively and positively charged clusters were directly measured by an atmospheric pressure interface time-of-flight (APi-TOF) mass spectrometer, as they initially formed from gas-phase NH3 and H2SO4, and then grew to larger clusters containing more than 50 molecules of NH3 and H2SO4, corresponding to mobility-equivalent diameters greater than 2 nm. Water molecules evaporate from these clusters during sampling and are not observed. We found that the composition of the NH3–H2SO4 clusters is primarily determined by the ratio of gas-phase concentrations [NH3] / [H2SO4], as well as by temperature. Pure binary H2O–H2SO4 clusters (observed as clusters of only H2SO4) only form at [NH3] / [H2SO4] < 0.1 to 1. For larger values of [NH3] / [H2SO4], the composition of NH3–H2SO4 clusters was characterized by the number of NH3 molecules m added for each added H2SO4 molecule n (Δm/Δ n), where n is in the range 4–18 (negatively charged clusters) or 1–17 (positively charged clusters). For negatively charged clusters, Δ m/Δn saturated between 1 and 1.4 for [NH3] / [H2SO4] > 10. Positively charged clusters grew on average by Δm/Δn = 1.05 and were only observed at sufficiently high [NH3] / [H2SO4]. The H2SO4 molecules of these clusters are partially neutralized by NH3, in close resemblance to the acid–base bindings of ammonium bisulfate. Supported by model simulations, we substantiate previous evidence for acid–base reactions being the essential mechanism behind the formation of these clusters under atmospheric conditions and up to sizes of at least 2 nm. Our results also suggest that electrically neutral NH3–H2SO4 clusters, unobservable in this study, have generally the same composition as ionic clusters for [NH3] / [H2SO4] > 10. We expect that NH3–H2SO4 clusters form and grow also mostly by Δm/Δn > 1 in the atmosphere's boundary layer, as [NH3] / [H2SO4] is mostly larger than 10. We compared our results from CLOUD with APi-TOF measurements of NH3–H2SO4 anion clusters during new-particle formation in the Finnish boreal forest. However, the exact role of NH3–H2SO4 clusters in boundary layer particle formation remains to be resolved.
The formation of particles from precursor vapors is an important source of atmospheric aerosol. Research at the Cosmics Leaving OUtdoor Droplets (CLOUD) facility at CERN tries to elucidate which vapors are responsible for this new particle formation, and how in detail it proceeds. Initial measurement campaigns at the CLOUD stainless-steel aerosol chamber focused on investigating particle formation from ammonia (NH3) and sulfuric acid (H2SO4). Experiments were conducted in the presence of water, ozone and sulfur dioxide. Contaminant trace gases were suppressed at the technological limit. For this study, we mapped out the compositions of small NH3-H2SO4 clusters over a wide range of atmospherically relevant environmental conditions. We covered [NH3] in the range from <2 to 1400 pptv, [H2SO4] from 3.3 × 106 to 1.4 × 109 cm−3, and a temperature range from −25 to +20 °C. Negatively and positively charged clusters were directly measured by an atmospheric pressure interface time-of-flight (APi-TOF) mass spectrometer, as they initially formed from gas-phase NH3 and H2SO4, and then grew to larger clusters containing more than 50 molecules of NH3 and H2SO4, corresponding to mobility-equivalent diameters greater than 2 nm. Water molecules evaporate from these clusters during sampling and are not observed. We found that the composition of the NH3-H2SO4 clusters is primarily determined by the ratio of gas-phase concentrations [NH3] / [H2SO4], as well as by temperature. Pure binary H2O-H2SO4 clusters (observed as clusters of only H2SO4) only form at [NH3] / [H2SO4]<0.1 to 1. For larger values of [NH3] / [H2SO4], the composition of NH3-H2SO4 clusters was characterized by the number of NH3 molecules m added for each added H2SO4 molecule n (Δm / Δn), where n is in the range 4–18 (negatively charged clusters) or 1–17 (positively charged clusters). For negatively charged clusters, Δm / Δn saturated between 1 and 1.4 for [NH3] / [H2SO4]>10. Positively charged clusters grew on average by Δm / Δn = 1.05 and were only observed at sufficiently high [NH3] / [H2SO4]. The H2SO4 molecules of these clusters are partially neutralized by NH3, in close resemblance to the acid-base bindings of ammonium bisulfate. Supported by model simulations, we substantiate previous evidence for acid-base reactions being the essential mechanism behind the formation of these clusters under atmospheric conditions and up to sizes of at least 2 nm. Our results also suggest that yet unobservable electrically neutral NH3-H2SO4 clusters grow by generally the same mechanism as ionic clusters, particularly for [NH3] / [H2SO4]>10. We expect that NH3-H2SO4 clusters form and grow also mostly by Δm / Δn>1 in the atmosphere's boundary layer, as [NH3] / [H2SO4] is mostly larger than 10. We compared our results from CLOUD with APi-TOF measurements of NH3-H2SO4 anion clusters during new particle formation in the Finnish boreal forest. However, the exact role of NH3-H2SO4 clusters in boundary layer particle formation remains to be resolved.
Die Geologie des Elsaß
(1940)
Tropical cyclones (TC) represent a substantial threat to life and property for Caribbean and adjacent populations. The prospective increase of TC magnitudes, expressed in the 15th chapter of the IPCC AR5 report, entails a rising probability of ecological and social disasters, which were tragically exemplified by several severe Caribbean TC strikes during the past 20 years. Modern IPCC-grade climate models, however, still lack the required spatial and temporal resolution to accurately consider the underlying boundary conditions that modulate long-time TC patterns beyond the Instrumental Era. It is thus necessary to provide a synoptic mechanistic understanding regarding the origin of such long-time patterns, in order to predict reliable changes of TC magnitude and frequency under future climate scenarios. Caribbean TC records are still rare and often lack the necessary continuity and resolution to overcome these limitations. Here, we report on an annually-resolved sedimentary archive from the bottom of the Great Blue Hole (Lighthouse Reef, Belize). The TC record encompasses 1885 years and extends all existing site-specific TC archives both in terms of resolution and duration. We identified a likely connection between long-term TC patterns and climate phenomena responses to Common Era climate variations and offer a conceptual and comparative view considering several involved tropospheric and oceanographic control mechanisms such as the El-Niño-Southern-Oscillation, the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. These basin-scaled climate modes exercise internal control on TC activity by modulating the thermodynamic environment (sea-surface temperature and vertical wind shear stress dynamics) for enhanced/suppressed TC formation both on millennial (primary) and multi-decadal (secondary) time scales. We interpret the beginning of the Medieval Warm Period (MWP) as an important time interval of the Common Era record and suspect that the southward migration of the intertropical convergence zone (ITCZ) caused, in combination with extensive hydro-climate changes, a shift in the tropical Atlantic TC regime. The TC activity in the south-western Caribbean changed in general from a stable and less active stage (100–900 CE) to a more active and variable state (1,100 CE-modern).
In situ single particle analysis of ice particle residuals (IPRs) and out-of-cloud aerosol particles was conducted by means of laser ablation mass spectrometry during the intensive INUIT-JFJ/CLACE campaign at the high alpine research station Jungfraujoch (3580 m a.s.l.) in January–February 2013. During the 4-week campaign more than 70 000 out-of-cloud aerosol particles and 595 IPRs were analyzed covering a particle size diameter range from 100 nm to 3 µm. The IPRs were sampled during 273 h while the station was covered by mixed-phase clouds at ambient temperatures between −27 and −6 °C. The identification of particle types is based on laboratory studies of different types of biological, mineral and anthropogenic aerosol particles. The outcome of these laboratory studies was characteristic marker peaks for each investigated particle type. These marker peaks were applied to the field data. In the sampled IPRs we identified a larger number fraction of primary aerosol particles, like soil dust (13 ± 5 %) and minerals (11 ± 5 %), in comparison to out-of-cloud aerosol particles (2.4 ± 0.4 and 0.4 ± 0.1 %, respectively). Additionally, anthropogenic aerosol particles, such as particles from industrial emissions and lead-containing particles, were found to be more abundant in the IPRs than in the out-of-cloud aerosol. In the out-of-cloud aerosol we identified a large fraction of aged particles (31 ± 5 %), including organic material and secondary inorganics, whereas this particle type was much less abundant (2.7 ± 1.3 %) in the IPRs. In a selected subset of the data where a direct comparison between out-of-cloud aerosol particles and IPRs in air masses with similar origin was possible, a pronounced enhancement of biological particles was found in the IPRs.
In-situ single particle analysis of ice particle residuals (IPR) and out-of-cloud aerosol particles was conducted by means of laser ablation mass spectrometry during the intensive INUIT-JFJ/CLACE campaign at the high alpine research station Jungfraujoch (3580 m a.s.l.) in January/February 2013. During the four week campaign more than 70000 out-of-cloud aerosol particles and 595 IPR were analyzed covering a particle size diameter range from 100 nm to 3 μm. The IPR were sampled during 273 hours while the station was covered by mixed-phase clouds at ambient temperatures between -27 °C and -6 °C. The identification of particle types is based on laboratory studies of different types of biological, mineral and anthropogenic aerosol particles. As outcome instrument specific marker peaks for the different investigated particle types were obtained and applied to the field data. The results show that the sampled IPR contain a larger relative amount of natural, primary aerosol, like soil dust (13 %) and minerals (11 %), in comparison to out-of-cloud aerosol particles (2 % and <1 %, respectively). Additionally, anthropogenic aerosol particles, like particles from industrial emissions and lead-containing particles, were found to be more abundant in the IPR than in the out-of-cloud aerosol. The out of-cloud aerosol contained a large fraction of aged particles (30 %, including organic material and secondary inorganics), whereas this particle type was much less abundant (3 %) in the IPR. In a selected subset of the data where a direct comparison between out-of-cloud aerosol particles and IPR in air masses with similar origin was possible, a pronounced enhancement of biological particles was found in the IPR.
This paper presents results from the "INUIT-JFJ/CLACE 2013" field campaign at the high alpine research station Jungfraujoch in January/February 2013. The chemical composition of ice particle residuals (IPR) in a size diameter range of 200–900 nm was measured in orographic, convective and non-convective clouds with a single particle mass spectrometer (ALABAMA) under ambient conditions characterized by temperatures between −28 and −4 °C and wind speed from 0.1 to 21 km h−1. Additionally, background aerosol particles in cloud free air were investigated. The IPR were sampled from mixed-phase clouds with two inlets which selectively extract small ice crystals in-cloud, namely the Counterflow Virtual Impactor (Ice-CVI) and the Ice Selective Inlet (ISI). The IPR as well as the aerosol particles were classified into seven different particle types: (1) black carbon, (2) organic carbon, (3) black carbon internally mixed with organic carbon, (4) minerals, (5) one particle group (termed "BioMinSal") that may contain biological particles, minerals, or salts, (6) industrial metals, and (7) lead containing particles. For any sampled particle population it was determined by means of single particle mass spectrometer how many of the analyzed particles belonged to each of these categories. Accordingly, between 20 and 30% of the IPR and roughly 42% of the background particles contained organic carbon. The measured fractions of minerals in the IPR composition varied from 6 to 33%, while the values for the "BioMinSal" group were between 15 and 29%. Four percent to 31% of the IPR contained organic carbon mixed with black carbon. Both inlets delivered similar results of the chemical composition and of the particle size distribution, although lead was found only in the IPR sampled by the Ice-CVI. The results show that the ice particle residual composition varies substantially between different cloud events, which indicates the influence of different meteorological conditions, such as origin of the air masses, temperature and wind speed.
A thick Middle and Late Pleistocene loess/palaeosol sequence is exposed at the gravel quarry Gaul located east of Weilbach in the southern foreland of the Taunus Mountains. The loess/palaeosol sequence correlates to the last three glacial cycles. Seven samples were dated by luminescence methods using an elevated temperature IRSL (post-IR IRSL) protocol for polymineral fine-grains to determine the deposition age of the sediment and to set up a more reliable chronological framework for these deposits. The fading corrected IR50 and the pIRIR225 age estimates show a good agreement for almost all samples. The fading corrected IRSL ages range from 23.7 ± 1.6 ka to >350 ka indicating that the oldest loess was deposited during marine isotope stage (MIS) 10 or earlier and that the humic-rich horizon (Weilbacher Humuszone) was developed during the late phase of MIS 7. Loess taken above the fCc horizon most likely accumulated during MIS 6 indicating that the remains of the palaeosol are not belonging to the last interglacial soil. The two uppermost samples indicate that the youngest loess accumulated during MIS 2 (Upper Würmian). Age estimates for the loess-palaeosol sequence of the gravel quarry Gaul/Weilbach could be obtained up to ~350 ka using the pIRIR225 from feldspar. Keywords: loess, luminescence dating, IRSL, fading, Weilbach, chronostratigraphy
We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a numerical weather prediction model at km-scale resolution. In addition, the influence of grid resolution (2.2 km and 1.1 km), topography filtering, and land surface datasets on the accuracy of the simulated valley winds is investigated. The simulations are evaluated against a comprehensive set of surface observations for an 18-day fair-weather summer period in July 2006. The episode is characterized by strong diurnal wind systems and the formation of shallow convection over the mountains, which transitions to precipitating convection in some areas. The near-surface winds (10 m above ground level) follow a typical diurnal pattern with strong daytime up-valley flow and weaker nighttime down-valley flow. At a 2.2 km resolution the valley winds are poorly simulated for most stations, while at a 1.1 km resolution the diurnal cycle of the valley winds is well represented in most large (e.g., Rhein valley at Chur and Rhone valley at Visp) and medium-sized valleys (e.g., Linth valley at Glarus). In the smaller valleys (e.g., Maggia valley at Cevio), the amplitude of the valley wind is still significantly underestimated, even at a 1.1 km resolution. Detailed sensitivity experiments show that the use of high-resolution land surface datasets, for both the soil characteristics as well as for the land cover, and reduced filtering of the topography are essential to achieve good performance at a 1.1 km resolution
Diurnal valley winds frequently form over complex topography, particularly under fair weather conditions, and have a significant impact on the local weather and climate. Since diurnal valley winds result from complex and multi-scale interactions, their representation in numerical weather prediction models is challenging. Better understanding of these local winds based on observations is crucial to improve the accuracy of the forecasts. This study investigates the diurnal evolution of the three-dimensional mean wind structure in a deep Alpine valley, the Rhone valley at Sion, using data from a radar wind profiler and a surface weather station operated continuously from 1 September 2016 to 17 July 2017. In particular, the wind profiler data was analyzed for a subset of days on which fair weather conditions allowed for the full development of thermally driven winds. A pronounced diurnal cycle of the wind speed, as well as a reversal of the wind direction twice per day is documented for altitudes up to about 2 km above ground level (AGL) in the warm season and less than 1 km AGL in winter. The diurnal pattern undergoes significant changes during the course of the year. Particularly during the warm-weather months of May through to September, a low-level wind maximum occurs, where mean maximum up-valley velocities of 8–10 m s−1 are found between 15–16 UTC at altitudes around 200 m AGL. In addition, during nighttime, a down-valley jet with maximum wind speeds of 4–8 m s−1 around 1 km AGL is found. A case study of a three-day period in September 2016 illustrates the occurrence of an elevated layer of cross-valley flow around 1–1.5 km AGL.
Highlights
• Solidification and cooling of an intruded dyke or sill within the middle or shallow crust generate stresses of order 200 MPa, which relax on time scales of a few years to million years.
• Stresses may exceed the brittle strength forming tensile fractures.
• Combined with the pressure gradient within the over-pressurized felsic melts, this explains the migration of felsic contact melt into shrinkage cracks (Sederholm-type veins).
Abstract
Rapid emplacement of a mafic dyke or sill at mid-crustal depth heats and possibly melts the felsic wall rock followed by solidification. Associated volume changes generate stresses, possibly enforcing brittle failure and melt migration. We model the evolution of melting, solidification, temperature, and stress including visco-elastic relaxation in 1D - dykes or -sills using realistic rock rheologies of the Weschnitz pluton (Odenwald). For deep emplacement (Case 1, 15.3 km) extensive contact melting of the wall rock occurs, for shallow emplacement (Case 2, 10 km) it is negligible. The stresses are zero at high melt fractions, but increase during solidification and cooling: The intrusion orthogonal stress is always zero. The intrusion parallel stress σ‖ within the intrusion is tensile (O(200 MPa)). It relaxes on a time scale between a few years (Case 1) and 0.6 m.y. (Case 2). Within the wall rock σ‖ is compressive during heating, but becomes tensile under solidification and cooling. Wall rock stresses relax on a time scale of months to 100 years. A Deborah number is defined based on viscous to thermal relaxation allowing generalization of our results. Adding lithostatic stresses, the total stresses of Case 1 remain below the brittle strength, while for Case 2 they may exceed it. Adding the lithostatic pressure to the melt pressure, the effective stresses exceed the brittle strength and intrusion orthogonal tensile fractures are predicted. Combined with the pressure gradient within the over-pressurized felsic melts generated in the wall rock, this explains the migration of felsic contact melt into shrinkage cracks of the mafic sill in the Weschnitz pluton.
Ecological networks are more sensitive to plant than to animal extinction under climate change
(2016)
Impacts of climate change on individual species are increasingly well documented, but we lack understanding of how these effects propagate through ecological communities. Here we combine species distribution models with ecological network analyses to test potential impacts of climate change on >700 plant and animal species in pollination and seed-dispersal networks from central Europe. We discover that animal species that interact with a low diversity of plant species have narrow climatic niches and are most vulnerable to climate change. In contrast, biotic specialization of plants is not related to climatic niche breadth and vulnerability. A simulation model incorporating different scenarios of species coextinction and capacities for partner switches shows that projected plant extinctions under climate change are more likely to trigger animal coextinctions than vice versa. This result demonstrates that impacts of climate change on biodiversity can be amplified via extinction cascades from plants to animals in ecological networks.
Here we present a comprehensive attempt to correlate aragonitic Na∕Ca ratios from Desmophyllum pertusum (formerly known as Lophelia pertusa), Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na∕Ca ratios using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements.
The results reveal no apparent correlation with salinity (30.1–40.57 g kg−1) but a significant inverse correlation with temperature (−0.31±0.04 mmolmol−1∘C−1). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na∕Mg ratios show a similar temperature sensitivity to Na∕Ca ratios, but the combination of two ratios appears to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na∕Ca and Na∕Mg, however, limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na∕Ca ratios for an open and semi-enclosed calcifying space based on temperature-sensitive Na- and Ca-pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na∕Ca ratio.
Here we present a comprehensive attempt to correlate aragonitic Na / Ca ratios from Lophelia pertusa, Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na / Ca content using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements. The results reveal no apparent correlation with salinity (30.1–40.57 g/kg) but a significant inverse correlation with temperature (−0.31 mmol/mol/°C). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na / Mg ratios show a similar temperature sensitivity to Na / Ca ratios, but the combination of two ratios appear to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na / Ca and Na / Mg however limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na / Ca ratios for an open and semi-enclosed calcifying space based on temperature sensitive Na and Ca pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na / Ca ratio.
Acesta excavata (Fabricius, 1779) is a slow growing bivalve from the Limidae family and is often found associated with cold-water coral reefs along the European continental margin. Here we present the compositional variability of frequently used proxy elemental ratios (Mg/ Ca, Sr/Ca, Na/Ca) measured by laser-ablation mass spectrometry (LA-ICP-MS) and com- pare it to in-situ recorded instrumental seawater parameters such as temperature and salin- ity. Shell Mg/Ca measured in the fibrous calcitic shell section was overall not correlated with seawater temperature or salinity; however, some samples show significant correlations with temperature with a sensitivity that was found to be unusually high in comparison to other marine organisms. Mg/Ca and Sr/Ca measured in the fibrous calcitic shell section display significant negative correlations with the linear extension rate of the shell, which indicates strong vital effects in these bivalves. Multiple linear regression analysis indicates that up to 79% of elemental variability is explicable with temperature and salinity as independent pre- dictor values. Yet, the overall results clearly show that the application of Element/Ca (E/Ca) ratios in these bivalves to reconstruct past changes in temperature and salinity is likely to be complicated due to strong vital effects and the effects of organic material embedded in the shell. Therefore, we suggest to apply additional techniques, such as clumped isotopes, in order to exactly determine and quantify the underlying vital effects and possibly account for these. We found differences in the chemical composition between the two calcitic shell lay- ers that are possibly explainable through differences of the crystal morphology. Sr/Ca ratios also appear to be partly controlled by the amount of magnesium, because the small magne- sium ions bend the crystal lattice which increases the space for strontium incorporation. Oxi- dative cleaning with H2O2 did not significantly change the Mg/Ca and Sr/Ca composition of the shell. Na/Ca ratios decreased after the oxidative cleaning, which is most likely a leaching effect and not caused by the removal of organic matter.
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.
Hyrrokkin sarcophaga is a parasitic foraminifer that is commonly found in cold-water coral reefs where it infests the file clam Acesta excavata and the scleractinian coral Lophelia pertusa. Here, we present measurements of the elemental and isotopic composition of this parasitic foraminifer for the first time, analyzed by inductively coupled optical emission spectrometry (ICP-OES), electron probe micro analysis (EPMA) and mass spectrometry (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 L. pertusa, which could be an indication that dissolved host carbonate material is utilised in shell calcification, given that the aragonite of L. pertusa 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 on20 L. pertusa, 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 foraminifers 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 foraminifer, these data also indicate that the geochemistry of this parasitic foraminifer 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.
An easy-to-use model to evaluate conductivities at high and middle latitudes in the height range 70–100 km is presented. It is based on electron density profiles obtained with the EISCAT VHF radar during 11 years and on the neutral atmospheric model MSIS95. The model uses solar zenith angle, geomagnetic activity and season as input parameters. It was mainly constructed to study the properties of Schumann resonances that depend on such conductivity profiles.
Kalkreiche Niedermoore
(2002)
Die Standorte der nährstoffarmen kalkreichen Sümpfe und Moore weisen oberflächennah anstehendes kalkreiches Grundwasser auf. Bei geringer Sauerstoffversorgung kommt es hier, bedingt durch schlechtere Humuszersetzung, zur Bildung kalkreicher Torfablagerungen. Teilweise werden die Standorte jedoch auch von sauerstoffreichem Quellwasser durchsickert. Hier findet eine wesentlich stärkere Zersetzung des Humus statt, wodurch keine Torfablagerungen entstehen und es lediglich zur Bildung von Kalksümpfen kommt. Diese Kalksümpfe weisen häufig eine mehr oder weniger ausgeprägte Kalktuffbildung auf. Bei Kontakt mit der Luft gibt das hier zutage tretende kalziumbikarbonatreiche Quellwasser Kohlendioxid ab. Dadurch wird Kalziumkarbonat frei. Die mit Wasser benetzten Teile der an diesen Standorten vorkommenden Pflanzen werden mit einer Kalkschicht überzogen, wodurch Kalktuff gebildet wird. Nur die kalkreichen Niedermoore in der subalpinen Stufe der Alpen sind weitgehend gehölzfrei (GÖRS 1974, in OBERDORFER 1977). In tieferen Lagen kommen auf diesen Standorten natürlicherweise Erlenbruchwälder aus dem Verband Alnion glutinosae vor. Durch Rodung des Waldes und anschließender jahrhundertelanger Beweidung oder gelegentlicher Mahd (SUCCOW u. JESCHKE 1986) entwickelten sich hier Kalkbinsensümpfe und Kalkkleinseggenriede.
Global impact models represent process-level understanding of how natural and human systems may be affected by climate change. Their projections are used in integrated assessments of climate change. Here we test, for the first time, systematically across many important systems, how well such impact models capture the impacts of extreme climate conditions. Using the 2003 European heat wave and drought as a historical analogue for comparable events in the future, we find that a majority of models underestimate the extremeness of impacts in important sectors such as agriculture, terrestrial ecosystems, and heat-related human mortality, while impacts on water resources and hydropower are overestimated in some river basins; and the spread across models is often large. This has important implications for economic assessments of climate change impacts that rely on these models. It also means that societal risks from future extreme events may be greater than previously thought.
Aquatic ecosystems are globally contaminated with microplastics (MP). However, comparative data on MP levels in freshwater systems is still scarce. Therefore, the aim of this study is to quantify MP abundance in water and sediment of the German river Elbe using visual, spectroscopic (Fourier-transform infrared spectroscopy) and thermo analytical (pyrolysis gas chromatography mass spectrometry) methods. Samples from eleven German sites along the German part of the Elbe were collected, both in the water and sediment phase, in order to better understand MP sinks and transport mechanisms. MP concentrations differed between the water and sediment phase. Sediment concentrations (mean: 3,350,000 particles m−3, 125–5000 μm MP) were in average 600,000-fold higher than water concentrations (mean: 5.57 particles m−3, 150–5000 μm MP). The abundance varied between the sampling sites: In sediments, the abundance decreased in the course of the river while in water samples no such clear trend was observed. This may be explained by a barrage retaining sediments and limiting tidal influence in the upstream parts of the river. Particle shape differed site-specifically with one site having exceptionally high quantities of spheres, most probably due to industrial emissions of PS-DVB resin beads. Suspended MP consisted predominantly of polyethylene and polypropylene whereas sediments contained a higher diversity of polymer types. Determined MP concentrations correspond well to previous results from other European rivers. In a global context, MP levels in the Elbe relate to the lower (water) to middle section (sediment) of the global range of MP concentrations determined for rivers worldwide. This highlights that elevated MP levels are not only found in single countries or continents, but that MP pollution is an issue of global concern.
Aus der Redaktion
(2020)
Droughts are anticipated to intensify in many parts of the world due to climate change. However, the issue of drought definition, namely the diversity of drought indices, makes it difficult to compare drought assessments. This issue is widely known, but its relative importance has never been quantitatively evaluated in comparison to other sources of uncertainty. Here, encompassing three drought categories (meteorological, agricultural, and hydrological droughts) with four temporal scales of interest, we evaluated changes in the drought frequency using multi-model and multi-scenario simulations to identify areas where the definition issue could result in pronounced uncertainties and to what extent. We investigated the disagreement in the signs of changes between drought definitions and decomposed the variance into four main factors: drought definitions, greenhouse gas concentration scenarios, global climate models, and global water models, as well as their interactions. The results show that models were the primary sources of variance over 82% of the global land area. On the other hand, the drought definition was the dominant source of variance in the remaining 17%, especially in parts of northern high-latitudes. Our results highlight specific regions where differences in drought definitions result in a large spread among projections, including areas showing opposite signs of significant changes. At a global scale, 7% of the variance resulted independently from the definition issue, and that value increased to 44% when 1st and 2nd order interactions were considered. The quantitative results suggest that by clarifying hydrological processes or sectors of interest, one could avoid these uncertainties in drought assessments to obtain a clearer picture of future drought change.
Le phénomène du cuirassement reste une "curiosité" et une énigme pour le pays de la zone intertropicale. A cause de ses caractéristiques lithologiques et structurales assez particulières de ses roches (roches riches en éléments ferromagnésiens), le Burkina Faso est une véritable zone de prédilection des cuirasses. Malgré l'effort des chercheurs pour élucider le phénomène du cuirassement, force est de reconnaître que de nos jours, certains points d'ombre subsistent toujours; ce qui invite à une analyse plus poussée ... Loin de négliger les problèmes liés à la reconstitution de la génèse de la cuirasse, nous proposons ici une analyse assez originale des cuirasses sur la base des connaissances déjà acquises et de nos multiples observations sur le terrain burkinabé.
Land surface and hydrologic models (LSMs/HMs) are used at diverse spatial resolutions ranging from catchment-scale (1–10 km) to global-scale (over 50 km) applications. Applying the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the chosen resolution, i.e., fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent hydrologic parameter fields. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB, and WaterGAP models demonstrate the pitfalls of deficient parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge's 1982 statement on the unsolved problem of parameterization in these models remains true. Based on a review of existing parameter regionalization techniques, we postulate that the multiscale parameter regionalization (MPR) technique offers a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. Herein, we develop a general model protocol to describe how MPR can be applied to a particular model and present an example application using the PCR-GLOBWB model. Finally, we discuss potential advantages and limitations of MPR in obtaining the seamless prediction of hydrological fluxes and states across spatial scales.
Land surface and hydrologic models (LSM/HM) are used at diverse spatial resolutions ranging from 1-10 km in catchment-scale applications to over 50 km in global-scale applications. Application of the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the model resolution and fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent and realistic parameter fields for land surface geophysical properties. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB and WaterGAP models are conducted to demonstrate the pitfalls of poor parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge’s 1982 statement on the unsolved problem of parameterization in these models remains true. We provide a short review of existing parameter regionalization techniques and discuss a method for obtaining seamless hydrological predictions of water fluxes and states across multiple spatial resolutions. The multiscale parameter regionalization (MPR) technique is a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. A general model protocol is presented to describe how MPR can be applied to a specific model, with an example of this application using the PCR-GLOBWB model. Applying MPR to PCR-GLOBWB substantially improves the flux-matching condition. Estimation of evapotranspiration without MPR at 5 arcmin and 30 arcmin spatial resolutions for the Rhine river basin results in a difference of approximately 29%. Applying MPR reduce this difference to 9%. For total soil water, the differences without and with MPR are 25% and 7%, respectively.
During the SHIVA (Stratospheric Ozone: Halogen Impacts in a Varying Atmosphere) project an extensive dataset of all halogen species relevant for the atmospheric budget of total organic bromine has been collected in the West Pacific region using the FALCON aircraft operated by the German Aerospace agency DLR (Deutsches Zentrum für Luft- und Raumfahrt) covering a vertical range from the planetary boundary layer up to the ceiling altitude of the aircraft of 13 km. In total, more than 700 measurements were performed with the newly developed fully-automated in-situ instrument GHOST-MS (Gas cHromatograph for the Observation of Tracers – coupled with a Mass Spectrometer) by the Goethe University of Frankfurt (GUF) and with the onboard whole-air sampler WASP with subsequent ground based state-of-the-art GC/MS analysis by the University of East Anglia (UEA). Both instruments yield good agreement for all major (CHBr3 and CH2Br2) and minor (CHBrCl, CHBrCl2 and CHBr2Cl) VSLS (very short-lived substances), at least at the level of their 2 σ measurement uncertainties. In contrast to the suggestion that the Western Pacific could be a major source region for VSLS (Pyle et al., 2011), we found only slightly enhanced mixing ratios of brominated halogen source gases relative to the levels reported in Montzka et al. (2011) for other tropical regions. A budget for total organic bromine, including all four halons,CH3Br and the VSLS, is derived for the upper troposphere, the input region for the TTL and thus also for the stratosphere, compiled from the SHIVA dataset. With exception of the two minor VSLS CHBrCl2 and CHBr2Cl, excellent agreement with the values reported in Montzka et al. (2011) is found, while being slightly higher than previous studies from our group based on balloon-borne measurements.
Latitudinal and bathymetrical species richness patterns in the NW Pacific and adjacent Arctic Ocean
(2019)
Global scale analyses have recently revealed that the latitudinal gradient in marine species richness is bimodal, peaking at low-mid latitudes but with a dip at the equator; and that marine species richness decreases with depth in many taxa. However, these overall and independently studied patterns may conceal regional differences that help support or qualify the causes in these gradients. Here, we analysed both latitudinal and depth gradients of species richness in the NW Pacific and its adjacent Arctic Ocean. We analysed 324,916 distribution records of 17,414 species from 0 to 10,900 m depth, latitude 0 to 90°N, and longitude 100 to 180°N. Species richness per c. 50 000 km2 hexagonal cells was calculated as alpha (local average), gamma (regional total) and ES50 (estimated species for 50 records) per latitudinal band and depth interval. We found that average ES50 and gamma species richness decreased per 5° latitudinal bands and 100 m depth intervals. However, average ES50 per hexagon showed that the highest species richness peaked around depth 2,000 m where the highest total number of species recorded. Most (83%) species occurred in shallow depths (0 to 500 m). The area around Bohol Island in the Philippines had the highest alpha species richness (more than 8,000 species per 50,000 km2). Both alpha and gamma diversity trends increased from the equator to latitude 10°N, then further decreased, but reached another peak at higher latitudes. The latitudes 60–70°N had the lowest gamma and alpha diversity where there is almost no ocean in our study area. Model selection on Generalized Additive Models (GAMs) showed that the combined effects of all environmental predictors produced the best model driving species richness in both shallow and deep sea. The results thus support recent hypotheses that biodiversity, while highest in the tropics and coastal depths, is decreasing at the equator and decreases with depth below ~2000 m. While we do find the declines of species richness with latitude and depth that reflect temperature gradients, local scale richness proved poorly correlated with many environmental variables. This demonstrates that while regional scale patterns in species richness may be related to temperature, that local scale richness depends on a greater variety of variables.
We report measurements of the deuterium content of molecular hydrogen (H2) obtained from a suite of air samples that were collected during a stratospheric balloon flight between 12 and 33 km at 40º N in October 2002. Strong deuterium enrichments of up to 400 permil versus Vienna Standard Mean Ocean Water (VSMOW) are observed, while the H2 mixing ratio remains virtually constant. Thus, as hydrogen is processed through the H2 reservoir in the stratosphere, deuterium is accumulated in H2 . Using box model calculations we investigated the effects of H2 sources and sinks on the stratospheric enrichments. Results show that considerable isotope enrichments in the production of H2 from CH4 must take place, i.e., deuterium is transferred preferentially to H2 during the CH4 oxidation sequence. This supports recent conclusions from tropospheric H2 isotope measurements which show that H2 produced photochemically from CH4 and non-methane hydrocarbons must be enriched in deuterium to balance the tropospheric hydrogen isotope budget. In the absence of further data on isotope fractionations in the individual reaction steps of the CH4 oxidation sequence, this effect cannot be investigated further at present. Our measurements imply that molecular hydrogen has to be taken into account when the hydrogen isotope budget in the stratosphere is investigated.
The isotopic composition of methane in the stratosphere : high-altitude balloon sample measurements
(2011)
The isotopic composition of stratospheric methane has been determined on a large suite of air samples from stratospheric balloon flights covering subtropical to polar latitudes and a time period of 16 yr. 154 samples were analyzed for δC and 119 samples for δD, increasing the previously published dataset for balloon borne samples by an order of magnitude, and more than doubling the total available stratospheric data (including aircraft samples) published to date. The samples also cover a large range in mixing ratio from tropospheric values near 1800 ppb down to only 250 ppb, and the strong isotope fractionation processes accordingly increase the isotopic composition up to δ13C=−14‰ and δD= +190‰, the largest enrichments observed for atmospheric CH4 so far. When analyzing and comparing kinetic isotope effects (KIEs) derived from single balloon profiles, it is necessary to take into account the residence time in the stratosphere in combination with the observed mixing ratio and isotope trends in the troposphere, and the range of isotope values covered by the individual profile. Temporal isotope trends can also be determined in the stratosphere and compare reasonably well with the tropospheric trends. The effects of chemical and dynamical processes on the isotopic composition of CH4 in the stratosphere are discussed in detail. Different ways to interpret the data in terms of the relative fractions of the three important sink mechanisms (reaction with OH, O(1D)) and Cl, respectively), and their limitations, are investigated. The classical approach of using global mean KIE values can be strongly biased when profiles with different minimum mixing ratios are compared. Approaches for more local KIE investigations are suggested. It is shown that any approach for a formal sink partitioning from the measured data severely underestimates the fraction removed by OH, which is likely due to the insensitivity of the measurements to the kinetic fractionation in the lower stratosphere. Attempts can be made to correct for the lower stratospheric sink bias, but full quantitative interpretation of the CH4 isotope data in terms of the three sink reactions requires a global model.
The isotopic composition of methane in the stratosphere : high-altitude balloon sample measurements
(2011)
The isotopic composition of stratospheric methane has been determined on a large suite of air samples from stratospheric balloon flights covering subtropical to polar latitudes and a time period of 16 yr. 154 samples were analyzed for δ13C and 119 samples for δD, increasing the previously published dataset for balloon borne samples by an order of magnitude, and more than doubling the total available stratospheric data (including aircraft samples) published to date. The samples also cover a large range in mixing ratio from tropospheric values near 1800 ppb down to only 250 ppb, and the strong isotope fractionation processes accordingly increase the isotopic composition up to δ13C = −14‰ and δD = +190‰, the largest enrichments observed for atmospheric CH4 so far. When analyzing and comparing kinetic isotope effects (KIEs) derived from single balloon profiles, it is necessary to take into account the residence time in the stratosphere in combination with the observed mixing ratio and isotope trends in the troposphere, and the range of isotope values covered by the individual profile. The isotopic composition of CH4 in the stratosphere is affected by both chemical and dynamical processes. This severely hampers interpretation of the data in terms of the relative fractions of the three important sink mechanisms (reaction with OH, O(1D) and Cl). It is shown that a formal sink partitioning using the measured data severely underestimates the fraction removed by OH, which is likely due to the insensitivity of the measurements to the kinetic fractionation in the lower stratosphere. Full quantitative interpretation of the CH4 isotope data in terms of the three sink reactions requires a global model.
The Tarim River Basin, located in Xinjiang, NW China, is the largest endorheic river basin of China and one of the largest in whole Central Asia. Due to the extremely arid climate with an annual precipitation of less than 100 mm, the water supply along the Aksu and Tarim River solely depends on river water. This applies for anthropogenic activities (e.g. agriculture) as well as for the natural ecosystems so that both compete for water. The on-going increase of water consumption by agriculture and other human activities in this region has been enhancing the competition for water between human needs and nature. Against this background, 11 German and 6 Chinese universities and research institutes formed the consortium SuMaRiO (www.sumario.de), which aims at gaining a holistic picture of the availability of water resources in the Tarim River Basin and the impacts on anthropogenic activities and natural ecosystems caused by the water distribution within the Tarim River Basin. The discharge of the Aksu River, which is the major tributary to the Tarim, has been increasing over the past 6 decades due to enhanced glacier melt. Alone from 1989 to 2011, the area under agriculture more than doubled. Thereby, cotton became the major crop and there was a shift from small-scale farming to large-scale intensive farming. The major natural ecosystems along the Aksu and Tarim River are riparian ecosystems: Riparian (Tugai) forests, shrub vegetation, reed beds, and other grassland. Within the SuMaRiO Cluster the focus was laid on the Tugai forests, with Populus euphratica as dominant tree, because the most productive and species-rich natural ecosystems can be found among those forests. On sites with groundwater distance of less than 7.5 m the annual increments correlated with river runoffs of the previous year. But, the further downstream along the Tarim River, the more the natural river dynamics ceased, which impacts on the recruitment of Populus euphratica. Household surveys revealed that there is a considerable willingness to pay for conservation of those riparian forests with the mitigation of dust and sandstorms considered as the most important ecosystem service. This interdisciplinary project will result in a decision support tool (DST), build on the participation of regional stakeholders and models based on results and field experiments. This DST finally shall assist stakeholders in balancing the water competition acknowledging the major external effects of any water allocation.
The Tarim River basin, located in Xinjiang, NW China, is the largest endorheic river basin in China and one of the largest in all of Central Asia. Due to the extremely arid climate, with an annual precipitation of less than 100 mm, the water supply along the Aksu and Tarim rivers solely depends on river water. This is linked to anthropogenic activities (e.g., agriculture) and natural and semi-natural ecosystems as both compete for water. The ongoing increase in water consumption by agriculture and other human activities in this region has been enhancing the competition for water between human needs and nature. Against this background, 11 German and 6 Chinese universities and research institutes have formed the consortium SuMaRiO (Sustainable Management of River Oases along the Tarim River; http://www.sumario.de), which aims to create a holistic picture of the availability of water resources in the Tarim River basin and the impacts on anthropogenic activities and natural ecosystems caused by the water distribution within the Tarim River basin. On the basis of the results from field studies and modeling approaches as well as from suggestions by the relevant regional stakeholders, a decision support tool (DST) will be implemented that will then assist stakeholders in balancing the competition for water, acknowledging the major external effects of water allocation to agriculture and to natural ecosystems. This consortium was formed in 2011 and is funded by the German Federal Ministry of Education and Research. As the data collection phase was finished this year, the paper presented here brings together the results from the fields from the disciplines of climate modeling, cryology, hydrology, agricultural sciences, ecology, geoinformatics, and social sciences in order to present a comprehensive picture of the effects of different water availability schemes on anthropogenic activities and natural ecosystems along the Tarim River. The second objective is to present the project structure of the whole consortium, the current status of work (i.e., major new results and findings), explain the foundation of the decision support tool as a key product of this project, and conclude with application recommendations for the region. The discharge of the Aksu River, which is the major tributary of the Tarim, has been increasing over the past 6 decades. From 1989 to 2011, agricultural area more than doubled: cotton became the major crop and there was a shift from small-scale to large-scale intensive farming. The ongoing increase in irrigated agricultural land leads to the increased threat of salinization and soil degradation caused by increased evapotranspiration. Aside from agricultural land, the major natural and semi-natural ecosystems are riparian (Tugai) forests, shrub vegetation, reed beds, and other grassland, as well as urban and peri-urban vegetation. Within the SuMaRiO cluster, focus has been set on the Tugai forests, with Populus euphratica as the dominant tree species, because these forests belong to the most productive and species-rich natural ecosystems of the Tarim River basin. At sites close to the groundwater, the annual stem diameter increments of Populus euphratica correlated with the river runoffs of the previous year. However, the natural river dynamics cease along the downstream course and thus hamper the recruitment of Populus euphratica. A study on the willingness to pay for the conservation of the natural ecosystems was conducted to estimate the concern of the people in the region and in China's capital. These household surveys revealed that there is a considerable willingness to pay for conservation of the natural ecosystems, with mitigation of dust and sandstorms considered the most important ecosystem service. Stakeholder dialogues contributed to creating a scientific basis for a sustainable management in the future.
Tectonics and geothermal gradients from subduction to collision in the NW Variscan Iberian Massif
(2022)
The earliest tectonometamorphic record of tectonic slices incorporated to the base of an orogen holds the key to understand how an orogen is built. The tectonic pile of the NW Iberian section of the Variscan Orogen includes tectonic slices separated by crustal-scale thrusts. The earliest tectonometamorphic record in the uppermost parautochthon is calculated at 11–14 kbar and 450–500°C (P-T gradient about 13°C/km), suggesting a subduction-related metamorphic recrystallization at lower pressure than the overlying Lower Allochthon. Early conditions calculated in the autochthon (9–10 kbar and 425–450°C; 16°C/km) point to a relatively ‘cold’ collisional setting. Higher thermal gradients obtained from some sections of the autochthon (11–12 kbar and 700–725°C; 21°C/km) and the Lower Parautochthon (7.5 kbar and 550–700°C; 24–31°C/km), correspond to more advanced and ‘hot’ stages of collision. New U–Pb monazite geochronology indicates a 318–311 Ma age for the final formation of HT domes in the region. We propose the rapid decrease in P-T gradient (from <10 to 16°C/km) documents a fail to sustain further burial along a regular subduction zone. We consider the subsequent increase in the geothermal gradient (from 16 to 31°C/km) as the culmination of previous crustal accretion and the onset of crustal underthrusting and later processes in a collisional stage. We propose these switches in the early tectonometamorphic record of individual tectonic slices as potential markers to track the transition from subduction to collision in collisional orogens.
As part of the CLACE-6 campaign we performed size-resolved CCN measurements for a~supersaturation range of S = 0.079 % to 0.66% at the high-alpine research station Jungfraujoch, Switzerland, in March~2007. The derived effective hygroscopicity parameter κ describing the influence of particle composition on CCN activity was on average 0.23–0.30 for Aitken (50–100 nm) and 0.32–0.43 for accumulation mode particles (100–200 nm). The campaign average value of κ = 0.3 is similar to the average value of κ for other continental locations. When air masses came from southeasterly directions crossing the Po Valley in Italy, particles were much more hygroscopic (κ ≈ 0.42) due to large sulfate mass fractions. The κ values obtained at S = 0.079 % exhibited a good negative correlation with the organic mass fractions derived from PM1 aerosol mass spectrometer (AMS) measurements. Applying a simple mixing rule the organic and inorganic mass fractions observed by the AMS could be used to reproduce the temporal fluctuations of the hygroscopicity of accumulation mode particles quite well.
We show how during a cloud event the aerosol particles were activated as cloud droplets and then removed from the air by precipitation leaving behind only a small amount of accumulation mode particles consisting mainly of weakly CCN-active particles, most likely externally mixed unprocessed soot particles.
During the campaign we had the opportunity to directly compare two DMT CCN counters for a certain time. The total CCN concentration (NCCN,tot) obtained by the two instruments at equal supersaturations agreed well for both possible operating modes: detecting NCCN,tot directly by sampling the polydisperse aerosol with the CCNC, or indirectly by combining size-resolved measurements of the activated fraction with parallel measurements of the particle size distribution (e.g., by SMPS). However, some supersaturation setpoints differed between the two CCNCs by as much as 20% after applying the instrument calibrations, which resulted in differences of the corresponding NCCN,tot of up to 50%. This emphasizes that it is extremely important to carefully calibrate the supersaturation of the instrument, especially at low S.
Ternary aerosol nucleation experiments were conducted in the CLOUD chamber at CERN in order to investigate the influence of ions on new particle formation. Neutral and ion-induced nucleation experiments, i.e., with and without the presence of ions, were carried out under precisely controlled conditions. The sulphuric acid concentration was measured with a Chemical Ionization Mass Spectrometer (CIMS) during the new particle formation experiments. The added ternary trace gases were ammonia (NH3), dimethylamine (DMA, C2H7N) or oxidised products of pinanediol (PD, C10H18O2). When pinanediol was introduced into the chamber, an increase in the mass spectrometric signal used to determine the sulphuric acid concentration (m/z 97, i.e., HSO4−) was observed due to ions from the CLOUD chamber. The enhancement was only observed during ion-induced nucleation measurements by using either galactic cosmic rays (GCR) or the proton synchrotron (PS) pion beam for the ion generation, respectively. The ion effect typically involved an increase in the apparent sulphuric acid concentration by a factor of ~2 to 3 and was qualitatively verified by the ion measurements by an Atmospheric Pressure interface-Time Of Flight (APi-TOF) mass spectrometer. By applying a high voltage (HV) clearing field inside the CLOUD chamber the ion effect on the CIMS measurement was completely eliminated since, under these conditions, small ions are swept from the chamber in about one second. In order to exclude the ion effect and to provide corrected sulphuric acid concentrations during the GCR and PS beam nucleation experiments, a parameterisation was derived that utilizes the trace gas concentrations and the UV light intensity as input parameters. Atmospheric sulphuric acid measurements with a CIMS showed an insignificant ion effect.
Ternary aerosol nucleation experiments were conducted in the CLOUD chamber at CERN in order to investigate the influence of ions on new particle formation. Neutral and ion-induced nucleation experiments, i.e. without and with the presence of ions, respectively, were carried out under precisely controlled conditions. The sulfuric acid concentration was measured with a chemical ionisation mass spectrometer (CIMS) during the new particle formation experiments. The added ternary trace gases were ammonia (NH3), dimethylamine (DMA, C2H7N) or oxidised products of pinanediol (PD, C10H18O2). When pinanediol was introduced into the chamber, an increase in the mass spectrometric signal used to determine the sulfuric acid concentration (m/z 97, i.e. HSO4−) was observed due to ions from the CLOUD chamber. The enhancement was only observed during ion-induced nucleation measurements by using either galactic cosmic rays (GCRs) or the proton synchrotron (PS) pion beam for the ion generation, respectively. The ion effect typically involved an increase in the apparent sulfuric acid concentration by a factor of ~ 2 to 3 and was qualitatively verified by the ion measurements with an atmospheric-pressure interface-time of flight (APi-TOF) mass spectrometer. By applying a high-voltage (HV) clearing field inside the CLOUD chamber, the ion effect on the CIMS measurement was completely eliminated since, under these conditions, small ions are swept from the chamber in about 1 s. In order to exclude the ion effect and to provide corrected sulfuric acid concentrations during the GCR and PS beam nucleation experiments, a parameterisation was derived that utilises the trace gas concentrations and the UV light intensity as input parameters. Atmospheric sulfuric acid measurements with a CIMS showed an insignificant ion effect.
Balloon-borne measurements of CFC-11 (on flights of the DIRAC in situ gas chromatograph and the DESCARTES grab sampler), ClO and O3 were made during the 1999/2000 winter as part of the SOLVE-THESEO 2000 campaign. Here we present the CFC-11 data from nine flights and compare them first with data from other instruments which flew during the campaign and then with the vertical distributions calculated by the SLIMCAT 3-D CTM. We calculate ozone loss inside the Arctic vortex between late January and early March using the relation between CFC-11 and O3 measured on the flights, the peak ozone loss (1200 ppbv) occurs in the 440–470 K region in early March in reasonable agreement with other published empirical estimates. There is also a good agreement between ozone losses derived from three independent balloon tracer data sets used here. The magnitude and vertical distribution of the loss derived from the measurements is in good agreement with the loss calculated from SLIMCAT over Kiruna for the same days.
Balloon-borne measurements of CFC11 (from the DIRAC in situ gas chromatograph and the DESCARTES grab sampler), ClO and O3 were made during the 1999/2000 Arctic winter as part of the SOLVE-THESEO 2000 campaign, based in Kiruna (Sweden). Here we present the CFC11 data from nine flights and compare them first with data from other instruments which flew during the campaign and then with the vertical distributions calculated by the SLIMCAT 3D CTM. We calculate ozone loss inside the Arctic vortex between late January and early March using the relation between CFC11 and O3 measured on the flights. The peak ozone loss (~1200ppbv) occurs in the 440-470K region in early March in reasonable agreement with other published empirical estimates. There is also a good agreement between ozone losses derived from three balloon tracer data sets used here. The magnitude and vertical distribution of the loss derived from the measurements is in good agreement with the loss calculated from SLIMCAT over Kiruna for the same days.
Dating of extensive alluvial fan surfaces and fluvial features in the hyperarid core of the Atacama Desert, Chile, using cosmogenic nuclides provides unrivalled insights about the onset and variability of aridity. The predominantly hyperarid conditions help to preserve the traces of episodic climatic and/or slow tectonic change. Utilizing single clast exposure dating with cosmogenic 10Be and 21Ne, we determine the termination of episodes of enhanced fluvial erosion and deposition occurring at ~19, ~14, ~9.5 Ma; large scale fluvial modification of the landscape had ceased by ~2–3 Ma. The presence of clasts that record pre-Miocene exposure ages (~28 Ma and ~34 Ma) require stagnant landscape development during the Oligocene. Our data implies an early onset of (hyper-) aridity in the core region of the Atacama Desert, interrupted by wetter but probably still arid periods. The apparent conflict with interpretation that favour a later onset of (hyper-) aridity can be reconciled when the climatic gradients within the Atacama Desert are considered.
Seasonal forecasting systems still have difficulties predicting temperature over continental regions, while their performance is better over some maritime regions. On the other hand, the land surface is a substantial source of (sub-)seasonal predictability. A crucial land surface component in focus here is the snow cover, which stores water and modulates the surface radiation balance. This paper’s goal is to attribute snow cover seasonal forecasting biases and lack of skill to either initialization or parameterization errors. For this purpose, we compare the snow representation in five seasonal forecasting systems (from DWD, ECMWF, Météo-France, CMCC, and ECCC) and their performances in predicting snow and 2-m temperature over a Siberian region against ERA5 reanalysis and station data. Although all systems use similar atmospheric and land initialization approaches and data, their snow and temperature biases differ in sign and amplitude. Too-large initial snow biases persist over the forecast period, delaying and prolonging the melting phase. The simplest snow scheme (used in DWD’s system) shows too-early and fast melting in spring. However, systems including multi-layer snow schemes (Météo-France and CMCC) do not necessarily perform better. Both initialization and parameterization are causes of snow biases, but, depending on the system, one can be more dominant.
Background: Replicate population pairs that diverge in response to similar selective regimes allow for an investigation of (a) whether phenotypic traits diverge in a similar and predictable fashion, (b) whether there is gradual variation in phenotypic divergence reflecting variation in the strength of natural selection among populations, (c) whether the extent of this divergence is correlated between multiple character suites (i.e., concerted evolution), and (d) whether gradual variation in phenotypic divergence predicts the degree of reproductive isolation, pointing towards a role for adaptation as a driver of (ecological) speciation. Here, we use poeciliid fishes of the genera Gambusia and Poecilia that have repeatedly evolved extremophile lineages able to tolerate high and sustained levels of toxic hydrogen sulfide (H2S) to answer these questions.
Results: We investigated evolutionary divergence in response to H2S in Gambusia spp. (and to a lesser extent Poecilia spp.) using a multivariate approach considering the interplay of life history, body shape, and population genetics (nuclear miscrosatellites to infer population genetic differentiation as a proxy for reproductive isolation). We uncovered both shared and unique patterns of evolution: most extremophile Gambusia predictably evolved larger heads and offspring size, matching a priori predictions for adaptation to sulfidic waters, while variation in adult life histories was idiosyncratic. When investigating patterns for both genera (Gambusia and Poecilia), we found that divergence in offspring-related life histories and body shape were positively correlated across populations, but evidence for individual-level associations between the two character suites was limited, suggesting that genetic linkage, developmental interdependencies, or pleiotropic effects do not explain patterns of concerted evolution. We further found that phenotypic divergence was positively correlated with both environmental H2S-concentration and neutral genetic differentiation (a proxy for gene flow).
Conclusions: Our results suggest that higher toxicity exerts stronger selection, and that divergent selection appears to constrain gene flow, supporting a scenario of ecological speciation. Nonetheless, progress toward ecological speciation was variable, partially reflecting variation in the strength of divergent selection, highlighting the complexity of selective regimes even in natural systems that are seemingly governed by a single, strong selective agent.
The complexity of atmospherical processes has always yielded a multitude of ways of knowing about the weather. What has been lacking in the historiography of meteorology so far is a way to formulate differences between forms of knowledge in a way that does not privilege modern scientific structures, but focuses instead on the epistemological category of causality. Using causality as ground of comparison for different knowledge claims, I shall argue, may enable researchers to investigate meteorological knowledge across time periods, perhaps even geographical regions, in a more symmetrical manner. This review demonstrates this approach as a means to organize a large set of historical meteorological writings from German countries between 1750 and 1850. Three distinct forms of knowledge (Semiotics, Physics, and Organics of the weather) during that time and in that region are suggested and will be described. While a bibliography with a national perspective from the 1880s was the basis for the selection of historical sources, such a setup proved awkward even to contemporaries. In addition, the bibliography came with a number of biases and shortcomings that will be critically reviewed.
Contribution of sulfuric acid and oxidized organic compounds to particle formation and growth
(2012)
Lack of knowledge about the mechanisms underlying new particle formation and their subsequent growth is one of the main causes for the large uncertainty in estimating the radiative forcing of atmospheric aerosols in global models. We performed chamber experiments designed to study the contributions of sulfuric acid and organic vapors to formation and to the early growth of nucleated particles, respectively. Distinct experiments in the presence of two different organic precursors (1,3,5-trimethylbenzene and α-pinene) showed the ability of these compounds to reproduce the formation rates observed in the low troposphere. These results were obtained measuring the sulfuric acid concentrations with two Chemical Ionization Mass Spectrometers confirming the results of a previous study which modeled the sulfuric acid concentrations in presence of 1,3,5-trimethylbenzene.
New analysis methods were applied to the data collected with a Condensation Particle Counter battery and a Scanning Mobility Particle Sizer, allowing the assessment of the size resolved growth rates of freshly nucleated particles. The effect of organic vapors on particle growth was investigated by means of the growth rate enhancement factor (Γ), defined as the ratio between the measured growth rate in the presence of α-pinene and the kinetically limited growth rate of the sulfuric acid and water system. The observed Γ values indicate that the growth is dominated by organic compounds already at particle diameters of 2 nm. Both the absolute growth rates and Γ showed a strong dependence on particle size supporting the nano-Köhler theory. Moreover, the separation of the contributions from sulfuric acid and organic compounds to particles growth reveals that the organic contribution seems to be enhanced by the sulfuric acid concentration. The size resolved growth analysis finally indicates that both condensation of oxidized organic compounds and reactive uptake contribute to particle growth.
Contribution of sulfuric acid and oxidized organic compounds to particle formation and growth
(2012)
Lack of knowledge about the mechanisms underlying new particle formation and their subsequent growth is one of the main causes for the large uncertainty in estimating the radiative forcing of atmospheric aerosols in global models. We performed chamber experiments designed to study the contributions of sulfuric acid and organic vapors to the formation and early growth of nucleated particles. Distinct experiments in the presence of two different organic precursors (1,3,5-trimethylbenzene and α-pinene) showed the ability of these compounds to reproduce the formation rates observed in the low troposphere. These results were obtained measuring the sulfuric acid concentrations with two chemical ionization mass spectrometers confirming the results of a previous study which modeled the sulfuric acid concentrations in presence of 1,3,5-trimethylbenzene.
New analysis methods were applied to the data collected with a condensation particle counter battery and a scanning mobility particle sizer, allowing the assessment of the size resolved growth rates of freshly nucleated particles. The effect of organic vapors on particle growth was investigated by means of the growth rate enhancement factor (Γ), defined as the ratio between the measured growth rate in the presence of α-pinene and the kinetically limited growth rate of the sulfuric acid and water system. The observed Γ values indicate that the growth is already dominated by organic compounds at particle diameters of 2 nm. Both the absolute growth rates and Γ showed a strong dependence on particle size, supporting the nano-Köhler theory. Moreover, the separation of the contributions from sulfuric acid and organic compounds to particle growth reveals that the organic contribution seems to be enhanced by the sulfuric acid concentration. Finally, the size resolved growth analysis indicates that both condensation of oxidized organic compounds and reactive uptake contribute to particle growth.
Unter dem Begriff „Stillgewässer“ wird eine fast nicht überschaubare Vielzahl von
stehenden Gewässern zusammengefaßt, die das gesamte Spektrum von der kurzlebigen
Regenwasserpfütze bis hin zu großen Binnenseen wie Bodensee oder Müritz umfassen. Bei
dieser Spannbreite wundert es nicht, daß Stillgewässer Lebensräume mit sehr
unterschiedlichen Eigenschaften sind und dementsprechend auch sehr unterschiedliche
Lebensgemeinschaften beherbergen.
Constraining the architecture of complex 3D volcanic plumbing systems within active rifts, and their impact on rift processes, is critical for examining the interplay between faulting, magmatism and magmatic fluids in developing rift segments. The Natron basin of the East African Rift System provides an ideal location to study these processes, owing to its recent magmatic-tectonic activity and ongoing active carbonatite volcanism at Oldoinyo Lengai. Here, we report seismicity and fault plane solutions from a 10-month temporary seismic network spanning Oldoinyo Lengai, Naibor Soito volcanic field and Gelai volcano. We locate 6827 earthquakes with ML -0.85 to 3.6, which are related to previous and ongoing magmatic and volcanic activity in the region, as well as regional tectonic extension. We observe seismicity down to ~17 km depth north and south of Oldoinyo Lengai and shallow seismicity (3 - 10 km) beneath Gelai, including two swarms. The deepest seismicity (~down to 20 km) occurs above a previously imaged magma body below Naibor Soito. These seismicity patterns reveal a detailed image of a complex volcanic plumbing system, supporting potential lateral and vertical connections between shallow- and deep-seated magmas, where fluid and melt transport to the surface is facilitated by intrusion of dikes and sills. Focal mechanisms vary spatially. T-axis trends reveal dominantly WNW-ESE extension near Gelai, while strike-slip mechanisms and a radial trend in P-axes are observed in the vicinity of Oldoinyo Lengai. These data support local variations in the state of stress, resulting from a combination of volcanic edifice loading and magma-driven stress changes imposed on a regional extensional stress field. Our results indicate that the southern Natron basin is a segmented rift system, in which fluids preferentially percolate vertically and laterally in a region where strain transfers from a border fault to a developing magmatic rift segment.
Highlights
• We show the first observations of seismo-acoustic tremor at Oldoinyo Lengai, the world's only active carbonatite volcano.
• We observe significant changes in seismic and acoustic tremor properties and their correlation in one year of data collection.
• Using satellite-based thermal data, we identify different volcanic processes (degassing, lava pond dynamics and spattering).
Abstract
We analyze volcanic tremor from Oldoinyo Lengai, Tanzania, which is currently the only active volcano on Earth producing carbonatitic lavas. Here, we use data from the recent SEISVOL deployment and focus on a co-located seismic and infrasound station about 200 m below the summit. We show the very first observations of seismo-acoustic tremor caused by carbonatitic eruptions. This seismo-acoustic tremor is highly variable throughout the ∼one year of data which we characterize by analyzing its seismic amplitude, duration, recurrence, dominant seismic frequency and harmonics. Frequency gliding occurs frequently and over short (minutes to hours) to long time scales (hours to days) and likely reflects different time-dependent mechanisms, such as evenly-spaced repeating events with a change in inter-event times, crater dynamics that alter resonators, and dike intrusions. Seismic and acoustic wavefields correlate well for stronger eruptive sequences but are only partially coherent which suggests that high-frequency seismic tremor (up to 25 Hz) may be caused by the low viscosity of the carbonatitic melt and not by ground-coupled airwaves. In addition, the comparison between seismic-acoustic and satellite InfraRed thermal data allows us to infer different volcanic activity styles which partially alternate throughout the year: intrusive activity and the construction of hornitos, degassing, activity from a lava pond, and varying styles of extrusive activity, in particular spattering. Our study provides important insights into the eruption dynamics of this peculiar volcano which suggests shallow melt storage within the crater floor.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
In global hydrological models, groundwater (GW) is typically represented by a bucket-like linear groundwater reservoir. Reservoir models, however, (1) can only simulate GW discharge to surface water (SW) bodies but not recharge from SW to GW, (2) provide no information on the location of the GW table, and (3) assume that there is no GW flow among grid cells. This may lead, for example, to an underestimation of groundwater resources in semiarid areas where GW is often replenished by SW or to an underestimation of evapotranspiration where the GW table is close to the land surface. To overcome these limitations, it is necessary to replace the reservoir model in global hydrological models with a hydraulic head gradient-based GW flow model.
We present G3M, a new global gradient-based GW model with a spatial resolution of 5′ (arcminutes), which is to be integrated into the 0.5∘ WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables in-memory coupling to WGHM while keeping overall runtime relatively low, which allows sensitivity analyses, calibration, and data assimilation. This paper presents the G3M concept and model design decisions that are specific to the large grid size required for a global-scale model. Model results under steady-state naturalized conditions, i.e., neglecting GW abstractions, are shown. Simulated hydraulic heads show better agreement to observations around the world compared to the model output of de Graaf et al. (2015). Locations of simulated SW recharge to GW are found, as is expected, in dry and mountainous regions but areal extent of SW recharge may be underestimated. Globally, GW discharge to rivers is by far the dominant flow component such that lateral GW flows only become a large fraction of total diffuse and focused recharge in the case of losing rivers, some mountainous areas, and some areas with very low GW recharge. A strong sensitivity of simulated hydraulic heads to the spatial resolution of the model and the related choice of the water table elevation of surface water bodies was found. We suggest to investigate how global-scale groundwater modeling at 5′ spatial resolution can benefit from more highly resolved land surface elevation data.
To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of Abstract. To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of capillary rise on evapotranspiration by global hydrological models (GHMs), it is necessary to replace the bucket-like linear GW reservoir model typical for hydrological models with a fully integrated gradient-based GW flow model. Linear reservoir models can only simulate GW discharge to SW bodies, provide no information on the location of the GW table and assume that there is no GW flow among grid cells. A gradient-based GW model simulates not only GW storage but also hydraulic head, which together with information on SW table elevation enables the quantification of water flows from GW to SW and vice versa. In addition, hydraulic heads are the basis for calculating lateral GW flow among grid cells and capillary rise.
G³M is a new global gradient-based GW model with a spatial resolution of 5' that will replace the current linear GW reservoir in the 0.5° WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables inmemory coupling to WGHM while keeping overall runtime relatively low, allowing sensitivity analyses and data assimilation. This paper presents the G³M concept and specific model design decisions together with results under steady-state naturalized conditions, i.e. neglecting GW abstractions. Cell-specific conductances of river beds, which govern GW-SW interaction, were determined based on the 30'' steady-state water table computed by Fan et al. (2013). Together with an appropriate choice for the effective elevation of the SW table within each grid cell, this enables a reasonable simulation of drainage from GW to SW such that, in contrast to the GW model of de Graaf et al. (2015, 2017), no additional drainage based on externally provided values for GW storage above the floodplain is required in G³M. Comparison of simulated hydraulic heads to observations around the world shows better agreement than de Graaf et al. (2015). In addition, G³M output is compared to the output of two established macro-scale models for the Central Valley, California, and the continental United States, respectively. As expected, depth to GW table is highest in mountainous and lowest in flat regions. A first analysis of losing and gaining rivers and lakes/wetlands indicates that GW discharge to rivers is by far the dominant flow, draining diffuse GW recharge, such that lateral flows only become a large fraction of total diffuse and focused recharge in case of losing rivers and some areas with very low GW recharge. G³M does not represent losing rivers in some dry regions. This study presents the first steps towards replacing the linear GW reservoir model in a GHM while improving on recent efforts, demonstrating the feasibility of the approach and the robustness of the newly developed framework.
In global hydrological models, groundwater storages and flows are generally simulated by linear reservoir models. Recently, the first global gradient-based groundwater models were developed in order to improve the representation of groundwater-surface water interactions, capillary rise, lateral flows and human water use impacts. However, the reliability of model outputs is limited by a lack of data as well as model assumptions required due to the necessarily coarse spatial resolution. The impact of data quality is presented by showing the sensitivity of a groundwater model to changes in the only available global hydraulic conductivity data-set. To better understand the sensitivity of model output to uncertain spatially distributed parameter inputs, we present the first application of a global sensitivity method for a global-scale groundwater model using nearly 2000 steady-state model runs of the global gradient-based groundwater model G3M. By applying the Morris method in a novel domain decomposition approach that identifies global hydrological response units, spatially distributed parameter sensitivities are determined for a computationally expensive model. Results indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, though parameter sensitivities vary regionally. For large areas of the globe, rivers are simulated to be either losing or gaining, depending on the parameter combination, indicating a high uncertainty of simulating the direction of flow between the two compartments. Mountainous and dry regions show a high variance in simulated head due to numerical difficulties of the model, limiting the reliability of computed sensitivities in these regions. This instability is likely caused by the uncertainty in surface water body elevation. We conclude that maps of spatially distributed sensitivities can help to understand complex behaviour of models that incorporate data with varying spatial uncertainties. The findings support the selection of possible calibration parameters and help to anticipate challenges for a transient coupling of the model.
Mémoire géologique sur le Bassin d'Amiens, et en partienlier sur les Cantons Littoraux de la Somme
(1834)
The CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) passenger aircraft observatory performed in situ measurements at 10–12 km altitude in the South Asian summer monsoon anticyclone between June and September 2008. These measurements enable us to investigate this atmospheric region (which so far has mostly been observed from satellites) using the broad suite of trace gases and aerosol particles measured by CARIBIC. Elevated levels of a variety of atmospheric pollutants (e.g. carbon monoxide, total reactive nitrogen oxides, aerosol particles, and several volatile organic compounds) were recorded. The measurements provide detailed information about the chemical composition of air in different parts of the monsoon anticyclone, particularly of ozone precursors. While covering a range of 3500 km inside the monsoon anticyclone, CARIBIC observations show remarkable consistency, i.e. with distinct latitudinal patterns of trace gases during the entire monsoon period.
Using the CARIBIC trace gas and aerosol particle measurements in combination with the Lagrangian particle dispersion model FLEXPART, we investigated the characteristics of monsoon outflow and the chemical evolution of air masses during transport. The trajectory calculations indicate that these air masses originated mainly from South Asia and mainland Southeast Asia. Estimated photochemical ages of the air were found to agree well with transport times from a source region east of 90–95° E. The photochemical ages of the air in the southern part of the monsoon anticyclone were systematically younger (less than 7 days) and the air masses were mostly in an ozone-forming chemical mode. In its northern part the air masses were older (up to 13 days) and had unclear ozone formation or destruction potential. Based on analysis of forward trajectories, several receptor regions were identified. In addition to predominantly westward transport, we found evidence for efficient transport (within 10 days) to the Pacific and North America, particularly during June and September, and also of cross-tropopause exchange, which was strongest during June and July. Westward transport to Africa and further to the Mediterranean was the main pathway during July.
The CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) passenger aircraft observatory performed in situ measurements at 10–12 km altitude in the South Asian summer monsoon anticyclone between June and September 2008. These measurements enable us to investigate this atmospheric region, which so far has mostly been observed from satellites, using the broad suite of trace gases and aerosols measured by CARIBIC. Elevated levels of a range of atmospheric pollutants were recorded e.g. carbon monoxide, total reactive nitrogen oxides, aerosol particles and several volatile organic compounds. The measurements provide detailed information about the chemical composition of air in different parts of the monsoon anticyclone, particularly of ozone precursors. While covering a range of 3500 km inside the monsoon anticyclone, CARIBIC observations show remarkable consistency, i.e. with regular latitudinal patterns of trace gases during the entire monsoon period. Trajectory calculations indicate that these air masses originated mainly from South Asia and Mainland Southeast Asia.
Using the CARIBIC trace gas and aerosol measurements in combination with the Lagrangian particle dispersion model FLEXPART we investigated the characteristics of monsoon outflow and the chemical evolution of air masses during transport. Estimated photochemical ages of the air were found to agree well with transport times from a source region east of 95 ° E. The photochemical ages of the air in the southern part of the monsoon anticyclone were consistently younger (less than 7 days) and the air masses mostly in an ozone forming chemical regime. In its northern part the air masses were older (up to 13 days) and had unclear ozone formation or destruction potential. Based on analysis of forward trajectories several receptor regions were identified. In addition to predominantly westward transport, we found evidence for efficient transport (within 10 days) to the Pacific and North America, particularly during June and September, and also of cross-tropopause exchange, which was strongest during June and July. Westward transport to Africa and further to the Mediterranean was the main pathway during July.
Highlights
• New fumarole and thermal water data for Askja and Kverkfjöll volcanoes, Iceland.
• Data compared to modelled compositions and fluxes of magmatic gas.
• Fumarole compositions compatible with origin of CO2 and S from degassing intrusions.
• Intrusive magmatic fluxes sufficient to sustain hydrothermal fluxes of CO2 and S in Iceland
• Magma degassing insignificant/minor source of H2O and Cl to Icelandic hydrothermal fluids
Abstract
Mantle volatiles are transported to Earth's crust and surface by basaltic volcanism. During subaerial eruptions, vast amounts of carbon, sulfur and halogens can be released to the atmosphere during a short time-interval, with impacts ranging in scale from the local environment to the global climate. By contrast, passive volatile release at the surface originating from magmatic intrusions is characterized by much lower flux, yet may outsize eruptive volatile quantities over long timescales. Volcanic hydrothermal systems (VHSs) act as conduits for such volatile release from degassing intrusions and can be used to gauge the contribution of intrusive magmatism to global volatile cycles. Here, we present new compositional and isotopic (δD and δ18O-H2O, 3He/4He, δ13C-CO2, Δ33S-δ34S-H2S and SO4) data for thermal waters and fumarole gases from the Askja and Kverkfjöll volcanoes in central Iceland. We use the data together with magma degassing modelling and mass balance calculations to constrain the sources of volatiles in VHSs and to assess the role of intrusive magmatism to the volcanic volatile emission budgets in Iceland.
The CO2/ΣS (10−30), 3He/4He (8.3–10.5 RA; 3He/4He relative to air), δ13C-CO2 (−4.1 to −0.2 ‰) and Δ33S-δ34S-H2S (−0.031 to 0.003 ‰ and −1.5 to +3.6‰) values in high-gas flux fumaroles (CO2 > 10 mmol/mol) are consistent with an intrusive magmatic origin for CO2 and S at Askja and Kverkfjöll. We demonstrate that deep (0.5–5 kbar, equivalent to ∼2–18 km crustal depth) decompression degassing of basaltic intrusions in Iceland results in CO2 and S fluxes of 330–5060 and 6–210 kt/yr, respectively, which is sufficient to account for the estimated CO2 flux of Icelandic VHSs (3365–6730 kt/yr), but not the VHS S flux (220–440 kt/yr). Secondary, crystallization-driven degassing from maturing intrusions and leaching of crustal rocks are suggested as additional sources of S. Only a minor proportion of the mantle flux of Cl is channeled via VHSs whereas the H2O flux remains poorly constrained, because magmatic signals in Icelandic VHSs are masked by a dominant shallow groundwater component of meteoric water origin. These results suggest that the bulk of the mantle CO2 and S flux to the atmosphere in Iceland is supplied by intrusive, not eruptive magmatism, and is largely vented via hydrothermal fields.
Assessing the uncertainties of simulation results of ecological models is becoming increasingly important, specifically if these models are used to estimate greenhouse gas emissions on site to regional/national levels. Four general sources of uncertainty effect the outcome of process-based models: (i) uncertainty of information used to initialise and drive the model, (ii) uncertainty of model parameters describing specific ecosystem processes, (iii) uncertainty of the model structure, and (iv) accurateness of measurements (e.g., soil-atmosphere greenhouse gas exchange) which are used for model testing and development.
The aim of our study was to assess the simulation uncertainty of the process-based biogeochemical model LandscapeDNDC. For this we set up a Bayesian framework using a Markov Chain Monte Carlo (MCMC) method, to estimate the joint model parameter distribution. Data for model testing, parameter estimation and uncertainty assessment were taken from observations of soil fluxes of nitrous oxide (N2O), nitric oxide (NO) and carbon dioxide (CO2) as observed over a 10 yr period at the spruce site of the Höglwald Forest, Germany. By running four independent Markov Chains in parallel with identical properties (except for the parameter start values), an objective criteria for chain convergence developed by Gelman et al. (2003) could be used.
Our approach shows that by means of the joint parameter distribution, we were able not only to limit the parameter space and specify the probability of parameter values, but also to assess the complex dependencies among model parameters used for simulating soil C and N trace gas emissions. This helped to improve the understanding of the behaviour of the complex LandscapeDNDC model while simulating soil C and N turnover processes and associated C and N soil-atmosphere exchange. In a final step the parameter distribution of the most sensitive parameters determining soil-atmosphere C and N exchange were used to obtain the parameter-induced uncertainty of simulated N2O, NO and CO2 emissions. These were compared to observational data of an calibration set (6 yr) and an independent validation set of 4 yr. The comparison showed that most of the annual observed trace gas emissions were in the range of simulated values and were predicted with a high certainty (Root-mean-squared error (RMSE) NO: 2.4 to 18.95 g N ha−1 d−1, N2O: 0.14 to 21.12 g N ha−1 d−1, CO2: 5.4 to 11.9 kg C ha−1 d−1). However, LandscapeDNDC simulations were sometimes still limited to accurately predict observed seasonal variations in fluxes.
Assessing the uncertainties of simulation results of ecological models is becoming of increasing importance, specifically if these models are used to estimate greenhouse gas emissions at site to regional/national levels. Four general sources of uncertainty effect the outcome of process-based models: (i) uncertainty of information used to initialise and drive the model, (ii) uncertainty of model parameters describing specific ecosystem processes, (iii) uncertainty of the model structure and (iv) accurateness of measurements (e.g. soil-atmosphere greenhouse gas exchange) which are used for model testing and development.
The aim of our study was to assess the simulation uncertainty of the process-based biogeochemical model LandscapeDNDC. For this we set up a Bayesian framework using a Markov Chain Monte Carlo (MCMC) method, to estimate the joint model parameter distribution. Data for model testing, parameter estimation and uncertainty assessment were taken from observations of soil fluxes of nitrous oxide (N2O), nitric oxide (NO), and carbon dioxide (CO2) as observed over a 10 yr period at the spruce site of the Höglwald Forest, Germany. By running four independent Markov Chains in parallel with identical properties (except for the parameter start values), an objective criteria for chain convergence developed by Gelman et al. (2003) could be used.
Our approach showed that by means of the joined parameter distribution, we were able not only to limit the parameter space and specify the probability of parameter values, but also to assess the complex dependencies among model parameters used for simulating soil C and N trace gas emissions. This helped to improve the understanding of the behaviour of the complex LandscapeDNDC model while simulating soil C and N turnover processes and associated C and N soil-atmosphere exchange.
In a final step the parameter distribution of the most sensitive parameters determining soil-atmosphere C and N exchange were used to obtain the parameter-induced uncertainty of simulated N2O, NO and CO2 emissions. These were compared to observational data of the calibration set (6 yr) and an independent validation set of 4 yr.
The comparison showed that most of the annual observed trace gas emissions were in the range of simulated values and were predicted with a high certainty (Residual mean squared error (RMSE) NO: 2.5 to 21.3 g N ha−1 d−1, N2O: 0.2 to 21.4 g N ha−1 d−1, CO2: 5.8 to 12.6 kg C ha−1 d−1). However, LandscapeDNDC simulations were sometimes limited to accurately predict observed seasonal variations in fluxes.
Living on the edge: environmental variability of a shallow late Holocene cold-water coral mound
(2022)
Similar to their tropical counterparts, cold-water corals (CWCs) are able to build large three-dimensional reef structures. These unique ecosystems are at risk due to ongoing climate change. In particular, ocean warming, ocean acidification and changes in the hydrological cycle may jeopardize the existence of CWCs. In order to predict how CWCs and their reefs or mounds will develop in the near future one important strategy is to study past fossil CWC mounds and especially shallow CWC ecosystems as they experience a greater environmental variability compared to other deep-water CWC ecosystems. We present results from a CWC mound off southern Norway. A sediment core drilled from this relatively shallow (~ 100 m) CWC mound exposes in full detail hydrographical changes during the late Holocene, which were crucial for mound build-up. We applied computed tomography, 230Th/U dating, and foraminiferal geochemical proxy reconstructions of bottom-water-temperature (Mg/Ca-based BWT), δ18O for seawater density, and the combination of both to infer salinity changes. Our results demonstrate that the CWC mound formed in the late Holocene between 4 kiloannum (ka) and 1.5 ka with an average aggradation rate of 104 cm/kiloyears (kyr), which is significantly lower than other Holocene Norwegian mounds. The reconstructed BWTMg/Ca and seawater density exhibit large variations throughout the entire period of mound formation, but are strikingly similar to modern in situ observations in the nearby Tisler Reef. We argue that BWT does not exert a primary control on CWC mound formation. Instead, strong salinity and seawater density variation throughout the entire mound sequence appears to be controlled by the interplay between the Atlantic Water (AW) inflow and the overlying, outflowing Baltic-Sea water. CWC growth and mound formation in the NE Skagerrak was supported by strong current flow, oxygen replenishment, the presence of a strong boundary layer and larval dispersal through the AW, but possibly inhibited by the influence of fresh Baltic Water during the late Holocene. Our study therefore highlights that modern shallow Norwegian CWC reefs may be particularly endangered due to changes in water-column stratification associated with increasing net precipitation caused by climate change.
The Fire Modeling Intercomparison Project (FireMIP), phase 1: experimental and analytical protocols
(2016)
The important role of fire in regulating vegetation community composition and contributions to emissions of greenhouse gases and aerosols make it a critical component of dynamic global vegetation models and Earth system models. Over two decades of development, a wide variety of model structures and mechanisms have been designed and incorporated into global fire models, which have been linked to different vegetation models. However, there has not yet been a systematic examination of how these different strategies contribute to model performance. Here we describe the structure of the first phase of the Fire Model Intercomparison Project (FireMIP), which for the first time seeks to systematically compare a number of models. By combining a standardized set of input data and model experiments with a rigorous comparison of model outputs to each other and to observations, we will improve the understanding of what drives vegetation fire, how it can best be simulated, and what new or improved observational data could allow better constraints on model behavior. Here we introduce the fire models used in the first phase of FireMIP, the simulation protocols applied, and the benchmarking system used to evaluate the models. The works published in this journal are distributed under the Creative Commons Attribution 3.0 License. This license does not affect the Crown copy-right work, which is re-usable under the Open Government Licence (OGL). The Creative Commons Attribution 3.0 License and the OGL are interoperable and do not conflict with, reduce, or limit each other.
The important role of fire in regulating vegetation community composition and contributions to emissions of greenhouse gases and aerosols make it a critical component of dynamic global vegetation models and Earth system models. Over 2 decades of development, a wide variety of model structures and mechanisms have been designed and incorporated into global fire models, which have been linked to different vegetation models. However, there has not yet been a systematic examination of how these different strategies contribute to model performance. Here we describe the structure of the first phase of the Fire Model Intercomparison Project (FireMIP), which for the first time seeks to systematically compare a number of models. By combining a standardized set of input data and model experiments with a rigorous comparison of model outputs to each other and to observations, we will improve the understanding of what drives vegetation fire, how it can best be simulated, and what new or improved observational data could allow better constraints on model behavior. In this paper, we introduce the fire models used in the first phase of FireMIP, the simulation protocols applied, and the benchmarking system used to evaluate the models. We have also created supplementary tables that describe, in thorough mathematical detail, the structure of each model.
Highlights
• We present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series.
• The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios.
• The model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33.
Abstract
High-rate Global Navigation Satellite System (HR-GNSS) data can be highly useful for earthquake analysis as it provides continuous high-frequency measurements of ground motion. This data can be used to analyze diverse parameters related to the seismic source and to assess the potential of an earthquake to prompt strong motions at certain distances and even generate tsunamis. In this work, we present the first results of a deep learning model based on a convolutional neural network for earthquake magnitude estimation, using HR-GNSS displacement time series. The influence of different dataset configurations, such as station numbers, epicentral distances, signal duration, and earthquake size, were analyzed to figure out how the model can be adapted to various scenarios. We explored the potential of the model for global application and compared its performance using both synthetic and real data from different seismogenic regions. The performance of our model at this stage was satisfactory in estimating earthquake magnitude from synthetic data with 0.07 ≤ RMS ≤ 0.11. Comparable results were observed in tests using synthetic data from a different region than the training data, with RMS ≤ 0.15. Furthermore, the model was tested using real data from different regions and magnitudes, resulting in the best cases with 0.09 ≤ RMS ≤ 0.33, provided that the data from a particular group of stations had similar epicentral distance constraints to those used during the model training. The robustness of the DL model can be improved to work independently from the window size of the time series and the number of stations, enabling faster estimation by the model using only near-field data. Overall, this study provides insights for the development of future DL approaches for earthquake magnitude estimation with HR-GNSS data, emphasizing the importance of proper handling and careful data selection for further model improvements.
We evaluate the influence of a forest parametrization on the simulation of the boundary layer flow over moderate complex terrain in the context of the Perdigão 2017 field campaign. The numerical simulations are performed using the Weather Research and Forecasting model in large eddy simulation mode (WRF-LES). The short-term, high-resolution (40 m horizontal grid spacing) and long-term (200 m horizontal grid spacing) WRF-LES are evaluated for an integration time of 12 h and 1.5 months, respectively, with and without forest parameterization. The short-term simulations focus on low-level jet events over the valley, while the long-term simulations cover the whole intensive observation period (IOP) of the field campaign. The results are validated using lidar and meteorological tower observations. The mean diurnal cycle during the IOP shows a significant improvement of the along-valley wind speed and the wind direction when using the forest parametrization. However, the drag imposed by the parametrization results in an underestimation of the cross-valley wind speed, which can be attributed to a poor representation of the land surface characteristics. The evaluation of the high-resolution WRF-LES shows a positive influence of the forest parametrization on the simulated winds in the first 500 m above the surface.
We evaluate the influence of a forest parametrization on the simulation of the boundary layer flow over moderate complex terrain in the context of the Perdigão 2017 field campaign. The numerical simulations are performed using the Weather Research and Forecasting model in large eddy simulation mode (WRF-LES). The short-term, high-resolution (40 m horizontal grid spacing) and long-term (200 m horizontal grid spacing) WRF-LES are evaluated for an integration time of 12 h and 1.5 months, respectively, with and without forest parameterization. The short-term simulations focus on low-level jet events over the valley, while the long-term simulations cover the whole intensive observation period (IOP) of the field campaign. The results are validated using lidar and meteorological tower observations. The mean diurnal cycle during the IOP shows a significant improvement of the along-valley wind speed and the wind direction when using the forest parametrization. However, the drag imposed by the parametrization results in an underestimation of the cross-valley wind speed, which can be attributed to a poor representation of the land surface characteristics. The evaluation of the high-resolution WRF-LES shows a positive influence of the forest parametrization on the simulated winds in the first 500 m above the surface.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations as well as hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a one-year period and full seasonal cycle (March 2014 - February 2015). The presented measurements provide a climatology of CCN properties for a characteristic central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The observed mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), elevated values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter κ exhibits remarkably little temporal variability: no pronounced diurnal cycles, weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modelling purposes, we compare different approaches of predicting CCN number concentration and present a novel parameterization, which allows accurate CCN predictions based on a small set of input data.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a 1-year period and full seasonal cycle (March 2014–February 2015). The measurements provide a climatology of CCN properties characteristic of a remote central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), higher values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter, κ, exhibits remarkably little temporal variability: no pronounced diurnal cycles, only weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modeling purposes, we compare different approaches of predicting CCN number concentration and present a novel parametrization, which allows accurate CCN predictions based on a small set of input data.
Size-resolved measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a full seasonal cycle (Mar 2014–Feb 2015). In a companion part 1 paper, we presented an in-depth CCN characterization based on annually as well as seasonally averaged time intervals and discuss different parametrization strategies to represent the Amazonian CCN cycling in modelling studies (M. Pöhlker et al., 2016b). The present part 2 study analyzes the aerosol and CCN variability in original time resolution and, thus, resolves aerosol advection and transformation for the following case studies, which represent the most characteristic states of the Amazonian atmosphere:
1. Near-pristine (NP) conditions, defined as the absence of detectable black carbon (< 0.01 µg m−3), showed their highest occurrence (up to 30 %) in the wet season (i.e., Mar–May). On average, the NP episodes are characterized by a bimodal aerosol size distribution (strong Aitken mode: DAit = 70 nm, NAit = ~ 200 cm−3 vs. weaker accumulation mode: Dacc = 170 nm, Nacc = ~ 60 cm−3), a mostly organic particle composition, and relatively low hygroscopicity levels (κAit = 0.12 vs. κacc = 0.18). The NP CCN efficiency spectrum shows that the CCN population is sensitive to changes in supersaturation (S) over a wide S range.
2. Long-range transport (LRT) conditions frequently mix Saharan dust, African combustion smoke, and sea spray aerosols into the Amazonian wet season atmosphere. The LRT episodes (i.e., Feb–Apr) are characterized by an accumulation mode dominated size distribution (DAit = 80 nm, NAit = 120 cm−3 vs. Dacc = 180 nm, Nacc = 300 cm−3), a clearly increased abundance of dust and salt compounds, and relatively high hygroscopicity levels (κAit = 0.18, κacc = 0.34). The LRT CCN efficiency spectrum shows that the CCN population is highly sensitive to changes in S in the low S regime.
3. Biomass burning (BB) conditions dominate the Amazonian dry season. A selected characteristic BB episode shows a very strong accumulation mode (DAit = 70 nm, NAit = ~ 140 cm−3 vs. Dacc = 170 nm, Nacc = ~ 3400 cm−3), particles with very high organic fractions (> 90 %), and correspondingly low hygroscopicity levels (κAit = 0.14, κacc = 0.17). The BB CCN efficiency spectrum shows that the CCN population is highly sensitive to changes in S in the low S regime.
4. Mixed pollution conditions show the superposition of African (i.e., volcanic) and Amazonian (i.e., biomass burning) aerosol emissions during the dry season. The African aerosols showed a broad monomodal distribution (D = 130 nm, N = ~ 1300 cm−3), with very high sulfate fractions (20 %), and correspondingly high hygroscopicity (κAit = 0.14, κacc = 0.22). This was superimposed by fresh smoke from nearby fires with one strong mode (D = 113 nm, Nacc = ~ 2800 cm−3), an organic-dominated aerosol, and sharply decreased hygroscopicity (κAit = 0.10, κacc = 0.20). These conditions underline the rapidly changing pollution regimes with clear impacts on the aerosol and CCN properties.
Overall, this study provides detailed insights into the CCN cycling in relation to aerosol-cloud interaction in the vulnerable and climate-relevant Amazon region. The detailed analysis of aerosol and CCN key properties and particularly the extracted CCN efficiency spectra with the associated fit parameters provide a basis for an in-depth analysis of aerosol-cloud interaction in the Amazon and beyond.
Two aircraft instruments for the measurement of total odd nitrogen (NOy) were compared side by side aboard a Learjet A35 in April 2003 during a campaign of the AFO2000 project SPURT (Spurengastransport in der Tropopausenregion). The instruments albeit employing the same measurement principle (gold converter and chemiluminescence) had different inlet configurations. The ECO-Physics instrument operated by ETH-Zürich in SPURT had the gold converter mounted outside the aircraft, whereas the instrument operated by FZ-Jülich in the European project MOZAIC III (Measurements of ozone, water vapour, carbon monoxide and nitrogen oxides aboard Airbus A340 in-service aircraft) employed a Rosemount probe with 80 cm of FEP-tubing connecting the inlet to the gold converter. The NOy concentrations during the flight ranged between 0.3 and 3 ppb. The two data sets were compared in a blind fashion and each team followed its normal operating procedures. On average, the measurements agreed within 7%, i.e. within the combined uncertainty of the two instruments. This puts an upper limit on potential losses of HNO3 in the Rosemount inlet of the MOZAIC instrument. Larger transient deviations were observed during periods after calibrations and when the aircraft entered the stratosphere. The time lag of the MOZAIC instrument observed in these instances is in accordance with the time constant of the MOZAIC inlet line determined in the laboratory for HNO3.
Two aircraft instruments for the measurement of total odd nitrogen (NOy) were compared side by side aboard a Learjet A35 in April 2003 during a campaign of the AFO2000 project SPURT (Spurengastransport in der Tropopausenregion). The instruments albeit employing the same measurement principle (gold converter and chemiluminescence) had different inlet configurations. The ECO-Physics instrument operated by ETH-Zürich in SPURT had the gold converter mounted outside the aircraft, whereas the instrument operated by FZ-Jülich in the European project MOZAIC III (Measurements of ozone, water vapour, carbon monoxide and nitrogen oxides aboard Airbus A340 in-service aircraft) employed a Rosemount probe with 80 cm of FEP-tubing connecting the inlet to the gold converter. The NOy concentrations during the flight ranged between 0.3 and 3 ppb. The two data sets were compared in a blind fashion and each team followed its normal operating procedures. On average, the measurements agreed within 6%, i.e. within the combined uncertainty of the two instruments. This puts an upper limit on potential losses of HNO3 in the Rosemount inlet of the MOZAIC instrument. Larger transient deviations were observed during periods after calibrations and when the aircraft entered the stratosphere. The time lag of the MOZAIC instrument observed in these instances is in accordance with the time constant of the MOZAIC inlet line determined in the laboratory for HNO3.
Convective rain cell properties and the resulting precipitation scaling in a warm-temperate climate
(2022)
Convective precipitation events have been shown to intensify at rates exceeding the Clausius–Clapeyron rate (CC rate) of ca. 7% K−1 under current climate conditions. In this study, we relate atmospheric variables (low-level dew point temperature, convective available potential energy, and vertical wind shear), which are regarded as ingredients for severe deep convection, to properties of convective rain cells (cell area, maximum precipitation intensity, lifetime, precipitation sum, and cell speed). The rain cell properties are obtained from a rain gauge-adjusted radar dataset in a mid-latitude region, which is characterized by a temperate climate with warm summers (Germany). Different Lagrangian cell properties scale with dew point temperature at varying rates. While the maximum precipitation intensity of cells scales consistently at the CC rate, the area and precipitation sum per cell scale at varying rates above the CC rate. We show that this super-CC scaling is caused by a covarying increase of convective available potential energy with dew point temperature. Wind shear increases the precipitation sum per cell mainly by increasing the spatial cell extent. From a Eulerian point of view, this increase is partly compensated by a higher cell velocity, which leads to Eulerian precipitation scaling rates close to and slightly above the CC rate. Thus, Eulerian scaling rates of convective precipitation are modulated by convective available potential energy and vertical wind shear, making it unlikely that present scaling rates can be applied to future climate conditions. Furthermore, we show that cells that cause heavy precipitation at fixed locations occur at low vertical wind shear and, thus, move relatively slowly compared to typical cells.
Convective shower characteristics simulated with the convection-permitting climate model COSMO-CLM
(2019)
This paper evaluates convective precipitation as simulated by the convection-permitting climate model (CPM) Consortium for Small-Scale Modeling in climate mode (COSMO-CLM) (with 2.8 km grid-spacing) over Germany in the period 2001–2015. Characteristics of simulated convective precipitation objects like lifetime, area, mean intensity, and total precipitation are compared to characteristics observed by weather radar. For this purpose, a tracking algorithm was applied to simulated and observed precipitation with 5-min temporal resolution. The total amount of convective precipitation is well simulated, with a small overestimation of 2%. However, the simulation underestimates convective activity, represented by the number of convective objects, by 33%. This underestimation is especially pronounced in the lowlands of Northern Germany, whereas the simulation matches observations well in the mountainous areas of Southern Germany. The underestimation of activity is compensated by an overestimation of the simulated lifetime of convective objects. The observed mean intensity, maximum intensity, and area of precipitation objects increase with their lifetime showing the spectrum of convective storms ranging from short-living single-cell storms to long-living organized convection like supercells or squall lines. The CPM is capable of reproducing the lifetime dependence of these characteristics but shows a weaker increase in mean intensity with lifetime resulting in an especially pronounced underestimation (up to 25%) of mean precipitation intensity of long-living, extreme events. This limitation of the CPM is not identifiable by classical evaluation techniques using rain gauges. The simulation can reproduce the general increase of the highest percentiles of cell area, total precipitation, and mean intensity with temperature but fails to reproduce the increase of lifetime. The scaling rates of mean intensity and total precipitation resemble observed rates only in parts of the temperature range. The results suggest that the evaluation of coarse-grained (e.g., hourly) precipitation fields is insufficient for revealing challenges in convection-permitting simulations.
Extreme convective precipitation is expected to increase with global warming. However, the rate of increase and the understanding of contributing processes remain highly uncertain. We investigated characteristics of convective rain cells like area, intensity, and lifetime as simulated by a convection-permitting climate model in the area of Germany under historical (1976–2005) and future (end-of-century, RCP8.5 scenario) conditions. To this end, a tracking algorithm was applied to 5-min precipitation output. While the number of convective cells is virtually similar under historical and future conditions, there are more intense and larger cells in the future. This yields an increase in hourly precipitation extremes, although mean precipitation decreases. The relative change in the frequency distributions of area, intensity, and precipitation sum per cell is highest for the most extreme percentiles, suggesting that extreme events intensify the most. Furthermore, we investigated the temperature and moisture scaling of cell characteristics. The temperature scaling drops off at high temperatures, with a shift in drop-off towards higher temperatures in the future, allowing for higher peak values. In contrast, dew point temperature scaling shows consistent rates across the whole dew point range. Cell characteristics scale at varying rates, either below (mean intensity), at about (maximum intensity and area), or above (precipitation sum) the Clausius–Clapeyron rate. Thus, the widely investigated extreme precipitation scaling at fixed locations is a complex product of the scaling of different cell characteristics. The dew point scaling rates and absolute values of the scaling curves in historical and future conditions are closest for the highest percentiles. Therefore, near-surface humidity provides a good predictor for the upper limit of for example, maximum intensity and total precipitation of individual convective cells. However, the frequency distribution of the number of cells depending on dew point temperature changes in the future, preventing statistical inference of extreme precipitation from near-surface humidity.
Das Cranium eines fossilen Hominiden des Formenkreises Homo sapiens sapiens wurde relativ-geologisch sowie absolut durch Radiokohlenstoff und Aminosäuren auf ungefähr 31 000 Jahre B.P. datiert. Andere absolute sowie relative Daten an Mollusken und Mammutzähnen in überlagernden jüngeren Straten datieren auf 18 000 — 21000 und 16 000 Jahren B.P. Geomorphologische und geophysikalische Datierungen stimmen somit gut überein. Er ist der älteste datierte und früheste Bewohner Zentraleuropas, der dem Homo sapiens sapiens angehört.
The frequency of extreme events has changed, having a direct impact on human lives. Regional climate models help us to predict these regional climate changes. This work presents an atmosphere–ocean coupled regional climate system model (RCSM; with the atmospheric component COSMO-CLM and the ocean component NEMO) over the European domain, including three marginal seas: the Mediterranean, North, and Baltic Sea. To test the model, we evaluate a simulation of more than 100 years (1900–2009) with a spatial grid resolution of about 25 km. The simulation was nested into a coupled global simulation with the model MPI-ESM in a low-resolution configuration, whose ocean temperature and salinity were nudged to the ocean–ice component of the MPI-ESM forced with the NOAA 20th Century Reanalysis (20CR). The evaluation shows the robustness of the RCSM and discusses the added value by the coupled marginal seas over an atmosphere-only simulation. The coupled system is stable for the complete 20th century and provides a better representation of extreme temperatures compared to the atmosphere-only model. The produced long-term dataset will help us to better understand the processes leading to meteorological and climate extremes.
Parabens and sorbic acid are commonly used as food preservatives due to their antimicrobial effect. However, their use in foods for infants and young children is not permitted in the European Union. Previous studies found these compounds in some gel-filled baby teethers, whereby parabens, which are well-known as endocrine disruptors, were identified in the polymer-based chewing surface consisting of ethylene-vinyl acetate (EVA). To assess the exposure of infants and young children to these products, the application of parabens in teethers should be thoroughly investigated. Therefore, the present study aimed to apply a representative migration test procedure combined with an accurate analytical method to examine gel-filled baby teethers without elaborate sample preparation, high costs, and long processing times. Accordingly, solid-phase extraction (SPE), in combination with a stable isotope dilution assay (SIDA) and subsequent gas chromatography–mass spectrometry (GC–MS) for analysis of methyl-, ethyl-, and n-propylparaben (MeP, EtP, and n-PrP), was found to be well-suited, with recoveries ranging from 93 to 99%. The study compared the release of these parabens from intact teether surfaces into water and saliva simulant under real-life conditions, with total amounts of detected parabens found to be in the range of 101–162 µg 100 mL−1 and 57–148 µg 100 mL−1, respectively. Furthermore, as a worst-case scenario, the release into water was examined using a long-term migration study.
Often in climate system studies, linear and symmetric statistical measures are applied to quantify interactions among subsystems or variables. However, they do not allow identification of the driving and responding subsystems. Therefore, in this study, we aimed to apply asymmetric measures from information theory: the axiomatically proposed transfer entropy and the first principle-based information flow to detect and quantify climate interactions. As their estimations are challenging, we initially tested nonparametric estimators like transfer entropy (TE)-binning, TE-kernel, and TE k-nearest neighbor and parametric estimators like TE-linear and information flow (IF)-linear with idealized two-dimensional test cases along with their sensitivity on sample size. Thereafter, we experimentally applied these methods to the Lorenz-96 model and to two real climate phenomena, i.e., (1) the Indo-Pacific Ocean coupling and (2) North Atlantic Oscillation (NAO)–European air temperature coupling. As expected, the linear estimators work for linear systems but fail for strongly nonlinear systems. The TE-kernel and TE k-nearest neighbor estimators are reliable for linear and nonlinear systems. Nevertheless, the nonparametric methods are sensitive to parameter selection and sample size. Thus, this work proposes a composite use of the TE-kernel and TE k-nearest neighbor estimators along with parameter testing for consistent results. The revealed information exchange in Lorenz-96 is dominated by the slow subsystem component. For real climate phenomena, expected bidirectional information exchange between the Indian and Pacific SSTs was detected. Furthermore, expected information exchange from NAO to European air temperature was detected, but also unexpected reversal information exchange. The latter might hint to a hidden process driving both the NAO and European temperatures. Hence, the limitations, availability of time series length and the system at hand must be taken into account before drawing any conclusions from TE and IF-linear estimations.
In 1998 the German Universities of Kassel and Giessen organised a workshop on water and solute transport in large drainage basins. The workshop focused on analysing and summarising the state of research, existing problems and perspectives in this research area. It was the second of a series of annual workshops since 1997 that became an important discussion forum for the German-speaking research community in the field of hydrological modelling. Now the 11th Workshop on Large-scale Hydrological Modelling referred to the same questions as posed in 1998 in order to evaluate the developments and advances of the last ten years. Based on keynote presentations, the workshop focused on discussion in working groups where also posters were presented. This volume of "Advances in Geosciences" comprises seven papers referring to the poster contributions. At the end of the volume, an overview paper summarises the outcome of the workshop presentations and discussions (Doll et al.). ...
Drought is understood as both a lack of water (i.e., a deficit compared to demand) and a temporal anomaly in one or more components of the hydrological cycle. Most drought indices, however, only consider the anomaly aspect, i.e., how unusual the condition is. In this paper, we present two drought hazard indices that reflect both the deficit and anomaly aspects. The soil moisture deficit anomaly index, SMDAI, is based on the drought severity index, DSI (Cammalleri et al., 2016), but is computed in a more straightforward way that does not require the definition of a mapping function. We propose a new indicator of drought hazard for water supply from rivers, the streamflow deficit anomaly index, QDAI, which takes into account the surface water demand of humans and freshwater biota. Both indices are computed and analyzed at the global scale, with a spatial resolution of roughly 50 km, for the period 1981–2010, using monthly time series of variables computed by the global water resources and the model WaterGAP 2.2d. We found that the SMDAI and QDAI values are broadly similar to values of purely anomaly-based indices. However, the deficit anomaly indices provide more differentiated spatial and temporal patterns that help to distinguish the degree and nature of the actual drought hazard to vegetation health or the water supply. QDAI can be made relevant for stakeholders with different perceptions about the importance of ecosystem protection, by adapting the approach for computing the amount of water that is required to remain in the river for the well-being of the river ecosystem. Both deficit anomaly indices are well suited for inclusion in local or global drought risk studies.
Drought is understood as both a lack of water (i.e., a deficit as compared to some requirement) and an anomaly in the condition of one or more components of the hydrological cycle. Most drought indices, however, only consider the anomaly aspect, i.e., how unusual the condition is. In this paper, we present two drought hazard indices that reflect both the deficit and anomaly aspects. The soil moisture deficit anomaly index, SMDAI, is based on the drought severity index, DSI, but is computed in a more straightforward way that does not require the definition of a mapping function. We propose a new indicator of drought hazard for water supply from rivers, the streamflow deficit anomaly index, QDAI, which takes into account the surface water demand of humans and freshwater biota. Both indices are computed and analyzed at the global scale, with a spatial resolution of roughly 50 km, for the period 1981-2010, using monthly time series of variables computed by the global water resources and the model WaterGAP2.2d. We found that the SMDAI and QDAI values are broadly similar to values of purely anomaly-based indices. However, the deficit anomaly indices provide more differentiated, spatial and temporal patterns that help to distinguish the degree of the actual drought hazard to vegetation health or the water supply. QDAI can be made relevant for stakeholders with different perceptions about the importance of ecosystem protection, by adapting the approach for computing the amount of water that is required to remain in the river for the well being of the river ecosystem. Both deficit anomaly indices are well suited for inclusion in local or global drought risk studies.
An accelerating Brewer-Dobson circulation (BDC) is a robust signal of climate change in model predictions but has been questioned by trace gas observations. We analyze stratospheric mean age of air and the full age spectrum as measures for the BDC and its trend. Age of air is calculated with the Chemical Lagrangian Model of the Stratosphere (CLaMS) driven by ERA-Interim, JRA-55 and MERRA-2 reanalysis data to assess the robustness of the representation of the BDC in current generation meteorological reanalyses. We find that climatological mean age significantly depends on the reanalysis, with JRA-55 showing the youngest and MERRA-2 the oldest mean age. Consideration of the age spectrum indicates that the older age for MERRA-2 is related to a stronger spectrum tail, likely related to weaker tropical upwelling and stronger recirculation. Seasonality of stratospheric transport is robustly represented in reanalyses, with similar mean age variations and age spectrum peaks. Long-term changes over 1989–2015 turn out to be similar for the reanalyses with mainly decreasing mean age accompanied by a shift of the age spectrum peak towards shorter transit times, resembling the forced response in climate model simulations to increasing greenhouse gas concentrations. For the shorter periods 1989–2001 and 2002–2015 age of air changes are less robust. Only ERA-Interim shows the hemispheric dipole pattern in age changes during 2002–2015 as viewed by recent satellite observations. Consequently, the representation of decadal variability of the BDC in current generation reanalyses appears less robust and a major uncertainty of modelling the BDC.
This paper investigates the global stratospheric Brewer–Dobson circulation (BDC) in the ERA5 meteorological reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMWF). The analysis is based on simulations of stratospheric mean age of air, including the full age spectrum, with the Lagrangian transport model CLaMS (Chemical Lagrangian Model of the Stratosphere), driven by reanalysis winds and total diabatic heating rates. ERA5-based results are compared to results based on the preceding ERA-Interim reanalysis. Our results show a significantly slower BDC for ERA5 than for ERA-Interim, manifesting in weaker diabatic heating rates and higher age of air. In the tropical lower stratosphere, heating rates are 30 %–40 % weaker in ERA5, likely correcting a bias in ERA-Interim. At 20 km and in the Northern Hemisphere (NH) stratosphere, ERA5 age values are around the upper margin of the uncertainty range from historical tracer observations, indicating a somewhat slow–biased BDC. The age trend in ERA5 over the 1989–2018 period is negative throughout the stratosphere, as climate models predict in response to global warming. However, the age decrease is not linear but steplike, potentially caused by multi-annual variability or changes in the observations included in the assimilation. During the 2002–2012 period, the ERA5 age shows a similar hemispheric dipole trend pattern as ERA-Interim, with age increasing in the NH and decreasing in the Southern Hemisphere (SH). Shifts in the age spectrum peak and residual circulation transit times indicate that reanalysis differences in age are likely caused by differences in the residual circulation. In particular, the shallow BDC branch accelerates in both reanalyses, whereas the deep branch accelerates in ERA5 and decelerates in ERA-Interim.
This paper presents an analysis of the recent tropospheric molecular hydrogen (H2) budget with a particular focus on soil uptake and European surface emissions. A variational inversion scheme is combined with observations from the RAMCES and EUROHYDROS atmospheric networks, which include continuous measurements performed between mid-2006 and mid-2009. Net H2 surface flux, then deposition velocity and surface emissions and finally, deposition velocity, biomass burning, anthropogenic and N2 fixation-related emissions were simultaneously inverted in several scenarios. These scenarios have focused on the sensibility of the soil uptake value to different spatio-temporal distributions. The range of variations of these diverse inversion sets generate an estimate of the uncertainty for each term of the H2 budget. The net H2 flux per region (High Northern Hemisphere, Tropics and High Southern Hemisphere) varies between −8 and +8 Tg yr−1. The best inversion in terms of fit to the observations combines updated prior surface emissions and a soil deposition velocity map that is based on bottom-up and top-down estimations. Our estimate of global H2 soil uptake is −59±9 Tg yr−1. Forty per cent of this uptake is located in the High Northern Hemisphere and 55% is located in the Tropics. In terms of surface emissions, seasonality is mainly driven by biomass burning emissions. The inferred European anthropogenic emissions are consistent with independent H2 emissions estimated using a H2/CO mass ratio of 0.034 and CO emissions within the range of their respective uncertainties. Additional constraints, such as isotopic measurements would be needed to infer a more robust partition of H2 sources and sinks.
We present an estimation of the uptake coefficient (γ) and yield of nitryl chloride (ClNO2) (f) for the heterogeneous processing of dinitrogen pentoxide (N2O5) using simultaneous measurements of particle and trace gas composition at a semi-rural, non-coastal, mountain site in the summer of 2011. The yield of ClNO2 varied between (0.035 ± 0.027) and (1.38 ± 0.60) with a campaign average of (0.49 ± 0.35). The large variability in f reflects the highly variable chloride content of particles at the site. Uptake coefficients were also highly variable with minimum, maximum and average γ values of 0.004, 0.11 and 0.028 ± 0.029, respectively, with no significant correlation with particle composition, but a weak dependence on relative humidity. The uptake coefficients obtained are compared to existing parameterizations based on laboratory datasets and with other values obtained by analysis of field data.
We present an estimation of the uptake coefficient (γ) and yield of nitryl chloride (ClNO2) (f) for the heterogeneous processing of dinitrogen pentoxide (N2O5) using simultaneous measurements of particle and trace gas composition at a semi-rural, non-coastal, mountain site in the summer of 2011. The yield of ClNO2 varied between (0.035 ± 0.027) and (1.38 ± 0.60) with a campaign average of (0.49 ± 0.35). The large variability in f reflects the highly variable chloride content of particles at the site. Uptake coefficients were also highly variable with minimum, maximum and average γ values of 0.004, 0.11 and 0.028 ± 0.029, respectively, with no significant correlation with particle composition, but a weak dependence on relative humidity. The uptake coefficients obtained are compared to existing parameterisations based on laboratory datasets and with other values obtained by analysis of field data.
Fire is the primary disturbance factor in many terrestrial ecosystems. Wildfire alters vegetation structure and composition, affects carbon storage and biogeochemical cycling, and results in the release of climatically relevant trace gases including CO2, CO, CH4, NOx, and aerosols. One way of assessing the impacts of global wildfire on centennial to multi-millennial timescales is to use process-based fire models linked to dynamic global vegetation models (DGVMs). Here we present an update to the LPJ-DGVM and a new fire module based on SPITFIRE that includes several improvements to the way in which fire occurrence, behaviour, and the effects of fire on vegetation are simulated. The new LPJ-LMfire model includes explicit calculation of natural ignitions, the representation of multi-day burning and coalescence of fires, and the calculation of rates of spread in different vegetation types. We describe a new representation of anthropogenic biomass burning under preindustrial conditions that distinguishes the different relationships between humans and fire among hunter-gatherers, pastoralists, and farmers. We evaluate our model simulations against remote-sensing-based estimates of burned area at regional and global scale. While wildfire in much of the modern world is largely influenced by anthropogenic suppression and ignitions, in those parts of the world where natural fire is still the dominant process (e.g. in remote areas of the boreal forest and subarctic), our results demonstrate a significant improvement in simulated burned area over the original SPITFIRE. The new fire model we present here is particularly suited for the investigation of climate–human–fire relationships on multi-millennial timescales prior to the Industrial Revolution.
Fire is the primary disturbance factor in many terrestrial ecosystems. Wildfire alters vegetation structure and composition, affects carbon storage and biogeochemical cycling, and results in the release of climatically relevant trace gases, including CO2, CO, CH4, NO and aerosols. Assessing the impacts of global wildfire on centennial to multimillennial timescales requires the linkage of process-based fire modeling with vegetation modeling using Dynamic Global Vegetation Models (DGVMs). Here we present a new fire module, SPITFIRE-2, and an update to the LPJ-DGVM that includes major improvements to the way in which fire occurrence, behavior, and the effect of fire on vegetation is simulated. The new fire module includes explicit calculation of natural ignitions, the representation of multi-day burning and coalescence of fires and the calculation of rates of spread in different vegetation types, as well as a simple scheme to model crown fires. We describe a new representation of anthropogenic biomass burning under preindustrial conditions that distinguishes the way in which the relationship between humans and fire are different between hunter-gatherers, obligate pastoralists, and farmers. Where and when available, we evaluate our model simulations against remote-sensing based estimates of burned area. While wildfire in much of the modern world is largely influenced by anthropogenic suppression and ignitions, in those parts of the world where natural fire is still the dominant process, e.g. in remote areas of the boreal forest, our results demonstrate a significant improvement in simulated burned area over previous models. With its unique properties of being able to simulate preindustrialfire, the new module we present here is particularly well suited for the investigation of climate-human-fire relationships on multi-millennial timescales.