550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (671)
- Doctoral Thesis (116)
- Conference Proceeding (14)
- Working Paper (11)
- Book (6)
- Part of a Book (4)
- Master's Thesis (3)
- Preprint (2)
- Review (2)
- Bachelor Thesis (1)
Language
- English (831) (remove)
Has Fulltext
- yes (831)
Is part of the Bibliography
- no (831)
Keywords
- climate change (11)
- Climate change (7)
- COSMO-CLM (6)
- precipitation (6)
- Atmospheric chemistry (5)
- Palaeoclimate (5)
- Biogeochemistry (4)
- Geochemistry (4)
- Palaeoceanography (4)
- Bewässerung (3)
Institute
- Geowissenschaften (616)
- Geowissenschaften / Geographie (125)
- Biodiversität und Klima Forschungszentrum (BiK-F) (61)
- Senckenbergische Naturforschende Gesellschaft (45)
- Geographie (22)
- Biowissenschaften (21)
- Institut für Ökologie, Evolution und Diversität (6)
- Institut für sozial-ökologische Forschung (ISOE) (5)
- Frankfurt Institute for Advanced Studies (FIAS) (4)
- Extern (3)
Small-scale thermal upwellings under the northern East African Rift from S travel time tomography
(2016)
There is a long-standing debate over how many and what types of plumes underlie the East African Rift and whether they do or do not drive its extension and consequent magmatism and seismicity. Here we present a new tomographic study of relative teleseismic S and SKS residuals that expands the resolution from previous regional studies below the northern East African Rift to image structure from the surface to the base of the transition zone. The images reveal two low-velocity clusters, below Afar and west of the Main Ethiopian Rift, that extend throughout the upper mantle and comprise several smaller-scale (about 100 km diameter), low-velocity features. These structures support those of our recent P tomographic study below the region. The relative magnitude of S to P residuals is around 3.5, which is consistent with a predominantly thermal nature of the anomalies. The S and P velocity anomalies in the low-velocity clusters can be explained by similar excess temperatures in the range of 100–200°C, consistent with temperatures inferred from other seismic, geochemical, and petrological studies. Somewhat stronger VS anomalies below Afar than west of the Main Ethiopian Rift may include an expression of volatiles and/or melt in this region. These results, together with a comparison with previous larger-scale tomographic models, indicate that these structures are likely small-scale upwellings with mild excess temperatures, rising from a regional thermal boundary layer at the base of the upper mantle.
Metasomatic reaction zones between mafic and ultramafic rocks exhumed from subduction zones provide a window into mass-transfer processes at high pressure. However, accurate interpretation of the rock record requires distinguishing high-pressure metasomatic processes from inherited oceanic signatures prior to subduction. We integrated constraints from bulk-rock geochemical compositions and petrophysical properties, mineral chemistry, and thermodynamic modeling to understand the formation of reaction zones between juxtaposed metagabbro and serpentinite as exemplified by the Voltri Massif (Ligurian Alps, Italy). Distinct zones of variably metasomatized metagabbro are dominated by chlorite, amphibole, clinopyroxene, epidote, rutile, ilmenite, and titanite between serpentinite and eclogitic metagabbro. Whereas the precursor serpentinite and oxide gabbro formed and were likely already in contact in an oceanic setting, the reaction zones formed by diffusional Mg-metasomatism between the two rocks from prograde to peak, to retrograde conditions in a subduction zone. Metasomatism of mafic rocks by Mg-rich fluids that previously equilibrated with serpentinite could be widespread along the subduction interface, within the subducted slab, and the mantle wedge. Furthermore, the models predict that talc formation by Si-metasomatism of serpentinite in subduction zones is limited by pressure-dependent increase in the silica activity buffered by the serpentine-talc equilibrium. Elevated activities of aqueous Ca and Al species would also favor the formation of chlorite and garnet. Accordingly, unusual conditions or processes would be required to stabilize abundant talc at high P-T conditions. Alternatively, a different set of mineral assemblages, such as serpentine- or chlorite-rich rocks, may be controlling the coupling-decoupling transition of the plate interface.
Talc formation via silica-metasomatism of ultramafic rocks is believed to play key roles in subduction zone processes. Yet, the conditions of talc formation remain poorly constrained. We used thermodynamic reaction-path models to assess the formation of talc at the slab-mantle interface and show that it is restricted to a limited set of pressure–temperature conditions, protolith, and fluid compositions. In contrast, our models predict that chlorite formation is ubiquitous at conditions relevant to the slab-mantle interface of subduction zones. The scarcity of talc and abundance of chlorite is evident in the rock record of exhumed subduction zone terranes. Talc formation during Si-metasomatism may thus play a more limited role in volatile cycling, strain localization, and in controlling the decoupling-coupling transition of the plate interface. Conversely, the observed and predicted ubiquity of chlorite corroborates its prominent role in slab-mantle interface processes that previous studies attributed to talc.
Key Points:
Limited talc formation by Si-metasomatism of ultramafic rocks in subduction zones
Chlorite formation is likely pervasive at the slab-mantle interface
Preferential formation of chlorite has wide-ranging chemical and physical implications for subduction zone processes
Plain Language Summary: In subduction zones, talc can form during chemical reactions of mantle rocks with silica-enriched fluids at the interface between descending oceanic plates and the overriding mantle. Its formation and distribution in subduction zones are believed to affect the volatile budget, rheological properties, and the down-dip limit of the decoupling of the slab-mantle interface. Therefore, illuminating the conditions that facilitate talc formation at high pressure-temperature conditions is key in assessing its roles in fundamental subduction zone processes. Using thermodynamic reaction-path models, we show that the formation of talc at the slab-mantle interface is restricted to a limited set of environmental conditions, because its formation is highly sensitive to the compositions of the mantle rocks and reactant fluids. Contrary to common belief, talc is unlikely to form in high abundance in ultramafic rocks metasomatized by Si-rich slab-derived fluids. Rather, our models predict the ubiquitous formation of chlorite along with other silicate minerals during Si-metasomatism due to the competing effects from other dissolved components that favor their formation over talc. This study calls into question the importance of talc during Si-metasomatism in subduction zones but highlights the more predominant role of chlorite.
The mechanisms of transfer of crustal material from the subducting slab to the overlying mantle wedge are still debated. Mélange rocks, formed by mixing of sediments, oceanic crust, and ultramafics along the slab-mantle interface, are predicted to ascend as diapirs from the slab-top and transfer their compositional signatures to the source region of arc magmas. However, the compositions of melts that result from the interaction of mélanges with a peridotite wedge remain unknown. Here we present experimental evidence that melting of peridotite hybridized by mélanges produces melts that carry the major and trace element abundances observed in natural arc magmas. We propose that differences in nature and relative contributions of mélanges hybridizing the mantle produce a range of primary arc magmas, from tholeiitic to calc-alkaline. Thus, assimilation of mélanges into the wedge may play a key role in transferring subduction signatures from the slab to the source of arc magmas.
Terrestrial climate and ecosystem evolution during ‘Greenhouse Earth’ phases of the early Paleogene remain incompletely known. Particularly, paleobotanical records from high southern latitudes are giving only limited insights into the Paleocene and early Eocene vegetation of the region. Hence, data from continuous well-calibrated sequences are required to make progress with the reconstruction of terrestrial climate and ecosystem dynamics from the southern latitudes during the early Paleogene.
In order to elucidate the terrestrial conditions from the high southern latitudes during the early Paleogene, terrestrial palynology was applied in the present study to two well-dated deep-marine sediment cores located at the Australo-Antarctic region: (i) IODP Site U1356 (Wilkes Land margin, East Antarctica) and (ii) ODP Site 1172 (East Tasman Plateau, southwest Pacific Ocean). The studied sequence from IODP Site U1356 comprises mid-shelfal sediments from the early to middle Eocene (53.9 – 46 million years ago [Ma]). For the ODP Site 1172, the studied succession is characterized by sediments deposited in shallow marine environments of the middle Paleocene to the early Eocene (60.7 – 54.2 Ma).
Based on the obtained pollen and spores (sporomorphs) results from the studied sequences of Site U1356 and Site 1172, this study aims to: (1) decipher the terrestrial climate conditions along the Australo-Antarctic region from the middle Paleocene to the middle Eocene; (2) evaluate the structure, diversity and compositional patterns of forests that throve in the Australo-Antarctic region during the early Paleogene; (3) understand the response of forests from the high southern latitudes to the climate dynamics from the early Paleogene; (4) establish a connection between the generated terrestrial palynomorph data and published Sea Surface Temperatures (SSTs) from the same cores.
To decipher the terrestrial climatic conditions on the Australo-Antarctic region, this study relies on the nearest living relative (NLR) concept that assumes that fossil taxa have similar climate requirements as their modern counterparts. This approach was applied to the sporomorph results of Site U1356 and Site 1172, following mainly the bioclimatic analysis. With regard to the structure and diversity patterns of the vegetation from the same region, the present study presents combined qualitative (i.e., reconstruction of the vegetation based mainly on the habitats of the known living relatives) and quantitative (i.e., application of ordination techniques, rarefaction and diversity indices) analyses of the fossil sporomorphs results.
The overall results from the paleoclimatic and vegetation reconstruction approaches applied in the present study, indicate that temperate and paratropical forests during the early Paleogene throve under different climatic conditions on the Wilkes Land margin and on Tasmania, at paleolatitudes of ∼70°S and ∼65°S, respectively.
Specifically, the sporomorph results from Site U1356, suggest that a highly diverse forest similar to present-day forests from New Caledonia was thriving on Antarctica during the early Eocene (53.9 – 51.9 Ma). These forests were characterized by the presence of termophilous taxa that are restricted today to tropical and subtropical settings, notably Bombacoideae, Strasburgeria, Beauprea, Spathiphyllum, Anacolosa and Lygodium. In combination with MBT/CBT paleotemperature results, they provide strong evidence for near-tropical warmth at least in the coastal lowlands along the Wilkes Land margin. The coeval presence of frost tolerant taxa such as Nothofagus, Araucariaceae and Podocarpaceae during the early Eocene on the same record suggests that paratropical forests were thriving along the Wilkes Land margin. Due to the presence of this kind of vegetation, it is possible to suggest that forests in this region were subject to a climatic gradient related to differences in elevation and/or the proximity to the coastline.
By the middle Eocene, the paratropical forests that characterized the vegetation of the early Eocene on the Wilkes Land margin were replaced by low diversity temperate forests dominated by Nothofagus, and similar to present-day cool-temperate forests from New Zealand. The dominance of these forests and the absence of thermophilous elements together with the lower temperatures suggested by the MBT/CBT and the sporomorph-based temperatures indicate consistently cooler conditions during this time interval.
With regard to the sporomorph results of Site 1172, this study suggests that three vegetation types were thriving on Tasmania from the middle Paleocene to the early Eocene under different climatic conditions. During the middle to late Paleocene, warm-temperate forests dominated by Podocarpaceae and Araucariaceae were the prevailing vegetation on Tasmania. The dominance of these forests was interrupted by the transient predominance of cool-temperate forests dominated by Nothofagus and Araucariaceae across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma). This cool-temperate forest was characterized by a lack of frost-sensitive elements (i.e., palms and cycads) indicating cooler conditions with harsher winters on Tasmania during this time interval. By the early Eocene, and linked with the Paleocene Eocene Thermal Maximum (PETM), Paleocene temperate forests dominated by gymnosperms were replaced by paratropical rainforests with the remarkable presence of the tropical mangrove palm Nypa during the PETM and the earliest Eocene. The overall results from Site U1356 and Site 1172, provide a new assessment of the terrestrial climatic conditions in the Australo-Antarctic region for validating climate models and understanding the response of high-latitude terrestrial ecosystems to the climate dynamics of the early Paleogene on southern latitudes.
The climatic conditions in the higher latitudes during the early Paleogene were further unravelled by comparing the obtained terrestrial and marine results. The integration of the obtained sporomorph data with previously published TEX86-based SSTs from Site 1172 documents that the vegetation dynamics were closely linked with the temperature evolution from the Australo-Antarctic region. Moreover, the comparison of TEX86-based SSTs and sporomorph-based climatic estimations from Site 1172 suggests a warm-season bias of both calibrations of TEX86 (i.e., TEX86Hand TEX86H), when this proxy is applied to high southern latitudes records of the early Paleogene.
Global warming, changes in the hydrological cycle and enhanced marine primary productivity all have been invoked to have contributed to the occurrence of widespread ocean anoxia during the Cenomanian-Turonian Oceanic Anoxic Event (OAE2; ~ 94 Ma), but disentangling these factors on a regional scale has remained problematic. We generated palynological and organic geochemical records that allow the separation of these forcing factors in a core spanning the OAE2 from Wunstorf, Lower Saxony Basin (LSB; North Gemany), which exhibits cyclic black shale–marl alternations related to the orbital precession cycle.
Despite the widely varying depositional conditions complicating the interpretation of the obtained records, TEX86H indicates that sea-surface temperature (SST) evolution in the LSB during OAE2 resembles that of previously studied sites throughout the proto-North Atlantic. Cooling during the so-called Plenus Cold Event interrupted black shale deposition during the early stages of OAE2. However, TEX86 does not vary significantly across marl–black shale alternations, suggesting that temperature variations did not force the formation of the cyclic black shale horizons. Relative (i.e., with respect to marine palynomorphs) and absolute abundances of pollen and spores are elevated during phases of black shale deposition, indicative of enhanced precipitation and run-off. High abundances of cysts from inferred heterotrophic and euryhaline dinoflagellates supports high run-off, which likely introduced additional nutrients to the epicontinental shelf resulting in elevated marine primary productivity.
We conclude that orbitally-forced enhanced precipitation and run-off, in tandem with elevated marine primary productivity, were critical in cyclic black shale formation on the northwest European epicontinental shelf and potentially for other OAE2 sections in the proto-Atlantic and Western Interior Seaway at similar latitudes as well.
Reconstructing the early Paleogene climate dynamics of terrestrial settings in the high southern latitudes is important to assess the role of high-latitude physical and biogeochemical processes in the global climate system. However, whereas a number of high-quality Paleogene climate records has become available for the marine realm of the high southern latitudes over the recent past, the long-term evolution of coeval terrestrial climates and ecosystems is yet poorly known. We here explore the climate and vegetation dynamics on Tasmania from the middle Paleocene to the early Eocene (60.7–54.2 Ma) based on a sporomorph record from Ocean Drilling Program (ODP) Site 1172 on the East Tasman Plateau. Our results show that three distinctly different vegetation types thrived on Tasmania under a high-precipitation regime during the middle Paleocene to early Eocene, with each type representing different temperature conditions: (i) warm-temperate forests dominated by gymnosperms that were dominant during the middle and late Paleocene; (ii) cool-temperate forests dominated by southern beech (Nothofagus) and araucarians across the middle/late Paleocene transition interval (~59.5 to ~59.0 Ma); and (iii) paratropical forests rich in ferns that were established during and in the wake of the Paleocene–Eocene Thermal Maximum (PETM). The transient establishment of cool-temperate forests lacking any frost-sensitive elements (i.e., palms and cycads) across the middle/late Paleocene transition interval indicates markedly cooler conditions, with the occurrence of frosts in winter, on Tasmania during that time. The integration of our sporomorph data with previously published TEX86-based sea-surface temperatures from ODP Site 1172 documents that the vegetation dynamics on Tasmania were closely linked with the temperature evolution in the Tasman sector of the Southwest Pacific region. Moreover, the comparison of our season-specific climate estimates for the sporomorph assemblages from ODP Site 1172 with the TEX86L- and TEX86H-based temperature data suggests a warm-season bias of both calibrations for the early Paleogene of the high southern latitudes.
A new, two-channel instrument for simultaneous NO3 and N2O5 monitoring was used to make the first comprehensive set of nocturnal NOx measurements (NO, NO2, NO3 and N2O5) at the Taunus Observatory, a rural mountain site (Kleiner Feldberg) in South-western Germany. In May 2008, NO3 and N2O5 mixing ratios were well above the instrumental detection limit (a few ppt) on all nights of the campaign and were characterised by large variability resulting from inhomogeneously distributed sinks. The concentrations of NO3, N2O5 and NO2 were consistent with the equilibrium constant, K2, defining the rates of formation and thermal dissociation of N2O5. A steady-state lifetime analysis showed that nocturnal NOx losses were generally dominated by reaction of NO3 with volatile organic compounds in this forested region, with N2O5 uptake to aerosols of secondary importance. Analysis of a limited dataset obtained at high relative humidity indicated that the loss of N2O5 by reaction with water vapour is less efficient (> factor 3) than derived using laboratory kinetic data. The fraction of NOx present as NO3 and N2O5 reached ≈20% on some nights, with night-time losses of NOx competing with daytime losses.
A new, two-channel instrument for simultaneous NO3 and N2O5 monitoring was used to make the first comprehensive set of nocturnal NOx measurements (NO, NO2, NO3 and N2O5) at the Taunus Observatory, a rural mountain site (Kleiner Feldberg) in South-western Germany. In May 2008, NO3 and N2O5 mixing ratios were well above the instrumental detection limit (a few ppt) on all nights of the campaign and were characterised by large variability. The concentrations of NO3, N2O5 and NO2 were consistent with the equilibrium constant, K2, defining the rates of formation and thermal dissociation of N2O5. A steady-state lifetime analysis is consistent with the loss of nocturnal NOx being dominated by the reaction of NO3 with volatile organic compounds in this forested region, with N2O5 uptake to aerosols of secondary importance. Analysis of a limited dataset obtained at high relative humidity indicated that the loss of N2O5 by reaction with water vapour is less efficient (>factor 3) than derived using laboratory kinetic data. The fraction of NOx present as NO3 and N2O5 reached ~20% on some nights, with night-time losses of NOx competing with daytime losses.
We combine structural analysis of fractures with 22 U–Pb dates measured in fracture-filling carbonate cements from bed-parallel fibrous calcite veins (beef), conjugated veins and faults within the Vaca Muerta Formation along the Andean fold and thrust belt in the Neuquén Basin. The measured ages constrain accurately the relationships between overpressures caused by hydrocarbon generation and Andean compression as mechanisms for natural fracturing and vein formation. Two generations of fibres have been identified in beef. The first one, consists of dark fibres from the inner zones, which are perpendicular to bedding and contain abundant cone-in-cone structures and hydrocarbon inclusions. U–Pb dating of these fibres yielded Early to Late Cretaceous ages from 116.7 ± 17.7 to 78.8 ± 10.2 Ma. The second generation of fibres corresponds to the outer zones and consists of white fibres oblique to bedding, indicating growth during layer-parallel shortening. Bed-perpendicular veins cutting beef yielded Late Cretaceous-late Palaeocene dates from 72.8 ± 22.4 to 60.9 ± 10.4 Ma. Eocene ages from 52.0 ± 2.9 to 42.2 ± 18.9 Ma were measured in bed-parallel slip surfaces and reverse and strike-slip faults, whereas Miocene dates from 13.9 ± 2.6 to 6.2 ± 1.1 Ma were measured in E-W calcite veins. U–Pb dating of veins, structural analysis of fractures and subsidence curves, indicate that beef inner zones formed in the oil window during burial of the Neuquén basin, and that tectonic stresses could enhance their formation. Beef outer zones and bed-perpendicular veins formed during E-W Late Cretaceous-late Palaeocene layer-parallel shortening. Contrarily, late Palaeocene-late Eocene bed-parallel slip surfaces and faults and Miocene E-W veins formed during NE-SW and E-W syn-to post-folding deformation, respectively. In both cases, syn-to post-folding compression occurred synchronously with forelandward migration of magmatic activity attributed to flat subduction of the Pacific slab beneath the Andes.
Non-forest ecosystems, dominated by shrubs, grasses and herbaceous plants, provide ecosystem services including carbon sequestration and forage for grazing, and are highly sensitive to climatic changes. Yet these ecosystems are poorly represented in remotely sensed biomass products and are undersampled by in situ monitoring. Current global change threats emphasize the need for new tools to capture biomass change in non-forest ecosystems at appropriate scales. Here we developed and deployed a new protocol for photogrammetric height using unoccupied aerial vehicle (UAV) images to test its capability for delivering standardized measurements of biomass across a globally distributed field experiment. We assessed whether canopy height inferred from UAV photogrammetry allows the prediction of aboveground biomass (AGB) across low-stature plant species by conducting 38 photogrammetric surveys over 741 harvested plots to sample 50 species. We found mean canopy height was strongly predictive of AGB across species, with a median adjusted R2 of 0.87 (ranging from 0.46 to 0.99) and median prediction error from leave-one-out cross-validation of 3.9%. Biomass per-unit-of-height was similar within but different among, plant functional types. We found that photogrammetric reconstructions of canopy height were sensitive to wind speed but not sun elevation during surveys. We demonstrated that our photogrammetric approach produced generalizable measurements across growth forms and environmental settings and yielded accuracies as good as those obtained from in situ approaches. We demonstrate that using a standardized approach for UAV photogrammetry can deliver accurate AGB estimates across a wide range of dynamic and heterogeneous ecosystems. Many academic and land management institutions have the technical capacity to deploy these approaches over extents of 1–10 ha−1. Photogrammetric approaches could provide much-needed information required to calibrate and validate the vegetation models and satellite-derived biomass products that are essential to understand vulnerable and understudied non-forested ecosystems around the globe.
Two types of particles exist in the atmosphere, primary and secondary particles. While primary particles such as soot, mineral dust, sea salt particles or pollen are introduced directly as particles into the atmosphere, secondary particles are formed in the atmosphere by condensation of gases. The formation of such new aerosol particles takes place frequently and at a broad variety of atmospheric conditions and geographic locations. A considerable fraction of the atmospheric particles is formed by such nucleation processes. The newly formed particles may grow by condensation to sizes where they are large enough to act as cloud condensation nuclei and therefore may affect cloud properties. The fundamental processes of aerosol nucleation are described and typical atmospheric observations are discussed. Two recent studies are introduced that potentially change our current understanding of atmospheric nucleation substantially.
Number concentrations of total and non-volatile aerosol particles with size diameters >0.01 µm as well as particle size distributions (0.4–23 µm diameter) were measured in situ in the Arctic lower stratosphere (10–20.5 km altitude). The measurements were obtained during the campaigns European Polar Stratospheric Cloud and Lee Wave Experiment (EUPLEX) and Envisat-Arctic-Validation (EAV). The campaigns were based in Kiruna, Sweden, and took place from January to March 2003. Measurements were conducted onboard the Russian high-altitude research aircraft Geophysica using the low-pressure Condensation Nucleus Counter COPAS (COndensation PArticle Counter System) and a modified FSSP 300 (Forward Scattering Spectrometer Probe). Around 18–20 km altitude typical total particle number concentrations nt range at 10–20 cm−3 (ambient conditions). Correlations with the trace gases nitrous oxide (N2O) and trichlorofluoromethane (CFC-11) are discussed. Inside the polar vortex the total number of particles >0.01 µm increases with potential temperature while N2O is decreasing which indicates a source of particles in the above polar stratosphere or mesosphere. A separate channel of the COPAS instrument measures the fraction of aerosol particles non-volatile at 250°C. Inside the polar vortex a much higher fraction of particles contained non-volatile residues than outside the vortex (~24% outside vortex). This is most likely due to a strongly increased fraction of meteoritic material in the particles which is transported downward from the mesosphere inside the polar vortex. The high fraction of non-volatile residual particles gives therefore experimental evidence for downward transport of mesospheric air inside the polar vortex. It is also shown that the fraction of non-volatile residual particles serves directly as a suitable experimental vortex tracer. Nanometer-sized meteoritic smoke particles may also serve as nuclei for the condensation of gaseous sulfuric acid and water in the polar vortex and these additional particles may be responsible for the increase in the observed particle concentration at low N2O. The number concentrations of particles >0.4 µm measured with the FSSP decrease markedly inside the polar vortex with increasing potential temperature, also a consequence of subsidence of air from higher altitudes inside the vortex. Another focus of the analysis was put on the particle measurements in the lowermost stratosphere. For the total particle density relatively high number concentrations of several hundred particles per cm3 at altitudes below ~14 km were observed in several flights. To investigate the origin of these high number concentrations we conducted air mass trajectory calculations and compared the particle measurements with other trace gas observations. The high number concentrations of total particles in the lowermost stratosphere are probably caused by transport of originally tropospheric air from lower latitudes and are potentially influenced by recent particle nucleation.
Number concentrations of total and non-volatile aerosol particles with size diameters >0.01 μm as well as particle size distributions (0.4–23 μm diameter) were measured in situ in the Arctic lower stratosphere (10–20.5 km altitude). The measurements were obtained during the campaigns European Polar Stratospheric Cloud and Lee Wave Experiment (EUPLEX) and Envisat-Arctic-Validation (EAV). The campaigns were based in Kiruna, Sweden, and took place from January to March 2003. Measurements were conducted onboard the Russian high-altitude research aircraft Geophysica using the low-pressure Condensation Nucleus Counter COPAS (COndensation PArticle Counter System) and a modified FSSP 300 (Forward Scattering Spectrometer Probe). Around 18–20 km altitude typical total particle number concentrations nt range at 10–20 cm−3 (ambient conditions). Correlations with the trace gases nitrous oxide (N2O) and trichlorofluoromethane (CFC-11) are discussed. Inside the polar vortex the total number of particles >0.01 μm increases with potential temperature while N2O is decreasing which indicates a source of particles in the above polar stratosphere or mesosphere. A separate channel of the COPAS instrument measures the fraction of aerosol particles non-volatile at 250°C. Inside the polar vortex a much higher fraction of particles contained non-volatile residues than outside the vortex (~67% inside vortex, ~24% outside vortex). This is most likely due to a strongly increased fraction of meteoric material in the particles which is transported downward from the mesosphere inside the polar vortex. The high fraction of non-volatile residual particles gives therefore experimental evidence for downward transport of mesospheric air inside the polar vortex. It is also shown that the fraction of non-volatile residual particles serves directly as a suitable experimental vortex tracer. Nanometer-sized meteoric smoke particles may also serve as nuclei for the condensation of gaseous sulfuric acid and water in the polar vortex and these additional particles may be responsible for the increase in the observed particle concentration at low N2O. The number concentrations of particles >0.4 μm measured with the FSSP decrease markedly inside the polar vortex with increasing potential temperature, also a consequence of subsidence of air from higher altitudes inside the vortex. Another focus of the analysis was put on the particle measurements in the lowermost stratosphere. For the total particle density relatively high number concentrations of several hundred particles per cm3 at altitudes below ~14 km were observed in several flights. To investigate the origin of these high number concentrations we conducted air mass trajectory calculations and compared the particle measurements with other trace gas observations. The high number concentrations of total particles in the lowermost stratosphere are probably caused by transport of originally tropospheric air from lower latitudes and are potentially influenced by recent particle nucleation.
This study reports and discusses new radiometric ages, petrographical and volcanological observations and whole rock geochemical data of the rocks of the Rudnik Mts. volcano-intrusive complex. The complex hosts a Pb-Zn-Ag deposit and belongs to the Serbo-Macedonian metallogenetic belt. Two distinct igneous events are distignuished. The first occurred >30 Ma and was mainly characterized by extrusive and shallow intrusive dacites and andesites and was unrelated to mineralization. The second igneous event occurred <23 Ma and was highly heterogeneous in terms of volcanic products and petrographic varieties, but with predominance of quartzlatites. The dacite-andesites (first event) and the quartzlatites (second event) are geochemically similar and display a calc-alkaline affinity and highly incompatible element enriched patterns on spider diagrams, but the younger quartzlatites are richer in K2O, Rb and Ba and poorer in Sr. This is taken as evidence that mixing between an ultrapotassic lamprophyre/lamproite magma and an acid calc-alkaline (dacite-like) magma was essential petrogenetic processes during the second event. The proposed simplified volcanological model suggests that this mixing was responsible for triggering strongly explosive volcanic activity as well as for providing conditions for active hydrothermal and mineralization processes. The observed link between a specific magmatic phase and ore deposit formation can be a general phennomenon in the Balkans, and must be addressed by further and more advanced studies.
This article presents an environmental remote sensing application using a UAV that is specifically aimed at reducing the data gap between field scale and satellite scale in soil erosion monitoring in Morocco. A fixed-wing aircraft type Sirius I (MAVinci, Germany) equipped with a digital system camera (Panasonic) is employed. UAV surveys are conducted over different study sites with varying extents and flying heights in order to provide both very high resolution site-specific data and lower-resolution overviews, thus fully exploiting the large potential of the chosen UAV for multi-scale mapping purposes. Depending on the scale and area coverage, two different approaches for georeferencing are used, based on high-precision GCPs or the UAV’s log file with exterior orientation values respectively. The photogrammetric image processing enables the creation of Digital Terrain Models (DTMs) and ortho-image mosaics with very high resolution on a sub-decimetre level. The created data products were used for quantifying gully and badland erosion in 2D and 3D as well as for the analysis of the surrounding areas and landscape development for larger extents.
This article presents a multiscale approach for detecting and monitoring soil erosion phenomena (i.e. gully erosion) in the agro-industrial area around the city of Taroudannt, Souss basin, Morocco. The study area is characterized as semi-arid with an annual average precipitation of 200 mm. Water scarcity, high population dynamics and changing land use towards huge areas of irrigation farming present numerous threats to sustainability. The agro-industry produces citrus fruits and vegetables in monocropping, mainly for the European market. Badland areas strongly affected by gully erosion border the agricultural areas as well as residential areas. To counteract the significant loss of land, land-leveling measures are attempted to create space for plantations and greenhouses. In order to develop sustainable approaches to limit gully growth the detection and monitoring of gully systems is fundamental. Specific gully sites are monitored with unmanned aerial vehicle (UAV) taking small-format aerial photographs (SFAP). This enables extremely high-resolution analysis (SFAP resolution: 2-10 cm) of the actual size of the gully channels as well as a detailed continued surveillance of their growth. Transferring the methodology on a larger scale using Quickbird satellite data (resolution: 60 cm) leads to the possibility of a large-scale analysis of the whole area around the city of Taroudannt (Area extent: ca. 350 km²). The results will then reveal possible relationships of gully growth and agro-industrial management and may even illustrate further interdependencies. The main objective is the identification of areas with high gully-erosion risk due to non-sustainable land use and the development of mitigation strategies for the study area.
Climatology of morphology and cloud-radiative properties of marine low-level mixed-phase clouds
(2023)
Marine stratocumuli cover about 40 - 60% of the ocean surface. They self-organize into different morphological regimes. The two organized cellular regimes are called open and closed mesoscale-cellular convective (MCC) clouds. In mid-to-high latitudes, open and closed cells are the two most frequent types of MCC clouds. In particular, many MCC clouds consist of a mixture of vapor, liquid droplets, and ice particles, referred to as mixed-phase clouds (MPCs). Even for the same cloud fraction, the albedo of open cells is, on average, lower than that of closed MCC clouds. Cloud phase and morphology individually influence the cloud radiative effect. Thus, this thesis investigates the relationships between the cloud phase, MCC organization, cell size, and differences regarding the cloud-radiative effect.
This thesis focuses on space-borne retrievals to achieve extensive temporal and spatial coverage. The liDAR-raDAR (DARDAR) version 2 product collocates two active and one passive satellite: CloudSat, Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), and Moderate Resolution Imaging Spectroradiometer (MODIS). The cloud phase of DARDAR is vertically integrated to establish a single cloud phase at each data point. The MCC classification data set based on the liquid water path (LWP) of MODIS scenes is collocated with the DARDAR product to determine the MCC organization. Cell-size statistics of both MCC clouds are obtained using a marker-based image segmentation method on MODIS reflectance scenes. In addition, based on MODIS reflectance scenes, a convolutional neural network (CNN) is developed to classify open and closed MCC scenes to avoid missing mature MPCs with a low LWP.
The first part of this thesis explores the relationships between cloud phase, morphology, and cloud albedo in the Southern Ocean (SO). At a given cloud-top temperature (CTT), seasonal changes in the mixed-phase fraction, defined as the number of MPCs divided by the sum of MPC and supercooled liquid cloud (SLC) pixels, are stronger than the morphological changes. Therefore, external factors seem to influence these changes instead of morphology. The dependence of cloud phase on cloud-top height (CTH) is more substantial than on CTT in clouds with CTHs below 2.5 km. The previously observed acceleration of closed-to-open transition in MPCs, known as preconditioning, is not the primary driver of climatological cloud morphology statistics in the SO. The morphological differences in cloud albedo are more pronounced in SLCs than in MPCs. This change in albedo alters the cloud radiative effect in the SO by 21Wm−2 to 39Wm−2 depending onseason and cloud phase.
Open and closed MCC clouds exhibit larger equivalent cell diameters in the MPCs than in SLCs in austral summer, whereas, in austral winter, the SLCs are larger. The cell’s aspect ratio accounts for varying CTHs. Closed cells have smaller aspect ratios than open cells, so their cell diameter is smaller, independent of CTH. While the seasonal differences in closed cells are due to changes in CTH, the seasonal aspect ratio differences in open cells are mainly caused by MPCs. With increasing aspect ratios, the cloud albedo decreases in both open and closed MCC clouds, with the most substantial decrease in open MPCs clouds. This leads to cloud-radiative changes of 60 - 75Wm−2 in the SO, depending on cloud phase and aspect ratio.
The established CNN exhibits a good accuracy of 80.6%, with even higher accuracies in the Open (85.5%) and Closed (87.3%) categories. The global MCC climatology based on the CNN generally agrees well with previous MCC distributions. The most notable difference occurs in the Northern Hemisphere (NH) in boreal winter, with a higher occurrence frequency of closed and open MCC clouds. This might indicate missing MPCs in previous studies based on the LWP and some restricted to warm cloud scenes. Thus, the developed CNN seems to better represent the different morphologies in MPCs than in previous classifications.
In conclusion, this thesis shows that understanding the dependencies of cloud phase, cloud morphology, and cell size is important to enhance predictions of the cloud-radiative effect and thus, it is important to evaluate how cloud phase, cloud morphology, and cellsize change in a warming climate.
Community trait assembly in highly diverse tropical rainforests is still poorly understood. Based on more than a decade of field measurements in a biodiversity hotspot of southern Ecuador, we implemented plant trait variation and improved soil organic matter dynamics in a widely used dynamic vegetation model (the Lund-Potsdam-Jena General Ecosystem Simulator, LPJ-GUESS) to explore the main drivers of community assembly along an elevational gradient. In the model used here (LPJ-GUESS-NTD, where NTD stands for nutrient-trait dynamics), each plant individual can possess different trait combinations, and the community trait composition emerges via ecological sorting. Further model developments include plant growth limitation by phosphorous (P) and mycorrhizal nutrient uptake. The new model version reproduced the main observed community trait shift and related vegetation processes along the elevational gradient, but only if nutrient limitations to plant growth were activated. In turn, when traits were fixed, low productivity communities emerged due to reduced nutrient-use efficiency. Mycorrhizal nutrient uptake, when deactivated, reduced net primary production (NPP) by 61–72% along the gradient. Our results strongly suggest that the elevational temperature gradient drives community assembly and ecosystem functioning indirectly through its effect on soil nutrient dynamics and vegetation traits. This illustrates the importance of considering these processes to yield realistic model predictions.
Projections of future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed. Projections of change from the baseline period (1981–2010) to the future (2070–2099) from a number of different ecosystems and hydrological models were studied. The differences between projections from the two types of model were looked at globally and regionally. Typically, across different regions the ecosystem models tended to project larger increases and smaller decreases in runoff than the hydrological models. However, the differences varied both regionally and seasonally. Sensitivity experiments were also used to investigate the contributions of varying CO2 and allowing vegetation distribution to evolve on projected changes in runoff. In two out of four models which had data available from CO2 sensitivity experiments, allowing CO2 to vary was found to increase runoff more than keeping CO2 constant, while in two models runoff decreased. This suggests more uncertainty in runoff responses to elevated CO2 than previously considered. As CO2 effects on evapotranspiration via stomatal conductance and leaf-area index are more commonly included in ecosystems models than in hydrological models, this may partially explain some of the difference between model types. Keeping the vegetation distribution static in JULES runs had much less effect on runoff projections than varying CO2, but this may be more pronounced if looked at over a longer timescale as vegetation changes may take longer to reach a new state.
Comparing projections of future changes in runoff from hydrological and biome models in ISI-MIP
(2013)
Future changes in runoff can have important implications for water resources and flooding. In this study, runoff projections from ISI-MIP (Inter-sectoral Impact Model Intercomparison Project) simulations forced with HadGEM2-ES bias-corrected climate data under the Representative Concentration Pathway 8.5 have been analysed for differences between impact models. Projections of change from a baseline period (1981–2010) to the future (2070–2099) from 12 impacts models which contributed to the hydrological and biomes sectors of ISI-MIP were studied. The biome models differed from the hydrological models by the inclusion of CO2 impacts and most also included a dynamic vegetation distribution. The biome and hydrological models agreed on the sign of runoff change for most regions of the world. However, in West Africa, the hydrological models projected drying, and the biome models a moistening. The biome models tended to produce larger increases and smaller decreases in regionally averaged runoff than the hydrological models, although there is large inter-model spread. The timing of runoff change was similar, but there were differences in magnitude, particularly at peak runoff. The impact of vegetation distribution change was much smaller than the projected change over time, while elevated CO2 had an effect as large as the magnitude of change over time projected by some models in some regions. The effect of CO2 on runoff was not consistent across the models, with two models showing increases and two decreases. There was also more spread in projections from the runs with elevated CO2 than with constant CO2. The biome models which gave increased runoff from elevated CO2 were also those which differed most from the hydrological models. Spatially, regions with most difference between model types tended to be projected to have most effect from elevated CO2, and seasonal differences were also similar, so elevated CO2 can partly explain the differences between hydrological and biome model runoff change projections. Therefore, this shows that a range of impact models should be considered to give the full range of uncertainty in impacts studies.
The objective of the present doctoral thesis was to investigate the occurrence, distribution, and behaviour of six hydrophilic ethers: ethyl tert-butyl ether (ETBE), 1,4-dioxane, ethylene glycol dimethyl ether (monoglyme), diethylene glycol dimethyl ether (diglyme), triethylene glycol dimethyl ether (triglyme), and tetraethylene glycol dimethyl ether (tetraglyme) in surface-, waste-, ground- and drinking water samples. Solid phase extraction and gas chromatography/mass spectrometry were used to analyze the six hydrophilic ethers. Altogether more than 150 surface water samples, almost 100 of each groundwater and wastewater samples, and 10 raw and drinking water samples were analyzed during the research project.
Initially, the method was validated in order to simultaneously determine the analytes of interest in various aquatic environments. A solid phase extraction method that uses coconut charcoal (Resprep® activated coconut charcoal, Restek) or carbon molecular sieve material (SupelcleanTM Envi-CarbTM Plus, Supelco) for analyte absorption were found suitable for determination of ETBE, 1,4-dioxane, and glymes in surface-, drinking-, ground- and wastewater samples. Precision and accuracy of both methods was demonstrated for all analytes of interest. The recovery of target compounds from the ultrapure water spiked at 1.0 µg L−1 was between 86.8 % and 98.2 %, with relative standard deviation below 6 %. The samples spiked at 10.0 µg L−1 gave slightly higher recovery of 90.6 % to 112.2 % with a relative standard deviation below 3.4 % for each analyte. Detection and quantification limits in ultrapure water and surface waters were furthermore established. The limit of quantitation (LOQ) in ultrapure water ranged between 0.024 µg L−1 to 0.057 µg L−1 using Restek cartridges, and 0.030 µg L−1 to 0.069 µg L−1 using Supelco cartridges. In the surface water samples the calculated LOQ was 0.032 µg L−1 to 0.067µg L−1 using coconut charcoal material and 0.032 µg L−1 to 0.052 µg L−1 using the carbon molecular sieve material. Moreover, stability of the unpreserved and preserved water samples as well as the extracts was determined. Preservation of samples with sodium bisulfate (at 1 gram per Liter) resulted in much better stability of the ethers in water samples. Subsequently, 27 samples obtained from seven surface water bodies in Germany (Rivers Rhine, Lippe, Main, Oder, Rur, Schwarzbach and Wesel-Datteln Canal) were analyzed for the six hydrophilic ethers. ETBE was present in only two surface waters (Rhine River and Wesel-Datteln Canal) with concentrations close to the LOQ (up to 0.065 µg L−1). 1,4-Dioxane was detected in all of the water samples at concentrations reaching 1.93 µg L–1. Monoglyme was identified only in the Main and Rhine Rivers at the maximum concentration of 0.114 µg L–1 and 0.427 µg L–1, respectively. Very high concentrations (up to 1.73 µg L−1) of diglyme, triglyme, and tetraglyme were detected in the samples from the Oder River. These glymes were also detected in the Rhine River; however the concentrations did not exceed 0.200 µg L–1. Furthermore, tetraglyme was detected in the Main River at an average concentration of 0.409 µg L–1 (n = 6) and in one sample from the Rur River at 0.192 µg L–1.
Four sampling campaigns were conducted at the Oderbruch polder between October 2009 and May 2012, in order to study the behavior of the hydrophilic ethers and organophosphates during riverbank filtration and in the anoxic aquifer. Moreover the suitability of these target compounds was assessed for their use as groundwater organic tracers. At the time of each sampling campaign, concentrations of triglyme and tetraglyme in the Oder River were between 20–185 ng L–1 (n = 4) and 273¬–1576 ng L–1 (n = 4). Monoglyme, diglyme, and 1,4-dioxane were analyzed only during the two last sampling campaigns. At that time, the concentration of diglyme in Oder River was 65¬–94 ng L-1 (n = 2) and 1,4-dioxane 1610¬–3290 ng L–1 (n = 2). In the drainage ditch, following bank filtration, concentrations of ethers ranged between 1090 ng L–1 and 1467 ng L–1 for 1,4-dioxane, 23¬ng L–1 and 41 ng L–1 for diglyme, 37 ng L–1 and 149 ng L–1 for triglyme, and 496 ng L–1 and 1403 ng L–1 for tetraglyme. In the anoxic aquifer, 1,4-dioxane showed the greatest persistence during the groundwater passage. At the distance of 1150 m from the river and an estimated groundwater age of 41.9 years, a concentration above 200 ng L−1 was detected. A positive correlation was found for the inorganic tracer chloride (Cl−) with 1,4-dioxane and tetraglyme. Similarities in the behavior of Cl− and the organic compound suggested that 1,4-dioxane and tetraglyme are controlled by the same hydraulic process and therefore can be used as additional tracers to study the dynamics of the groundwater system. These results show that high concentrations of ethers are present in the surface water and are not removed during bank filtration processes. Moreover, the hydrophilic ethers persist in the anoxic aquifer and little or no degradation is expected, supporting, their possible application as organic tracers.
A separate sampling project was conducted for 1,4-dioxane that focused primarily on its fate in the aquatic environment. This study provided missing information on the extent of water pollution with 1,4-dioxane is Germany. Numerous waste-, surface-, ground- and drinking water samples were collected in order to determine the persistence of 1,4-dioxane in the aquatic environment. The occurrence of 1,4-dioxane was determined in wastewater samples from four municipal sewage treatment plants (STP). The influent and effluent samples were collected during weekly campaigns. The average influent concentrations in all four plants ranged from 262 ± 32 ng L−1 to 834 ± 480 ng L−1, whereas the average effluents concentrations were between 267 ± 35 ng L−1 and 62,260 ± 36,000 ng L−1. The source of increased 1,4-dioxane concentrations in one of the effluents was identified to originate from impurities in the methanol used in the postanoxic denitrification process. Spatial and temporal distribution of 1,4-dioxane in the river Main, Rhine, and Oder was also examined. Concentrations reaching 2,200 ng L−1 in the Oder River, and 860 ng L−1 in both Main and Rhine River were detected. The average load during the sampling was estimated to be 6.5 kg d−1 in the Main, 34.1 kg d−1 in the Oder, and 134.5 kg d−1 in the Rhine River. In all of the sampled rivers, concentrations of 1,4-dioxane increased with distance from the mouth of the river and were found to negatively correlate with the discharge of the river. In order to determine if 1,4-dioxane can reach drinking water supplies, samples from a Rhine River bank filtration site and potable water from two drinking water production facilities were analyzed for the presence of 1,4-dioxane in the raw water and finished potable water. The raw water (following bank filtration) contained 650 ng L−1 to 670 ng L−1 of 1,4-dioxane, whereas the concentration in the finished drinking water fell only to 600 ng L−1 and 490 ng L−1, respectively.
During the final project, investigations of the source identification of high glyme concentrations in the Oder River were carried out. During four sampling campaigns between January, 2012 and April, 2013, 50 samples from the Oder River in the Oderbruch region and Poland were collected. During the first two samplings in the Oderbruch polder, glymes were detected at concentration reaching 0.07 µg L-1 (diglyme), 0.54 µg L−1 (triglyme) and 1.73 µg L−1 (tetraglyme) in the Oder River. The extensive sampling campaign of the Oder River (about 500 km) in Poland helped to identify the area of possible glyme entry into the river. During that sampling the maximum concentrations of triglyme and tetraglyme were 0.46 µg L−1 and 2.21 µg L−1, respectively. A closer investigation of the identified area of pollution, helped to determine the possible sources of glymes in the Oder River. Hence, the final sampling focused on the Kaczawa River, a left tributary of the Oder River and Czarna Woda, a left tributary of Kaczawa River. Moreover, samples from an industrial wastewater treatment plant were collected. Samples from Czarna Woda stream and Kaczawa River contained even higher concentrations of diglyme, triglyme, and tetraglyme, reaching 5.18 µg L−1, 12.87 µg L−1 and 80.81 µg L−1, respectively. Finally, three water samples from a wastewater treatment plant receiving influents from a copper smelter were analyzed. Diglyme, triglyme, and tetraglyme were present at an average concentration of 569 µg L−1, 4300 µg L−1, and 65900 µg L−1, respectively in the wastewater. Further research helped to identify the source of the glymes in the wastewater. The gas desulfurization process – Solinox implemented in the nearby copper smelter uses glymes as physical absorption medium for sulfur dioxide.
Results of this doctoral research provide important information about the occurrence, distribution, and behavior of hydrophilic ethers: 1,4-dioxane, monoglyme, diglyme, triglyme, and tetraglyme in the aquatic environment. A method capable of analyzing a wide range of ether compounds: from a volatile ETBE to a high molecular weight tetraglyme was validated. 1,4-Dioxane and tetraglyme were found to be applicable as organic tracers, since they are not easily attenuated during bank filtration and the anoxic groundwater passage. The extent of water pollution with 1,4-dioxane was shown in waste-, surface-, ground-, and drinking waters. One source of extremely high concentrations of 1,4-dioxane in a municipal sewage treatment plant applying postanoxic denitrification was identified, however more information is needed on the entry of 1,4-dioxane into surface waters. Moreover, 1,4-dioxane was present in drinking water samples from river bank filtration, which demonstrates its persistence in the aquatic environment and its low degradation potential during bank filtration and subsequent water treatment. Furthermore, this was the first study that focused primarily on identifying sources of glymes in surface waters. Glymes find a widespread use in industrial sectors, hence establishing their origin in the surface water is difficult (as with 1,4-dioxane). In this work, a gas desulphurization process was identified to be a dominating source of glyme pollution in the Oder River.
Remote sensing data are essential for monitoring the Earth’s surface waters, especially since the amount of publicly available in-situ data is declining. Satellite altimetry provides valuable information on the water levels and variations of lakes, reservoirs and rivers. In combination with satellite imagery, the derived time series allow the monitoring of lake storage changes and river discharge. However, satellite altimetry is limited in terms of its spatial resolution due to its measurement geometry, only providing information in the nadir direction beneath the satellite’s orbit. In a case study in the Mississippi River Basin (MRB), this study investigates the potential and limitations of past and current satellite missions for the monitoring of basin-wide storage changes. For that purpose, an automated target detection is developed and the extracted lake surfaces are merged with the satellites’ tracks. This reveals that the current altimeter configuration misses about 80% of all lakes larger than 0.1 km2 in the MRB and 20% of lakes larger than 10 km2, corresponding to 30% and 7% of the total water area, respectively. Past altimetry configurations perform even more poorly. From the larger water bodies represented by a global hydrology model, at least 91% of targets and 98% of storage changes are captured by the current altimeter configuration. This will improve significantly with the launch of the planned Surface Water and Ocean Topography (SWOT) mission.
In this study, we construct a new monthly zonal mean carbon dioxide (CO2) distribution from the upper troposphere to the stratosphere over the 2000–2010 time period. This reconstructed CO2 product is based on a Lagrangian backward trajectory model driven by ERA-Interim reanalysis meteorology and tropospheric CO2 measurements. Comparisons of our CO2 product to extratropical in situ measurements from aircraft transects and balloon profiles show remarkably good agreement. The main features of the CO2 distribution include (1) relatively large mixing ratios in the tropical stratosphere; (2) seasonal variability in the extratropics, with relatively high mixing ratios in the summer and autumn hemisphere in the 15–20 km altitude layer; and (3) decreasing mixing ratios with increasing altitude from the upper troposphere to the middle stratosphere ( ∼ 35 km). These features are consistent with expected variability due to the transport of long-lived trace gases by the stratospheric Brewer–Dobson circulation. The method used here to construct this CO2 product is unique from other modelling efforts and should be useful for model and satellite validation in the upper troposphere and stratosphere as a prior for inversion modelling and to analyse features of stratosphere–troposphere exchange as well as the stratospheric circulation and its variability.
In this study, we aim to reconstruct a relevant and new database of monthly zonal mean distribution of carbon dioxide (CO2) at global scale extending from the upper-troposphere (UT) to stratosphere (S). This product can be used for model and satellite validation in the UT/S, as a prior for inversion modelling and mainly to analyse a plausible feature of the stratospherictropospheric exchange as well as the stratospheric circulation and its variability. To do so, we investigate the ability of a Lagrangian trajectory model guided by ERA-Interim reanalysis to construct the CO2 abundance in the UT/S. From 10 year backward trajectories and tropospheric observations of CO2, we reconstruct upper-tropospheric and stratospheric CO2 over the period 2000–2010. The inter-comparisons of the reconstructed CO2 with mid-latitude vertical profiles measured by balloon samples as well as quasi-horizontal air samples from ER-2 aircraft during SOLVE and CONTRAIL campaigns exhibit a remarkable agreement. That demonstrates the potential of Lagrangian model to reconstruct CO2 in the UT/S. The zonal mean distribution exhibits relatively large CO2 in the tropical stratosphere due to the seasonal variation of the tropical upwelling of Brewer-Dobson circulation. During winter and spring, the tropical pipe is relatively isolated but is less confined during summer and autumn so that high CO2 values are more readily transported out of the tropics to the mid- and high latitude stratosphere. The shape of the vertical profiles suggests that relatively high CO2 above 20 km altitude mainly enter the stratosphere through tropical upwelling. CO2 mixing ratio is relatively low in the polar and tropical regions above 25 km. On average the CO2 mixing ratio decreases with altitude by 6-8 ppmv from the UT to stratosphere (e.g. up to 35 km) and is nearly constant with altitude.
The CLOUD (Cosmics Leaving OUtdoor Droplets) experiment at CERN is studying the nucleation and growth of aerosol particles under atmospheric conditions, and their activation into cloud droplets. A key feature of the CLOUD experiment is precise control of the experimental parameters. Temperature uniformity and stability in the chamber are important since many of the processes under study are sensitive to temperature and also to contaminants that can be released from the stainless steel walls by upward temperature fluctuations. The air enclosed within the 3 m CLOUD chamber is equipped with several arrays (strings) of high precision, fast-response thermometers to measure its temperature. Here we present a study of the air temperature uniformity inside the CLOUD chamber under various experimental conditions. Measurements were performed under calibration conditions and run conditions, which are distinguished by the flow rate of fresh air and trace gases entering the chamber: 20 l/min and up to 210 l/min, respectively. During steady-state calibration runs between −70 °C and +20 °C, the air temperature uniformity is better than +/−0.06 °C in the radial direction and +/−0.1 °C in the vertical direction. Larger non-uniformities are present during experimental runs, depending on the temperature control of the make-up air and trace gases (since some trace gases require elevated temperatures until injection into the chamber). The temperature stability is a few times 0.01 °C over periods of several hours during either calibration or steady-state run conditions. During rapid adiabatic expansions to activate cloud droplets and ice particles, the chamber walls are up to 10 °C warmer than the enclosed air. This results in larger non-uniformities while the air returns to its equilibrium temperature with time constant of about 200 s.
Projections from coarse-grid global circulation models are not suitable for regional estimates of water balance or trends of extreme precipitation and temperature, especially not in complex terrain. Thus, downscaling of global to regionally resolved projections is necessary to provide input to integrated water resources management approaches for river basins like the Upper Danube River Basin (UDRB) and the Upper Brahmaputra River Basin (UBRB).
This paper discusses the application of the regional climate model COSMO-CLM as a dynamical downscaling tool. To provide accurate data the COSMO-CLM model output was post-processed by statistical means. This downscaling chain performs well in the baseline period 1971 to 2000. However, COSMO-CLM performs better in the UDRB than in the UBRB because of a longer application experience and a less complex climate in Europe.
Different climate change scenarios were downscaled for the time period 1960–2100. The projections show an increase of temperature in both basins and for all seasons. The values are generally higher in the UBRB with the highest values occurring in the region of the Tibetan Plateau. Annual precipitation shows no substantial change. However, seasonal amounts show clear trends, for instance an increasing amount of spring precipitation in the UDRB. Again, the largest trends for different precipitation statistics are projected in the region of the Tibetan Plateau. Here, the projections show up to 50% longer dry periods in the months June to September with a simultaneous increase of about 10% for the maximum amount of precipitation on five consecutive days. For the Assam region in India, the projections also show an increase of 25% in the number of consecutive dry days during the monsoon season leading to prolonged monsoon breaks.
Melting inside earth is a common phenomenon and can be observed in many different regions where melt travels through the mantle and crust to eventually reach the surface where it crystallizes to build large volcanic provinces, whole stratigraphic layers of flood basalts, or even the oceanic crust. Often, melt reaching the surface is a good source of information. It can be used to achieve a better understanding about processes taking place in deeper regions inside the mantle and it is therefore essential to fundamentally understand melting and melt percolation processes. In order to achieve a deeper understanding, the aim of this thesis is to investigate processes that are connected to melting by using numerical models.
The physical model used is a so called two-phase flow model which describes the ability of melt to percolate through a viscously deforming, partially molten matrix. A famous feature of two-phase flow are solitary porosity waves, which are waves of locally higher porosity ascending through a partially molten background, keeping its shape constant, driven by decompaction and compaction of the solid matrix in front and behind the wave.
The viscosity law for shear- and volume viscosity was strongly simplified in most previous studies that modeled solitary waves. Often the porosity dependency is underestimated or its influence on the volume viscosity is even neglected, leading to too high viscosities. In this work more realistic laws are used that strongly decrease for small melt fractions. Those laws are incorporated into a 2D Finite Difference mantle convection code with two-phase flow to study the ascent of solitary porosity waves.
The model results show that an initial Gaussian-shaped wave rapidly evolves into a solitary wave with a certain amplitude, traveling upwards with constant velocity. Even though strongly weaker viscosities are used, the effect on dispersion curves and wave shape are only minor as long as the background porosity is rather small. The results are still in agreement to semi-analytical solutions which neglect shear stresses in the melt segregation equation. Higher background porosities and wave amplitudes lead to significant decrease in phase velocity and wave width, as the viscosity is strongly effected. However, the models show that solitary waves are still a possible mechanism for more realistic matrix viscosities.
While the ascending of porosity waves are mostly described by the movement of fluid melt, partially molten regions inside Earth trigger upwelling of both, solid and fluid phases, which can be called diapirism. While diapirs can have a wide variety of wavelengths, porosity waves are restricted to a few times the compaction length. The size of a melt perturbation in terms of compaction length therefore describes whether material is transported by diapirism or porosity waves. In this thesis we study the transition from diapiric rise to solitary porosity waves by systematically changing the size of a porosity perturbation from 1.8 to 120 times the compaction length. In case of a perturbation of the size of a few times the compaction length a single porosity wave will emerge, either with a positive or negative vertical matrix flux and if melt is not allowed to move relative to the matrix a diapir will emerge. In between these physical end members a regime can be observed where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is low compared to the Stokes velocity and the swarm of solitary waves will ascend jointly as a diapir, slowly elongating due to a higher amplitude main solitary wave.
Solitary waves will always emerge from a melt perturbation as long as two-phase flow is enabled, but the time for a solitary wave to emerge increases non-linearly with the perturbation radius in terms of compaction length. In nature, in many cases this time might be too long for solitary waves to emerge.
Another important feature when it comes to two-phase flow is the transport of trace elements in melt. Incompatible elements prefer to go into the melt, which eventually enriches the area where it crystallizes again. In order to model this redistribution, the code FDCON was extended to allow for fully consistent transport of elements in melt, including melting, freezing and re-equilibration with time. A 2D model, a simple representation of a volcanic back arc, is set up to investigate the behavior of trace elements. The influence of retention number and re-equilibration time is examined. Lava-lamp like convection can be observed in the lower part of the model, producing melt, that eventually leads to enrichment in trace elements in the upper high-viscous layer. The total enrichment in this layer approaches an asymptotic value and a 0D model is introduced to recreate this behavior.
In partially molten regions inside the earth melt buoyancy may trigger upwelling of both solid and fluid phases, i.e. diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation controls whether a diapir or a porosity wave will emerge. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1 to 100 times the compaction length. If the perturbation is much larger than a regular solitary wave, its Stokes velocity is large and therefore faster than the segregating melt. Consequently, the fluid is not able to form a porosity wave and a diapir emerges. For small perturbations solitary waves emerge, either with a positive or negative vertical matrix velocity inside. In between the diapir and solitary wave regimes we observe a third regime of solitary wave induced focusing of melt. In these cases, diapirism is dominant but the fluid is still fast enough to locally build up small solitary waves which rise slightly faster than the diapir and form finger like structures at the front of the diapir. In our numerical simulations the width of these fingers is controlled by the compaction length or the grid size, whichever is larger. In cases where the compaction length becomes similar to or smaller than the grid size the finger-like leading solitary porosity waves are no more properly resolved, and too big and too fast waves may be the result. Therefore, one should be careful in large scale two-phase flow modelling with melt focusing especially when compaction length and grid size are of similar order.
In partially molten regions inside the Earth, melt buoyancy may trigger upwelling of both solid and fluid phases, i.e., diapirism. If the melt is allowed to move separately with respect to the matrix, melt perturbations may evolve into solitary porosity waves. While diapirs may form on a wide range of scales, porosity waves are restricted to sizes of a few times the compaction length. Thus, the size of a partially molten perturbation in terms of compaction length controls whether material is dominantly transported by porosity waves or by diapirism. We study the transition from diapiric rise to solitary porosity waves by solving the two-phase flow equations of conservation of mass and momentum in 2D with porosity-dependent matrix viscosity. We systematically vary the initial size of a porosity perturbation from 1.8 to 120 times the compaction length. If the perturbation is of the order of a few compaction lengths, a single solitary wave will emerge, either with a positive or negative vertical matrix flux. If melt is not allowed to move separately to the matrix a diapir will emerge. In between these end members we observe a regime where the partially molten perturbation will split up into numerous solitary waves, whose phase velocity is so low compared to the Stokes velocity that the whole swarm of waves will ascend jointly as a diapir, just slowly elongating due to a higher amplitude main solitary wave. Only if the melt is not allowed to move separately to the matrix will no solitary waves build up, but as soon as two-phase flow is enabled solitary waves will eventually emerge. The required time to build them up increases nonlinearly with the perturbation radius in terms of compaction length and might be too long to allow for them in nature in many cases.
Motions on planetary spatial scales in the atmosphere are governed by the planetary geostrophic equations. However, little attention has been paid to the interaction between the baroclinic and barotropic flows within the planetary geostrophic scaling. This is the focus of the present study, which utilizes planetary geostrophic equations for a Boussinesq fluid supplemented by a novel evolution equation for the barotropic flow. The latter is affected by meridional momentum flux due to baroclinic flow and drag by the surface wind. The barotropic wind, on the other hand, affects the baroclinic flow through buoyancy advection. Via a relaxation towards a prescribed buoyancy profile the model produces realistic major features of the zonally symmetric wind and temperature fields. We show that there is considerable cancellation between the barotropic and the baroclinic surface zonal mean zonal winds. Linear and nonlinear model responses to steady diabatic zonally asymmetric forcing are investigated, and the arising stationary waves are interpreted in terms of analytical solutions. We also study the problem of baroclinic instability on the sphere within the present model.
Balloon-borne stratospheric BrO measurements : comparison with Envisat/SCIAMACHY BrO limb profiles
(2005)
For the first time, results of all four existing stratospheric BrO profiling instruments, are presented and compared with reference to the SLIMCAT 3-dimensional chemical transport model (3-D CTM). Model calculations are used to infer a BrO profile validation set, measured by 3 different balloon sensors, for the new Envisat/SCIAMACHY (ENVIronment SATellite/SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) satellite instrument. The balloon observations include (a) balloon-borne in situ resonance fluorescence detection of BrO, (b) balloon-borne solar occultation DOAS measurements (Differential Optical Absorption Spectroscopy) of BrO in the UV, and (c) BrO profiling from the solar occultation SAOZ (Systeme d'Analyse par Observation Zenithale) balloon instrument. Since stratospheric BrO is subject to considerable diurnal variation and none of the measurements are performed close enough in time and space for a direct comparison, all balloon observations are considered with reference to outputs from the 3-D CTM. The referencing is performed by forward and backward air mass trajectory calculations to match the balloon with the satellite observations. The diurnal variation of BrO is considered by 1-D photochemical model calculation along the trajectories. The 1-D photochemical model is initialised with output data of the 3-D model with additional constraints on the vertical transport, the total amount and photochemistry of stratospheric bromine as given by the various balloon observations. Total [Bry]=(20.1±2.8)pptv obtained from DOAS BrO observations at mid-latitudes in 2003, serves as an upper limit of the comparison. Most of the balloon observations agree with the photochemical model predictions within their given error estimates. First retrieval exercises of BrO limb profiling from the SCIAMACHY satellite instrument agree to <±50% with the photochemically-corrected balloon observations, and tend to show less agreement below 20 km.
Balloon-borne stratospheric BrO measurements: comparison with Envisat/SCIAMACHY BrO limb profiles
(2006)
For the first time, results of all four existing stratospheric BrO profiling instruments, are presented and compared with reference to the SLIMCAT 3-dimensional chemical transport model (3-D CTM). Model calculations are used to infer a BrO profile validation set, measured by 3 different balloon sensors, for the new Envisat/SCIAMACHY (ENVIronment SATellite/SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY) satellite instrument. The balloon observations include (a) balloon-borne in situ resonance fluorescence detection of BrO, (b) balloon-borne solar occultation DOAS measurements (Differential Optical Absorption Spectroscopy) of BrO in the UV, and (c) BrO profiling from the solar occultation SAOZ (Systeme d'Analyse par Observation Zenithale) balloon instrument. Since stratospheric BrO is subject to considerable diurnal variation and none of the measurements are performed close enough in time and space for a direct comparison, all balloon observations are considered with reference to outputs from the 3-D CTM. The referencing is performed by forward and backward air mass trajectory calculations to match the balloon with the satellite observations. The diurnal variation of BrO is considered by 1-D photochemical model calculation along the trajectories. The 1-D photochemical model is initialised with output data of the 3-D model with additional constraints on the vertical transport, the total amount and photochemistry of stratospheric bromine as given by the various balloon observations. Total [Bry]=(20.1±2.8)pptv obtained from DOAS BrO observations at mid-latitudes in 2003, serves as an upper limit of the comparison. Most of the balloon observations agree with the photochemical model predictions within their given error estimates. First retrieval exercises of BrO limb profiling from the SCIAMACHY satellite instrument agree to <±50% with the photochemically-corrected balloon observations, and tend to show less agreement below 20 km.
Pollen-based climate reconstructions were performed on two high-resolution pollen marines cores from the Alboran and Aegean Seas in order to unravel the climatic variability in the coastal settings of the Mediterranean region between 15 000 and 4000 years BP (the Lateglacial, and early to mid-Holocene). The quantitative climate reconstructions for the Alboran and Aegean Sea records focus mainly on the reconstruction of the seasonality changes (temperatures and precipitation), a crucial parameter in the Mediterranean region. This study is based on a multi-method approach comprising 3 methods: the Modern Analogues Technique (MAT), the recent Non-Metric Multidimensional Scaling/Generalized Additive Model method (NMDS/GAM) and Partial Least Squares regression (PLS). The climate signal inferred from this comparative approach confirms that cold and dry conditions prevailed in the Mediterranean region during the Oldest and Younger Dryas periods, while temperate conditions prevailed during the Bølling/Allerød and the Holocene. Our records suggest a West/East gradient of decreasing precipitation across the Mediterranean region during the cooler Late-glacial and early Holocene periods, similar to present-day conditions. Winter precipitation was highest during warm intervals and lowest during cooling phases. Several short-lived cool intervals (i.e. Older Dryas, another oscillation after this one (GI-1c2), Gerzensee/Preboreal Oscillations, 8.2 ka event, Bond events) connected to the North Atlantic climate system are documented in the Alboran and Aegean Sea records indicating that the climate oscillations associated with the successive steps of the deglaciation in the North Atlantic area occurred in both the western and eastern Mediterranean regions. This observation confirms the presence of strong climatic linkages between the North Atlantic and Mediterranean regions.
Pollen-based climate reconstructions were performed on two high-resolution pollen – marines cores from the Alboran and Aegean Seas in order to unravel the climatic variability in the coastal settings of the Mediterranean region between 15 000 and 4000 cal yrs BP (the Lateglacial, and early to mid-Holocene). The quantitative climate reconstructions for the Alboran and Aegean Sea records focus mainly on the reconstruction of the seasonality changes (temperatures and precipitation), a crucial parameter in the Mediterranean region. This study is based on a multi-method approach comprising 3 methods: the Modern Analogues Technique (MAT), the recent Non-Metric Multidimensional Scaling/Generalized Additive Model method (NMDS/GAM) and Partial Least Squares regression (PLS). The climate signal inferred from this comparative approach confirms that cold and dry conditions prevailed in the Mediterranean region during the Heinrich event 1 and Younger Dryas periods, while temperate conditions prevailed during the Bølling/Allerød and the Holocene. Our records suggest a West/East gradient of decreasing precipitation across the Mediterranean region during the cooler Late-glacial and early Holocene periods, similar to present-day conditions. Winter precipitation was highest during warm intervals and lowest during cooling phases. Several short-lived cool intervals (i.e., Older Dryas, another oscillation after this one (GI-1c2), Gerzensee/Preboreal Oscillations, 8.2 ka event, Bond events) connected to the North Atlantic climate system are documented in the Alboran and Aegean Sea records indicating that the climate oscillations associated with the successive steps of the deglaciation in the North Atlantic area occurred in both the western and eastern Mediterranean regions. This observation confirms the presence of strong climatic linkages between the North Atlantic and Mediterranean regions.
With the Earth system being about to leave Holocene conditions and thus the known safe operating space for humanity, frameworks such as the Planetary Boundaries (PBs) and the Sustainable Development Goals (SDGs) provide quantitative metrics to guide sustainability transformations. In order to strive, not only for compliance with the PBs but also for societal well-being, some approaches attempt to combine both PBs and SDGs within a single assessment. We focus on two prominent examples, the “Doughnut” by Kate Raworth and the #SDGinPB project of the 2018 report to the Club of Rome, which are not only aimed at public outreach, but also at guiding policy-making. To meet these objectives, the approaches should possess a certain accuracy in determining the progress in achieving the SDGs and in complying with the PBs. We evaluate, by using a multi-indicator approach for comparison, whether both approaches’ limited set of indicators can still represent the SDGs’ complexity. This comparative approach estimates the progress in achieving SDGs, especially in the Global North, to be significantly lower. Based on these results and against the approaches’ purposes, we discuss their simplifications and at which point the results are no longer reliable. We conclude that global assessments can be an important factor in initiating transformative processes by stimulating public discourse, but that the actual implementation of these would require approaches with greater recognition of local particularities.
In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multimillion-year decrease (i.e., Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating hydrographic restriction via South Atlantic–Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological timescales (> 1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (> 6 %) and relatively low (∼ 3.5 %) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive intervals before and after OAE 1a. In all cases, these short-term intervals of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we show that fluctuations in weathering-derived nutrient input from the southern African continent, linked to changes in orbitally driven humidity and aridity, were the underlying drivers of repetitive episodes of enhanced organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short-term fluctuations in riverine nutrient fluxes. We explain this relationship using the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
Extensive black shale deposits formed in the Early Cretaceous South Atlantic, supporting the notion that this emerging ocean basin was a globally important site of organic carbon burial. The magnitude of organic carbon burial in marine basins is known to be controlled by various tectonic, oceanographic, hydrological, and climatic processes acting on different temporal and spatial scales, the nature and relative importance of which are poorly understood for the young South Atlantic. Here we present new bulk and molecular geochemical data from an Aptian–Albian sediment record recovered from the deep Cape Basin at Deep Sea Drilling Project (DSDP) Site 361, which we combine with general circulation model results to identify driving mechanisms of organic carbon burial. A multi-million year decrease (i.e. Early Aptian–Albian) in organic carbon burial, reflected in a lithological succession of black shale, gray shale, and red beds, was caused by increasing bottom water oxygenation due to abating tectonic restriction via South Atlantic-Southern Ocean gateways. These results emphasize basin evolution and ocean gateway development as a decisive primary control on enhanced organic carbon preservation in the Cape Basin at geological time scales (>1 Myr). The Early Aptian black shale sequence comprises alternations of shales with high (>5%) and relatively low (~3%) organic carbon content of marine sources, the former being deposited during the global Oceanic Anoxic Event (OAE) 1a, as well as during repetitive events before and after OAE 1a. In all cases, these short-term events of enhanced organic carbon burial coincided with strong influxes of sediments derived from the proximal African continent, indicating closely coupled climate–land–ocean interactions. Supported by our model results, we propose that fluctuations in weathering-derived nutrient input from the southern African continent, linked to fluctuations in pCO2 and/or orbitally driven humidity/aridity, were the underlying drivers of short-term organic carbon burial in the deep Cape Basin. These results suggest that deep marine environments of emerging ocean basins responded sensitively and directly to short term fluctuations in riverine nutrient fluxes. We explain this relationship by the lack of wide and mature continental shelf seas that could have acted as a barrier or filter for nutrient transfer from the continent into the deep ocean.
A hygroscopicity tandem differential mobility analyzer (HTDMA) was used to measure the water uptake (hygroscopicity) of secondary organic aerosol (SOA) formed during the chemical and photochemical oxidation of several organic precursors in a smog chamber. Electron ionization mass spectra of the non-refractory submicron aerosol were simultaneously determined with an aerosol mass spectrometer (AMS), and correlations between the two different signals were investigated. SOA hygroscopicity was found to strongly correlate with the relative abundance of the ion signal m/z 44 expressed as a fraction of total organic signal (f44). m/z 44 is due mostly to the ion fragment CO2+ for all types of SOA systems studied, and has been previously shown to strongly correlate with organic O/C for ambient and chamber OA. The analysis was also performed on ambient OA from two field experiments at the remote site Jungfraujoch, and the megacity Mexico City, where similar results were found. A simple empirical linear relation between the hygroscopicity of OA at subsaturated RH, as given by the hygroscopic growth factor (GF) or "ϰorg" parameter, and f44 was determined and is given by ϰorg = 2.2 × f44 − 0.13. This approximation can be further verified and refined as the database for AMS and HTDMA measurements is constantly being expanded around the world. The use of this approximation could introduce an important simplification in the parameterization of hygroscopicity of OA in atmospheric models, since f44 is correlated with the photochemical age of an air mass.
A hygroscopicity tandem differential mobility analyzer (HTDMA) was used to measure the water uptake (hygroscopicity) of secondary organic aerosol (SOA) formed during the chemical and photochemical oxidation of several organic precursors in a smog chamber. Electron ionization mass spectra of the non-refractory submicron aerosol were simultaneously determined with an aerosol mass spectrometer (AMS), and correlations between the two different signals were investigated. SOA hygroscopicity was found to strongly correlate with the relative abundance of the ion signal m/z 44 expressed as a fraction of total organic signal (f44). m/z 44 is due mostly to the ion fragment CO2+ for all types of SOA systems studied, and has been previously shown to strongly correlate with organic O/C for ambient and chamber OA. The analysis was also performed on ambient OA from two field experiments at the remote site Jungfraujoch, and the megacity Mexico City, where similar results were found. A simple empirical linear relation between the hygroscopicity of OA at subsaturated RH, as given by the hygroscopic growth factor (GF) or "κorg" parameter, and f44 was determined and is given by κorg=2.2×f44−0.13. This approximation can be further verified and refined as the database for AMS and HTDMA measurements is constantly being expanded around the world. The use of this approximation could introduce an important simplification in the parameterization of hygroscopicity of OA in atmospheric models, since f44 is correlated with the photochemical age of an air mass.
During a 4-week run in October–November 2006, a pilot experiment was performed at the CERN Proton Synchrotron in preparation for the CLOUD1 experiment, whose aim is to study the possible influence of cosmic rays on clouds. The purpose of the pilot experiment was firstly to carry out exploratory measurements of the effect of ionising particle radiation on aerosol formation from trace H2SO4 vapour and secondly to provide technical input for the CLOUD design. A total of 44 nucleation bursts were produced and recorded, with formation rates of particles above the 3 nm detection threshold of between 0.1 and 100 cm−3s−1, and growth rates between 2 and 37 nm h−1. The corresponding H2SO4 concentrations were typically around 106 cm−3 or less. The experimentally-measured formation rates and H2SO4 concentrations are comparable to those found in the atmosphere, supporting the idea that sulphuric acid is involved in the nucleation of atmospheric aerosols. However, sulphuric acid alone is not able to explain the observed rapid growth rates, which suggests the presence of additional trace vapours in the aerosol chamber, whose identity is unknown. By analysing the charged fraction, a few of the aerosol bursts appear to have a contribution from ion-induced nucleation and ion-ion recombination to form neutral clusters. Some indications were also found for the accelerator beam timing and intensity to influence the aerosol particle formation rate at the highest experimental SO2 concentrations of 6 ppb, although none was found at lower concentrations. Overall, the exploratory measurements provide suggestive evidence for ion-induced nucleation or ion-ion recombination as sources of aerosol particles. However in order to quantify the conditions under which ion processes become significant, improvements are needed in controlling the experimental variables and in the reproducibility of the experiments. Finally, concerning technical aspects, the most important lessons for the CLOUD design include the stringent requirement of internal cleanliness of the aerosol chamber, as well as maintenance of extremely stable temperatures (variations below 0.1°C).
During a 4-week run in October–November 2006, a pilot experiment was performed at the CERN Proton Synchrotron in preparation for the Cosmics Leaving OUtdoor Droplets (CLOUD) experiment, whose aim is to study the possible influence of cosmic rays on clouds. The purpose of the pilot experiment was firstly to carry out exploratory measurements of the effect of ionising particle radiation on aerosol formation from trace H2SO4 vapour and secondly to provide technical input for the CLOUD design. A total of 44 nucleation bursts were produced and recorded, with formation rates of particles above the 3 nm detection threshold of between 0.1 and 100 cm -3 s -1, and growth rates between 2 and 37 nm h -1. The corresponding H2O concentrations were typically around 106 cm -3 or less. The experimentally-measured formation rates and htwosofour concentrations are comparable to those found in the atmosphere, supporting the idea that sulphuric acid is involved in the nucleation of atmospheric aerosols. However, sulphuric acid alone is not able to explain the observed rapid growth rates, which suggests the presence of additional trace vapours in the aerosol chamber, whose identity is unknown. By analysing the charged fraction, a few of the aerosol bursts appear to have a contribution from ion-induced nucleation and ion-ion recombination to form neutral clusters. Some indications were also found for the accelerator beam timing and intensity to influence the aerosol particle formation rate at the highest experimental SO2 concentrations of 6 ppb, although none was found at lower concentrations. Overall, the exploratory measurements provide suggestive evidence for ion-induced nucleation or ion-ion recombination as sources of aerosol particles. However in order to quantify the conditions under which ion processes become significant, improvements are needed in controlling the experimental variables and in the reproducibility of the experiments. Finally, concerning technical aspects, the most important lessons for the CLOUD design include the stringent requirement of internal cleanliness of the aerosol chamber, as well as maintenance of extremely stable temperatures (variations below 0.1 °C)
The present work was devised to address the systematic analysis of samples from a range of Roman non-ferrous metal artefacts from different archaeological contexts and sites in the Roman provinces of Germania Superior. One of the focal points of this study is the provenancing of different lead objects from five important Roman settlements between 15 BC and the beginning of fourth century AD. For this purpose, measurements were made on lead and copper ore samples from the Siegerland, Eifel, Hunsrück and Lahn-Dill area in Germany and supplemented with data from the literature to create a data bank of lead isotope ratios of European deposits. Compositional analysis of lead objects by Electron Microprobe analysis showed that Romans were able to purify lead from ore up to 99%. Multi-Collector Inductively Coupled Plasma Mass-Spectrometry was used to determine the source of lead, which played an important role in nearly all aspects of Roman life. Lead isotope ratios were measured for ore samples from German deposits from the eastern side of the Rhine (Siegerland, Lahn-Dill, Ems) and the western side of the Rhine (Eifel, Hunsrück), which contained enough ore reserves to answer the increasing local demand and are believed to have been mined during the Roman period. This data together with those from Mediterranean ore deposits from the literature was used to establish a data bank. The Mediterranean ore deposits range from Cambrian (high 207Pb/206Pb) to tertiary (lower 207Pb/206Pb) values. In particular, the Cypriot deposits are younger, while the Spanish deposits fall either with the younger Sardic ores or close to the older Cypriot ores. The lead isotope ratios of most German ore deposits fall in between the 208Pb/206Pb vs. 207Pb/206Pb ratios of Sardinia and Cyprus, where the lead isotope signature of ore deposits from France and Britain are also found. Over 240 lead objects were measured from Wallendorf (second century BC to first century AD) Dangstetten (15-8 BC), Waldgirmes (AD 1-10), Mainz (AD 1-300), Martberg (first to fourth centuries AD) & Trier (third to fourth centuries AD). Comparing the lead isotope ratios of lead objects and those from German ores shows that the source of over 85 percent of objects are Eifel ore deposits, but the Roman’s had also imported lead from the Southern Massif Central and from Great Britain. A further topic of this work was the systematic study of the variation of copper isotope ratios in different copper minerals and the mechanisms, which controls copper isotope fractionation in ores deposits. For this purpose, copper isotope analyses were made by Multi-Collector Inductively Coupled Plasma Mass-Spectrometry from a series of hydrothermal copper sulphides and their alteration products. Copper and lead isotope ratios were measured in coexisting phases of chalcopyrite and malachite and also coexisting malachite and azurite. No significant fractionation was observed in malachite-azurite phases, but in chalcopyrite-malachite coexisting phases, malachite always shows a positive fractionation to heavier isotope values. Zhu et al. and Larson et al. showed that isotopic variations in copper principally reflect mass fractionation in response to low temperature processes rather than source heterogeneity. The low temperature ore formation processes are mostly represented by weathering of primary sulphide ores to produce secondary carbonate phases and therefore are usually observed on the surface of ore deposits, which were probably removed during the early Bronze Age. Using this concept, copper isotope ratios were measured in some Early Bronze Age copper alloys and Roman copper alloys. However, no large copper isotope fractionation has been observed. Lead and copper isotope ratios were measured on samples from the Kupferschiefer. Two profiles were investigated; 1) Sangerhausen, which was not directly influenced by the oxidizing brines of Rote Fäule and 2) Oberkatz, where both Rote Fäule-controlled and structure-controlled mineralization were observed. Results from maturation studies of organic matter suggest the maximum temperature affecting the Kupferschiefer did not exceed 130°C. delta-65-Cu ranges between -0.78-+0.58‰, shows a positive correlation with copper concentration. Maximum temperature in the Kupferschiefer profile from Oberkatz is supposed to be around 150°C. delta-65Cu in this profile ranges between -0.71-+0.68‰. The pattern of copper isotope fractionation and copper concentration is same as the for profile of Sangerhausen. Origina lead isotope ratios are strongly overprinted by high concentrations of uranium in bottom of both profiles causing more radiogenic lead.
Large-scale hydrological modelling has become increasingly wide-spread during the last decade. An annual workshop series on large-scale hydrological modelling has provided, since 1997, a forum to the German-speaking community for discussing recent developments and achievements in this research area. In this paper we present the findings from the 2007 workshop which focused on advances and visions in large-scale hydrological modelling. We identify the state of the art, difficulties and research perspectives with respect to the themes "sensitivity of model results", "integrated modelling" and "coupling of processes in hydrosphere, atmosphere and biosphere". Some achievements in large-scale hydrological modelling during the last ten years are presented together with a selection of remaining challenges for the future.
Quantification of spatially and temporally resolved water flows and water storage variations for all land areas of the globe is required to assess water resources, water scarcity and flood hazards, and to understand the Earth system. This quantification is done with the help of global hydrological models (GHMs). What are the challenges and prospects in the development and application of GHMs? Seven important challenges are presented. (1) Data scarcity makes quantification of human water use difficult even though significant progress has been achieved in the last decade. (2) Uncertainty of meteorological input data strongly affects model outputs. (3) The reaction of vegetation to changing climate and CO2 concentrations is uncertain and not taken into account in most GHMs that serve to estimate climate change impacts. (4) Reasons for discrepant responses of GHMs to changing climate have yet to be identified. (5) More accurate estimates of monthly time series of water availability and use are needed to provide good indicators of water scarcity. (6) Integration of gradient-based groundwater modelling into GHMs is necessary for a better simulation of groundwater–surface water interactions and capillary rise. (7) Detection and attribution of human interference with freshwater systems by using GHMs are constrained by data of insufficient quality but also GHM uncertainty itself. Regarding prospects for progress, we propose to decrease the uncertainty of GHM output by making better use of in situ and remotely sensed observations of output variables such as river discharge or total water storage variations by multi-criteria validation, calibration or data assimilation. Finally, we present an initiative that works towards the vision of hyperresolution global hydrological modelling where GHM outputs would be provided at a 1-km resolution with reasonable accuracy.
Long-term average groundwater recharge, which is equivalent to renewable groundwater resources, is the major limiting factor for the sustainable use of groundwater. Compared to surface water resources, groundwater resources are more protected from pollution, and their use is less restricted by seasonal and inter-annual flow variations. To support water management in a globalized world, it is necessary to estimate groundwater recharge at the global scale. Here, we present a best estimate of global-scale long-term average diffuse groundwater recharge (i.e. renewable groundwater resources) that has been calculated by the most recent version of the WaterGAP Global Hydrology Model WGHM (spatial resolution of 0.5° by 0.5°, daily time steps). The estimate was obtained using two state-of-the-art global data sets of gridded observed precipitation that we corrected for measurement errors, which also allowed to quantify the uncertainty due to these equally uncertain data sets. The standard WGHM groundwater recharge algorithm was modified for semi-arid and arid regions, based on independent estimates of diffuse groundwater recharge, which lead to an unbiased estimation of groundwater recharge in these regions. WGHM was tuned against observed long-term average river discharge at 1235 gauging stations by adjusting, individually for each basin, the partitioning of precipitation into evapotranspiration and total runoff. We estimate that global groundwater recharge was 12 666 km3/yr for the climate normal 1961–1990, i.e. 32% of total renewable water resources. In semi-arid and arid regions, mountainous regions, permafrost regions and in the Asian Monsoon region, groundwater recharge accounts for a lower fraction of total runoff, which makes these regions particularly vulnerable to seasonal and inter-annual precipitation variability and water pollution. Average per-capita renewable groundwater resources of countries vary between 8 m3/(capita yr) for Egypt to more than 1 million m3/(capita yr) for the Falkland Islands, the global average in the year 2000 being 2091 m3/(capita yr). Regarding the uncertainty of estimated groundwater resources due to the two precipitation data sets, deviation from the mean is 1.1% for the global value, and less than 1% for 50 out of the 165 countries considered, between 1 and 5% for 62, between 5 and 20% for 43 and between 20 and 80% for 10 countries. Deviations at the grid scale can be much larger, ranging between 0 and 186 mm/yr.
Long-term average groundwater recharge, which is equivalent to renewable groundwater resources, is the major limiting factor for the sustainable use of groundwater. Compared to surface water resources, groundwater resources are more protected from pollution, and their use is less restricted by seasonal and inter-annual flow variations. To support water management in a globalized world, it is necessary to estimate groundwater recharge at the global scale. Here, we present a best estimate of global-scale long-term average diffuse groundwater recharge (i.e. renewable groundwater resources) that has been calculated by the most recent version of the WaterGAP Global Hydrology Model WGHM (spatial resolution of 0.5° by 0.5°, daily time steps). The estimate was obtained using two state-of-the art global data sets of gridded observed precipitation that we corrected for measurement errors, which also allowed to quantify the uncertainty due to these equally uncertain data sets. The standard WGHM groundwater recharge algorithm was modified for semi-arid and arid regions, based on 15 independent estimates of diffuse groundwater recharge, which lead to an unbiased estimation of groundwater recharge in these regions. WGHM was tuned against observed long-term average river discharge at 1235 gauging stations by adjusting, individually for each basin, the partitioning of precipitation into evapotranspiration and total runoff. We estimate that global groundwater recharge was 12 666 km3/yr for the climate nor20 mal 1961–1990, i.e. 32% of total renewable water resources. In semi-arid and arid regions, mountainous regions, permafrost regions and in the Asian Monsoon region, groundwater recharge accounts for a lower fraction of total runoff, which makes these regions particularly vulnerable to seasonal and inter-annual precipitation variability and water pollution. Average per-capita renewable groundwater resources of countries vary 25 between 8m3/(capita yr) for Egypt to more than 1 million m3/(capita yr) for the Falkland Islands, the global average in the year 2000 being 2091m3/(capita yr). Regarding the uncertainty of estimated groundwater resources due to the two precipitation data sets, deviation from the mean is less than 1% for 50 out of the 165 countries considered, between 1 and 5% for 62, between 5 and 20% for 43 and between 20 and 80% for 10 countries. Deviations at the grid scale can be much larger, ranging between 0 and 186 mm/yr.
Global-scale information on natural river flows and anthropogenic river flow alterations is required to identify areas where aqueous ecosystems are expected to be strongly degraded. Such information can support the identification of environmental flow guidelines and a sustainable water management that balances the water demands of humans and ecosystems. This study presents the first global assessment of the anthropogenic alteration of river flow regimes, in particular of flow variability, by water withdrawals and dams/reservoirs. Six ecologically relevant flow indicators were quantified using an improved version of the global water model WaterGAP. WaterGAP simulated, with a spatial resolution of 0.5 degree, river discharge as affected by human water withdrawals and dams around the year 2000, as well as naturalized discharge without this type of human interference. Compared to naturalized conditions, long-term average global discharge into oceans and internal sinks has decreased by 2.7% due to water withdrawals, and by 0.8% due to dams. Mainly due to irrigation, long-term average river discharge and statistical low flow Q90 (monthly river discharge that is exceeded in 9 out of 10 months) have decreased by more than 10% on one sixth and one quarter of the global land area (excluding Antarctica and Greenland), respectively. Q90 has increased significantly on only 5% of the land area, downstream of reservoirs. Due to both water withdrawals and reservoirs, seasonal flow amplitude has decreased significantly on one sixth of the land area, while interannual variability has increased on one quarter of the land area mainly due to irrigation. It has decreased on only 8% of the land area, in areas downstream of reservoirs where consumptive water use is low. The impact of reservoirs is likely underestimated by our study as small reservoirs are not taken into account. Areas most affected by anthropogenic river flow alterations are the Western and Central USA, Mexico, the western coast of South America, the Mediterranean rim, Southern Africa, the semi-arid and arid countries of the Near East and Western Asia, Pakistan and India, Northern China and the Australian Murray-Darling Basin, as well as some Arctic rivers. Due to a large number of uncertainties related e.g. to the estimation of water use and reservoir operation rules, the analysis is expected to provide only first estimates of river flow alterations that should be refined in the future.
Global-scale information on natural river flows and anthropogenic river flow alterations is required to identify areas where aqueous ecosystems are expected to be strongly degraded. Such information can support the identification of environmental flow guidelines and a sustainable water management that balances the water demands of humans and ecosystems. This study presents the first global assessment of the anthropogenic alteration of river flow regimes by water withdrawals and dams, focusing in particular on the change of flow variability. Six ecologically relevant flow indicators were quantified using an improved version of the global water model WaterGAP. WaterGAP simulated, with a spatial resolution of 0.5 degree, river discharge as affected by human water withdrawals and dams, as well as naturalized discharge without this type of human interference. Mainly due to irrigation, long-term average river discharge and statistical low flow Q90 (monthly river discharge that is exceeded in 9 out of 10 months) have decreased by more than 10% on one sixth and one quarter of the global land area (excluding Antarctica and Greenland), respectively. Q90 has increased significantly on only 5% of the land area, downstream of reservoirs. Due to both water withdrawals and dams, seasonal flow amplitude has decreased significantly on one sixth of the land area, while interannual variability has increased on one quarter of the land area mainly due to irrigation. It has decreased on only 8% of the land area, in areas with little consumptive water use that are downstream of dams. Areas most affected by anthropogenic river flow alterations are the western and central USA, Mexico, the western coast of South America, the Mediterranean rim, Southern Africa, the semi-arid and arid countries of the Near East and Western Asia, Pakistan and India, Northern China and the Australian Murray-Darling Basin, as well as some Arctic rivers. Due to a large number of uncertainties related e.g. to the estimation of water use and reservoir operation rules, the analysis is expected to provide only first estimates of river flow alterations that should be refined in the future.
Groundwater recharge is the major limiting factor for the sustainable use of groundwater. To support water management in a globalized world, it is necessary to estimate, in a spatially resolved way, global-scale groundwater recharge. In this report, improved model estimates of diffuse groundwater recharge at the global-scale, with a spatial resolution of 0.5° by 0.5°, are presented. They are based on calculations of the global hydrological model WGHM (WaterGAP Global Hydrology Model) which, for semi-arid and arid areas of the globe, was tuned against independent point estimates of diffuse groundwater recharge. This has led to a decrease of estimated groundwater recharge under semi-arid and arid conditions as compared to the model results before tuning, and the new estimates are more similar to country level data on groundwater recharge. Using the improved model, the impact of climate change on groundwater recharge was simulated, applying two greenhouse gas emissions scenarios as interpreted by two different climate models.
The Earth's future depends on how we manage the manifold risks of climate change (CC). It is state-of-the-art to assume that risk reduction requires participatory management involving a broad range of stakeholders and scientists. However, there is still little knowledge about the optimal design of participatory climate change risk management processes (PRMPs), in particular with respect to considering the multitude of substantial uncertainties that are relevant for PRMPs. To support the many local to regional PRMPs that are necessary for a successful global-scale reduction of CC risks, we present a roadmap for designing such transdisciplinary knowledge integration processes. The roadmap suggests ways in which uncertainties can be comprehensively addressed within a PRMP. We discuss the concept of CC risks and their management and propose an uncertainty framework that distinguishes epistemic, ontological, and linguistic uncertainty as well as ambiguity. Uncertainties relevant for CC risk management are identified. Communicative and modeling methods that support social learning as well as the development of risk management strategies are proposed for each of six phases of a PRMP. Finally, we recommend how to evaluate PRMPs as such evaluations and their publication are paramount for achieving a reduction of CC risks.
Wetlands such as bogs, swamps, or freshwater marshes are hotspots of biodiversity. For 5.1 million km2 of inland wetlands, the dynamics of area and water storage, which strongly impact biodiversity and ecosystem services, were simulated using the global hydrological model WaterGAP. For the first time, the impacts of both human water use and man‐made reservoirs (WUR) and future climate change (CC) on wetlands around the globe were quantified. WUR impacts are concentrated in arid/semiarid regions, where WUR decreased mean wetland water storage by more than 5% on 8.2% of the mean wetland area during 1986–2005 (Am), with highest decreases in groundwater depletion area. Using output of three climate models, CC impacts on wetlands were quantified, distinguishing unavoidable impacts [i.e., at 2 °C global warming (GW)] from avoidable impacts (difference between 3 °C and 2 °C impacts). Even unavoidable CC impacts are projected to be much larger than WUR impacts, also in arid/semiarid regions. On most wetland area with reliable estimates, avoidable CC impacts are more than twice as large as unavoidable impacts. In case of 2 °C GW, half of Am is estimated to be unaffected by mean storage changes of more than 5%, but only one third in case of 3 °C GW. Temporal variability of water storage will increase for most wetlands. Wetlands in dry regions will be affected the most, particularly by water storage decreases in the dry season. Different from wealthier countries, low‐income countries will dominantly suffer from a decrease in wetland water storage due to CC.
River flow regimes, including long-term average flows, seasonality, low flows, high flows and other types of flow variability, play an important role for freshwater ecosystems. Thus, climate change affects freshwater ecosystems not only by increased temperatures but also by altered river flow regimes. However, with one exception, transferable quantitative relations between flow alterations and ecological responses have not yet been derived. While discharge decreases are generally considered to be detrimental for ecosystems, the effect of future discharge increases is unclear. As a first step towards a global-scale analysis of climate change impacts on freshwater ecosystems, we quantified the impact of climate change on five ecologically relevant river flow indicators, using the global water model WaterGAP 2.1g to simulate monthly time series of river discharge with a spatial resolution of 0.5 degrees. Four climate change scenarios based on two global climate models and two greenhouse gas emissions scenarios were evaluated. We compared the impact of climate change by the 2050s to the impact of water withdrawals and dams on natural flow regimes that had occurred by 2002. Climate change was computed to alter seasonal flow regimes significantly (i.e. by more than 10%) on 90% of the global land area (excluding Greenland and Antarctica), as compared to only one quarter of the land area that had suffered from significant seasonal flow regime alterations due to dams and water withdrawals. Due to climate change, the timing of the maximum mean monthly river discharge will be shifted by at least one month on one third on the global land area, more often towards earlier months (mainly due to earlier snowmelt). Dams and withdrawals had caused comparable shifts on less than 5% of the land area only. Long-term average annual river discharge is predicted to significantly increase on one half of the land area, and to significantly decrease on one quarter. Dams and withdrawals had led to significant decreases on one sixth of the land area, and nowhere to increases. Thus, by the 2050s, climate change may have impacted ecologically relevant river flow characteristics more strongly than dams and water withdrawals have up to now. The only exception refers to the decrease of the statistical low flow Q90, with significant decreases both by past water withdrawals and future climate change on one quarter of the land area. However, dam impacts are likely underestimated by our study. Considering long-term average river discharge, only a few regions, including Spain, Italy, Iraq, Southern India, Western China, the Australian Murray Darling Basin and the High Plains Aquifer in the USA, all of them with extensive irrigation, are expected to be less affected by climate change than by past anthropogenic flow alterations. In some of these regions, climate change will exacerbate the discharge reductions, while in others climate change provides opportunities for reducing past reductions. Emissions scenario B2 leads to only slightly reduced alterations of river flow regimes as compared to scenario A2 even though emissions are much smaller. The differences in alterations resulting from the two applied climate models are larger than those resulting from the two emissions scenarios. Based on general knowledge about ecosystem responses to flow alterations and data related to flow alterations by dams and water withdrawals, we expect that the computed climate change induced river flow alterations will impact freshwater ecosystems more strongly than past anthropogenic alterations.
River flow regimes, including long-term average flows, seasonality, low flows, high flows and other types of flow variability, play an important role for freshwater ecosystems. Thus, climate change affects freshwater ecosystems not only by increased temperatures but also by altered river flow regimes. However, with one exception, transferable quantitative relations between flow alterations and ecosystem responses have not yet been derived. While discharge decreases are generally considered to be detrimental for ecosystems, the effect of future discharge increases is unclear. As a first step towards a global-scale analysis of climate change impacts on freshwater ecosystems, we quantified the impact of climate change on five ecologically relevant river flow indicators, using the global water model WaterGAP 2.1g to simulate monthly time series of river discharge with a spatial resolution of 0.5 degrees. Four climate change scenarios based on two global climate modelsand two greenhouse gas emissions scenarios were evaluated.
We compared the impact of climate change by the 2050s to the impact of water withdrawals and dams on natural flow regimes that had occurred by 2002. Climate change was computed to alter seasonal flow regimes significantly (i.e. by more than 10%) on 90% of the global land area (excluding Greenland and Antarctica), as compared to only one quarter of the land area that had suffered from significant seasonal flow regime alterations due to dams and water withdrawals. Due to climate change, the timing of the maximum mean monthly river discharge will be shifted by at least one month on one third on the global land area, more often towards earlier months (mainly due to earlier snowmelt). Dams and withdrawals had caused comparable shifts on less than 5% of the land area only. Long-term average annual river discharge is predicted to significantly increase on one half of the land area, and to significantly decrease on one quarter. Dams and withdrawals had led to significant decreases on one sixth of the land area, and nowhere to increases.
Thus, by the 2050s, climate change will have impacted ecologically relevant river flow characteristics much more strongly than dams and water withdrawals have up to now. The only exception refers to the decrease of the statistical low flow Q90, with significant decreases both by past water withdrawals and future climate change on one quarter of the land area. Considering long-term average river discharge, only a few regions, including Spain, Italy, Iraq, Southern India, Western China, the Australian Murray Darling Basin and the High Plains Aquifer in the USA, all of them with extensive irrigation, are expected to be less affected by climate change than by past anthropogenic flow alterations. In some of these regions, climate change will exacerbate the discharge reduction. Emissions scenario B2 leads to only slightly reduced alterations of river flow regimes as compared to scenario A2 even though emissions are much smaller. The differences in alterations resulting from the two applied climate models are larger than those resulting from the two emissions scenarios. Based on general knowledge about ecosystem responses to flow alterations and data related to flow alterations by dams and water withdrawals, we expect that the computed climate change induced river flow alterations will impact freshwater ecosystems more strongly than past anthropogenic alterations.
U–Pb age spectra of detrital zircons related to the East European Platform could be traced in paragneiss through the whole Mid-German-Crystalline Zone (Variscides, Central Europe) from the Odenwald via the Spessart to the Ruhla crystalline forming an exotic unit between Armorica and Laurussia. The depositional ages of the paragneiss are defined by the youngest age of the detrital zircons and the oldest intrusion ages as Ordovician to Silurian. The Ediacaran dominated age spectrum of detrital zircons from the paragneiss of the East Odenwald suggests the latter to be derived from the shelf of the East European Platform (Baltica), which was influenced by the 1.5 Ga old detritus delivered from a giant intrusion (Mazury granitoid, Poland). The detrital zircon age spectrum of the lower Palaeozoic paragneiss of the East Odenwald and sandstone of the northern Holy Cross Mountains are identical. The pure Sveconorwegian spectrum of the lower Palaeozoic quartzite from the Spessart, (Kirchner and Albert Int J Earth Sci 2020) and the Ruhla (Zeh and Gerdes Gondwana Res 17:254–263, 2010) could be sourced from Bornholm and southern Sweden. A U–Pb age spectrum with 88% Palaeozoic detrital zircons from a volcano-sedimentary rock of the East Odenwald is interpreted to be derived from a Silurian magmatic arc (46%), which was probably generated during the drift of the Mid-German-Crystalline Zone micro-continent to the south. A tentative plate tectonic model of Mid-German-Crystalline Zone is presented taking into account (a) the East European Platform related age spectra of the detrital zircons (b) the Ordovician to Silurian depositional age of the metasediments (c) the Silurian and Early Devonian intrusion age of the plutonic and volcanic rocks and (d) the U–Pb ages of the Middle Devonian high-grade metamorphism. The East European Platform-related part of the Mid-German-Crystalline Zone is interpreted as a micro-continent, which drifted through the Rheic Ocean to the south and collided with the Saxothuringian (Armorican Terrane Assemblage) during the Early Devonian. Such large-scale tectonic transport from the northern continent to the southern continent is also known from the SW Iberia, where Laurussia-related metasediments of the Rheic suture zone are explained by a large scale tectonic escape (Braid et al. J Geol Soc Lond 168:383–392, 2011).
U–Pb age spectra of detrital zircons related to the East European Platform could be traced in paragneiss through the whole Mid-German-Crystalline Zone (Variscides, Central Europe) from the Odenwald via the Spessart to the Ruhla crystalline forming an exotic unit between Armorica and Laurussia. The depositional ages of the paragneiss are defined by the youngest age of the detrital zircons and the oldest intrusion ages as Ordovician to Silurian. The Ediacaran dominated age spectrum of detrital zircons from the paragneiss of the East Odenwald suggests the latter to be derived from the shelf of the East European Platform (Baltica), which was influenced by the 1.5 Ga old detritus delivered from a giant intrusion (Mazury granitoid, Poland). The detrital zircon age spectrum of the lower Palaeozoic paragneiss of the East Odenwald and sandstone of the northern Holy Cross Mountains are identical. The pure Sveconorwegian spectrum of the lower Palaeozoic quartzite from the Spessart, (Kirchner and Albert Int J Earth Sci 2020) and the Ruhla (Zeh and Gerdes Gondwana Res 17:254–263, 2010) could be sourced from Bornholm and southern Sweden. A U–Pb age spectrum with 88% Palaeozoic detrital zircons from a volcano-sedimentary rock of the East Odenwald is interpreted to be derived from a Silurian magmatic arc (46%), which was probably generated during the drift of the Mid-German-Crystalline Zone micro-continent to the south. A tentative plate tectonic model of Mid-German-Crystalline Zone is presented taking into account (a) the East European Platform related age spectra of the detrital zircons (b) the Ordovician to Silurian depositional age of the metasediments (c) the Silurian and Early Devonian intrusion age of the plutonic and volcanic rocks and (d) the U–Pb ages of the Middle Devonian high-grade metamorphism. The East European Platform-related part of the Mid-German-Crystalline Zone is interpreted as a micro-continent, which drifted through the Rheic Ocean to the south and collided with the Saxothuringian (Armorican Terrane Assemblage) during the Early Devonian. Such large-scale tectonic transport from the northern continent to the southern continent is also known from the SW Iberia, where Laurussia-related metasediments of the Rheic suture zone are explained by a large scale tectonic escape (Braid et al. J Geol Soc Lond 168:383–392, 2011).
Profiles of CFC-11 (CCl3F) and CFC-12 (CCl2F2) of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aboard the European satellite Envisat have been retrieved from versions MIPAS/4.61 to MIPAS/4.62 and MIPAS/5.02 to MIPAS/5.06 level-1b data using the scientific level-2 processor run by Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK) and Consejo Superior de Investigaciones Científicas (CSIC), Instituto de Astrofísica de Andalucía (IAA). These profiles have been compared to measurements taken by the balloon-borne cryosampler, Mark IV (MkIV) and MIPAS-Balloon (MIPAS-B), the airborne MIPAS-STRatospheric aircraft (MIPAS-STR), the satellite-borne Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) and the High Resolution Dynamic Limb Sounder (HIRDLS), as well as the ground-based Halocarbon and other Atmospheric Trace Species (HATS) network for the reduced spectral resolution period (RR: January 2005–April 2012) of MIPAS. ACE-FTS, MkIV and HATS also provide measurements during the high spectral resolution period (full resolution, FR: July 2002–March 2004) and were used to validate MIPAS CFC-11 and CFC-12 products during that time, as well as profiles from the Improved Limb Atmospheric Spectrometer, ILAS-II. In general, we find that MIPAS shows slightly higher values for CFC-11 at the lower end of the profiles (below ∼ 15 km) and in a comparison of HATS ground-based data and MIPAS measurements at 3 km below the tropopause. Differences range from approximately 10 to 50 pptv ( ∼ 5–20 %) during the RR period. In general, differences are slightly smaller for the FR period. An indication of a slight high bias at the lower end of the profile exists for CFC-12 as well, but this bias is far less pronounced than for CFC-11 and is not as obvious in the relative differences between MIPAS and any of the comparison instruments. Differences at the lower end of the profile (below ∼ 15 km) and in the comparison of HATS and MIPAS measurements taken at 3 km below the tropopause mainly stay within 10–50 pptv (corresponding to ∼ 2–10 % for CFC-12) for the RR and the FR period. Between ∼ 15 and 30 km, most comparisons agree within 10–20 pptv (10–20 %), apart from ILAS-II, which shows large differences above ∼ 17 km. Overall, relative differences are usually smaller for CFC-12 than for CFC-11. For both species – CFC-11 and CFC-12 – we find that differences at the lower end of the profile tend to be larger at higher latitudes than in tropical and subtropical regions. In addition, MIPAS profiles have a maximum in their mixing ratio around the tropopause, which is most obvious in tropical mean profiles. Comparisons of the standard deviation in a quiescent atmosphere (polar summer) show that only the CFC-12 FR error budget can fully explain the observed variability, while for the other products (CFC-11 FR and RR and CFC-12 RR) only two-thirds to three-quarters can be explained. Investigations regarding the temporal stability show very small negative drifts in MIPAS CFC-11 measurements. These instrument drifts vary between ∼ 1 and 3 % decade−1. For CFC-12, the drifts are also negative and close to zero up to ∼ 30 km. Above that altitude, larger drifts of up to ∼ 50 % decade−1 appear which are negative up to ∼ 35 km and positive, but of a similar magnitude, above.
Profiles of CFC-11 (CCl3F) and CFC-12 (CCl2F2) of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aboard the European satellite Envisat have been retrieved from versions MIPAS/4.61 to MIPAS/4.62 and MIPAS/5.02 to MIPAS/5.06 level-1b data using the scientific level-2 processor run by Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK) and Consejo Superior de Investigaciones Científicas (CSIC), Instituto de Astrofísica de Andalucía (IAA). These profiles have been compared to measurements taken by the balloon-borne cryosampler, Mark IV (MkIV) and MIPAS-Balloon (MIPAS-B), the airborne MIPAS-STRatospheric aircraft (MIPAS-STR), the satellite-borne Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) and the High Resolution Dynamic Limb Sounder (HIRDLS), as well as the ground-based Halocarbon and other Atmospheric Trace Species (HATS) network for the reduced spectral resolution period (RR: January 2005–April 2012) of MIPAS. ACE-FTS, MkIV and HATS also provide measurements during the high spectral resolution period (full resolution, FR: July 2002–March 2004) and were used to validate MIPAS CFC-11 and CFC-12 products during that time, as well as profiles from the Improved Limb Atmospheric Spectrometer, ILAS-II. In general, we find that MIPAS shows slightly higher values for CFC-11 at the lower end of the profiles (below ∼ 15 km) and in a comparison of HATS ground-based data and MIPAS measurements at 3 km below the tropopause. Differences range from approximately 10 to 50 pptv ( ∼ 5–20 %) during the RR period. In general, differences are slightly smaller for the FR period. An indication of a slight high bias at the lower end of the profile exists for CFC-12 as well, but this bias is far less pronounced than for CFC-11 and is not as obvious in the relative differences between MIPAS and any of the comparison instruments. Differences at the lower end of the profile (below ∼ 15 km) and in the comparison of HATS and MIPAS measurements taken at 3 km below the tropopause mainly stay within 10–50 pptv (corresponding to ∼ 2–10 % for CFC-12) for the RR and the FR period. Between ∼ 15 and 30 km, most comparisons agree within 10–20 pptv (10–20 %), apart from ILAS-II, which shows large differences above ∼ 17 km. Overall, relative differences are usually smaller for CFC-12 than for CFC-11. For both species – CFC-11 and CFC-12 – we find that differences at the lower end of the profile tend to be larger at higher latitudes than in tropical and subtropical regions. In addition, MIPAS profiles have a maximum in their mixing ratio around the tropopause, which is most obvious in tropical mean profiles. Comparisons of the standard deviation in a quiescent atmosphere (polar summer) show that only the CFC-12 FR error budget can fully explain the observed variability, while for the other products (CFC-11 FR and RR and CFC-12 RR) only two-thirds to three-quarters can be explained. Investigations regarding the temporal stability show very small negative drifts in MIPAS CFC-11 measurements. These instrument drifts vary between ∼ 1 and 3 % decade−1. For CFC-12, the drifts are also negative and close to zero up to ∼ 30 km. Above that altitude, larger drifts of up to ∼ 50 % decade−1 appear which are negative up to ∼ 35 km and positive, but of a similar magnitude, above.
The SAWNUC microphysical aerosol nucleation model is used to study the effect of reactor walls on the interpretation of nucleation experiments with respect to nucleation theory. This work shows that loss processes, such as wall losses, influence the interpretation of nucleation experiments, especially at low growth rates and short lifetime of freshly nucleated particles. In these cases the power dependency of the formation rates, determined at a certain particle size, with respect to H2SO4 does not correspond to the approximate number of H2SO4 molecules in the critical cluster as expected by the first nucleation theorem. Observed ∂log(J)/∂log([H2SO4]) therefore can vary widely for identical nucleation conditions but different sink terms.
The SAWNUC (Sulphuric Acid Water NUCleation) microphysical aerosol nucleation model is used to study the effect of reactor walls on the interpretation of nucleation experiments with respect to nucleation theory. This work shows that loss processes, such as wall losses, influence the interpretation of nucleation experiments, especially at low growth rates and short lifetimes of freshly nucleated particles. In these cases the power dependency of the formation rates, determined at a certain particle size, with respect to H2SO4 does not correspond to the approximate number of H2SO4 molecules in the critical cluster as expected by the first nucleation theorem. Observed ∂log(J)/∂log([H2SO4]) therefore can vary widely for identical nucleation conditions but different sink terms.
Binary nucleation of sulphuric acid-water particles is expected to be an important process in the free troposphere at low temperatures. SAWNUC (Sulphuric Acid Water Nucleation) is a model of binary nucleation that is based on laboratory measurements of the binding energies of sulphuric acid and water in charged and neutral clusters. Predictions of SAWNUC are compared for the first time comprehensively with experimental binary nucleation data from the CLOUD chamber at European Organization for Nuclear Research. The experimental measurements span a temperature range of 208–292 K, sulphuric acid concentrations from 1·106 to 1·109 cm−3, and distinguish between ion-induced and neutral nucleation. Good agreement, within a factor of 5, is found between the experimental and modeled formation rates for ion-induced nucleation at 278 K and below and for neutral nucleation at 208 and 223 K. Differences at warm temperatures are attributed to ammonia contamination which was indicated by the presence of ammonia-sulphuric acid clusters, detected by an Atmospheric Pressure Interface Time of Flight (APi-TOF) mass spectrometer. APi-TOF measurements of the sulphuric acid ion cluster distributions (math formula with i = 0, 1, ..., 10) show qualitative agreement with the SAWNUC ion cluster distributions. Remaining differences between the measured and modeled distributions are most likely due to fragmentation in the APi-TOF. The CLOUD results are in good agreement with previously measured cluster binding energies and show the SAWNUC model to be a good representation of ion-induced and neutral binary nucleation of sulphuric acid-water clusters in the middle and upper troposphere.
Positive plant–plant interactions are thought to drive vegetation patterns in harsh environments, such as semi-arid areas. According to the stress-gradient hypothesis (SGH), the role of positive interactions between species (facilitation) is expected to increase with harshness, predicting associated variation in species composition along environmental gradients. However, the relation between stress and facilitation along environmental gradients is debated. Furthermore, differentiating facilitative interactions from other underlying mechanisms, such as microtopographic heterogeneity, is not trivial. We analysed the spatial co-occurrence relationships of vascular plant species that form patchy vegetation in arid lapilli fields (tephra) from recent volcanic eruptions on La Palma, Canary Islands. Assuming a harshness gradient negatively correlated with elevation because the lower elevations are more arid and water availability is considered the most limiting resource, and that an outcome of facilitation is plants co-occurring in the same patch, from the SGH we expected a greater degree of co-occurrence at lower elevation. We tested this at both the species and the individual plant level. We analysed the species composition of 1277 shrubby vegetation patches at 64 different sampling points, ranging from the coast to around 700 m a.s.l. Patch morphology and microtopographic heterogeneity variables were also measured, to account for their potential effects on the species composition of patches. We used generalized linear models and generalized mixed-effects models to analyse species richness, number of individuals in patches and percentage of patches with positive co-occurrences, and a pairwise co-occurrence analysis combined with a graphical network analysis to reveal positive links between 13 of the species. We found that the percentage of patches with positive co-occurrences increased at higher elevations, in contrast to the predictions of the SGH, but in accordance with a refined stress-gradient hypothesis for arid sites, in which characteristics of the interacting species are incorporated.
Questions: Habitat islands are often characterized by the presence of more or less sharp boundaries to adjacent matrix habitats. However, knowledge on boundaries of natural habitat islands is scarce, especially regarding patterns of beta diversity and its two underlying components: species turnover and nestedness. We therefore aim to quantify the effects of fine-scaled and sharp boundaries of quartz islands (quartz gravel-covered soils) on the different components of plant beta diversity and how they are linked to different soil environmental drivers. Location: Knersvlakte, Western Cape, South Africa. Methods: We sampled plant species richness in 56 fine-scale transects of 6 m × 1 m plots across eight different boundary types (four quartz island to matrix, four between habitats on quartz islands). Soil depth and chemistry (pH, electrical conductivity) were analyzed for each 1 m2 plot. Differences in the two beta diversity components (turnover and nestedness) for each boundary type were tested by t tests. We used linear models to test relationships between species and environmental dissimilarity. Results: All boundary types showed high beta diversity. Species turnover was the prevailing component for six boundary types, the nestedness component was only important for two boundary types. We found a significant linear increase of species dissimilarity with increasing dissimilarity in soil pH and distinct plant communities for the habitat types, but no significant increase for electrical conductivity or soil depth. Conclusions: The spatial distinctiveness of the quartz islands leads to sharp boundaries, which result in high beta diversity, mainly through species turnover. This reflects the high levels of diversification and adaptation of the local plant communities. Nestedness occurred at two boundaries to the matrix, indicating that the latter does not necessarily represent an impermeable boundary for all species of the respective ecosystem. Studying diversity patterns across boundaries contributes to the question of applicability of island biogeography theory to habitat islands.
The Late Tertiary to Quaternary evolution of the Ntem interior delta in SW Cameroon shall be modelled. A step fault was formed along neotectonically remobilized Precambrian structures. Uncalibrated 14C-datations in this ‘sediment trap’ show Pleistocene to Holocene ages. Both within and below the interior delta pebbles and clasts which are cemented in an iron and manganese matrix were found. These ‘fanglomerates’ are used to discuss different processes of the younger evolution also concerning climatic fluctuations in the study area.
Inappropriate land management leads to soil loss with destruction of the land’s resource and sediment input into the receiving river. Part of the sediment budget of a catchment is the estimation of soil loss. In the Ruzizi catchment in the Eastern Democratic Republic of the Congo (DRC), only limited research has been conducted on soil loss mainly dealing with local observations on geomorphological forms or river load measurements; a regional quantification of soil loss is missing so far. Such quantifications can be calculated using the Universal Soil Loss Equation (USLE). It is composed of four factors: precipitation (R), soil (K), topography (LS), and vegetation cover (C). The factors can be calculated in different ways according to the characteristics of the study area. In this paper, different approaches for calculating the single factors are reviewed and validated with field work in two sub-catchments of Ruzizi River supplying the water for the reservoirs of Ruzizi I and II hydroelectric dams. It became obvious that the (R)USLE model provides the best results with revised R and LS factors. C factor calculations required to conduct a supervised classification using the Maximum Likelihood Procedure. Different C factor values were assigned to the land cover classes. The calculations resulted in a soil loss rate for the predominantly occurring Ferralsols and Leptosols of around 576 kt/yr in both catchments, when 2016 landcover and precipitation are used. This represents an area-normalized value of 40.4 t/ha/yr for Ruzizi I and 50.5 t/ha/yr for Ruzizi II due to different landcover in the two sub-catchments. The mean value for the whole study area is 47.8 t/ha/yr or even 27.1 t/ha/yr when considering land management techniques like terracing on the slopes (P factor). This work has shown that the (R)USLE model can serve as an easy to handle tool for soil loss quantification when comprehensive field work results are sparse. The model can be implemented in Geographic Information Systems (GIS) with free data; hence, a validation is crucial. It becomes apparent that the use of high resolution Sentinel 2a MSI data as the basis for C factor calculations is an appropriate method for considering heterogeneous Land Use Land Cover (LULC) patterns. To transfer the approach to other regions, the calculation of factor R needs to be modified
We present the first 3-D model of seismic P and S velocities in the crust and uppermost mantle beneath the Gulf of Aqaba and surrounding areas based on the results of passive travel time tomography. The tomographic inversion was performed based on travel time data from ∼ 9000 regional earthquakes provided by the Egyptian National Seismological Network (ENSN), and this was complemented with data from the International Seismological Centre (ISC). The resulting P and S velocity patterns were generally consistent with each other at all depths. Beneath the northern part of the Red Sea, we observed a strong high-velocity anomaly with abrupt limits that coincide with the coastal lines. This finding may indicate the oceanic nature of the crust in the Red Sea, and it does not support the concept of gradual stretching of the continental crust. According to our results, in the middle and lower crust, the seismic anomalies beneath the Gulf of Aqaba seem to delineate a sinistral shift (∼ 100 km) in the opposite flanks of the fault zone, which is consistent with other estimates of the left-lateral displacement in the southern part of the Dead Sea Transform fault. However, no displacement structures were visible in the uppermost lithospheric mantle.
One possible approach to study systematically the influence of the deformation regime on the geometry of geological structures like folds and boudins is analogue modelling. For a complete understanding of the resulting structures, consideration of the third dimension is required. This PhD study deals with scaled analogue modelling under constriction and plane-strain conditions to improve our knowledge of folding and boudinage of lower crustal rocks in space and time. Plasticine is an appropriate analogue material for rocks in the lower crust. Therefore, this material was used for the experiments. The macroscopic behaviour of most types of plasticine is quite similar to rocks undergoing strain-rate softening and strain hardening regardless of the different microscopic aspects of deformation. Therefore, if one is aware that the stress exponent and viscosity increase with increasing strain, the original plasticine types used with stress exponents ranging from 5.8 to 8.0 are adequate for modelling geologic structures. The same holds for plasticine/oil mixtures. Thus, plasticine and plasticine/oil mixtures can be used to model the viscous flow of different rock types in the lower crust. If climb-accommodated dislocation creep and associated steady-state flow is assumed for the natural rocks, the plasticine/oil mixtures should be used, which flow under steady-state conditions. Three different experimental studies of plane-strain coaxial deformation of stiff layers, with viscosity η2 and stress exponent n2, embedded in a weak matrix, with viscosity η1 and stress exponent n1, have been carried out. The undeformed samples (matrix plus layer) were cubes with an edge length of 12 cm. All experimental runs have been carried out at T = 25 ± 1°C and varying strain rates ė, ranging from 7.9 x 10 high -6 s high -1 to 1.7 x 10 high -2 s high -1, until a finite longitudinal strain of 30% – 40% was achieved. The first experimental study improved the understanding about the evolution of folds and boudins when the layer is oriented perpendicular to the Y-axis of the finite strain ellipsoid. The rock analogues used were Beck’s green plasticine (matrix) and Beck’s black plasticine (competent layer), both of which are strain-rate softening modelling materials with stress exponent n = ca. 8. The effective viscosity η of the matrix plasticine was changed by adding different amounts of oil to the original plasticine. At a strain rate ė of 10 high -3 s high -1 and a finite strain e of 10%, the effective viscosity of the matrix ranges from 1.2 x 10 high 6 to 7.2 x 10 high 6 Pa s. The effective viscosity of the competent layer has been determined as 4.2 x 10 high 7 Pa s. If the viscosity ratio is large (> ca. 20) and the initial thickness of the competent layer is small, both folds and boudins develop simultaneously. Although the growth rate of the folds seems to be higher than the growth rate of the boudins, the wavelength of both structures is approximately the same as is suggested by analytical solutions. A further unexpected, but characteristic, aspect of the deformed competent layer is a significant increase in thickness, which can be used to distinguish plane-strain folds and boudins from constrictional folds and boudins. In the second experimental study, the impact of varying strain rates on growing folds and boudins under plane strain have been investigated. The strain rates used range from 7.9 x 10 high -6 s high -1 to 1.7 x 10 high -2 s high -1. The stiff layer and matrix consist of non-linear viscous Kolb grey and Beck’s green plasticine, respectively, both of which are strain-rate softening modelling materials with power law exponents (n) and apparent viscosities (η) ranging from 6.5 to 7.9 and 8.5 x 10 high 6 to 7.2 x 10 high 6 Pa s, respectively. The effective viscosity (η) of the matrix plasticine was partly modified by adding oil to the original plasticine. At the strain rates used in the experiments the viscosity ratio between layer and matrix ranges between 3 and 10. Different runs have been carried out where the layer was oriented perpendicular to the principal strain axes (X>Y>Z). The results suggest a considerable influence of the strain rate on the geometry of the deformed stiff layer including its thickness. This holds for every type of layer orientation (S ┴ X, S ┴ Y, S ┴ Z). If the stiff layer is oriented perpendicular to the short axis Z of the finite strain ellipsoid, the number of the resulting boudins and the thickness of the stiff layer increase, whereas the length of boudins decreases with increasing strain rate. If the stiff layer is oriented perpendicular to the long axis, X, of the finite strain ellipsoid, enlargement of the strain rate results in increasing wavelength of folds, whereas the number of folds and the degree of thickening of the stiff layer decreased. If the stiff layer is oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid enlargement of the strain rate results in a decreasing number of boudins and folds associated with increasing wavelengths of both structures. The wavelength of folds is approximately half of the boudins wavelength. This is true for the case where folds and boudins develop simultaneously (S ┴ Y) and for cases where both structures develop independently (folds at S ┴ X and boudins at S ┴ Z). In the third experimental study, scaled analogue experiments have been carried out to demonstrate the growth of plane-strain folds and boudins through space and time. Previous 3D-studies are based only on finite deformation structures. Their results can therefore not be used to prove if both structures grew simultaneously or in sequence. Plane strain acted on a single stiff layer that was embedded in a weak matrix, with the layer oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid. Two different experimental runs have been carried out using computer tomography (CT) to analyse the results. The first run was carried out without interruption. During the second run, the deformation was stopped in each case at longitudinal strain increments of 10%. Every experiment was carried out at a temperature T of 25°C and a strain rate, ė, of ca. 4 x 10 high -3 s high -1 until a finite longitudinal strain of 40% was achieved with a viscosity contrast m of 18.6 between the non-linear viscous layer (Kolb brown plasticine) and the matrix (Beck’s green plasticine with 150 ml oil kg high -1). The apparent viscosity, η, and the stress exponent, n, for the layer at a strain rate ė = ca. 10 high -3 s high -1 and a finite strain e = 10% are 2.23 x 10 high 7 Pa s and n = 5.8 and for the matrix 1.2 x 10 high 6 Pa s and 10.5. These new data that result from incremental analogue modelling corroborate previous suggestions that folds and boudins are coeval structures in cases of plane-strain coaxial deformation with the stiff layer oriented perpendicular to the intermediate Y-axis of the finite strain ellipsoid. They will be of interest for all workers who are dealing with plane-strain boudins and folds, where the fold axes are parallel to the major axis (X) of the finite strain ellipsoid. As has been demonstrated by the first experimental study, coeval folding and boudinage under plane strain, with S ┴ Y, are associated with a significant increase in the thickness of the competent layer. The latter phenomenon does not occur in other cases of simultaneous folding and boudinage, such as bulk pure constriction. To study the impact of layer thickness on the geometry of folds and boudins under pure constriction, we carried out additional experiments using different types of plasticine for a stiff layer and a weaker matrix to model folding and boudinaging under pure constriction, with the initially planar layer oriented parallel to the Xaxis of the finite strain ellipsoid. The stiff layer and matrix consist of non-linear viscous Kolb brown and Beck’s green plasticine, respectively, both of which are strain-rate softening modelling materials. Six runs have been carried out using thicknesses of the stiff layer of 1, 2, 4, 6, 8 and 10 ± 0.2 mm. All experimental runs were carried out at a temperature T of 30 ± 2°C and a strain rate, ė, of ca. 1.1 x 10 high -4 s high -1 until a finite longitudinal strain of 40% was achieved with a viscosity contrast m of 3.1 between the stiff layer (Kolb brown plasticine) and the matrix (Beck’s green plasticine). The apparent viscosity, η, and the stress exponent, n, for the layer at a strain rate ė = ca. 10 high -3 s high -1 and a finite strain e = 10% are 2.23 x 10 high 7 Pa s and n = 5.8 and for the matrix 7.2 x 10 high 6 Pa s and 7.9. Our results suggest a considerable influence of the initial thickness of the stiff layer on the geometry of the deformed stiff layer. There is no evidence for folding in XY=XZ-sections if the initial thickness of the competent layer is larger than ca. 8 mm. If the initial thickness of the competent layer is set at ca. 10 ± 0.2 mm, both folds and boudins develop simultaneously. However, the growth rate of the boudins seems to be higher than the growth rate of the folds. A further expected, but characteristic, aspect of the deformed competent layer is no change in thickness of the competent layer, which can be used to distinguish plane-strain folds and boudins from constrictional folds and boudins. The model results are important for the analysis and interpretation of deformation structures in rheologically stratified rocks undergoing dislocation creep under bulk constriction. Tectonic settings where constrictional folds and boudins may develop simultaneously are stems of salt diapirs, subduction zones or thermal plumes. To make (paleo) viscosimetric statements possible, the rheological data of the different plasticine types were related to the geometrical data. When comparing the normalized dominant wavelength Wd obtained from the deformed layer of the models with the theoretical dominant wavelength (Ld) calculated using the Smith equation (1977, 1979), the latter probably also holds when folding and boudinage develop simultaneously (S ┴ Y) and when boudins develop independently (S ┴ Z), but can obviously not be applied at very low viscosity ratios as is indicated by the low-strain-rate experiments.
During SPURT (Spurenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region) we performed measurements of a wide range of trace gases with different lifetimes and sink/source characteristics in the northern hemispheric upper troposphere (UT) and lowermost stratosphere (LMS). A large number of in-situ instruments were deployed on board a Learjet 35A, flying at altitudes up to 13.7 km, at times reaching to nearly 380 K potential temperature. Eight measurement campaigns (consisting of a total of 36 flights), distributed over all seasons and typically covering latitudes between 35° N and 75° N in the European longitude sector (10° W–20° E), were performed. Here we present an overview of the project, describing the instrumentation, the encountered meteorological situations during the campaigns and the data set available from SPURT. Measurements were obtained for N2O, CH4, CO, CO2, CFC12, H2, SF6, NO, NOy, O3 and H2O. We illustrate the strength of this new data set by showing mean distributions of the mixing ratios of selected trace gases, using a potential temperature – equivalent latitude coordinate system. The observations reveal that the LMS is most stratospheric in character during spring, with the highest mixing ratios of O3 and NOy and the lowest mixing ratios of N2O and SF6. The lowest mixing ratios of NOy and O3 are observed during autumn, together with the highest mixing ratios of N2O and SF6 indicating a strong tropospheric influence. For H2O, however, the maximum concentrations in the LMS are found during summer, suggesting unique (temperature- and convection-controlled) conditions for this molecule during transport across the tropopause. The SPURT data set is presently the most accurate and complete data set for many trace species in the LMS, and its main value is the simultaneous measurement of a suite of trace gases having different lifetimes and physical-chemical histories. It is thus very well suited for studies of atmospheric transport, for model validation, and for investigations of seasonal changes in the UT/LMS, as demonstrated in accompanying and elsewhere published studies.
During SPURT (Spurenstofftransport in der Tropopausenregion, trace gas transport in the tropopause region) we performed measurements of a wide range of trace gases with different lifetimes and sink/source characteristics in the northern hemispheric upper troposphere (UT) and lowermost stratosphere (LMS). A large number of in-situ instruments were deployed on board a Learjet 35A, flying at altitudes up to 13.7 km, at times reaching to nearly 380 K potential temperature. Eight measurement campaigns (consisting of a total of 36 flights), distributed over all seasons and typically covering latitudes between 35° N and 75° N in the European longitude sector (10° W–20° E), were performed. Here we present an overview of the project, describing the instrumentation, the encountered meteorological situations during the campaigns and the data set available from SPURT. Measurements were obtained for N2O, CH4, CO, CO2, CFC12, H2, SF6, NO, NOy, O3 and H2O. We illustrate the strength of this new data set by showing mean distributions of the mixing ratios of selected trace gases, using a potential temperature – equivalent latitude coordinate system. The observations reveal that the LMS is most stratospheric in character during spring, with the highest mixing ratios of O3 and NOy and the lowest mixing ratios of N2O and SF6. The lowest mixing ratios of NOy and O3 are observed during autumn, together with the highest mixing ratios of N2O and SF6 indicating a strong tropospheric influence. For H2O, however, the maximum concentrations in the LMS are found during summer, suggesting unique (temperature- and convection-controlled) conditions for this molecule during transport across the tropopause. The SPURT data set is presently the most accurate and complete data set for many trace species in the LMS, and its main value is the simultaneous measurement of a suite of trace gases having different lifetimes and physical-chemical histories. It is thus very well suited for studies of atmospheric transport, for model validation, and for investigations of seasonal changes in the UT/LMS, as demonstrated in accompanying and elsewhere published studies.
Chlorine and bromine atoms can lead to catalytic destruction of ozone in the stratosphere. Therefore the use and production of ozone depleting substances (ODS) containing chlorine and bromine is regulated by the Montreal Protocol to protect the ozone layer. Equivalent Effective Stratospheric Chlorine (EESC) has been adapted as an appropriate metric to describe the combined effects of chlorine and bromine released from halocarbons on stratospheric ozone. Here we revisit the concept of calculating EESC. We derive a new formulation of EESC based on an advanced concept of ODS propagation into the stratosphere and reactive halogen release. A new transit time distribution is introduced in which the age spectrum for an inert tracer is weighted with the release function for inorganic halogen from the source gases. This distribution is termed the “release time distribution”. The improved formulation shows that EESC levels in the year 1980 for the mid latitude lower stratosphere were significantly lower than previously calculated. 1980 marks the year commonly defined as the onset of anthropogenic ozone depletion in the stratosphere. Assuming that the EESC value must return to the same level in order for ozone to fully recover, we show that it will take more than 10 years longer than currently assumed in this region of the stratosphere. Based on the improved formulation, EESC level at mid-latitudes will reach this landmark only in 2060. We also present a range of sensitivity studies to investigate the effect of changes and uncertainties in the fractional release factors and in the assumptions on the shape of the release time distributions. We conclude that, under the assumptions that all other atmospheric parameters like stratospheric dynamics and chemistry are unchanged, the recovery of mid latitude stratospheric ozone would be expected to be delayed by about a 10 years, in a similar way as EESC.
Chlorine and bromine atoms lead to catalytic depletion of ozone in the stratosphere. Therefore the use and production of ozone-depleting substances (ODSs) containing chlorine and bromine is regulated by the Montreal Protocol to protect the ozone layer. Equivalent effective stratospheric chlorine (EESC) has been adopted as an appropriate metric to describe the combined effects of chlorine and bromine released from halocarbons on stratospheric ozone. Here we revisit the concept of calculating EESC. We derive a refined formulation of EESC based on an advanced concept of ODS propagation into the stratosphere and reactive halogen release. A new transit time distribution is introduced in which the age spectrum for an inert tracer is weighted with the release function for inorganic halogen from the source gases. This distribution is termed the release time distribution. We show that a much better agreement with inorganic halogen loading from the chemistry transport model TOMCAT is achieved compared with using the current formulation. The refined formulation shows EESC levels in the year 1980 for the mid-latitude lower stratosphere, which are significantly lower than previously calculated. The year 1980 is commonly used as a benchmark to which EESC must return in order to reach significant progress towards halogen and ozone recovery. Assuming that – under otherwise unchanged conditions – the EESC value must return to the same level in order for ozone to fully recover, we show that it will take more than 10 years longer than estimated in this region of the stratosphere with the current method for calculation of EESC. We also present a range of sensitivity studies to investigate the effect of changes and uncertainties in the fractional release factors and in the assumptions on the shape of the release time distributions. We further discuss the value of EESC as a proxy for future evolution of inorganic halogen loading under changing atmospheric dynamics using simulations from the EMAC model. We show that while the expected changes in stratospheric transport lead to significant differences between EESC and modelled inorganic halogen loading at constant mean age, EESC is a reasonable proxy for modelled inorganic halogen on a constant pressure level.
MIPAS-Envisat is a satellite-borne sensor which measured vertical profiles of a wide range of trace gases from 2002 to 2012 using IR emission spectroscopy. We present geophysical validation of the MIPAS-Envisat operational retrieval (version 6.0) of N2O, CH4, CFC-12, and CFC-11 by the European Space Agency (ESA). The geophysical validation data are derived from measurements of samples collected by a cryogenic whole air sampler flown to altitudes of up to 34 km by means of large scientific balloons. In order to increase the number of coincidences between the satellite and the balloon observations, we applied a trajectory matching technique. The results are presented for different time periods due to a change in the spectroscopic resolution of MIPAS in early 2005. Retrieval results for N2O, CH4, and CFC-12 show partly good agreement for some altitude regions, which differs for the periods with different spectroscopic resolution. The more recent low spectroscopic resolution data above 20 km altitude show agreement with the combined uncertainties, while there is a tendency of the earlier high spectral resolution data set to underestimate these species above 25 km. The earlier high spectral resolution data show a significant overestimation of the mixing ratios for N2O, CH4, and CFC-12 below 20 km. These differences need to be considered when using these data. The CFC-11 results from the operation retrieval version 6.0 cannot be recommended for scientific studies due to a systematic overestimation of the CFC-11 mixing ratios at all altitudes.
MIPAS-Envisat is a satellite-borne sensor which measured vertical profiles of a wide range of trace gases from 2002 to 2012 using IR emission spectroscopy. We present geophysical validation of the MIPAS-Envisat operational retrieval (version 6.0) of N2O, CH4, CFC-12, and CFC-11 by the European Space Agency (ESA). The geophysical validation data are derived from measurements of samples collected by a cryogenic whole air sampler flown to altitudes of up to 34 km by means of large scientific balloons. In order to increase the number of coincidences between the satellite and the balloon observations, we applied a trajectory matching technique. The results are presented for different time periods due to a change in the spectroscopic resolution of MIPAS in early 2005. Retrieval results for N2O, CH4, and CFC-12 show partly good agreement for some altitude regions, which differs for the periods with different spectroscopic resolution. The more recent low spectroscopic resolution data above 20 km altitude show agreement with the combined uncertainties, while there is a tendency of the earlier high spectral resolution data set to underestimate these species above 25 km. The earlier high spectral resolution data show a significant overestimation of the mixing ratios for N2O, CH4, and CFC-12 below 20 km. These differences need to be considered when using these data. The CFC-11 results from the operation retrieval version 6.0 cannot be recommended for scientific studies due to a systematic overestimation of the CFC-11 mixing ratios at all altitudes.
Mean age of stratospheric air can be derived from observations of sufficiently long-lived trace gases with approximately linear trends in the troposphere. Mean age can serve as a tracer to investigate stratospheric transport and long-term changes in the strength of the overturning Brewer–Dobson circulation of the stratosphere. For this purpose, a low-cost method is required in order to allow for regular observations up to altitudes of about 30 km. Despite the desired low costs, high precision and accuracy are required in order to determine mean age. We present balloon-borne AirCore observations from two midlatitude sites: Timmins in Ontario/Canada and Lindenberg in Germany. During the Timmins campaign, five AirCores sampled air in parallel with a large stratospheric balloon and were analysed for CO2, CH4 and partly CO. We show that there is good agreement between the different AirCores (better than 0.1 %), especially when vertical gradients are small. The measurements from Lindenberg were performed using small low-cost balloons and yielded very comparable results. We have used the observations to extend our long-term data set of mean age observations at Northern Hemisphere midlatitudes. The time series now covers more than 40 years and shows a small, statistically non-significant positive trend of 0.15 ± 0.18 years decade−1. This trend is slightly smaller than the previous estimate of 0.24 ± 0.22 years decade−1 which was based on observations up to the year 2006. These observations are still in contrast to strong negative trends of mean age as derived from some model calculations.
Mean age of stratospheric air can be derived from observations of sufficiently long lived trace gases with approximately linear trends in the troposphere. Mean age can serve as a tracer to investigate stratospheric transport and long term changes in the strength of the overturning Brewer-Dobson circulation of the stratosphere. For this purpose, a low-cost method is required in order to allow for regular observations up to altitudes of about 30 km. Despite the desired low costs, high precision and accuracy are required in order to allow determination of mean age. We present balloon borne AirCore observations from two mid latitude sites: Timmins in Ontario/Canada and Lindenberg in Germany. During the Timmins campaign five AirCores sampled air in parallel from a large stratospheric balloon and were analysed for CO2, CH4 and partly CO. We show that there is good agreement between the different AirCores (better than 0.1 %) especially when vertical gradients are small. The measurements from Lindenberg were performed using small low-cost balloons and yielded very comparable results. We have used the observations to extend our long term data set of mean age observations at Northern Hemi-sphere mid latitudes. The time series now covers more than 40 years and shows a small, statis-tically not significant positive trend of 0.15 ± 0.18 years/decade. This trend is slightly smaller than the previous estimate of 0.24 ± 0.22 years/decade which was based on observations up to the year 2006. These observations are still in contrast to strong negative trends of mean age as derived from some model calculations.
On the observation of mesospheric air inside the arctic stratospheric polar vortex in early 2003
(2005)
During several balloon flights inside the Arctic polar vortex in early 2003, unusual trace gas distributions were observed, which indicate a strong influence of mesospheric air in the stratosphere. The tuneable diode laser (TDL) instrument SPIRALE (Spectroscopie InFrarouge par Absorption de Lasers Embarqués) measured unusually high CO values (up to 600 ppb) on 27 January at about 30 km altitude. The cryosampler BONBON sampled air masses with very high molecular Hydrogen, extremely low SF6 and enhanced CO values on 6 March at about 25 km altitude. Finally, the MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) Fourier Transform Infra-Red (FTIR) spectrometer showed NOy values which are significantly higher than NOy* (the NOy derived from a correlation between N2O and NOy under undisturbed conditions), on 21 and 22 March in a layer centred at 22 km altitude. Thus, the mesospheric air seems to have been present in a layer descending from about 30 km in late January to 25 km altitude in early March and about 22 km altitude on 20 March. We present corroborating evidence from a model study using the KASIMA (KArlsruhe Simulation model of the Middle Atmosphere) model that also shows a layer of mesospheric air, which descended into the stratosphere in November and early December 2002, before the minor warming which occurred in late December 2002 lead to a descent of upper stratospheric air, cutting of a layer in which mesospheric air is present. This layer then descended inside the vortex over the course of the winter. The same feature is found in trajectory calculations, based on a large number of trajectories started in the vicinity of the observations on 6 March. Based on the difference between the mean age derived from SF6 (which has an irreversible mesospheric loss) and from CO2 (whose mesospheric loss is much smaller and reversible) we estimate that the fraction of mesospheric air in the layer observed on 6 March, must have been somewhere between 35% and 100%.
During several balloon flights inside the Arctic polar vortex in early 2003, unusual trace gas distributions were observed, which indicate a strong influence of mesospheric air in the stratosphere. The tuneable diode laser (TDL) instrument SPIRALE (Spectroscopie InFrarouge par Absorption de Lasers Embarqués) measured unusually high CO values (up to 600 ppb) on 27 January at about 30 km altitude. The cryosampler BONBON sampled air masses with very high molecular Hydrogen, extremely low SF6 and enhanced CO values on 6 March at about 25 km altitude. Finally, the MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) Fourier Transform Infra-Red (FTIR) spectrometer showed NOy values which are significantly higher than NOy* (the NOy derived from a correlation between N2O and NOy under undisturbed conditions), on 21 and 22 March in a layer centred at 22 km altitude. Thus, the mesospheric air seems to have been present in a layer descending from about 30 km in late January to 25 km altitude in early March and about 22 km altitude on 20 March. We present corroborating evidence from a model study using the KASIMA (KArlsruhe Simulation model of the Middle Atmosphere) model that also shows a layer of mesospheric air, which descended into the stratosphere in November and early December 2002, before the minor warming which occurred in late December 2002 lead to a descent of upper stratospheric air, cutting of a layer in which mesospheric air is present. This layer then descended inside the vortex over the course of the winter. The same feature is found in trajectory calculations, based on a large number of trajectories started in the vicinity of the observations on 6 March. Based on the difference between the mean age derived from SF6 (which has an irreversible mesospheric loss) and from CO2 (whose mesospheric loss is much smaller and reversible) we estimate that the fraction of mesospheric air in the layer observed on 6 March, must have been somewhere between 35% and 100%.
Abstract:
The mid Miocene represents an important target for paleoclimatic study because the atmospheric CO2 concentration ranged from near modern values to ∼800 ppm, while a large, dynamic Antarctic ice sheet was likely to have been present throughout much of this interval. In this special issue, Modestou et al. (2020) (doi.org/10.1029/2020PA003927) reconstruct deep ocean warmth based on the clumped isotopic composition of benthic foraminifera, a technique that allows the ice volume and thermal components of the benthic oxygen isotope stack to be separated. These data reveal a very warm deep ocean while simultaneously suggesting that continental ice volume may, at times, have been greater than today. Here, I review these results in the context of recent developments in geochemical proxies and ice sheet modeling, and explore how the presence of a large Miocene ice sheet could be reconciled with CO2 at least as high as present. More broadly, I argue that many of the 'paradoxes' that pepper the paleoclimate literature result as much from our imperfect understanding of the proxies, as from our understanding of the climate system. Robust proxies with a well-understood mechanistic basis, as employed by Modestou et al. (2020), as well as advances in model-data comparability usher in a new era of palaeoclimate research; an exciting future of untangling Earth's myriad past climate states awaits.
Plain Language Summary:
Reconstructing climate variation in Earth's geologic past informs us of the broad features of warm climates, which is relevant to preparing for climate change over the coming centuries. Moreover, these data can be compared to state-of-the-art climate models, which provides a test of the degree to which our models can reproduce warm climate states. A paper recently published in this journal applies a new method in order to reconstruct the temperature of the deep ocean in the middle Miocene (between 17 and 12 million years ago), when the atmospheric CO2 concentration was naturally similar to or higher than it is today. Coupled with decades of previous study, these exciting results depict an unfamiliar world characterized by a warm deep ocean, and yet a large ice sheet was present on Antarctica. Both models and data agree that the Antarctic ice sheet in the Miocene was highly responsive to changes in the atmospheric CO2 concentration, a clear cause of concern in the context of ongoing anthropogenic climate change.
The turnover time of terrestrial ecosystem carbon is an emergent ecosystem property that quantifies the strength of land surface on the global carbon cycle–climate feedback. However, observation- and modeling-based estimates of carbon turnover and its response to climate are still characterized by large uncertainties. In this study, by assessing the apparent whole ecosystem carbon turnover times (τ) as the ratio between carbon stocks and fluxes, we provide an update of this ecosystem level diagnostic and its associated uncertainties in high spatial resolution (0.083∘) using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble of data, we estimated the global median τ to be 43+7−7 yr (median+difference to percentile 75−difference to percentile 25) when the full soil is considered, in contrast to limiting it to 1 m depth. Only considering the top 1 m of soil carbon in circumpolar regions (assuming maximum active layer depth is up to 1 m) yields a global median τ of 37+3−6 yr, which is longer than the previous estimates of 23+7−4 yr (Carvalhais et al., 2014). We show that the difference is mostly attributed to changes in global Csoil estimates. Csoil accounts for approximately 84 % of the total uncertainty in global τ estimates; GPP also contributes significantly (15 %), whereas Cveg contributes only marginally (less than 1 %) to the total uncertainty. The high uncertainty in Csoil is reflected in the large range across state-of-the-art data products, in which full-depth Csoil spans between 3362 and 4792 PgC. The uncertainty is especially high in circumpolar regions with an uncertainty of 50 % and a low spatial correlation between the different datasets (0.2<r<0.5) when compared to other regions (0.6<r<0.8). These uncertainties cast a shadow on current global estimates of τ in circumpolar regions, for which further geographical representativeness and clarification on variations in Csoil with soil depth are needed. Different GPP estimates contribute significantly to the uncertainties of τ mainly in semiarid and arid regions, whereas Cveg causes the uncertainties of τ in the subtropics and tropics. In spite of the large uncertainties, our findings reveal that the latitudinal gradients of τ are consistent across different datasets and soil depths. The current results show a strong ensemble agreement on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30–60∘ N) than in the subtropical and tropical zones (30∘ S–30∘ N). Additionally, while the strength of the τ–precipitation correlation was dependent on the Csoil data source, the latitudinal gradients also agree among different ensemble members. Overall, and despite the large variation in τ, we identified robust features in the spatial patterns of τ that emerge beyond the differences stemming from the data-driven estimates of Csoil, Cveg and GPP. These robust patterns, and associated uncertainties, can be used to infer τ–climate relationships and for constraining contemporaneous behavior of Earth system models (ESMs), which could contribute to uncertainty reductions in future projections of the carbon cycle–climate feedback. The dataset of τ is openly available at https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
The turnover time of terrestrial carbon (τ) controls the global carbon cycle – climate feedback and, yet, is poorly simulated by the current Earth System Models (ESMs). In this study, by assessing apparent carbon turnover time as the ratio between carbon stocks and fluxes, we provide a new, updated ensemble of diagnostic terrestrial carbon turnover times and associated uncertainties on a global scale using multiple, state-of-the-art, observation-based datasets of soil organic carbon stock (Csoil), vegetation biomass (Cveg) and gross primary productivity (GPP). Using this new ensemble, we estimated the global average τ to be 42$% &' years when the full soil depth is considered, longer than the previous estimates of 23$) &* years. Only considering the top 1 m (assuming maximum active layer depth is up to 1 meter) of soil carbon in circumpolar regions yields a global τ of 35$) &' years. Csoil in circumpolar regions account for two thirds of the total uncertainty in global τ estimates, whereas Csoil in non-circumpolar contributes merely 9.38%. GPP (2.25%) and Cveg (0.05%) contribute even less to the total uncertainty. Therefore, the high uncertainty in Csoil is the main factor behind the uncertainty in global τ, as reflected in the larger range of full-depth Csoil (3152-4372 PgC). The uncertainty is especially high in circumpolar regions with a behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI: 10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019) uncertainty of 50% and the spatial correlations among different datasets are also low compared to other regions. Overall, we argue that current global datasets do not support robust estimates of τ globally, for which we need clarification on variations of Csoil with soil depth and stronger estimates of Csoil in circumpolar regions. Despite the large variation in both magnitude and spatial patterns of τ, we identified robust features in the spatial patterns of τ that emerge regardless of soil depth and differences in data sources of Csoil, Cveg and GPP. Our findings show that the latitudinal gradients of τ are consistent across different datasets and soil depth. Furthermore, there is a strong consensus on the negative correlation between τ and temperature along latitude that is stronger in temperate zones (30ºN-60ºN) than in subtropical and tropical zones (30ºS30ºN). The identified robust patterns can be used to infer the response of τ to climate and for constraining contemporaneous behaviour of ESMs which could contribute to uncertainty reductions in future projections of the carbon cycle - climate feedback. The dataset of the terrestrial turnover time ensemble (DOI:10.17871/bgitau.201911) is openly available from the data portal: https://doi.org/10.17871/bgitau.201911 (Fan et al., 2019).
The crystal structure of the high temperature phase of anilinium bromide, C6H5NH3⊕Br⊖ , was studied by X-ray and neutron diffraction at T = 343 K. The refinement supports disordered positions of the -NH3⊕ group. A split-atom model is proposed which includes disorder of the benzene ring. The thermal parameters, hydrogen bond distances, and other experimental data (NMR, NQR, inelastic neutron scattering) are in accordance with this model.
Artificial drainage of agricultural land, for example with ditches or drainage tubes, is used to avoid water logging and to manage high groundwater tables. Among other impacts it influences the nutrient balances by increasing leaching losses and by decreasing denitrification. To simulate terrestrial transport of nitrogen on the global scale, a digital global map of artificially drained agricultural areas was developed. The map depicts the percentage of each 5’ by 5’ grid cell that is equipped for artificial drainage. Information on artificial drainage in countries or sub-national units was mainly derived from international inventories. Distribution to grid cells was based, for most countries, on the "Global Croplands Dataset" of Ramankutty et al. (1998) and the "Digital Global Map of Irrigation Areas" of Siebert et al. (2005). For some European countries the CORINE land cover dataset was used instead of the both datasets mentioned above. Maps with outlines of artificially drained areas were available for 6 countries. The global drainage area on the map is 167 Mio hectares. For only 11 out of the 116 countries with information on artificial drainage areas, sub-national information could be taken into account. Due to this coarse spatial resolution of the data sources, we recommended to use the map of artificially drained areas only for continental to global scale assessments. This documentation describes the dataset, the data sources and the map generation, and it discusses the data uncertainty.
We have used the SLIMCAT 3-D off-line chemical transport model (CTM) to quantify the Arctic chemical ozone loss in the year 2002/2003 and compare it with similar calculations for the winters 1999/2000 and 2003/2004. Recent changes to the CTM have improved the model's ability to reproduce polar chemical and dynamical processes. The updated CTM uses σ-θ as a vertical coordinate which allows it to extend down to the surface. The CTM has a detailed stratospheric chemistry scheme and now includes a simple NAT-based denitrification scheme in the stratosphere.
In the model runs presented here the model was forced by ECMWF ERA40 and operational analyses. The model used 24 levels extending from the surface to ~55 km and a horizontal resolution of either 7.5°×7.5° or 2.8°×2.8°. Two different radiation schemes, MIDRAD and the CCM scheme, were used to diagnose the vertical motion in the stratosphere. Based on tracer observations from balloons and aircraft, the more sophisticated CCM scheme gives a better representation of the vertical transport in this model which includes the troposphere. The higher resolution model generally produces larger chemical O3 depletion, which agrees better with observations.
The CTM results show that very early chemical ozone loss occurred in December 2002 due to extremely low temperatures and early chlorine activation in the lower stratosphere. Thus, chemical loss in this winter started earlier than in the other two winters studied here. In 2002/2003 the local polar ozone loss in the lower stratosphere was ~40% before the stratospheric final warming. Larger ozone loss occurred in the cold year 1999/2000 which had a persistently cold and stable vortex during most of the winter. For this winter the current model, at a resolution of 2.8°×2.8°, can reproduce the observed loss of over 70% locally. In the warm and more disturbed winter 2003/2004 the chemical O3 loss was generally much smaller, except above 620 K where large losses occurred due to a period of very low minimum temperatures at these altitudes.
We have used the SLIMCAT 3-D off-line chemical transport model (CTM) to quantify the Arctic chemical ozone loss in the year 2002/2003 and compare it with similar calculations for the winters 1999/2000 and 2003/2004. Recent changes to the CTM have improved the model's ability to reproduce polar chemical and dynamical processes. The updated CTM uses σ-θ as a vertical coordinate which allows it to extend down to the surface. The CTM has a detailed stratospheric chemistry scheme and now includes a simple NAT-based denitrification scheme in the stratosphere.
In the model runs presented here the model was forced by ECMWF ERA40 and operational analyses. The model used 24 levels extending from the surface to ~55km and a horizontal resolution of either 7.5° x 7.5° or 2.8° x 2.8°. Two different radiation schemes, MIDRAD and the CCM scheme, were used to diagnose the vertical motion in the stratosphere. Based on tracer observations from balloons and aircraft, the more sophisticated CCM scheme gives a better representation of the vertical transport in this model which includes the troposphere. The higher resolution model generally produces larger chemical O3 depletion, which agrees better with observations.
The CTM results show that very early chemical ozone loss occurred in December 2002 due to extremely low temperatures and early chlorine activation in the lower stratosphere. Thus, chemical loss in this winter started earlier than in the other two winters studied here. In 2002/2003 the local polar ozone loss in the lower stratosphere was ~40% before the stratospheric final warming. Larger ozone loss occurred in the cold year 1999/2000 which had a persistently cold and stable vortex during most of the winter. For this winter the current model, at a resolution of 2.8° x 2.8°, can reproduce the observed loss of over 70% locally. In the warm and more disturbed winter 2003/2004 the chemical O3 loss was generally much smaller, except above 620K where large losses occurred due to a period of very low minimum temperatures at these altitudes.
The analysis of charcoal fragments in peat and lake sediments is the most widely used approach to reconstruct past biomass burning. With a few exceptions, this method typically relies on the quantification of the total charcoal content of the sediment. To enhance charcoal analyses for the reconstruction of past fire regimes, and to make the method more relevant to studies of both plant evolution and fire management, more information must be extracted from charcoal particles. Here, I burned in the laboratory seven fuel types comprising 17 species from boreal Siberia, and build on published schemes to develop morphometric and finer diagnostic classifications of the experimentally charred particles. As most of the species used in this study are common to Northern Hemisphere forests and peatlands, these results can be directly applicable over a broad geographical scale. Results show that the effect of temperature on charcoal production is fuel dependent. Graminoids and Sphagnum, and wood (trunk) lose the most mass at low burn temperatures, whereas heathland shrub leaves, brown moss, and ferns retain the most mass at high burn temperatures. In contrast to the wood of trunk, the wood of twigs retained their mass at intermediate temperature. This suggests that species with low mass retention at hotter burning temperatures might be underrepresented in the fossil charcoal record. Charred particle aspect ratio (L/W) appeared to be the strongest indicator of the fuel type burnt. Graminoid charcoals are more elongate than those of all other fuel types, leaf charcoals are the shortest and bulkiest, and twig and wood charcoals are intermediate. Finer diagnostic features were the most useful in distinguishing between wood, graminoid, and leaf particles, but further distinctions within these fuel types are difficult. High-aspect-ratio particles dominated by graminoid and Sphagnum morphologies are robust indicators of cooler surface fires. Contrastingly, abundant wood and leaf morphologies and low-aspect-ratio particles likely indicate higher-temperature fires. However, the overlapping morphologies of leaves and wood from trees and shrubs make it hard to distinguish between high-intensity surface fires combusting living shrubs and dead wood and leaves or high-intensity crown fires combusting living trees. Despite these limitations, the combined use of charred-particle aspect ratios and fuel morphotypes can aid in more robustly interpreting changes in fuel source and fire type, thereby substantially refining histories of past wildfires. Further fields of investigation to improve the interpretation of the fossil charcoal records will require: i) More in-depth knowledge of plant anatomy for a better determination of fuel sources; ii) Relate the proportion of particular charcoal morphotypes to the quantity of biomass; iii) Link the chemical composition of fuels, combustion temperature, and charcoal production. The advanced use of image-recognition software to collect data on other charcoal features could also aid in extracting fire temperatures as well as a change in particles morphology and morphometry during particles transportation.
The analysis of charcoal fragments in peat and lake sediments is the most widely used approach to reconstruct past biomass burning. With a few exceptions, this method typically relies on the quantification of the total charcoal content of the sediment. To enhance charcoal analyses for the reconstruction of past fire regimes and make the method more relevant to studies of both plant evolution and fire management, the extraction of more information from charcoal particles is critical. Here, I used a muffle oven to burn seven fuel types comprising 17 species from boreal Siberia (near Teguldet village), which are also commonly found in the Northern Hemisphere, and built on published schemes to develop morphometric and finer diagnostic classifications of the experimentally charred particles. I then combined these results with those from fossil charcoal from a peat core taken from the same location (Ulukh-Chayakh mire) in order to demonstrate the relevance of these experiments to the fossil charcoal records. Results show that graminoids, Sphagnum, and wood (trunk) lose the most mass at low burn temperatures (<300 ∘C), whereas heathland shrub leaves, brown moss, and ferns lose the most mass at high burn temperatures. This suggests that species with low mass retention in high-temperature fires are likely to be under-represented in the fossil charcoal record. The charcoal particle aspect ratio appeared to be the strongest indicator of the fuel type burnt. Graminoid charcoal particles are the most elongate (6.7–11.5), with a threshold above 6 that may be indicative of wetland graminoids; leaves are the shortest and bulkiest (2.1–3.5); and twigs and wood are intermediate (2.0–5.2). Further, the use of fine diagnostic features was more successful in separating wood, graminoids, and leaves, but it was difficult to further differentiate these fuel types due to overlapping features. High-aspect-ratio particles, dominated by graminoid and Sphagnum morphologies, may be robust indicators of low-temperature surface fires, whereas abundant wood and leaf morphologies as well as low-aspect-ratio particles are indicative of higher-temperature fires. However, the overlapping morphologies of leaves and wood from trees and shrubs make it hard to distinguish between high-intensity surface fires, combusting living shrubs and dead wood and leaves, and high-intensity crown fires that have burnt living trees. Distinct particle shape may also influence charcoal transportation, with elongated particles (graminoids) potentially having a more heterogeneous distribution and being deposited farther away from the origin of fire than the rounder, polygonal leaf particles. Despite these limitations, the combined use of charred-particle aspect ratios and fuel morphotypes can aid in the more robust interpretation of fuel source and fire-type changes. Lastly, I highlight the further investigations needed to refine the histories of past wildfires.
Wildfire is the most common disturbance type in boreal forests and can trigger significant changes in forest composition. Waterlogging in peatlands determines the degree of tree cover and the depth of the burning horizon associated with wildfires. However, interactions between peatland moisture, vegetation composition and flammability, and fire regime in forested peatland in Eurasia remain largely unexplored, despite their huge extent in boreal regions. To address this knowledge gap, we reconstructed the Holocene fire regime, vegetation composition and peatland hydrology at two sites in Western Siberia near Tomsk Oblast, Russia. The palaeoecological records originate from forested peatland areas in predominantly light taiga (Pinus-Betula) with increase in dark taiga communities (Pinus sibirica, Picea obovata, Abies sibirica) towards the east. We found that the past water level fluctuated between 8 and 30 cm below the peat surface. Wet peatland conditions promoted broadleaf trees (Betula), whereas dry peatland conditions favoured conifers and a greater forest density (dark-to-light-taiga ratio). The frequency and severity of fire increased with a declining water table that enhanced fuel dryness and flammability and at an intermediate forest density. We found that the probability of intensification in fire severity increased when the water
level declined below 20 cm suggesting a tipping point in peatland hydrology at which wildfire regime intensifies. On a Holocene scale, we found two scenarios of moisture-vegetation-fire interactions. In the first, severe fires were recorded 45 between 7.5 and 4.5 ka BP with lower water level and an increased proportion of dark taiga and fire avoiders (Pinus sibirica at Rybanya and Abies sibirica at Ulukh Chayakh) mixed into the dominantly light taiga and fire-resister community of Pinus
sylvestris. The second occurred over the last 1.5 ka and was associated with fluctuating water tables, a declining abundance of fire avoiders, and an expansion of fire invaders (Betula). These findings suggest that frequent high-severity fires can lead to compositional and structural changes in forests when trees fail to reach reproductive maturity between fire events or where extensive forest gaps limit seed dispersal. This study also shows prolonged periods of synchronous fire activity across the sites, particularly during the early to mid-Holocene, suggesting a regional imprint of centennial to millennial-scale Holocene climate
variability on wildfire activity. Increasing human presence in the region of the Ulukh-Chayakh Mire near Teguldet over the last four centuries drastically enhanced ignitions compared to natural background levels. Frequent warm and dry spells predicted for the future in Siberia by climate change scenarios will enhance peatland drying and may convey a competitive advantage to conifer taxa. However, dry conditions, particularly a water table decline below the threshold of 20 cm, will probably exacerbate the frequency and severity of wildfire, disrupt conifers’ successional pathway and accelerate shifts towards more fire-adapted broadleaf tree cover. Furthermore, climate-disturbance-fire feedbacks will accelerate changes in the carbon balance of forested boreal peatlands and affect their overall future resilience to climate change.
Wildfire is the most common disturbance type in boreal forests and can trigger significant changes in forest composition. Waterlogging in peatlands determines the degree of tree cover and the depth of the burnt horizon associated with wildfires. However, interactions between peatland moisture, vegetation composition and flammability, and fire regime in forest and forested peatland in Eurasia remain largely unexplored, despite their huge extent in boreal regions. To address this knowledge gap, we reconstructed the Holocene fire regime, vegetation composition, and peatland hydrology at two sites located in predominantly light taiga (Pinus sylvestris Betula) with interspersed dark taiga communities (Pinus sibirica, Picea obovata, Abies sibirica) in western Siberia in the Tomsk Oblast, Russia. We found marked shifts in past water levels over the Holocene. The probability of fire occurrence and the intensification of fire frequency and severity increased at times of low water table (drier conditions), enhanced fuel dryness, and an intermediate dark-to-light taiga ratio. High water level, and thus wet peat surface conditions, prevented fires from spreading on peatland and surrounding forests. Deciduous trees (i.e. Betula) and Sphagnum were more abundant under wetter peatland conditions, and conifers and denser forests were more prevalent under drier peatland conditions. On a Holocene scale, severe fires were recorded between 7.5 and 4.5 ka with an increased proportion of dark taiga and fire avoiders (Pinus sibirica at Rybnaya and Abies sibirica at Ulukh–Chayakh) in a predominantly light taiga and fire-resister community characterised by Pinus sylvestris and lower local water level. Severe fires also occurred over the last 1.5 kyr and were associated with a declining abundance of dark taiga and fire avoiders, an expansion of fire invaders (Betula), and fluctuating water tables. These findings suggest that frequent, high-severity fires can lead to compositional and structural changes in forests when trees fail to reach reproductive maturity between fire events or where extensive forest gaps limit seed dispersal. This study also shows prolonged periods of synchronous fire activity across the sites, particularly during the early to mid-Holocene, suggesting a regional imprint of centennial- to millennial-scale Holocene climate variability on wildfire activity. Humans may have affected vegetation and fire from the Neolithic; however, increasing human presence in the region, particularly at the Ulukh–Chayakh Mire over the last 4 centuries, drastically enhanced ignitions compared to natural background levels. Frequent warm and dry spells predicted by climate change scenarios for Siberia in the future will enhance peatland drying and may convey a competitive advantage to conifer taxa. However, dry conditions will probably exacerbate the frequency and severity of wildfire, disrupt conifers' successional pathway, and accelerate shifts towards deciduous broadleaf tree cover. Furthermore, climate–disturbance–fire feedbacks will accelerate changes in the carbon balance of boreal peatlands and affect their overall future resilience to climate change.
Abiotic formation of n-alkane hydrocarbons has been postulated to occur within Earth's crust. Apparent evidence was primarily based on uncommon carbon and hydrogen isotope distribution patterns that set methane and its higher chain homologues apart from biotic isotopic compositions associated with microbial production and closed system thermal degradation of organic matter. Here, we present the first global investigation of the carbon and hydrogen isotopic compositions of n-alkanes in volcanic-hydrothermal fluids hosted by basaltic, andesitic, trachytic and rhyolitic rocks. We show that the bulk isotopic compositions of these gases follow trends that are characteristic of high temperature, open system degradation of organic matter. In sediment-free systems, organic matter is supplied by surface waters (seawater, meteoric water) circulating through the reservoir rocks. Our data set strongly implies that thermal degradation of organic matter is able to satisfy isotopic criteria previously classified as being indicative of abiogenesis. Further considering the ubiquitous presence of surface waters in Earth’s crust, abiotic hydrocarbon occurrences might have been significantly overestimated.
Global modelling of continental water storage changes : sensitivity to different climate data sets
(2007)
Since 2002, the GRACE satellite mission provides estimates of the Earth's dynamic gravity field with unprecedented accuracy. Differences between monthly gravity fields contain a clear hydrological signal due to continental water storage changes. In order to evaluate GRACE results, the state-of-the-art WaterGAP Global Hydrological Model (WGHM) is applied to calculate terrestrial water storage changes on a global scale. WGHM is driven by different climate data sets to analyse especially the influence of different precipitation data on calculated water storage. The data sets used are the CRU TS 2.1 climate data set, the GPCC Full Data Product for precipitation and data from the ECMWF integrated forecast system. A simple approach for precipitation correction is introduced. WGHM results are then compared with GRACE data. The use of different precipitation data sets leads to considerable differences in computed water storage change for a large number of river basins. Comparing model results with GRACE observations shows a good spatial correlation and also a good agreement in phase. However, seasonal variations of water storage as derived from GRACE tend to be significantly larger than those computed by WGHM, regardless of which climate data set is used.
We present the analysis of the impact of convection on the composition of the tropical tropopause layer region (TTL) in West-Africa during the AMMA-SCOUT campaign. Geophysica M55 aircraft observations of water vapor, ozone, aerosol and CO2 during August 2006 show perturbed values at altitudes ranging from 14 km to 17 km (above the main convective outflow) and satellite data indicates that air detrainment is likely to have originated from convective cloud east of the flights. Simulations of the BOLAM mesoscale model, nudged with infrared radiance temperatures, are used to estimate the convective impact in the upper troposphere and to assess the fraction of air processed by convection. The analysis shows that BOLAM correctly reproduces the location and the vertical structure of convective outflow. Model-aided analysis indicates that convection can influence the composition of the upper troposphere above the level of main outflow for an event of deep convection close to the observation site. Model analysis also shows that deep convection occurring in the entire Sahelian transect (up to 2000 km E of the measurement area) has a non negligible role in determining TTL composition.
We present the analysis of the impact of convection on the composition of the tropical tropopause layer region (TTL) in West-Africa during the AMMA-SCOUT campaign. Geophysica M55 aircraft observations of water vapor, ozone, aerosol and CO2 show perturbed values at altitudes ranging from 14 km to 17 km (above the main convective outflow) and satellite data indicates that air detrainment is likely originated from convective cloud east of the flight. Simulations of the BOLAM mesoscale model, nudged with infrared radiance temperatures, are used to estimate the convective impact in the upper troposphere and to assess the fraction of air processed by convection. The analysis shows that BOLAM correctly reproduces the location and the vertical structure of convective outflow. Model-aided analysis indicates that in the outflow of a large convective system, deep convection can largely modify chemical composition and aerosol distribution up to the tropical tropopause. Model analysis also shows that, on average, deep convection occurring in the entire Sahelian transect (up to 2000 km E of the measurement area) has a non negligible role in determining TTL composition.
This study aims to assess the skill of regional climate models (RCMs) at reproducing the climatology of Mediterranean cyclones. Seven RCMs are considered, five of which were also coupled with an oceanic model. All simulations were forced at the lateral boundaries by the ERA-Interim reanalysis for a common 20-year period (1989–2008). Six different cyclone tracking methods have been applied to all twelve RCM simulations and to the ERA-Interim reanalysis in order to assess the RCMs from the perspective of different cyclone definitions. All RCMs reproduce the main areas of high cyclone occurrence in the region south of the Alps, in the Adriatic, Ionian and Aegean Seas, as well as in the areas close to Cyprus and to Atlas mountains. The RCMs tend to underestimate intense cyclone occurrences over the Mediterranean Sea and reproduce 24–40 % of these systems, as identified in the reanalysis. The use of grid nudging in one of the RCMs is shown to be beneficial, reproducing about 60 % of the intense cyclones and keeping a better track of the seasonal cycle of intense cyclogenesis. Finally, the most intense cyclones tend to be similarly reproduced in coupled and uncoupled model simulations, suggesting that modeling atmosphere–ocean coupled processes has only a weak impact on the climatology and intensity of Mediterranean cyclones.
This paper is a contribution to the special issue on Med-CORDEX, an international coordinated initiative dedicated to the multi-component regional climate modelling (atmosphere, ocean, land surface, river) of the Mediterranean under the umbrella of HyMeX, CORDEX, and Med-CLIVAR and coordinated by Samuel Somot, Paolo Ruti, Erika Coppola, Gianmaria Sannino, Bodo Ahrens, and Gabriel Jordà.
Abrupt climate changes of the last deglaciation detected in a Western Mediterranean forest record
(2010)
Abrupt changes in Western Mediterranean climate during the last deglaciation (20 to 6 cal ka BP) are detected in marine core MD95-2043 (Alboran Sea) through the investigation of high-resolution pollen data and pollen-based climate reconstructions by the modern analogue technique (MAT) for annual precipitation (Pann) and mean temperatures of the coldest and warmest months (MTCO and MTWA). Changes in temperate Mediterranean forest development and composition and MAT reconstructions indicate major climatic shifts with parallel temperature and precipitation changes at the onsets of Heinrich stadial 1 (equivalent to the Oldest Dryas), the Bölling-Allerød (BA), and the Younger Dryas (YD). Multi-centennial-scale oscillations in forest development occurred throughout the BA, YD, and early Holocene. Shifts in vegetation composition and (Pann reconstructions indicate that forest declines occurred during dry, and generally cool, episodes centred at 14.0, 13.3, 12.9, 11.8, 10.7, 10.1, 9.2, 8.3 and 7.4 cal ka BP. The forest record also suggests multiple, low-amplitude Preboreal (PB) climate oscillations, and a marked increase in moisture availability for forest development at the end of the PB at 10.6 cal ka BP. Dry atmospheric conditions in the Western Mediterranean occurred in phase with Lateglacial events of high-latitude cooling including GI-1d (Older Dryas), GI-1b (Intra-Allerød Cold Period) and GS-1 (YD), and during Holocene events associated with high-latitude cooling, meltwater pulses and N. Atlantic ice-rafting. A possible climatic mechanism for the recurrence of dry intervals and an opposed regional precipitation pattern with respect to Western-central Europe relates to the dynamics of the westerlies and the prevalence of atmospheric blocking highs. Comparison of radiocarbon and ice-core ages for well-defined climatic transitions in the forest record suggests possible enhancement of marine reservoir ages in the Alboran Sea by 200 years (surface water age 600 years) during the Lateglacial.
Abrupt climate changes of the last deglaciation detected in a western Mediterranean forest record
(2009)
Abrupt changes in Western Mediterranean climate during the last deglaciation (20 to 6 cal ka BP) are detected in marine core MD95-2043 (Alboran Sea) through the investigation of high-resolution pollen data and pollen-based climate reconstructions by the modern analogue technique (MAT) for annual precipitation (Pann) and mean temperatures of the coldest and warmest months (MTCO and MTWA). Changes in temperate Mediterranean forest development and composition and MAT reconstructions indicate major climatic shifts with parallel temperature and precipitation changes at the onsets of Heinrich stadial 1 (equivalent to the Oldest Dryas), the Bölling-Allerød (BA), and the Younger Dryas (YD). Multi-centennial-scale oscillations in forest development occurred throughout the BA, YD, and early Holocene. Shifts in vegetation composition and (Pann reconstructions indicate that forest declines occurred during dry, and generally cool, episodes centred at 14.0, 13.3, 12.9, 11.8, 10.7, 10.1, 9.2, 8.3 and 7.4 cal ka BP. The forest record also suggests multiple, low-amplitude Preboreal (PB) climate oscillations, and a marked increase in moisture availability for forest development at the end of the PB at 10.6 cal ka BP. Dry atmospheric conditions in the Western Mediterranean occurred in phase with Lateglacial events of high-latitude cooling including GI-1d (Older Dryas), GI-1b (Intra-Allerød Cold Period) and GS-1 (YD), and during Holocene events associated with high-latitude cooling, meltwater pulses and N. Atlantic ice-rafting. A possible climatic mechanism for the recurrence of dry intervals and an opposed regional precipitation pattern with respect to Western-central Europe relates to the dynamics of the westerlies and the prevalence of atmospheric blocking highs. Comparison of radiocarbon and ice-core ages for well-defined climatic transitions in the forest record suggests possible enhancement of marine reservoir ages in the Alboran Sea by 200 years (surface water age 600 years) during the Lateglacial.
Holocene climate was characterised by variability on multi-centennial to multi-decadal time scales. In central Europe, these fluctuations were most pronounced during winter. Here we present a record of past winter climate variability for the last 10.8 ka based on four speleothems from Bunker Cave, western Germany. Due to its central European location, the cave site is particularly well suited to record changes in precipitation and temperature in response to changes in the North Atlantic realm. We present high-resolution records of δ18O, δ13C values and Mg/Ca ratios. Changes in the Mg/Ca ratio are attributed to past meteoric precipitation variability. The stable C isotope composition of the speleothems most likely reflects changes in vegetation and precipitation, and variations in the δ18O signal are interpreted as variations in meteoric precipitation and temperature. We found cold and dry periods between 8 and 7 ka, 6.5 and 5.5 ka, 4 and 3 ka as well as between 0.7 and 0.2 ka. The proxy signals in the Bunker Cave stalagmites compare well with other isotope records and, thus, seem representative for central European Holocene climate variability. The prominent 8.2 ka event and the Little Ice Age cold events are both recorded in the Bunker Cave record. However, these events show a contrasting relationship between climate and δ18O, which is explained by different causes underlying the two climate anomalies. Whereas the Little Ice Age is attributed to a pronounced negative phase of the North Atlantic Oscillation, the 8.2 ka event was triggered by cooler conditions in the North Atlantic due to a slowdown of the thermohaline circulation.