550 Geowissenschaften
Refine
Year of publication
Document Type
- Article (896) (remove)
Has Fulltext
- yes (896)
Keywords
- Climate change (7)
- climate change (7)
- COSMO-CLM (6)
- Palaeoclimate (6)
- Atmospheric chemistry (5)
- Biogeochemistry (5)
- Palaeoceanography (5)
- precipitation (5)
- loess (4)
- Clausius–Clapeyron scaling (3)
Institute
- Geowissenschaften (546)
- Geowissenschaften / Geographie (100)
- Biodiversität und Klima Forschungszentrum (BiK-F) (59)
- Senckenbergische Naturforschende Gesellschaft (51)
- Extern (22)
- Biowissenschaften (21)
- Geographie (15)
- Präsidium (12)
- Institut für Ökologie, Evolution und Diversität (9)
- Physik (6)
Ice-nucleating particle concentrations of the past: insights from a 600-year-old Greenland ice core
(2020)
Ice-nucleating particles (INPs) affect the microphysics in cloud and precipitation processes. Hence, they modulate the radiative properties of clouds. However, atmospheric INP concentrations of the past are basically unknown. Here, we present INP measurements from an ice core in Greenland, which dates back to the year 1370. In total 135 samples were analyzed with the FRIDGE droplet freezing assay in the temperature range from −14 to −35 ∘C. The sampling frequency was set to 1 in 10 years from 1370 to 1960. From 1960 to 1990 the frequency was increased to one sample per year. Additionally, a few special events were probed, including volcanic episodes. The typical time coverage of a sample was on the order of a few months. Historical atmospheric INP concentrations were estimated with a conversion factor, which depends on the snow accumulation rate of the ice core, particle dry deposition velocity, and wet scavenging ratio. Typical atmospheric INP concentrations were on the order of 0.1 L−1 at −25 ∘C. The INP variability was found to be about 1–2 orders of magnitude. Yet, the short-term variability from samples over a seasonal cycle was considerably lower. INP concentrations were significantly correlated to some chemical tracers derived from continuous-flow analysis (CFA) and ion chromatography (IC) over a broad range of nucleation temperatures. The highest correlation coefficients were found for the particle concentration (spherical diameter dp > 1.2 µm). The correlation is higher for a time period of seasonal samples, where INP concentrations follow a clear annual pattern, highlighting the importance of the annual dust input in Greenland from East Asian deserts during spring. Scanning electron microscopy (SEM) analysis of selected samples found mineral dust to be the dominant particle fraction, verifying their significance as INPs. Overall, the concentrations compare reasonably well to present-day INP concentrations, albeit they are on the lower side. However, we found that the INP concentration at medium supercooled temperatures differed before and after 1960. Average INP concentrations at −23, −24, −25, −26, and −28 ∘C were significantly higher (and more variable) in the modern-day period, which could indicate a potential anthropogenic impact, e.g., from land-use change.
In this study, we use simulations from seven global vegetation models to provide the first multi‐model estimate of fire impacts on global tree cover and the carbon cycle under current climate and anthropogenic land use conditions, averaged for the years 2001–2012. Fire globally reduces the tree covered area and vegetation carbon storage by 10%. Regionally, the effects are much stronger, up to 20% for certain latitudinal bands, and 17% in savanna regions. Global fire effects on total carbon storage and carbon turnover times are lower with the effect on gross primary productivity (GPP) close to 0. We find the strongest impacts of fire in savanna regions. Climatic conditions in regions with the highest burned area differ from regions with highest absolute fire impact, which are characterized by higher precipitation. Our estimates of fire‐induced vegetation change are lower than previous studies. We attribute these differences to different definitions of vegetation change and effects of anthropogenic land use, which were not considered in previous studies and decreases the impact of fire on tree cover. Accounting for fires significantly improves the spatial patterns of simulated tree cover, which demonstrates the need to represent fire in dynamic vegetation models. Based upon comparisons between models and observations, process understanding and representation in models, we assess a higher confidence in the fire impact on tree cover and vegetation carbon compared to GPP, total carbon storage and turnover times. We have higher confidence in the spatial patterns compared to the global totals of the simulated fire impact. As we used an ensemble of state‐of‐the‐art fire models, including effects of land use and the ensemble median or mean compares better to observational datasets than any individual model, we consider the here presented results to be the current best estimate of global fire effects on ecosystems.
Aim: Plant life‐forms characterize key morphological strategies that enable large‐scale comparisons of plant communities. This study applies Raunkiær's plant life‐form concept that was developed for temperate climate to a subtropical island flora, in parts, dominated by summer aridity. We quantify how plant life‐form patterns as well as patterns of important plant functional traits (PFTs) relate to important climate and topographic characteristics.
Location: La Palma, Canary Islands.
Taxon: Flora of La Palma.
Methods: We assigned each native plant species a plant life‐form, that is, phanerophyte, chamaephyte, hemicryptophyte, geophyte and therophyte, as well as PFTs (succulence and N‐fixer). We used stacked species distribution models to assess occurrence probability for each species using the Atlantis database (500 m × 500 m grid). We related richness and percentage values for each plant life‐form and PFT to climate and topography.
Results: Plant life‐forms and PFTs showed a clear pattern within geographic but also climate space, while topography had a minor effect. Phanerophytes mainly contributed to the flora in humid areas. Chamaephytes and hemicryptophytes most strongly contributed to the summit scrub flora and, to some degree, also to the arid coastal regions. Geophytes and therophytes were mainly found in dry coastal regions. N‐fixers contributed mainly to warm‐arid and cool‐arid regions, while succulent species were mainly found in arid coastal regions.
Main conclusions: Raunkiær's plant life‐form concept can be comprehensively transferred to a subtropical island flora by adapting to local unfavourable growing conditions, that is, aridity. Using the strong environmental gradients offered by our study island, we identify substantial climate‐driven variation in patterns of plant life‐forms and PFTs that might be used for large‐scale comparisons in macroecological studies. The growth strategies reflected in Raunkiær's plant life‐forms suggest differences in species establishment and coexistence dynamics within different parts of the island's climate space.
The complexity of atmospherical processes has always yielded a multitude of ways of knowing about the weather. What has been lacking in the historiography of meteorology so far is a way to formulate differences between forms of knowledge in a way that does not privilege modern scientific structures, but focuses instead on the epistemological category of causality. Using causality as ground of comparison for different knowledge claims, I shall argue, may enable researchers to investigate meteorological knowledge across time periods, perhaps even geographical regions, in a more symmetrical manner. This review demonstrates this approach as a means to organize a large set of historical meteorological writings from German countries between 1750 and 1850. Three distinct forms of knowledge (Semiotics, Physics, and Organics of the weather) during that time and in that region are suggested and will be described. While a bibliography with a national perspective from the 1880s was the basis for the selection of historical sources, such a setup proved awkward even to contemporaries. In addition, the bibliography came with a number of biases and shortcomings that will be critically reviewed.
We apply seismic full waveform inversion to SH‐ and Love‐wave data for investigating the near‐surface lithology at an archaeological site. We evaluate the resolution of the applied full waveform inversion algorithm through ground truthing in the form of an excavation and sediment core studies. Thereby, we investigate the benefits of full waveform inversion in comparison with other established methods of near‐surface prospecting in terms of resolution capabilities and interpretation security. The study is performed in a presumed harbour area of the ancient Thracian city of Ainos. The exemplary target is the source of a linear magnetic anomaly oriented perpendicular to the coast, which was found in a previous magnetic gradiometry survey, suggesting a mole. The SH‐wave full waveform inversion recovered a subsurface SH‐wave velocity model with submeter resolution showing lateral and vertical velocity variation between 40 and 150 m/s. To tame the non‐linearity of the full waveform inversion, a sequential inversion of frequency bands has to be combined with time‐windowing in order to separate the Love wave from the reflected SH wavefield. We compare the full waveform inversion results with multichannel analysis of surface waves, standard seismic reflection imaging, electrical resistivity tomography and electromagnetic induction. It turns out that the respective depth sections are correlated to a certain degree with the full waveform inversion results. However, the structural resolution of the other geophysical methods is significantly lower than for the full waveform inversion. An exception is the reflection seismic imaging, which shows the same resolution as full waveform inversion but can only be interpreted together with the full waveform inversion–based velocity model. An archaeological excavation as well as coring data allows ground truthing and a direct understanding of the geophysical structures. The results show that the target was a sort of near‐surface trench of about 3–4 m width and 0.8 m to 1.0 m depth, filled with silty sediment, which differs from the layered surrounding in colour and composition. The ground truthing revealed that only SH‐wave full waveform inversion and seismic reflection imaging could image the trench and sediment structure with satisfying lateral and depth resolution. We emphasize that the velocity distribution from SH‐wave full waveform inversion agrees closely with the excavated subsurface structures, and that the discovered changes in seismic velocity correlate with changes in the sand content in the respective sediment facies sequences. The study demonstrated that SH‐wave full waveform inversion is capable to image structural and lithological changes in the near subsurface at scales as low as 0.5 m, thus providing the high resolution needed for archaeological and geoarchaeological prospection.
Opportunities and challenges for paleoaltimetry in "small" orogens: insights from the European Alps
(2020)
Many stable isotope paleoaltimetry studies have focused on paleoelevation reconstructions of orogenic plateaus such as the Tibetan or Andean Plateaus. We address the opportunities and challenges of applying stable isotope paleoaltimetry to “smaller” orogens. We do this using a high‐resolution isotope tracking general circulation model (ECHAM5‐wiso) and explore the precipitation δ18O (δ18Op) signal of Cenozoic paleoclimate and topographic change in the European Alps. Results predict a maximum δ18Op change of 4–5‰ (relative to present day) during topographic development of the Alps. This signal of topographic change has the same magnitude as changes in δ18Op values resulting from Pliocene and Last Glacial Maximum global climatic change. Despite the similar magnitude of the isotopic signals resulting from topographic and paleoclimate changes, their spatial patterns across central Europe differ. Our results suggest that an integration of paleoclimate modeling, multiproxy approaches, and low‐elevation reference proxy records distal from an orogen improve topographic reconstructions.
Global‐scale gradient‐based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady‐state hydraulic heads simulated by two versions of the global groundwater model G3M, at spatial resolutions of 5 arc‐minutes (9 km) and 30 arc‐seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low‐resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global‐scale groundwater models will not improve the simulation of the global freshwater system.
Barriereinseln und Außensände, geformt durch eine Kombination aus Wind, Wellen, Strömung und Küstenlängstransport gelten als morphologisch hoch aktive Küstenbereiche und variieren häufig in Ursprung, Genese und Entwicklung. Sie besitzen durch ihre dissipative Wirkungsweise eine bedeutende Schutzfunktion für rückwärtige Inseln, Halligen und Festlandbereiche und bilden vor der Küste Nordfrieslands die westliche Außengrenze des Wattenmeeres.
Ziel der nachfolgenden Studie ist es, anhand hochauflösender Georadarmessungen und sedimentologischer Daten aus Bohrungen die Landschaftsgeschichte an der Westküste Amrums sowie die Entstehung und interne sedimentäre Architektur des der Insel westlich angelagerten Kniepsandes zu untersuchen und die Prozesse, die zur Genese geführt haben, durch Datierungen zeitlich einzuordnen. Auf Grundlage der gewonnenen Daten wurden zwei stratigraphische Modelle entwickelt, welche die geologisch-geomorphologischen Prozesse und Sedimentationsbedingungen im Westküstenvorfeld erklären.
Ein Modell zeigt die Landschaftsentwicklung an der Westküste Amrums und beschreibt die Sedimentationsbedingungen, die im Vorfeld der Westküste herrschten, bevor der Kniepsand an die Insel heranwanderte. Auf der Landoberfläche des ertrinkenden saaleeiszeitlichen Geestkerns wurden zu Beginn der Flandrischen Transgression feinkörnige Sedimente eines Mischund Schlickwatts abgelagert. Es ist davon auszugehen, dass der damals noch weit vor der Küste Amrums liegende Kniepsand eine Barriere bildete und so an der heute hochenergetischen Westküste für strömungsberuhigte Sedimentationsbedingungen sorgte. Durch Erosion am Geestkern bildeten sich in unterschiedlichen Höhenpositionen fossile Kliffs, die dem damaligen Meeresspiegelstand entsprechen.
Ein weiteres Modell beschreibt den Andockmechanismus des Kniepsandes an die Insel Amrum. Durch die Anlagerung des ehemaligen Außensandes und den damit einhergehenden Sedimentinput wurden die Bedingungen für eine großflächige Dünenbildung geschaffen.
Als erste Ergebnisse der derzeit laufenden paläomagnetischen Untersuchungen im Rhein-Main-Gebiet werden die Befunde von einigen stratigraphisch bedeutenden Profilen mitgeteilt. Während die Remanenzwerte aus altpleistozänen Tonen der Kelsterbacher Terrasse mit einer Ausnahme revers sind, konnte im Autobahneinschnitt Abenheim (NW Worms) und im Lößprofil der Ziegeleigrube Bad Soden am Taunus der Ubergang von reverser zu normaler Magnetisierung gefunden werden. Die Grenze zwischen MATUYAMA- und BRUNHES-Epoche liegt im Profil Bad Soden unter dem 6., das Jaramillo-Event unter dem 7. fossilen Bt-Horizont.
Die über dem letztinterglazialen Boden (Parabraunerde bzw. pseudovergleyte Parabraunerde) folgenden würmzeitlichen äolischen Sedimente können durch kennzeichnende Bodenhorizonte in 3 Abschnitte (Alt-, Mittel- und Jungwürm) gegliedert werden. Charakteristisch für das Altwürm sind Humuszonen, für das Mittelwürm neben einigen „Naßböden" braune Verlehmungszonen bis zu einer Mächtigkeit von 1,1 m und für das Jungwürm mehrere schwach ausgebildete, geringmächtige Verbraunungszonen und „Naßböden". Den wichtigsten Leithorizont des Jungwürms bildet das Kärlicher Tuffbändchen, das in jüngster Zeit auch in Nordhessen aufgefunden wurde. Abschließend wird das in Hessen aufgestellte Gliederungsschema mit den Würmlöß-Gliederungen in anderen Teilen Europas verglichen und eine Parallelisierung versucht.
Der Auffassung, das Substrat der Lockerbraunerde auf dem Oberwald des Vogelsberges sei ausschließlich ein äolisches Sediment der Jüngeren Tundrenzeit, wird widersprochen. Verschiedene Befunde, vor allem 14C-Datierungen, lassen vielmehr den Schluß zu, daß große Teile des Substrates vieler Lockerbraunerden holozäne anthropogene Kolluvien sind. Bekräftigt wird dagegen der Befund gleichen (jungtundrenzeitlichen) Alters für Lockerbraunerde-Substrat und Deckschutt (Decksediment).
Aus dem Gebiet der weichselzeitlichen Vereisung in Polen werden allgemein verbreitete geringmächtige periglaziale Deckschichten beschrieben. Sie zeigen in der Regel eine äolische Beeinflussung und unterscheiden sich dadurch vom Liegenden. Es handelt sich um spätglaziale Bildungen, wie sie aus dem Jungmoränengebiet der DDR seit langem bekannt sind. Ähnliche Substrate wurden auch im nördlichen Alpenvorland gefunden.
Im Lößprofil der Ziegelei Glos bei Lisieux ist eine Abfolge fossiler Boden aufgeschlossen, die sich gut mit der hessischen Lößstratigraphie verbinden läßt. Über dem letztinterglazialen Boden liegen nach einer Diskordanz der Lohner Boden, der E1-, der E2- und der E4-Naß-boden. Letzterer wird häufig als Äquivalent des „Sol de Kesselt" angesehen. Somit zeigt sich auch für dieses Gebiet, daß dieser Boden stratigraphisch nicht dem Lohner Boden entsprechen kann.
Un tuf volcanique a ete mis en evidence dans trois profils de loess de Hesbaye (Rocourt et Lixhe) et du Limbourg neerlandais (Nagelbeek). Ce tuf correspond probablement ä l'Eltviller Tuff connu en Hesse (SEMMEL 1967). En consequence, cette decouverte permet de comparer directement une partie de la sequence des loess du Pleistocene superieur propre ä chaque region
In der Monte Cavallo-Gruppe fand Verf. in Stauseesedimenten fossile Holzstücke, für die die 14 C-Datierung ein Alter von 29 350 ± 460 Jahren vor 1950 n. Chr. ergab. Die schluffigen Ablagerungen, die in etwa 900 m, fast 80 m über dem heutigen Flußbett des T. Caltea aufgeschlossen sind, ruhen auf Schottern und werden von Moräne überlagert.
Aufgrund seiner stratigraphischen Lage und seines 14C-datierten Alters kann das Holz (Picea abies bzw. Larix) einem zeitlichen Äquivalent des Paudorf-Interstadials zugeordnet werden. Dieser Fund von Großresten ermöglicht somit eine erste absolute Datierung des Paudorf-Interstadials für die Südalpen und Oberitalien.
Das Cranium eines fossilen Hominiden des Formenkreises Homo sapiens sapiens wurde relativ-geologisch sowie absolut durch Radiokohlenstoff und Aminosäuren auf ungefähr 31 000 Jahre B.P. datiert. Andere absolute sowie relative Daten an Mollusken und Mammutzähnen in überlagernden jüngeren Straten datieren auf 18 000 — 21000 und 16 000 Jahren B.P. Geomorphologische und geophysikalische Datierungen stimmen somit gut überein. Er ist der älteste datierte und früheste Bewohner Zentraleuropas, der dem Homo sapiens sapiens angehört.
Aus Franken wird die Entwicklung quartärer Hohlformen beschrieben, deren Rekonstruktion mit Hilfe lößstratigraphischer Methoden (fossile Böden, Tuffbänder, Umlagerungszonen etc.) möglich ist. Bei vielen Formen zeigt sich, daß sie bereits größere Vorläuferformen präwürmzeitlichen Alters hatten. Die Entwicklung während des Würms läßt sich an manchen Beispielen in besonders instruktiver Weise verfolgen. Zu Beginn des Würms, im unteren Mittelwürm und im unteren Jungwürm dominierte zeitweise die Abtragung und Verlagerung. Im oberen Mittelwürm sowie im oberen Jungwürm herrschte äolische Lößsedimentation vor. Diese Ergebnisse stimmen gut mit den bereits aus anderen mitteleuropäischen Lößgebieten bekannten Befunden überein. Mit dem Trockental-System von Helmstadt wird die Entwicklung von Hohlformen beschrieben, deren Anlage bis in das ältere Pleistozän zurückreicht.
In vorliegender Untersuchung wurde der Rißlöß zwischen der 1. und 2. fossilen Parabraunerde anhand schwächerer Bodenbildungen und eingeschalteter Abtragungsphasen zu gliedern versucht. Im jüngeren Riß herrschte starke Lößsedimentation vor, wobei es in mindestens 6 kaltfeuchten Abschnitten zur Ausbildung schwacher periglazialer Naßböden kam. Die Naßbodenserie wurde als Bruchköbeler Böden (B) bezeichnet. Im jüngsten Rißlöß ist wenige dm unter dem Eemboden als tephrochronologischer Leithorizont der Krifteler Tuff (vgl. SEMMEL 1968) eingeschaltet. Den mittleren Profilbereich im Rißlöß zeichnen feuchtere Klimaabschnitte mit starken Verschwemmungsphasen aus, die in den meisten Profilen zu erheblichen Diskordanzen geführt haben. An der Basis der wenigen kompletten Rißlöß-Profile treten vorwiegend in Hessen über der zumeist gekappten 2. fossilen Parabraunerde maximal zwei Schwarzerden auf, die von SEMMEL (1968) als Weilbacher Humuszonen bezeichnet werden. Unmittelbar über diesen Schwarzerden folgt die Ostheimer Zone, eine Fließerde aus aufgearbeitetem Solumaterial der liegenden Böden. Insgesamt zeigt die aus den Rißböden rekonstruierte Klimaabfolge — neben geringfügigen Abweichungen — überraschende Parallelen zur paläopedologisch-klimatischen Gliederung der Würmkaltzeit.
Crustal structure at the western end of the North Anatolian Fault Zone from deep seismic sounding
(2001)
The first deep seismic sounding experiment in Northwestern Anatolia was carried out in October 1991 as part of the "German - Turkish Project on Earthquake Prediction Research" in the Mudurnu area of the North Anatolian Fault Zone. The experiment was a joint enterprise by the Institute of Meteorology and Geophysics of Frankfurt University, the Earthquake Research Institute (ERI) in Ankara, and the Turkish Oil Company (TPAO). Two orthogonal profiles, each 120 km in length with a crossing point near Akyazi, were covered in succession by 30 short period tape recording seismograph stations with 2 km station spacing. 12 shots, with charge sizes between 100 and 250 kg, were fired and 342 seismograms out of 360 were used for evaluation. By coincidence an M b = 4.5 earthquake located below Imroz Island was also recorded and provided additional information on Moho and the sub-Moho velocity. A ray tracing method orginally developed by Weber (1986) was used for travel time inversion. From a compilation of all data two generalized crustal models were derived, one with velocity gradients within the layers and one with constant layer velocities. The latter consists of a sediment cover of about 2 km with V p » 3.6 km/s, an upper crystalline crust down to 13 km with V p » 5.9 km/s, a middle crust down to 25 km depth with V p » 6.5 km/s, a lower crust down to 39 km Moho depth with V p » 7.0 km/s and V p » 8.05 km/s below the Moho. The structure of the individual profiles differs slightly. The thickest sediment cover is reached in the Izmit-Sapanca-trough and in the Akyazi basin. Of particular interest is a step of about 4 km in the lower crust near Lake Sapanca and probably an even larger one in the Moho (derived from the Imroz earthquake data). After the catastrophic earthquake of Izmit on 17 August 1999 this significant heterogeneity in crustal structure appears in a new light with regard to the possible cause of the Izmit earthquake. Heterogeneities in structure are frequently also heterogeneities in strength and stress that impede or even lock rupture. The Izmit earthquake is discussed in relation to a large stepover or jog at the North Anatolian Fault.
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 paper provides new geochemical data focusing on valuable elements in the coal, parting, and floor samples in the No. 5 coal seam of the Taiyuan Formation from the Wujiawan mine, Datong coalfield, northern China. The minerals mainly consist of kaolinite, calcite, and pyrite, as well as trace amounts of quartz and illite. The No. 5 coal is enriched in Li, Ga, high field strength elements (HFSEs), and rare earth elements and yttrium (REY) when compared with world hard coals. Of particular interest is the high average concentration of Li (67.66 μg/g), which is around seven times higher than the value for world hard coals. Lithium, Ga, and HFSEs have strong inorganic affinities, whereas REY have organic affinities. The main carrier of Li, Ga, and HFSEs is aluminosilicate minerals, while REY appear to occur with organophosphorus. These HFSEs are enriched, both in the parting and in the adjacent coal samples. This suggests that these elements are likely to leach out during the diagenetic process. The distribution patterns of REY, along with the ratio of Al2O3/TiO2 and the figure of Zr/TiO2 vs. Nb/Y are suggestive of their derivation from felsic parent material. In the northern and eastern part of the Datong coalfield, there are several regions where the Li content is higher than the mineable grade, in particular in the northern Datong coalfield where there is a mine with an Li content of 294.6 μg/g. This is significantly higher than the mineable grade. Therefore, there is a potential for financially viable recovery of Li in these coals of the Datong coalfield.
Tropical cyclones (TC) represent a substantial threat to life and property for Caribbean and adjacent populations. The prospective increase of TC magnitudes, expressed in the 15th chapter of the IPCC AR5 report, entails a rising probability of ecological and social disasters, which were tragically exemplified by several severe Caribbean TC strikes during the past 20 years. Modern IPCC-grade climate models, however, still lack the required spatial and temporal resolution to accurately consider the underlying boundary conditions that modulate long-time TC patterns beyond the Instrumental Era. It is thus necessary to provide a synoptic mechanistic understanding regarding the origin of such long-time patterns, in order to predict reliable changes of TC magnitude and frequency under future climate scenarios. Caribbean TC records are still rare and often lack the necessary continuity and resolution to overcome these limitations. Here, we report on an annually-resolved sedimentary archive from the bottom of the Great Blue Hole (Lighthouse Reef, Belize). The TC record encompasses 1885 years and extends all existing site-specific TC archives both in terms of resolution and duration. We identified a likely connection between long-term TC patterns and climate phenomena responses to Common Era climate variations and offer a conceptual and comparative view considering several involved tropospheric and oceanographic control mechanisms such as the El-Niño-Southern-Oscillation, the North Atlantic Oscillation and the Atlantic Multidecadal Oscillation. These basin-scaled climate modes exercise internal control on TC activity by modulating the thermodynamic environment (sea-surface temperature and vertical wind shear stress dynamics) for enhanced/suppressed TC formation both on millennial (primary) and multi-decadal (secondary) time scales. We interpret the beginning of the Medieval Warm Period (MWP) as an important time interval of the Common Era record and suspect that the southward migration of the intertropical convergence zone (ITCZ) caused, in combination with extensive hydro-climate changes, a shift in the tropical Atlantic TC regime. The TC activity in the south-western Caribbean changed in general from a stable and less active stage (100–900 CE) to a more active and variable state (1,100 CE-modern).
The Early Permian coal is of great value in the Tengxian Coalfield, Shandon Province, Eastern China. This work deals with the new data focusing on mineralogical characteristics in the Early Permian Shanxi Formation No. 3 coal from the Jinyuan Mine. The Jinyuan coal is a low ash and highly volatile A bituminous coal. Minerals in the No. 3 coal mainly comprise of kaolinite, ankerite, illite, calcite, siderite, and quartz, with varying compositions of trace amounts of pyrite, jarosite, bassanite, anatase, and rutile. According to mineral assemblage in the coal plies, three Types (A to C) can be identified in the No. 3 coal. The dominant minerals in Type A are poorly-ordered kaolinite, illite, quartz, pyrite, and jarosite. Type B is mainly composed of well-ordered kaolinite, illite, siderite, ankerite, and calcite. Type C, with just one sample (JY-3-7c), which contains high proportions of calcite (54%) and ankerite (34%). Terrigenous minerals are elevated in coal plies that typically have relatively high contents of ash yield. The formation of syngenetic pyrite was generally due to seawater, while the sulphate minerals (jarosite and coquimbite) were derived from the oxidation of pyrite. Epigenetic vein-like or fracture-fillings carbonate minerals (ankerite, calcite, and siderite), kaolinite, and pyrite, as well as authigenic quartz were derived from the influx of hydrothermal fluids during different periods, from the authigenic to epigenetic. The paragonite in the coal may have been formed by the precipitated from Na-rich hydrothermal fluids. No effects of magmatic intrusion on mineralogy were investigated in this research.
We present novel measurements of five short-lived brominated source gases (CH2Br2, CHBr3, CH2ClBr, CHCl2Br and CHClBr2) obtained using a gas chromatograph-mass spectrometer system on board the High Altitude and Long Range Research Aircraft (HALO). The instrument is extremely sensitive due to the use of chemical ionisation, allowing detection limits in the lower parts per quadrillion (10-15) range. Data from three campaigns using the HALO aircraft are presented, where the Upper Troposphere/Lower Stratosphere (UTLS) of the Northern Hemisphere mid to high latitudes were sampled during winter and during late summer to early fall. We show that an observed decrease with altitude in the stratosphere is consistent with the relative lifetimes of the different compounds. Distributions of the five source gases and total organic bromine just below the tropopause shows an increase in mixing ratio with latitude, in particular during polar winter. This increase in mixing ratio is explained by increasing lifetimes at higher latitudes during winter. As the mixing ratio at the extratropical tropopause are generally higher than those derived for the tropical tropopause, extratropical troposphere-to-stratosphere transport will result in elevated levels of organic bromine in comparison to air transported over the tropical tropopause. The observations are compared to model estimates using different emission scenarios. A scenario which has emissions most strongly concentrated to low latitudes cannot reproduce the observed latitudinal distributions and will tend to overestimate bromine input through the tropical tropopause from CH2Br2 and CHBr3. Consequently, the scenario also overestimates the amount of brominated organic gases in the stratosphere. The two scenarios with the highest overall emissions of CH2Br2 tend to overestimate mixing ratios at the tropical tropopause but are in much better agreement with extratropical tropopause values, showing that not only total emissions but also latitudinal distributions in the emissions are of importance. While an increase in tropopause values with latitude is reproduced with all emission scenarios during winter, the simulated extratropical tropopause values are on average lower than the observations during late summer to fall. We show that a good knowledge of the latitudinal distribution of tropopause mixing ratios and of the fractional contributions of tropical and extratropical air is needed to derive stratospheric inorganic bromine in the lowermost stratosphere from observations. Depending on the underlying emission scenario, differences of a factor 2 in reactive bromine derived from observations and model outputs are found for the lowermost stratosphere, based on source gas injection. We conclude that a good representation of the contributions of different source regions is required in models for a robust assessment of the role of short-lived halogen source gases on ozone depletion in the UTLS.
Background: Point of care devices for performing targeted coagulation substitution in bleeding patients have become increasingly important in recent years. New on the market is the Quantra® from HemoSonics (LC, Charlottesville, VA, US). It uses sonorheometry, a sonic estimation of elasticity via resonance (SEER), a novel ultrasound-based technology that measures viscoelastic properties of whole blood. Several studies have already shown the comparability with devices already established on the market such as the ROTEM® (TEM International GmbH, Munich, Germany).
Objective: In contrast to existing studies, the planned study will be the first prospective interventional study using the new Quantra® system in a cardiac surgical patient cohort. The aim is to investigate the non-inferiority between an already existing coagulation algorithm, based on ROTEM®/Multiplate®, and a new algorithm based on the Quantra®, for the treatment of coagulopathic cardiac surgical patients.
Methods: The study is divided into two phases. In an initial observation phase, whole blood samples of 20 patients will be analyzed using both ROTEM®/Multiplate® and Quantra® obtained at three defined points of time (prior to surgery, after completion of cardiopulmonary bypass, on arrival in the intensive care unit). The obtained threshold values will be used to create an algorithm for hemotherapy. In a second intervention phase, the new algorithm will be tested against an algorithm used routineously for years at our department for non-inferiority.
Results: The main objective of the examination is the cumulative loss of blood within 24 hours after surgery. Statistical calculations based on literature and in-house data suggest that the new algorithm is not inferior if the difference in cumulative blood loss is < 150ml/24 h.
Conclusions: Because of the comparability of the Quantra® sonorheometry system with ROTEM® rotational thromboelastometric measurement methods, the existing hemotherapy treatment algorithm can be adapted to the Quantra device with a proof of non-inferiority. Clinical Trial: International Registered Report Identifier (IRRID): clinicaltrials.gov: NCT03902275
Inclusions of breyite (previously known as walstromite-structured CaSiO3) in diamond are usually interpreted as retrogressed CaSiO3 perovskite trapped in the transition zone or the lower mantle. However, the thermodynamic stability field of breyite does not preclude its crystallization together with diamond under upper-mantle conditions (6–10 GPa). The possibility of breyite forming in subducted sedimentary material through the reaction CaCO3 + SiO2 = CaSiO3 + C + O2 was experimentally evaluated in the CaO–SiO2–C–O2 ± H2O system at 6–10 GPa, 900–1500 ∘C and oxygen fugacity 0.5–1.0 log units below the Fe–FeO (IW) buffer. One experimental series was conducted in the anhydrous subsystem and aimed at determining the melting temperature of the aragonite–coesite (or stishovite) assemblage. It was found that melting occurs at a lower temperature (∼1500 ∘C) than the decarbonation reaction, which indicates that breyite cannot be formed from aragonite and silica under anhydrous conditions and an oxygen fugacity above IW – 1. In the second experimental series, we investigated partial melting of an aragonite–coesite mixture under hydrous conditions at the same pressures and redox conditions. The melting temperature in the presence of water decreased strongly (to 900–1200 ∘C), and the melt had a hydrous silicate composition. The reduction of melt resulted in graphite crystallization in equilibrium with titanite-structured CaSi2O5 and breyite at ∼1000 ∘C. The maximum pressure of possible breyite formation is limited by the reaction CaSiO3 + SiO2 = CaSi2O5 at ∼8 GPa. Based on the experimental results, it is concluded that breyite inclusions found in natural diamond may be formed from an aragonite–coesite assemblage or carbonate melt at 6–8 GPa via reduction at high water activity.
The accurate knowledge of the groundwater storage variation (ΔGWS) is essential for reliable water resource assessment, particularly in arid and semi-arid environments (e.g., Australia, the North China Plain (NCP)) where water storage is significantly affected by human activities and spatiotemporal climate variations. The large-scale ΔGWS can be simulated from a land surface model (LSM), but the high model uncertainty is a major drawback that reduces the reliability of the estimates. The evaluation of the model estimate is then very important to assess its accuracy. To improve the model performance, the terrestrial water storage variation derived from the Gravity Recovery And Climate Experiment (GRACE) satellite mission is commonly assimilated into LSMs to enhance the accuracy of the ΔGWS estimate. This study assimilates GRACE data into the PCRaster Global Water Balance (PCR-GLOBWB) model. The GRACE data assimilation (DA) is developed based on the three-dimensional ensemble Kalman smoother (EnKS 3D), which considers the statistical correlation of all extents (spatial, temporal, vertical) in the DA process. The ΔGWS estimates from GRACE DA and four LSM simulations (PCR-GLOBWB, the Community Atmosphere Biosphere Land Exchange (CABLE), the Water Global Assessment and Prognosis Global Hydrology Model (WGHM), and World-Wide Water (W3)) are validated against the in situ groundwater data. The evaluation is conducted in terms of temporal correlation, seasonality, long-term trend, and detection of groundwater depletion. The GRACE DA estimate shows a significant improvement in all measures, notably the correlation coefficients (respect to the in situ data) are always higher than the values obtained from model simulations alone (e.g., ~0.15 greater in Australia, and ~0.1 greater in the NCP). GRACE DA also improves the estimation of groundwater depletion that the models cannot accurately capture due to the incorrect information of the groundwater demand (in, e.g., PCR-GLOBWB, WGHM) or the unavailability of a groundwater consumption routine (in, e.g., CABLE, W3). In addition, this study conducts the inter-comparison between four model simulations and reveals that PCR-GLOBWB and CABLE provide a more accurate ΔGWS estimate in Australia (subject to the calibrated parameter) while PCR-GLOBWB and WGHM are more accurate in the NCP (subject to the inclusion of anthropogenic factors). The analysis can be used to declare the status of the ΔGWS estimate, as well as itemize the possible improvements of the future model development.
Stability of halocarbons in whole air samples from the upper
troposphere and lowermost stratosphere
(2019)
Measurements of halogenated hydrocarbons of ambient air frequently rely on canister sampling followed by offline laboratory analysis. This allows for a large number of compounds to be analysed under stable conditions, maximising measurement precision. However, individual compounds might be affected during storage of canister samples. In order to assess halocarbon stability in whole air samples from the upper troposphere and lowermost stratosphere, we performed stability tests using the air sampling unit High REsolution Sampler (HIRES) which is part of the CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) instrument package. HIRES holds 88 light-weight stainless steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES sampling unit was first deployed in 2010, but has up to now not been used for regular halocarbon analysis with exception of chloromethane. The sample collection unit was tested for sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless steel canisters during storage of up to five weeks and on the influence of ozone, since flights take place in the upper troposphere and lowermost stratosphere with ozone mixing ratios of up to several hundred ppbV. Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to five weeks and were unaltered by ozone being present during pressurization. Some compounds such as for example dichloromethane, trichloromethane and tetrachloroethene started to decrease in the canisters after a storage time of more than two weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. Few compounds such as for example tetrachloromethane and tribromomethane were found to be not stable in the HIRES stainless-steel canisters independent of ozone levels. Also growth was observed during storage, namely for HFC-152a and HFC-23.
Measurements of halogenated trace gases in ambient air frequently rely on canister sampling followed by offline laboratory analysis. This allows for a large number of compounds to be analysed under stable conditions, maximizing measurement precision. However, individual compounds might be affected during the sampling and storage of canister samples. In order to assess halocarbon stability in whole-air samples from the upper troposphere and lowermost stratosphere, we performed stability tests using the high-resolution sampler (HIRES) air sampling unit, which is part of the Civil Aircraft for the Regular Investigation of the atmosphere Based on an Instrument Container (CARIBIC) instrument package. The HIRES unit holds 88 lightweight stainless-steel cylinders that are pressurized in flight to 4.5 bar using metal bellows pumps. The HIRES unit was first deployed in 2010 but has up to now not been used for regular halocarbon analysis with the exception of chloromethane analysis. The sample collection unit was tested for the sampling and storage effects of 28 halogenated compounds. The focus was on compound stability in the stainless-steel canisters during storage of up to 5 weeks and on the influence of ozone, since flights take place in the upper troposphere and lowermost stratosphere with ozone mixing ratios of up to several hundred parts per billion by volume (ppbv). Most of the investigated (hydro)chlorofluorocarbons and long-lived hydrofluorocarbons were found to be stable over a storage time of up to 5 weeks and were unaltered by ozone being present during pressurization. Some compounds such as dichloromethane, trichloromethane, and tetrachloroethene started to decrease in the canisters after a storage time of more than 2 weeks or exhibited lowered mixing ratios in samples pressurized with ozone present. A few compounds such as tetrachloromethane and tribromomethane were found to be unstable in the HIRES stainless-steel canisters independent of ozone levels. Furthermore, growth was observed during storage for some species, namely for HFC-152a, HFC-23, and Halon 1301.
A list of authors and their affiliations appears at the end of the paper New-particle formation is a major contributor to urban smog, but how it occurs in cities is often puzzling. If the growth rates of urban particles are similar to those found in cleaner environments (1–10 nanometres per hour), then existing understanding suggests that new urban particles should be rapidly scavenged by the high concentration of pre-existing particles. Here we show, through experiments performed under atmospheric conditions in the CLOUD chamber at CERN, that below about +5 degrees Celsius, nitric acid and ammonia vapours can condense onto freshly nucleated particles as small as a few nanometres in diameter. Moreover, when it is cold enough (below −15 degrees Celsius), nitric acid and ammonia can nucleate directly through an acid–base stabilization mechanism to form ammonium nitrate particles. Given that these vapours are often one thousand times more abundant than sulfuric acid, the resulting particle growth rates can be extremely high, reaching well above 100 nanometres per hour. However, these high growth rates require the gas-particle ammonium nitrate system to be out of equilibrium in order to sustain gas-phase supersaturations. In view of the strong temperature dependence that we measure for the gas-phase supersaturations, we expect such transient conditions to occur in inhomogeneous urban settings, especially in wintertime, driven by vertical mixing and by strong local sources such as traffic. Even though rapid growth from nitric acid and ammonia condensation may last for only a few minutes, it is nonetheless fast enough to shepherd freshly nucleated particles through the smallest size range where they are most vulnerable to scavenging loss, thus greatly increasing their survival probability. We also expect nitric acid and ammonia nucleation and rapid growth to be important in the relatively clean and cold upper free troposphere, where ammonia can be convected from the continental boundary layer and nitric acid is abundant from electrical storms.
Diamond formation in the Earth has been extensively discussed in recent years on the basis of geochemical analysis of natural materials, high-pressure experimental studies, or theoretical aspects. Here, we demonstrate experimentally for the first time, the spontaneous crystallization of diamond from CH4-rich fluids at pressure, temperature and redox conditions approximating those of the deeper parts of the cratonic lithospheric mantle (5-7 GPa) without using diamond seed crystals or carbides. In these experiments the fluid phase is nearly pure methane, even though the oxygen fugacity was significantly above metal saturation. We propose several previously unidentified mechanisms that may promote diamond formation under such conditions and which may also have implications for the origin of sublithospheric diamonds. These include the hydroxylation of silicate minerals like olivine and pyroxene, H2 incorporation into these phases and the "etching" of graphite by H2 and CH4 and reprecipitation as diamond. This study also serves as a demonstration of our new high-pressure experimental technique for obtaining reduced fluids, which is not only relevant for diamond synthesis, but also for investigating the metasomatic origins of diamond in the upper mantle, which has further implications for the deep carbon cycle.
The frequency of extreme events has changed, having a direct impact on human lives. Regional climate models help us to predict these regional climate changes. This work presents an atmosphere–ocean coupled regional climate system model (RCSM; with the atmospheric component COSMO-CLM and the ocean component NEMO) over the European domain, including three marginal seas: the Mediterranean, North, and Baltic Sea. To test the model, we evaluate a simulation of more than 100 years (1900–2009) with a spatial grid resolution of about 25 km. The simulation was nested into a coupled global simulation with the model MPI-ESM in a low-resolution configuration, whose ocean temperature and salinity were nudged to the ocean–ice component of the MPI-ESM forced with the NOAA 20th Century Reanalysis (20CR). The evaluation shows the robustness of the RCSM and discusses the added value by the coupled marginal seas over an atmosphere-only simulation. The coupled system is stable for the complete 20th century and provides a better representation of extreme temperatures compared to the atmosphere-only model. The produced long-term dataset will help us to better understand the processes leading to meteorological and climate extremes.
The most frequently used boundary-layer turbulence parameterization in numerical weather prediction (NWP) models are turbulence kinetic energy (TKE) based-based schemes. However, these parameterizations suffer from a potential weakness, namely the strong dependence on an ad-hoc quantity, the so-called turbulence length scale. The physical interpretation of the turbulence length scale is difficult and hence it cannot be directly related to measurements or large eddy simulation (LES) data. Consequently, formulations for the turbulence length scale in basically all TKE schemes are based on simplified assumptions and are model-dependent. A good reference for the independent evaluation of the turbulence length scale expression for NWP modeling is missing. Here we propose a new turbulence length scale diagnostic which can be used in the gray zone of turbulence without modifying the underlying TKE turbulence scheme. The new diagnostic is based on the TKE budget: The core idea is to encapsulate the sum of the molecular dissipation and the cross-scale TKE transfer into an effective dissipation, and associate it with the new turbulence length scale. This effective dissipation can then be calculated as a residuum in the TKE budget equation (for horizontal sub-domains of different sizes) using LES data. Estimation of the scale dependence of the diagnosed turbulence length scale using this novel method is presented for several idealized cases.
Metasomatic evolution of coesite-bearing diamondiferous eclogite from the Udachnaya Kimberlite
(2020)
A coesite-bearing diamondiferous eclogite from the Udachnaya kimberlite (Daldyn field, Siberian craton) has been studied to trace its complex evolution recorded in rock-forming and minor mineral constituents. The eclogite sample is composed of rock-forming omphacite (60 vol%), garnet (35 vol%) and quartz/coesite (5 vol%) and contains intergranular euhedral zoned olivine crystals, up to 200 µm long, coexisting with phlogopite, orthopyroxene, clinopyroxene (secondary), K-feldspar, plagioclase, spinel, sodalite and djerfisherite. Garnet grains are zoned, with a relatively homogeneous core and a more magnesian overgrowth rim. The rim zones further differ from the core in having higher Zr/Y (6 times that in the cores), ascribed to interaction with, or precipitation from, a kimberlite-related melt. Judging by pressure-temperature estimates (~1200 °C; 6.2 GPa), the xenolith originated at depths of ~180–200 km at the base of the continental lithosphere. The spatial coexistence of olivine, orthopyroxene and coesite/quartz with K-Na-Cl minerals in the xenolith indicates that eclogite reacted with a deep-seated kimberlite melt. However, Fe-rich olivine, orthopyroxene and low-pressure minerals (sodalite and djerfisherite) likely result from metasomatic reaction at shallower depths during transport of the eclogite by the erupting kimberlite melt. Our results demonstrate that a mixed eclogitic-peridotitic paragenesis, reported previously from inclusions in diamond, can form by interaction of eclogite and a kimberlite-related melt.
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
Accretion rates of Holocene tropical coral reefs in three areas in the Atlantic, Pacific and Indian Oceans have been quantified in 79 dated core sections in 34 reef cores from Belize, the Maldives and French Polynesia. Holocene vertical reef accretion rate averages 5.05 m/kyr and has decreased during the past 10 kyr. Accretion rates in branched and massive coral facies are statistically similar. Reef accretion rate is positively correlated with the rate of sea‐level rise, that is the degree of creation of accommodation space, and with climate as expressed in a Holocene sea surface temperature anomaly. Accommodation space is also created by subsidence, but at a rate one to two orders of magnitude lower than that created by glacio‐eustasy (0.04 to 0.16 m/kyr). Lagoonal background sedimentation in adjacent reef lagoons averages 0.89 m/kyr as measured in 72 dated core sections in 28 cores. Lagoonal carbonate sedimentation on top of underlying mangrove peat usually starts after a considerable hiatus of ca 3 kyr on average. The lagoonal background sedimentation rate increased during the Holocene, probably due to deepening. The differences between vertical reef accretion and lagoonal background sedimentation rates are a major factor in the production of the widely known saucer shapes typical of tropical reefs and carbonate platforms, that is the creation of unfilled accommodation space. Reef core recovery, used as a proxy for reef consolidation, and core depth exhibit a statistically negative correlation based on data from 326 core barrels. Recovery and marine cement abundance (average volume 8.6%) also decrease from windward to leeward core positions. These observations are presumably a result of both a decrease in the rate of sea‐level rise that is the increase in time available for submarine cementation during the Holocene and the amount of flushing of reef interstices by marine waters.
Diurnal valley winds frequently form over complex topography, particularly under fair weather conditions, and have a significant impact on the local weather and climate. Since diurnal valley winds result from complex and multi-scale interactions, their representation in numerical weather prediction models is challenging. Better understanding of these local winds based on observations is crucial to improve the accuracy of the forecasts. This study investigates the diurnal evolution of the three-dimensional mean wind structure in a deep Alpine valley, the Rhone valley at Sion, using data from a radar wind profiler and a surface weather station operated continuously from 1 September 2016 to 17 July 2017. In particular, the wind profiler data was analyzed for a subset of days on which fair weather conditions allowed for the full development of thermally driven winds. A pronounced diurnal cycle of the wind speed, as well as a reversal of the wind direction twice per day is documented for altitudes up to about 2 km above ground level (AGL) in the warm season and less than 1 km AGL in winter. The diurnal pattern undergoes significant changes during the course of the year. Particularly during the warm-weather months of May through to September, a low-level wind maximum occurs, where mean maximum up-valley velocities of 8–10 m s−1 are found between 15–16 UTC at altitudes around 200 m AGL. In addition, during nighttime, a down-valley jet with maximum wind speeds of 4–8 m s−1 around 1 km AGL is found. A case study of a three-day period in September 2016 illustrates the occurrence of an elevated layer of cross-valley flow around 1–1.5 km AGL.
Under the Kigali Amendment to the Montreal Protocol, new controls are being implemented to reduce emissions of HFC-23 (CHF3), a by-product during the manufacture of HCFC-22 (CHClF2). Starting in 2015, China and India, who dominate global HCFC-22 production (75% in 2017), set out ambitious programs to reduce HFC-23 emissions. Here, we estimate that these measures should have seen global emissions drop by 87% between 2014 and 2017. Instead, atmospheric observations show that emissions have increased and in 2018 were higher than at any point in history (15.9 ± 0.9 Gg yr−1). Given the magnitude of the discrepancy between expected and observation-inferred emissions, it is likely that the reported reductions have not fully materialized or there may be substantial unreported production of HCFC-22, resulting in unaccounted-for HFC-23 by-product emissions. The difference between reported and observation-inferred estimates suggests that an additional ~309 Tg CO2-equivalent emissions were added to the atmosphere between 2015 and 2017.
The endemic argan tree (Argania spinosa) populations in South Morocco are highly degraded due to their use as a biomass resource in dry years and illegal firewood extraction. The intensification and expansion of agricultural land lead to a retreat of the wooded area, while the remaining argan open woodlands are often overgrazed. Thus, canopy-covered areas decrease while areas without vegetation cover between the argan trees increase. In total, 36 rainfall simulation experiments as well as 60 infiltration measurements were conducted to investigate the potential difference between tree-covered areas and bare intertree areas. In addition, 60 soil samples were taken under the trees and in the intertree areas parallel to the contour lines. Significant differences using a t-test were found between tree and intertree areas for the studied parameters Ksat, Kh, pH, electric conductivity, percolation stability, total C-content, total N-content, K-content, Na-content, and Mg-content. Surface runoff and soil losses were not as conclusive but showed similar trends. The results showed that argan trees influence the soil underneath significantly, while the soil in intertree areas is less protected and more degraded. It is therefore reasonable to assume further degradation of the soil when intertree areas extend further due to lack of rejuvenation of argan trees.
Invasive plant species are increasingly altering species composition and the functioning of ecosystems from a local to a global scale. The grass species Pennisetum setaceum has recently raised concerns as an invader on different archipelagos worldwide. Among these affected archipelagos are the Canary Islands, which are a hotspot of endemism. Consequently, conservation managers and stakeholders are interested in the potential spreading of this species in the archipelago. We identify the current extent of the suitable habitat for P. setaceum on the island of La Palma to assess how it affects island ecosystems, protected areas (PAs), and endemic plant species richness. We recorded in situ occurrences of P. setaceum from 2010 to 2018 and compiled additional ones from databases at a 500 m × 500 m resolution. To assess the current suitable habitat and possible distribution patterns of P. setaceum on the island, we built an ensemble model. We projected habitat suitability for island ecosystems and PAs and identified risks for total as well as endemic plant species richness. The suitable habitat for P. setaceum is calculated to cover 34.7% of the surface of La Palma. In open ecosystems at low to mid elevations, where native ecosystems are already under pressure by land use and human activities, the spread of the invader will likely lead to additional threats to endemic plant species. Forest ecosystems (e.g., broadleaved evergreen and coniferous forests) are not likely to be affected by the spread of P. setaceum because of its heliophilous nature. Our projection of suitable habitat of P. setaceum within ecosystems and PAs on La Palma supports conservationists and policymakers in prioritizing management and control measures and acts as an example for the potential threat of this graminoid invader on other islands.
Often in climate system studies, linear and symmetric statistical measures are applied to quantify interactions among subsystems or variables. However, they do not allow identification of the driving and responding subsystems. Therefore, in this study, we aimed to apply asymmetric measures from information theory: the axiomatically proposed transfer entropy and the first principle-based information flow to detect and quantify climate interactions. As their estimations are challenging, we initially tested nonparametric estimators like transfer entropy (TE)-binning, TE-kernel, and TE k-nearest neighbor and parametric estimators like TE-linear and information flow (IF)-linear with idealized two-dimensional test cases along with their sensitivity on sample size. Thereafter, we experimentally applied these methods to the Lorenz-96 model and to two real climate phenomena, i.e., (1) the Indo-Pacific Ocean coupling and (2) North Atlantic Oscillation (NAO)–European air temperature coupling. As expected, the linear estimators work for linear systems but fail for strongly nonlinear systems. The TE-kernel and TE k-nearest neighbor estimators are reliable for linear and nonlinear systems. Nevertheless, the nonparametric methods are sensitive to parameter selection and sample size. Thus, this work proposes a composite use of the TE-kernel and TE k-nearest neighbor estimators along with parameter testing for consistent results. The revealed information exchange in Lorenz-96 is dominated by the slow subsystem component. For real climate phenomena, expected bidirectional information exchange between the Indian and Pacific SSTs was detected. Furthermore, expected information exchange from NAO to European air temperature was detected, but also unexpected reversal information exchange. The latter might hint to a hidden process driving both the NAO and European temperatures. Hence, the limitations, availability of time series length and the system at hand must be taken into account before drawing any conclusions from TE and IF-linear estimations.
Convective shower characteristics simulated with the convection-permitting climate model COSMO-CLM
(2019)
This paper evaluates convective precipitation as simulated by the convection-permitting climate model (CPM) Consortium for Small-Scale Modeling in climate mode (COSMO-CLM) (with 2.8 km grid-spacing) over Germany in the period 2001–2015. Characteristics of simulated convective precipitation objects like lifetime, area, mean intensity, and total precipitation are compared to characteristics observed by weather radar. For this purpose, a tracking algorithm was applied to simulated and observed precipitation with 5-min temporal resolution. The total amount of convective precipitation is well simulated, with a small overestimation of 2%. However, the simulation underestimates convective activity, represented by the number of convective objects, by 33%. This underestimation is especially pronounced in the lowlands of Northern Germany, whereas the simulation matches observations well in the mountainous areas of Southern Germany. The underestimation of activity is compensated by an overestimation of the simulated lifetime of convective objects. The observed mean intensity, maximum intensity, and area of precipitation objects increase with their lifetime showing the spectrum of convective storms ranging from short-living single-cell storms to long-living organized convection like supercells or squall lines. The CPM is capable of reproducing the lifetime dependence of these characteristics but shows a weaker increase in mean intensity with lifetime resulting in an especially pronounced underestimation (up to 25%) of mean precipitation intensity of long-living, extreme events. This limitation of the CPM is not identifiable by classical evaluation techniques using rain gauges. The simulation can reproduce the general increase of the highest percentiles of cell area, total precipitation, and mean intensity with temperature but fails to reproduce the increase of lifetime. The scaling rates of mean intensity and total precipitation resemble observed rates only in parts of the temperature range. The results suggest that the evaluation of coarse-grained (e.g., hourly) precipitation fields is insufficient for revealing challenges in convection-permitting simulations.
The sources and critical enrichment processes for granite related tin ores are still not well understood. The Erzgebirge represents one of the classical regions for tin mineralization. We investigated the four largest plutons from the Western Erzgebirge (Germany) for the geochemistry of bulk rocks and autocrystic zircons and relate this information to their intrusion ages. The source rocks of the Variscan granites were identified as high-grade metamorphic rocks based on the comparison of Hf-O isotope data on zircons, the abundance of xenocrystic zircon ages as well as Nd and Hf model ages. Among these rocks, restite is the most likely candidate for later Variscan melts. Based on the evolution with time, we could reconstruct enrichment factors for tin and tungsten starting from the protoliths (575 Ma) that were later converted to high-grade metamorphic rocks (340 Ma) and served as sources for the older biotite granites (323–318 Ma) and the tin granites (315–314 Ma). This evolution involved a continuous enrichment of both tin and tungsten with an enrichment factor of ~15 for tin and ~7 for tungsten compared to the upper continental crust (UCC). Ore level concentrations (>10–100 times enrichment) were achieved only in the greisen bodies and dykes by subsequent hydrothermal processes.
Video and image data are regularly used in the field of benthic ecology to document biodiversity. However, their use is subject to a number of challenges, principally the identification of taxa within the images without associated physical specimens. The challenge of applying traditional taxonomic keys to the identification of fauna from images has led to the development of personal, group, or institution level reference image catalogues of operational taxonomic units (OTUs) or morphospecies. Lack of standardisation among these reference catalogues has led to problems with observer bias and the inability to combine datasets across studies. In addition, lack of a common reference standard is stifling efforts in the application of artificial intelligence to taxon identification. Using the North Atlantic deep sea as a case study, we propose a database structure to facilitate standardisation of morphospecies image catalogues between research groups and support future use in multiple front-end applications. We also propose a framework for coordination of international efforts to develop reference guides for the identification of marine species from images. The proposed structure maps to the Darwin Core standard to allow integration with existing databases. We suggest a management framework where high-level taxonomic groups are curated by a regional team, consisting of both end users and taxonomic experts. We identify a mechanism by which overall quality of data within a common reference guide could be raised over the next decade. Finally, we discuss the role of a common reference standard in advancing marine ecology and supporting sustainable use of this ecosystem.
Here we present a comprehensive attempt to correlate aragonitic Na∕Ca ratios from Desmophyllum pertusum (formerly known as Lophelia pertusa), Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na∕Ca ratios using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements.
The results reveal no apparent correlation with salinity (30.1–40.57 g kg−1) but a significant inverse correlation with temperature (−0.31±0.04 mmolmol−1∘C−1). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na∕Mg ratios show a similar temperature sensitivity to Na∕Ca ratios, but the combination of two ratios appears to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na∕Ca and Na∕Mg, however, limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na∕Ca ratios for an open and semi-enclosed calcifying space based on temperature-sensitive Na- and Ca-pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na∕Ca ratio.
The Arctic Svalbard Archipelago hosts the world’s northernmost cold-water ‘carbonate factories’ thriving here despite of presumably unfavourable environmental conditions and extreme seasonality. Two contrasting sites of intense biogenic carbonate production, the rhodolith beds in Mosselbukta in the north of the archipelago and the barnacle-mollusc dominated carbonate sediments accumulating in the strong hydrodynamic regime of the Bjørnøy-Banken south of Spitsbergen, were the targets of the RV Maria S. Merian cruise 55 in June 2016. By integrating data from physical oceanography, marine biology, and marine geology, the present contribution characterises the environmental setting and biosedimentary dynamics of these two polar carbonate factories. Repetitive CTD profiling in concert with autonomous temperature/salinity loggers on a long-term settlement platform identified spatiotemporal patterns in the involved Atlantic and Polar water masses, whereas short-term deployments of a lander revealed fluctuations of environmental variables in the rhodolith beds in Mosselbukta and at same depth (46 m) at Bjørnøy-Banken. At both sites, dissolved inorganic nutrients in the water column were found depleted (except for elevated ammonium concentrations) and show an overall increase in concentration and N:P ratios toward deeper waters. This indicates that a recycling system was fuelling primary production after the phytoplankton spring bloom at the time of sampling in June 2016. Accordingly, oxygen levels were found elevated and carbon dioxide concentrations (pCO2) markedly reduced, on average only half the expected equilibrium values. Backed up by seawater stable carbon and oxygen isotope signatures, this is interpreted as an effect of limited air-sea gas exchange during seasonal ice cover in combination with a boost in community photosynthesis during the spring phytoplankton bloom. The observed trends are enhanced by the onset of rhodophyte photosynthesis in the rhodolith beds during the polar day upon retreat of sea-ice. Potential adverse effects of ocean acidification on the local calcifier community are thus predicted to be seasonally buffered by the marked drop in pCO2 during the phase of sea-ice cover and spring phyto-plankton bloom, but this effect will diminish should the seasonal sea-ice formation continue to decline. Among the 25 macrobenthos taxa identified from images captured by the lander’s camera system, all but three species were calcifiers contributing to the carbonate production. Biodiversity was found to be much higher in Mosselbukta (21 taxa) compared to Bjørnøy-Banken (8 taxa), which is considered as a result of enhanced habitat diversity provided in the rhodolith beds by the bioengineering crustose alga Lithothamnion glaciale. Filter-feeding activity of selected key species did reveal group-specific but no common activity patterns. Biotic disturbance of the filtering activity was common, in contrast to abiotic factors, with hermit crabs representing the primary trigger. Motion tracking of rhodoliths revealed a high frequency of dislocation, triggered not by abiotic factors but by the activity of benthic invertebrates, in particular echinoids ploughing below or moving over the rhodoliths. The echinoid Strongylocentrotus sp. is the most abundant component of the associated fauna, thereby considerably contributing both to carbonate production and to grazing bioerosion. Together, these results portray a high degree of seasonal as well as short-term dynamics in environmental conditions that despite many similarities support distinctly different communities and biodiversity patterns in the calcifying macrobenthos at the two studied polar carbonate factories.
The exchange of heat, momentum, and mass in the atmosphere over mountainous terrain is controlled by synoptic-scale dynamics, thermally driven mesoscale circulations, and turbulence. This article reviews the key challenges relevant to the understanding of exchange processes in the mountain boundary layer and outlines possible research priorities for the future. The review describes the limitations of the experimental study of turbulent exchange over complex terrain, the impact of slope and valley breezes on the structure of the convective boundary layer, and the role of intermittent mixing and wave–turbulence interaction in the stable boundary layer. The interplay between exchange processes at different spatial scales is discussed in depth, emphasizing the role of elevated and ground-based stable layers in controlling multi-scale interactions in the atmosphere over and near mountains. Implications of the current understanding of exchange processes over mountains towards the improvement of numerical weather prediction and climate models are discussed, considering in particular the representation of surface boundary conditions, the parameterization of sub-grid-scale exchange, and the development of stochastic perturbation schemes.
Throughout mankind’s history, the need to secure and protect the home settlement was an essential one. This holds especially true for the city of Ainos (modern Enez) in Turkish Thrace. Due to its continuous settlement history since the 7th/6th century BC, several different types of city walls were built—sometimes even on top of each other—several of which have been preserved over time. To decipher the construction style, the course and the age of a buried city wall segment in the southern part of the former city, a geoscientific multi-proxy approach including magnetic gradiometer (MG) and electrical resistivity tomography (ERT) measurements in combination with granulometrical, sedimentological and microfaunistical investigations on sediment cores was applied. We were able to (1) present reasonable arguments for its Hellenistic age; (2) reveal the course of this wall segment and extrapolate it further north into a less studied area; and (3) demonstrate that in this near-coastal area, the former swampy terrain had been consolidated for constructing the wall. Our multi-proxy approach serves as a valuable example for investigating buried structures in archaeological contexts, avoiding a less-economical, time-consuming, or even forbidden excavation.
Atmospheric aerosols and their effect on clouds are thought to be important for anthropogenic radiative forcing of the climate, yet remain poorly understood1. Globally, around half of cloud condensation nuclei originate from nucleation of atmospheric vapours2. It is thought that sulfuric acid is essential to initiate most particle formation in the atmosphere3,4, and that ions have a relatively minor role5. Some laboratory studies, however, have reported organic particle formation without the intentional addition of sulfuric acid, although contamination could not be excluded6,7. Here we present evidence for the formation of aerosol particles from highly oxidized biogenic vapours in the absence of sulfuric acid in a large chamber under atmospheric conditions. The highly oxygenated molecules (HOMs) are produced by ozonolysis of α-pinene. We find that ions from Galactic cosmic rays increase the nucleation rate by one to two orders of magnitude compared with neutral nucleation. Our experimental findings are supported by quantum chemical calculations of the cluster binding energies of representative HOMs. Ion-induced nucleation of pure organic particles constitutes a potentially widespread source of aerosol particles in terrestrial environments with low sulfuric acid pollution.
We present a study characterizing aerosol particles resulting from a skyscraper blasting. High mass concentrations with a maximum of 844.9 μg m-3 were present for a short time period of approximately 15 minutes. They result in a day mean of 32.6 μg m-3 compared to a 27.6 μg m-3 background not exceeding the 50 μg m-3 EU maximum permissive value. The increase in particle number concentration was less pronounced with a maximum concentration of 6.9 ⋅ 104 cm-3 compared to the local background value of 1.8 ⋅ 104 cm-3. The size-resolved number concentration shows a single mode of ultrafine particles at approximately 93 nm. The spatial distribution of deposited dust was investigated with Bergerhoff glass collection vessels, showing a decrease with distance. In the deposited dust samples the concentrations of twelve metals was determined, non of them exceeded the regional background concentrations significantly. The chemical composition of individual particles emitted by the demolition was studied by Scanning Electron Microscopy. They were mainly concrete and steel particles, with 60% calcium carbonates, 19% calcium sulfates, 19% silicates and 2% steel. In energy-dispersive X-Ray Spectroscopy, no fibers like asbestos were observed. Using a broad spectrum of instruments and methods, we obtain comprehensive characterization of the particles emitted by the demolition.
A twentieth century-long coupled atmosphere-ocean regional climate simulation with COSMO-CLM (Consortium for Small-Scale Modeling, Climate Limited-area Model) and NEMO (Nucleus for European Modelling of the Ocean) is studied here to evaluate the added value of coupled marginal seas over continental regions. The interactive coupling of the marginal seas, namely the Mediterranean, the North and the Baltic Seas, to the atmosphere in the European region gives a comprehensive modelling system. It is expected to be able to describe the climatological features of this geographically complex area even more precisely than an atmosphere-only climate model. The investigated variables are precipitation and 2 m temperature. Sensitivity studies are used to assess the impact of SST (sea surface temperature) changes over land areas. The different SST values affect the continental precipitation more than the 2 m temperature. The simulated variables are compared to the CRU (Climatic Research Unit) observational data, and also to the HOAPS/GPCC (Hamburg Ocean Atmosphere Parameters and Fluxes from Satellite Data, Global Precipitation Climatology Centre) data. In the coupled simulation, added skill is found primarily during winter over the eastern part of Europe. Our analysis shows that, over this region, the coupled system is dryer than the uncoupled system, both in terms of precipitation and soil moisture, which means a decrease in the bias of the system. Thus, the coupling improves the simulation of precipitation over the eastern part of Europe, due to cooler SST values and in consequence, drier soil.
Historically, the expansion of soy plantations has been a major driver of land-use/cover change (LUCC) in Brazil. While a series of recent public actions and supply-chain commitments reportedly curbed the replacement of forests by soy, the expansion of the agricultural commodity still poses a considerable threat to the Amazonian and Cerrado biomes. Identification of areas under high risk of soy expansion is thus paramount to assist conservation efforts in the region. We mapped the areas suitable for undergoing transition to soy plantations in the Legal Amazon with a machine-learning approach adopted from the ecological modeling literature. Simulated soy expansion for the year 2014 exhibited favorable validation scores compared to other LUCC models. We then used our model to simulate how potential future infrastructure improvements would affect the 2014 probabilities of soy occurrence in the region. In addition to the 2.3 Mha of planted soy in the Legal Amazon in 2014, our model identified another 14.7 Mha with high probability of soy conversion in the region given the infrastructure conditions at that time. Out of those, pastures and forests represented 9.8 and 0.4 Mha, respectively. Under the new infrastructure scenarios simulated, the Legal Amazonian area under high risk of soy conversion increased by up to 2.1 Mha (14.6%). These changes led to up to 11.4 and 51.4% increases in the high-risk of conversion areas of pastures and forests, respectively. If conversion occurs in the identified high-risk areas, at least 4.8 Pg of CO2 could be released into the atmosphere, a value that represents 10 times the total CO2 emissions of Brazil in 2014. Our results highlight the importance of targeting conservation policies and enforcement actions, including the Soy Moratorium, to mitigate future forest cover loss associated with infrastructure improvements in the region.
We evaluate the near-surface representation of thermally driven winds in the Swiss Alps in a numerical weather prediction model at km-scale resolution. In addition, the influence of grid resolution (2.2 km and 1.1 km), topography filtering, and land surface datasets on the accuracy of the simulated valley winds is investigated. The simulations are evaluated against a comprehensive set of surface observations for an 18-day fair-weather summer period in July 2006. The episode is characterized by strong diurnal wind systems and the formation of shallow convection over the mountains, which transitions to precipitating convection in some areas. The near-surface winds (10 m above ground level) follow a typical diurnal pattern with strong daytime up-valley flow and weaker nighttime down-valley flow. At a 2.2 km resolution the valley winds are poorly simulated for most stations, while at a 1.1 km resolution the diurnal cycle of the valley winds is well represented in most large (e.g., Rhein valley at Chur and Rhone valley at Visp) and medium-sized valleys (e.g., Linth valley at Glarus). In the smaller valleys (e.g., Maggia valley at Cevio), the amplitude of the valley wind is still significantly underestimated, even at a 1.1 km resolution. Detailed sensitivity experiments show that the use of high-resolution land surface datasets, for both the soil characteristics as well as for the land cover, and reduced filtering of the topography are essential to achieve good performance at a 1.1 km resolution
The synthesis of polynitrogen compounds is of fundamental importance due to their potential as environmentally-friendly high energy density materials. Attesting to the intrinsic difficulties related to their formation, only three polynitrogen ions, bulk stabilized as salts, are known. Here, magnesium and molecular nitrogen are compressed to about 50 GPa and laser-heated, producing two chemically simple salts of polynitrogen anions, MgN4 and Mg2N4. Single-crystal X-ray diffraction reveals infinite anionic polythiazyl-like 1D N-N chains in the crystal structure of MgN4 and cis-tetranitrogen N44− units in the two isosymmetric polymorphs of Mg2N4. The cis-tetranitrogen units are found to be recoverable at atmospheric pressure. Our results respond to the quest for polynitrogen entities stable at ambient conditions, reveal the potential of employing high pressures in their synthesis and enrich the nitrogen chemistry through the discovery of other nitrogen species, which provides further possibilities to design improved polynitrogen arrangements.
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.
To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of Abstract. To quantify water flows between groundwater (GW) and surface water (SW) as well as the impact of capillary rise on evapotranspiration by global hydrological models (GHMs), it is necessary to replace the bucket-like linear GW reservoir model typical for hydrological models with a fully integrated gradient-based GW flow model. Linear reservoir models can only simulate GW discharge to SW bodies, provide no information on the location of the GW table and assume that there is no GW flow among grid cells. A gradient-based GW model simulates not only GW storage but also hydraulic head, which together with information on SW table elevation enables the quantification of water flows from GW to SW and vice versa. In addition, hydraulic heads are the basis for calculating lateral GW flow among grid cells and capillary rise.
G³M is a new global gradient-based GW model with a spatial resolution of 5' that will replace the current linear GW reservoir in the 0.5° WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables inmemory coupling to WGHM while keeping overall runtime relatively low, allowing sensitivity analyses and data assimilation. This paper presents the G³M concept and specific model design decisions together with results under steady-state naturalized conditions, i.e. neglecting GW abstractions. Cell-specific conductances of river beds, which govern GW-SW interaction, were determined based on the 30'' steady-state water table computed by Fan et al. (2013). Together with an appropriate choice for the effective elevation of the SW table within each grid cell, this enables a reasonable simulation of drainage from GW to SW such that, in contrast to the GW model of de Graaf et al. (2015, 2017), no additional drainage based on externally provided values for GW storage above the floodplain is required in G³M. Comparison of simulated hydraulic heads to observations around the world shows better agreement than de Graaf et al. (2015). In addition, G³M output is compared to the output of two established macro-scale models for the Central Valley, California, and the continental United States, respectively. As expected, depth to GW table is highest in mountainous and lowest in flat regions. A first analysis of losing and gaining rivers and lakes/wetlands indicates that GW discharge to rivers is by far the dominant flow, draining diffuse GW recharge, such that lateral flows only become a large fraction of total diffuse and focused recharge in case of losing rivers and some areas with very low GW recharge. G³M does not represent losing rivers in some dry regions. This study presents the first steps towards replacing the linear GW reservoir model in a GHM while improving on recent efforts, demonstrating the feasibility of the approach and the robustness of the newly developed framework.
In global hydrological models, groundwater (GW) is typically represented by a bucket-like linear groundwater reservoir. Reservoir models, however, (1) can only simulate GW discharge to surface water (SW) bodies but not recharge from SW to GW, (2) provide no information on the location of the GW table, and (3) assume that there is no GW flow among grid cells. This may lead, for example, to an underestimation of groundwater resources in semiarid areas where GW is often replenished by SW or to an underestimation of evapotranspiration where the GW table is close to the land surface. To overcome these limitations, it is necessary to replace the reservoir model in global hydrological models with a hydraulic head gradient-based GW flow model.
We present G3M, a new global gradient-based GW model with a spatial resolution of 5′ (arcminutes), which is to be integrated into the 0.5∘ WaterGAP Global Hydrology Model (WGHM). The newly developed model framework enables in-memory coupling to WGHM while keeping overall runtime relatively low, which allows sensitivity analyses, calibration, and data assimilation. This paper presents the G3M concept and model design decisions that are specific to the large grid size required for a global-scale model. Model results under steady-state naturalized conditions, i.e., neglecting GW abstractions, are shown. Simulated hydraulic heads show better agreement to observations around the world compared to the model output of de Graaf et al. (2015). Locations of simulated SW recharge to GW are found, as is expected, in dry and mountainous regions but areal extent of SW recharge may be underestimated. Globally, GW discharge to rivers is by far the dominant flow component such that lateral GW flows only become a large fraction of total diffuse and focused recharge in the case of losing rivers, some mountainous areas, and some areas with very low GW recharge. A strong sensitivity of simulated hydraulic heads to the spatial resolution of the model and the related choice of the water table elevation of surface water bodies was found. We suggest to investigate how global-scale groundwater modeling at 5′ spatial resolution can benefit from more highly resolved land surface elevation data.
Latitudinal and bathymetrical species richness patterns in the NW Pacific and adjacent Arctic Ocean
(2019)
Global scale analyses have recently revealed that the latitudinal gradient in marine species richness is bimodal, peaking at low-mid latitudes but with a dip at the equator; and that marine species richness decreases with depth in many taxa. However, these overall and independently studied patterns may conceal regional differences that help support or qualify the causes in these gradients. Here, we analysed both latitudinal and depth gradients of species richness in the NW Pacific and its adjacent Arctic Ocean. We analysed 324,916 distribution records of 17,414 species from 0 to 10,900 m depth, latitude 0 to 90°N, and longitude 100 to 180°N. Species richness per c. 50 000 km2 hexagonal cells was calculated as alpha (local average), gamma (regional total) and ES50 (estimated species for 50 records) per latitudinal band and depth interval. We found that average ES50 and gamma species richness decreased per 5° latitudinal bands and 100 m depth intervals. However, average ES50 per hexagon showed that the highest species richness peaked around depth 2,000 m where the highest total number of species recorded. Most (83%) species occurred in shallow depths (0 to 500 m). The area around Bohol Island in the Philippines had the highest alpha species richness (more than 8,000 species per 50,000 km2). Both alpha and gamma diversity trends increased from the equator to latitude 10°N, then further decreased, but reached another peak at higher latitudes. The latitudes 60–70°N had the lowest gamma and alpha diversity where there is almost no ocean in our study area. Model selection on Generalized Additive Models (GAMs) showed that the combined effects of all environmental predictors produced the best model driving species richness in both shallow and deep sea. The results thus support recent hypotheses that biodiversity, while highest in the tropics and coastal depths, is decreasing at the equator and decreases with depth below ~2000 m. While we do find the declines of species richness with latitude and depth that reflect temperature gradients, local scale richness proved poorly correlated with many environmental variables. This demonstrates that while regional scale patterns in species richness may be related to temperature, that local scale richness depends on a greater variety of variables.
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of mean age of air values from long-lived trace gases. However, this provides more insight into general transport strength and less into its mechanism. Deriving complete transit time distributions (age spectra) is desirable, but their deduction from direct measurements is difficult and so far primarily achieved by assumptions about dynamics and spectra themselves. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases. For a full description of transport seasonality the formulation includes an imposed seasonal cycle to gain multimodal spectra. The EMAC model simulation used for a proof of concept features an idealized dataset of 40 radioactive trace gases with different chemical lifetimes as well as 40 chemically inert pulsed trace gases to calculate pulse age spectra. Annual and seasonal mean inverse spectra are compared to pulse spectra including first and second moments as well as the ratio between them to assess the performance on these time scales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years mean age of air. The imposed seasonal cycle emerges as a reliable tool to include transport seasonality in the age spectra. Below 1.5 years mean age of air, tropospheric influence intensifies and breaks the assumption of single entry through the tropical tropopause, leading to inaccurate spectra in particular in the northern hemisphere. The imposed seasonal cycle wrongly prescribes seasonal entry in this lower region and does not lead to a better agreement between inverse and pulse age spectra without further improvement. As the inverse method aims for future implementation on in situ observational data, possible critical factors for this purpose are delineated finally.
Deriving stratospheric age of air spectra using an idealized set of chemically active trace gases
(2019)
Analysis of stratospheric transport from an observational point of view is frequently realized by evaluation of the mean age of air values from long-lived trace gases. However, this provides more insight into general transport strength and less into its mechanism. Deriving complete transit time distributions (age spectra) is desirable, but their deduction from direct measurements is difficult. It is so far primarily based on model work. This paper introduces a modified version of an inverse method to infer age spectra from mixing ratios of short-lived trace gases and investigates its basic principle in an idealized model simulation. For a full description of transport seasonality the method includes an imposed seasonal cycle to gain multimodal spectra. An ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulation is utilized for a general proof of concept of the method and features an idealized dataset of 40 radioactive trace gases with different chemical lifetimes as well as 40 chemically inert pulsed trace gases to calculate pulse age spectra. It is assessed whether the modified inverse method in combination with the seasonal cycle can provide matching age spectra when chemistry is well-known. Annual and seasonal mean inverse spectra are compared to pulse spectra including first and second moments as well as the ratio between them to assess the performance on these timescales. Results indicate that the modified inverse age spectra match the annual and seasonal pulse age spectra well on global scale beyond 1.5 years of mean age of air. The imposed seasonal cycle emerges as a reliable tool to include transport seasonality in the age spectra. Below 1.5 years of mean age of air, tropospheric influence intensifies and breaks the assumption of single entry through the tropical tropopause, leading to inaccurate spectra, in particular in the Northern Hemisphere. The imposed seasonal cycle wrongly prescribes seasonal entry in this lower region and does not lead to a better agreement between inverse and pulse age spectra without further improvement. Tests with a focus on future application to observational data imply that subsets of trace gases with 5 to 10 species are sufficient for deriving well-matching age spectra. These subsets can also compensate for an average uncertainty of up to ±20 % in the knowledge of chemical lifetime if a deviation of circa ±10 % in modal age and amplitude of the resulting spectra is tolerated.
Understanding new particle formation and growth is important because of the strong impact of these processes on climate and air quality. Measurements to elucidate the main new particle formation mechanisms are essential; however, these mechanisms have to be implemented in models to estimate their impact on the regional and global scale. Parameterizations are computationally cheap ways of implementing nucleation schemes in models but they have their limitations, as they do not necessarily include all relevant parameters. Process models using sophisticated nucleation schemes can be useful for the generation of look-up tables in large scale models or for the analysis of individual new particle formation events. In addition, some other important properties can be derived from a process model that implicitly calculates the evolution of the full aerosol size distribution, e.g., the particle growth rates. Within this study, a model (SANTIAGO, Sulfuric acid Ammonia NucleaTIon And GrOwth model) is constructed that simulates new particle formation starting from the monomer of sulfuric acid up to a particle size of several hundred nanometers. The smallest sulfuric acid clusters containing one to four acid molecules and varying amount of base (ammonia) are allowed to evaporate in the model, whereas growth beyond the pentamer (5 sulfuric acid molecules) is assumed to be entirely collision-controlled. The main goal of the present study is to derive appropriate thermodynamic data needed to calculate the cluster evaporation rates as a function of temperature. These data are derived numerically from CLOUD (Cosmics Leaving OUtdoor Droplets) chamber new particle formation rates for neutral sulfuric acid-water-ammonia nucleation at temperatures between 208 K and 292 K. The numeric methods include an optimization scheme to derive the best estimates for the thermodynamic data (dH and dS) and a Monte Carlo method to derive their probability density functions. The derived data are compared to literature values. Using different data sets for dH and dS in SANTIAGO detailed comparison between model results and measured CLOUD new particle formation rates is discussed.
Understanding new particle formation and growth is important because of the strong impact of these processes on climate and air quality. Measurements to elucidate the main new particle formation mechanisms are essential; however, these mechanisms have to be implemented in models to estimate their impact on the regional and global scale. Parameterizations are computationally cheap ways of implementing nucleation schemes in models, but they have their limitations, as they do not necessarily include all relevant parameters. Process models using sophisticated nucleation schemes can be useful for the generation of look-up tables in large-scale models or for the analysis of individual new particle formation events. In addition, some other important properties can be derived from a process model that implicitly calculates the evolution of the full aerosol size distribution, e.g., the particle growth rates. Within this study, a model (SANTIAGO – Sulfuric acid Ammonia NucleaTIon And GrOwth model) is constructed that simulates new particle formation starting from the monomer of sulfuric acid up to a particle size of several hundred nanometers. The smallest sulfuric acid clusters containing one to four acid molecules and a varying amount of base (ammonia) are allowed to evaporate in the model, whereas growth beyond the pentamer (five sulfuric acid molecules) is assumed to be entirely collision-controlled. The main goal of the present study is to derive appropriate thermodynamic data needed to calculate the cluster evaporation rates as a function of temperature. These data are derived numerically from CLOUD (Cosmics Leaving OUtdoor Droplets) chamber new particle formation rates for neutral sulfuric acid–water–ammonia nucleation at temperatures between 208 and 292 K. The numeric methods include an optimization scheme to derive the best estimates for the thermodynamic data (dH and dS) and a Monte Carlo method to derive their probability density functions. The derived data are compared to literature values. Using different data sets for dH and dS in SANTIAGO detailed comparison between model results and measured CLOUD new particle formation rates is discussed.
Within the world’s oceans, regionally distinct ecological niches develop due to differences in water temperature, nutrients, food availability, predation and light intensity. This results in differences in the vertical dispersion of planktonic foraminifera on the global scale. Understanding the controls on these modern-day distributions is important when using these organisms for paleoceanographic reconstructions. As such, this study constrains modern depth habitats for the northern equatorial Indian Ocean, for 14 planktonic foraminiferal species (G. ruber, G. elongatus, G. pyramidalis, G. rubescens, T. sacculifer, G. siphonifera, G. glutinata, N. dutertrei, G. bulloides, G. ungulata, P. obliquiloculata, G. menardii, G. hexagonus, G. scitula) using stable isotopic signatures (δ18O and δ13C) and Mg/Ca ratios. We evaluate two aspects of inferred depth habitats: (1) the significance of the apparent calcification depth (ACD) calculation method/equations and (2) regional species-specific ACD controls. Through a comparison with five global, (sub)tropical studies we found the choice of applied equation and δ18Osw significant and an important consideration when comparing with the published literature. The ACDs of the surface mixed layer and thermocline species show a tight clustering between 73–109 m water depth coinciding with the deep chlorophyll maximum (DCM). Furthermore, the ACDs for the sub-thermocline species are positioned relative to secondary peaks in the local primary production. We surmise that food source plays a key role in the relative living depths for the majority of the investigated planktonic foraminifera within this oligotrophic environment of the Maldives and elsewhere in the tropical oceans.
Parabens and sorbic acid are commonly used as food preservatives due to their antimicrobial effect. However, their use in foods for infants and young children is not permitted in the European Union. Previous studies found these compounds in some gel-filled baby teethers, whereby parabens, which are well-known as endocrine disruptors, were identified in the polymer-based chewing surface consisting of ethylene-vinyl acetate (EVA). To assess the exposure of infants and young children to these products, the application of parabens in teethers should be thoroughly investigated. Therefore, the present study aimed to apply a representative migration test procedure combined with an accurate analytical method to examine gel-filled baby teethers without elaborate sample preparation, high costs, and long processing times. Accordingly, solid-phase extraction (SPE), in combination with a stable isotope dilution assay (SIDA) and subsequent gas chromatography–mass spectrometry (GC–MS) for analysis of methyl-, ethyl-, and n-propylparaben (MeP, EtP, and n-PrP), was found to be well-suited, with recoveries ranging from 93 to 99%. The study compared the release of these parabens from intact teether surfaces into water and saliva simulant under real-life conditions, with total amounts of detected parabens found to be in the range of 101–162 µg 100 mL−1 and 57–148 µg 100 mL−1, respectively. Furthermore, as a worst-case scenario, the release into water was examined using a long-term migration study.
About half of present-day cloud condensation nuclei originate from atmospheric nucleation, frequently appearing as a burst of new particles near midday1. Atmospheric observations show that the growth rate of new particles often accelerates when the diameter of the particles is between one and ten nanometres2,3. In this critical size range, new particles are most likely to be lost by coagulation with pre-existing particles4, thereby failing to form new cloud condensation nuclei that are typically 50 to 100 nanometres across. Sulfuric acid vapour is often involved in nucleation but is too scarce to explain most subsequent growth5,6, leaving organic vapours as the most plausible alternative, at least in the planetary boundary layer7,8,9,10. Although recent studies11,12,13 predict that low-volatility organic vapours contribute during initial growth, direct evidence has been lacking. The accelerating growth may result from increased photolytic production of condensable organic species in the afternoon2, and the presence of a possible Kelvin (curvature) effect, which inhibits organic vapour condensation on the smallest particles (the nano-Köhler theory)2,14, has so far remained ambiguous. Here we present experiments performed in a large chamber under atmospheric conditions that investigate the role of organic vapours in the initial growth of nucleated organic particles in the absence of inorganic acids and bases such as sulfuric acid or ammonia and amines, respectively. Using data from the same set of experiments, it has been shown15 that organic vapours alone can drive nucleation. We focus on the growth of nucleated particles and find that the organic vapours that drive initial growth have extremely low volatilities (saturation concentration less than 10−4.5 micrograms per cubic metre). As the particles increase in size and the Kelvin barrier falls, subsequent growth is primarily due to more abundant organic vapours of slightly higher volatility (saturation concentrations of 10−4.5 to 10−0.5 micrograms per cubic metre). We present a particle growth model that quantitatively reproduces our measurements. Furthermore, we implement a parameterization of the first steps of growth in a global aerosol model and find that concentrations of atmospheric cloud concentration nuclei can change substantially in response, that is, by up to 50 per cent in comparison with previously assumed growth rate parameterizations.
New geochemical data from the Malawi Rift (Chiwondo Beds, Karonga Basin) fill a major spatial gap in our knowledge of hominin adaptations on a continental scale. Oxygen (δ18O), carbon (δ13C), and clumped (Δ47) isotope data on paleosols, hominins, and selected fauna elucidate an unexpected diversity in the Pleistocene hominin diet in the various habitats of the East African Rift System (EARS). Food sources of early Homo and Paranthropus thriving in relatively cool and wet wooded savanna ecosystems along the western shore of paleolake Malawi contained a large fraction of C3 plant material. Complementary water consumption reconstructions suggest that ca. 2.4 Ma, early Homo (Homo rudolfensis) and Paranthropus (Paranthropus boisei) remained rather stationary near freshwater sources along the lake margins. Time-equivalent Paranthropus aethiopicus from the Eastern Rift further north in the EARS consumed a higher fraction of C4 resources, an adaptation that grew more pronounced with increasing openness of the savanna setting after 2 Ma, while Homo maintained a high versatility. However, southern African Paranthropus robustus had, similar to the Malawi Rift individuals, C3-dominated feeding strategies throughout the Early Pleistocene. Collectively, the stable isotope and faunal data presented here document that early Homo and Paranthropus were dietary opportunists and able to cope with a wide range of paleohabitats, which clearly demonstrates their high behavioral flexibility in the African Early Pleistocene.
The growth of freshly formed aerosol particles can be the bottleneck in their survival to cloud condensation nuclei. It is therefore crucial to understand how particles grow in the atmosphere. Insufficient experimental data has impeded a profound understanding of nano-particle growth under atmospheric conditions. Here we study nano-particle growth in the CLOUD (Cosmics Leaving OUtdoors Droplets) chamber, starting from the formation of molecular clusters. We present measured growth rates at sub-3 nm sizes with different atmospherically relevant concentrations of sulphuric acid, water, ammonia and dimethylamine. We find that atmospheric ions and small acid-base clusters, which are not generally accounted for in the measurement of sulphuric acid vapour, can participate in the growth process, leading to enhanced growth rates. The availability of compounds capable of stabilizing sulphuric acid clusters governs the magnitude of these effects and thus the exact growth mechanism. We bring these observations into a coherent framework and discuss their significance in the atmosphere.
Here we present a comprehensive attempt to correlate aragonitic Na / Ca ratios from Lophelia pertusa, Madrepora oculata and a caryophylliid cold-water coral (CWC) species with different seawater parameters such as temperature, salinity and pH. Living CWC specimens were collected from 16 different locations and analyzed for their Na / Ca content using solution-based inductively coupled plasma-optical emission spectrometry (ICP-OES) measurements. The results reveal no apparent correlation with salinity (30.1–40.57 g/kg) but a significant inverse correlation with temperature (−0.31 mmol/mol/°C). Other marine aragonitic organisms such as Mytilus edulis (inner aragonitic shell portion) and Porites sp. exhibit similar results highlighting the consistency of the calculated CWC regressions. Corresponding Na / Mg ratios show a similar temperature sensitivity to Na / Ca ratios, but the combination of two ratios appear to reduce the impact of vital effects and domain-dependent geochemical variation. The high degree of scatter and elemental heterogeneities between the different skeletal features in both Na / Ca and Na / Mg however limit the use of these ratios as a proxy and/or make a high number of samples necessary. Additionally, we explore two models to explain the observed temperature sensitivity of Na / Ca ratios for an open and semi-enclosed calcifying space based on temperature sensitive Na and Ca pumping enzymes and transport proteins that change the composition of the calcifying fluid and consequently the skeletal Na / Ca ratio.
An accelerating Brewer-Dobson circulation (BDC) is a robust signal of climate change in model predictions but has been questioned by trace gas observations. We analyze stratospheric mean age of air and the full age spectrum as measures for the BDC and its trend. Age of air is calculated with the Chemical Lagrangian Model of the Stratosphere (CLaMS) driven by ERA-Interim, JRA-55 and MERRA-2 reanalysis data to assess the robustness of the representation of the BDC in current generation meteorological reanalyses. We find that climatological mean age significantly depends on the reanalysis, with JRA-55 showing the youngest and MERRA-2 the oldest mean age. Consideration of the age spectrum indicates that the older age for MERRA-2 is related to a stronger spectrum tail, likely related to weaker tropical upwelling and stronger recirculation. Seasonality of stratospheric transport is robustly represented in reanalyses, with similar mean age variations and age spectrum peaks. Long-term changes over 1989–2015 turn out to be similar for the reanalyses with mainly decreasing mean age accompanied by a shift of the age spectrum peak towards shorter transit times, resembling the forced response in climate model simulations to increasing greenhouse gas concentrations. For the shorter periods 1989–2001 and 2002–2015 age of air changes are less robust. Only ERA-Interim shows the hemispheric dipole pattern in age changes during 2002–2015 as viewed by recent satellite observations. Consequently, the representation of decadal variability of the BDC in current generation reanalyses appears less robust and a major uncertainty of modelling the BDC.
The Paleocene-Eocene Thermal Maximum (PETM) offers insight into massive short-term carbon cycle perturbations that caused significant warming during a high-pCO2 world, affecting both marine and terrestrial ecosystems. PETM records from the marine-terrestrial interface (e.g. estuarine swamps and mire deposits) are, therefore, of great interest as their present-day counterparts are highly vulnerable to future climate and sea level change. Here, we assess paleoenvironmental changes of mid-latitudinal Late Paleocene-Early Eocene peat mire records along the paleo-North Sea coast. We provide carbon isotope data of bulk organic matter (δ13CTOC), organic carbon content (%TOC), and palynological data from an extensive peat mire deposited at a mid-latitudinal (ca. 41 °N) coastal site (Schöningen, Germany). The δ13CTOC data show a carbon isotope excursion (CIE) of −1.7 ‰ coeval with a conspicuous Apectodinium acme, calling for the presence of the PETM in this coastal section. Due to the exceptionally large stratigraphic thickness of the PETM at Schöningen (10 m of section) we established a detailed palynological record that indicates only minor changes in paleovegetation leading to and during the PETM. Instead, paleovegetation changes mostly follow natural successions in response to changes along the marine-terrestrial interface. Compared to other available peat mire records (Cobham, UK; Vasterival, France) it appears that wetland deposits around the Paleogene North Sea have a typical CIE magnitude of ca. −1.3 ‰ in δ13CTOC. Moreover, the Schöningen record shares major characteristics with the Cobham Lignite, including evidence for increased fire activity prior to the PETM, minor PETM-related plant species changes, a reduced CIE in δ13CTOC, and drowning of the mire (marine ingressions) during much of the PETM. This suggests that paleoenvironmental conditions during the Late Paleocene-Early Eocene, including the PETM, consistently affected major segments of the paleo-North Sea coast.
In global hydrological models, groundwater storages and flows are generally simulated by linear reservoir models. Recently, the first global gradient-based groundwater models were developed in order to improve the representation of groundwater-surface water interactions, capillary rise, lateral flows and human water use impacts. However, the reliability of model outputs is limited by a lack of data as well as model assumptions required due to the necessarily coarse spatial resolution. The impact of data quality is presented by showing the sensitivity of a groundwater model to changes in the only available global hydraulic conductivity data-set. To better understand the sensitivity of model output to uncertain spatially distributed parameter inputs, we present the first application of a global sensitivity method for a global-scale groundwater model using nearly 2000 steady-state model runs of the global gradient-based groundwater model G3M. By applying the Morris method in a novel domain decomposition approach that identifies global hydrological response units, spatially distributed parameter sensitivities are determined for a computationally expensive model. Results indicate that globally simulated hydraulic heads are equally sensitive to hydraulic conductivity, groundwater recharge and surface water body elevation, though parameter sensitivities vary regionally. For large areas of the globe, rivers are simulated to be either losing or gaining, depending on the parameter combination, indicating a high uncertainty of simulating the direction of flow between the two compartments. Mountainous and dry regions show a high variance in simulated head due to numerical difficulties of the model, limiting the reliability of computed sensitivities in these regions. This instability is likely caused by the uncertainty in surface water body elevation. We conclude that maps of spatially distributed sensitivities can help to understand complex behaviour of models that incorporate data with varying spatial uncertainties. The findings support the selection of possible calibration parameters and help to anticipate challenges for a transient coupling of the model.
Global impact models represent process-level understanding of how natural and human systems may be affected by climate change. Their projections are used in integrated assessments of climate change. Here we test, for the first time, systematically across many important systems, how well such impact models capture the impacts of extreme climate conditions. Using the 2003 European heat wave and drought as a historical analogue for comparable events in the future, we find that a majority of models underestimate the extremeness of impacts in important sectors such as agriculture, terrestrial ecosystems, and heat-related human mortality, while impacts on water resources and hydropower are overestimated in some river basins; and the spread across models is often large. This has important implications for economic assessments of climate change impacts that rely on these models. It also means that societal risks from future extreme events may be greater than previously thought.
In late 2013, a whole air flask collection program started at the Taunus Observatory (TO) in central Germany. Being a rural site in close vicinity to the densely populated Rhein-Main area, Taunus Observatory allows to assess local and regional emissions. Owed to its altitude of 825 m, the site also regularly experiences background conditions, especially when air masses approach from north-westerly directions. With a large footprint area mainly covering central Europe north of the Alps, halocarbon measurements at the site have the potential to improve the data base for estimation of regional and total European halogenated greenhouse gas emissions. Flask samples are collected weekly for offline analysis using a GC-MS system employing a quadrupole as well as a time-of-flight mass spectrometer. As background reference, additional samples are collected approximately bi-weekly at the Mace Head Atmospheric Research Station (MHD) when air masses approach from the site’s clean air sector. Thus the TO time series can be linked to the in-situ AGAGE measurements and the NOAA flask sampling program at MHD. An iterative baseline identification procedure separates polluted samples from baseline data. While there is good agreement of baseline mixing ratios between TO and MHD, with a larger variability of mixing ratios at the continental site, measurements at TO are regularly influenced by elevated halocarbon mixing ratios. Here, first time series are presented for CFC-11, CFC-12, HCFC-22, HFC-134a, HFC-227ea, HFC-245fa, and dichloromethane. While atmospheric mixing ratios of the CFCs decrease, they increase for the HCFC and the HFCs. Small unexpected differences between CFC-11 and CFC-12 are found with regard to the occurrence of high mixing ratio events and seasonality, although production and use of both compounds are strictly regulated by the Montreal Protocol, and therefore a similar decrease of atmospheric mixing ratios should occur. Dichloromethane, a solvent about which recently concerns have risen regarding its growing influence on stratospheric ozone depletion, does not show a significant trend with regard to both, baseline mixing ratios and the occurrence of pollution events at Taunus Observatory for the time period covered, indicating stable emissions in the regions that influence the site. An analysis of HYSPLIT trajectories reveals differences in halocarbon mixing ranges depending on air mass origin.
In late 2013, a whole air flask collection programme was started at Taunus Observatory (TO) in central Germany. Being a rural site in close proximity to the Rhine–Main area, Taunus Observatory allows assessment of emissions from a densely populated region. Owing to its altitude of 825 m, the site also regularly experiences background conditions, especially when air masses approach from north-westerly directions. With a large footprint area mainly covering central Europe north of the Alps, halocarbon measurements at the site have the potential to improve the database for estimation of regional and total European halogenated greenhouse gas emissions. Flask samples are collected weekly for offline analysis using a GC/MS system simultaneously employing a quadrupole as well as a time-of-flight mass spectrometer. As background reference, additional samples are collected approximately once every 2 weeks at the Mace Head Atmospheric Research Station (MHD) when air masses approach from the site's clean air sector. Thus the time series at TO can be linked to the in situ AGAGE measurements and the NOAA flask sampling programme at MHD. An iterative baseline identification procedure separates polluted samples from baseline data. While there is good agreement of baseline mixing ratios between TO and MHD, with a larger variability of mixing ratios at the continental site, measurements at TO are regularly influenced by elevated halocarbon mixing ratios. Here, first time series are presented for CFC-11, CFC-12, HCFC-22, HFC-134a, HFC-227ea, HFC-245fa, and dichloromethane. While atmospheric mixing ratios of the chlorofluorocarbons (CFCs) decrease, they increase for the hydrochlorofluorocarbons (HCFCs) and the hydrofluorocarbons (HFCs). Small unexpected differences between CFC-11 and CFC-12 are found with regard to frequency and relative enhancement of high mixing ratio events and seasonality, although production and use of both compounds are strictly regulated by the Montreal Protocol, and therefore a similar decrease in atmospheric mixing ratios should occur. Dichloromethane, a solvent about which recently concerns have been raised regarding its growing influence on stratospheric ozone depletion, does not show a significant trend with regard to both baseline mixing ratios and the occurrence of pollution events at Taunus Observatory for the time period covered, indicating stable emissions in the regions that influence the site. An analysis of trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model reveals differences in halocarbon mixing ranges depending on air mass origin.
A new method for size resolved chemical analysis of nucleation mode aerosol particles (size range from ~10 to ~30 nm) is presented. The Thermal Desorption Differential Mobility Analyzer (TD-DMA) uses an online, discontinuous principle. The particles are charged, a specific size is selected by differential mobility analysis and they are collected on a filament by electrostatic precipitation. Subsequently, the sampled mass is evaporated in a clean carrier gas and analyzed by a chemical ionization mass spectrometer. Gas phase measurements are performed with the same mass spectrometer during the sampling of particles. The characterization shows reproducible results, with a particle size resolution of 1.19 and the transmission efficiency for 15 nm particles being slightly above 50 %. The signal from the evaporation of a test substance can be detected starting from 0.01 ng and shows a linear response in the mass spectrometer. Instrument operation in the range of pg/m3 is demonstrated by an example measurement of 15 nm particles produced by nucleation from dimethylamine, sulfuric acid and water.
A new method for size-resolved chemical analysis of nucleation mode aerosol particles (size range from ∼10 to ∼30 nm) is presented. The Thermal Desorption Differential Mobility Analyzer (TD-DMA) uses an online, discontinuous principle. The particles are charged, a specific size is selected by differential mobility analysis and they are collected on a filament by electrostatic precipitation. Subsequently, the sampled mass is evaporated in a clean carrier gas and analyzed by a chemical ionization mass spectrometer. Gas-phase measurements are performed with the same mass spectrometer during the sampling of particles. The characterization shows reproducible results, with a particle size resolution of 1.19 and the transmission efficiency for 15 nm particles being slightly above 50 %. The signal from the evaporation of a test substance can be detected starting from 0.01 ng and shows a linear response in the mass spectrometer. Instrument operation in the range of pg m−3 is demonstrated by an example measurement of 15 nm particles produced by nucleation from dimethylamine, sulfuric acid and water.
Increasing atmospheric CO2 stimulates photosynthesis which can increase net primary production (NPP), but at longer timescales may not necessarily increase plant biomass. Here we analyse the four decade-long CO2-enrichment experiments in woody ecosystems that measured total NPP and biomass. CO2 enrichment increased biomass increment by 1.05 ± 0.26 kg C m−2 over a full decade, a 29.1 ± 11.7% stimulation of biomass gain in these early-secondary-succession temperate ecosystems. This response is predictable by combining the CO2 response of NPP (0.16 ± 0.03 kg C m−2 y−1) and the CO2-independent, linear slope between biomass increment and cumulative NPP (0.55 ± 0.17). An ensemble of terrestrial ecosystem models fail to predict both terms correctly. Allocation to wood was a driver of across-site, and across-model, response variability and together with CO2-independence of biomass retention highlights the value of understanding drivers of wood allocation under ambient conditions to correctly interpret and predict CO2 responses.
Formation of new aerosol particles from trace gases is a major source of cloud condensation nuclei (CCN) in the global atmosphere, with potentially large effects on cloud optical properties and Earth’s radiative balance. Controlled laboratory experiments have resolved, in detail, the different nucleation pathways likely responsible for atmospheric new particle formation, yet very little is known from field studies about the molecular steps and compounds involved in different regions of the atmosphere. The scarcity of primary particle sources makes secondary aerosol formation particularly important in the Antarctic atmosphere. Here, we report on the observation of ion-induced nucleation of sulfuric acid and ammonia—a process experimentally investigated by the CERN CLOUD experiment—as a major source of secondary aerosol particles over coastal Antarctica. We further show that measured high sulfuric acid concentrations, exceeding 107 molecules cm−3, are sufficient to explain the observed new particle growth rates. Our findings show that ion-induced nucleation is the dominant particle formation mechanism, implying that galactic cosmic radiation plays a key role in new particle formation in the pristine Antarctic atmosphere.
During the Late Cretaceous and early Cenozoic the Earth experienced prolonged climatic cooling most likely caused by decreasing volcanic activity and atmospheric CO2 levels. However, the causes and mechanisms of subsequent major global warming culminating in the late Paleocene to Eocene greenhouse climate remain enigmatic. We present deep and intermediate water Nd-isotope records from the North and South Atlantic to decipher the control of the opening Atlantic Ocean on ocean circulation and its linkages to the evolution of global climate. The marked convergence of Nd-isotope signatures 59 million years ago indicates a major intensification of deep-water exchange between the North and South Atlantic, which coincided with the turning point of deep-water temperatures towards early Paleogene warming. We propose that this intensification of Atlantic overturning circulation in concert with increased atmospheric CO2 from continental rifting marked a climatic tipping point contributing to a more efficient distribution of heat over the planet.
The multi-valence nature of vanadium means that its geochemical behaviour will be ƒO2-dependent, so that its concentration or V/Sc (or V/Ga), can serve as proxies for oxidation state in mantle peridotites. Compared to Fe3+/Fe2+-based equilibria, such trace elements may be less sensitive to metasomatic processes. To investigate these systematics, we have measured V, Sc, Ga and Fe3+ contents in clinopyroxene from well-characterised spinel peridotite xenoliths from the Massif Central, France. These samples were metasomatised by a variety of agents with different oxidation states.V contents can be modified by metasomatic interactions, and other geochemically similar elements including Sc and Ga can also be added, removed or remain constant. A link between V/Sc and Fe3+-Fe2+ equilibria is apparent. Partial removal of V is caused by different metasomatic agents; the common factor is that all agents were significantly more oxidised than the initial ambient mantle peridotite. This extraction can be understood by a decreasing partition coefficient for V for ΔlogƒO2 > ~FMQ-2. Considering that mineral/melt partitioning of V decreases similarly for all peridotite minerals, the bulk-rock V/Sc will also change during relatively oxidising metasomatic interactions and mirror the results obtained for clinopyroxene.
Nucleation and growth of aerosol particles from atmospheric vapors constitutes a major source of global cloud condensation nuclei (CCN). The fraction of newly formed particles that reaches CCN sizes is highly sensitive to particle growth rates, especially for particle sizes <10 nm, where coagulation losses to larger aerosol particles are greatest. Recent results show that some oxidation products from biogenic volatile organic compounds are major contributors to particle formation and initial growth. However, whether oxidized organics contribute to particle growth over the broad span of tropospheric temperatures remains an open question, and quantitative mass balance for organic growth has yet to be demonstrated at any temperature. Here, in experiments performed under atmospheric conditions in the Cosmics Leaving Outdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN), we show that rapid growth of organic particles occurs over the range from −25 ∘C to 25 ∘C. The lower extent of autoxidation at reduced temperatures is compensated by the decreased volatility of all oxidized molecules. This is confirmed by particle-phase composition measurements, showing enhanced uptake of relatively less oxygenated products at cold temperatures. We can reproduce the measured growth rates using an aerosol growth model based entirely on the experimentally measured gas-phase spectra of oxidized organic molecules obtained from two complementary mass spectrometers. We show that the growth rates are sensitive to particle curvature, explaining widespread atmospheric observations that particle growth rates increase in the single-digit-nanometer size range. Our results demonstrate that organic vapors can contribute to particle growth over a wide range of tropospheric temperatures from molecular cluster sizes onward.
Subduction zone magmas are more oxidised on eruption than those at mid-ocean ridges. This is attributed either to oxidising components, derived from subducted lithosphere (slab) and added to the mantle wedge, or to oxidation processes occurring during magma ascent via differentiation. Here we provide direct evidence for contributions of oxidising slab agents to melts trapped in the sub-arc mantle. Measurements of sulfur (S) valence state in sub-arc mantle peridotites identify sulfate, both as crystalline anhydrite (CaSO4) and dissolved SO42− in spinel-hosted glass (formerly melt) inclusions. Copper-rich sulfide precipitates in the inclusions and increased Fe3+/∑Fe in spinel record a S6+–Fe2+ redox coupling during melt percolation through the sub-arc mantle. Sulfate-rich glass inclusions exhibit high U/Th, Pb/Ce, Sr/Nd and δ34S (+ 7 to + 11‰), indicating the involvement of dehydration products of serpentinised slab rocks in their parental melt sources. These observations provide a link between liberated slab components and oxidised arc magmas.
Dating of extensive alluvial fan surfaces and fluvial features in the hyperarid core of the Atacama Desert, Chile, using cosmogenic nuclides provides unrivalled insights about the onset and variability of aridity. The predominantly hyperarid conditions help to preserve the traces of episodic climatic and/or slow tectonic change. Utilizing single clast exposure dating with cosmogenic 10Be and 21Ne, we determine the termination of episodes of enhanced fluvial erosion and deposition occurring at ~19, ~14, ~9.5 Ma; large scale fluvial modification of the landscape had ceased by ~2–3 Ma. The presence of clasts that record pre-Miocene exposure ages (~28 Ma and ~34 Ma) require stagnant landscape development during the Oligocene. Our data implies an early onset of (hyper-) aridity in the core region of the Atacama Desert, interrupted by wetter but probably still arid periods. The apparent conflict with interpretation that favour a later onset of (hyper-) aridity can be reconciled when the climatic gradients within the Atacama Desert are considered.
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.
During the Holocene, North American ice sheet collapse and rapid sea-level rise reconnected the Black Sea with the global ocean. Rapid meltwater releases into the North Atlantic and associated climate change arguably slowed the pace of Neolithisation across southeastern Europe, originally hypothesized as a catastrophic flooding that fueled culturally-widespread deluge myths. However, we currently lack an independent record linking the timing of meltwater events, sea-level rise and environmental change with the timing of Neolithisation in southeastern Europe. Here, we present a sea surface salinity record from the Northern Aegean Sea indicative of two meltwater events at ~8.4 and ~7.6 kiloyears that can be directly linked to rapid declines in the establishment of Neolithic sites in southeast Europe. The meltwater events point to an increased outflow of low salinity water from the Black Sea driven by rapid sea level rise >1.4 m following freshwater outbursts from Lake Agassiz and the final decay of the Laurentide ice sheet. Our results shed new light on the link between catastrophic sea-level rise and the Neolithisation of southeastern Europe, and present a historical example of how coastal populations could have been impacted by future rapid sea-level rise.
A recent CLOUD (Cosmics Leaving OUtdoor Droplets) chamber study showed that sulfuric acid and dimethylamine produce new aerosols very efficiently, and yield particle formation rates that are compatible with boundary layer observations. These previously published new particle formation (NPF) rates are re-analyzed in the present study with an advanced method. The results show that the NPF rates at 1.7 nm are more than a factor of 10 faster than previously published due to earlier approximations in correcting particle measurements made at larger detection threshold. The revised NPF rates agree almost perfectly with calculated rates from a kinetic aerosol model at different sizes (1.7 nm and 4.3 nm mobility diameter). In addition, modeled and measured size distributions show good agreement over a wide range (up to ca. 30 nm). Furthermore, the aerosol model is modified such that evaporation rates for some clusters can be taken into account; these evaporation rates were previously published from a flow tube study. Using this model, the findings from the present study and the flow tube experiment can be brought into good agreement. This confirms that nucleation proceeds at rates that are compatible with collision-controlled (a.k.a. kinetically-controlled) new particle formation for the conditions during the CLOUD7 experiment (278 K, 38% RH, sulfuric acid concentration between 1×106 and 3×107 cm-3 and dimethylamine mixing ratio of ~40 pptv). Finally, the simulation of atmospheric new particle formation reveals that even tiny mixing ratios of dimethylamine (0.1 pptv) yield NPF rates that could explain significant boundary layer particle formation. This highlights the need for improved speciation and quantification techniques for atmospheric gas-phase amine measurements.
A recent CLOUD (Cosmics Leaving OUtdoor Droplets) chamber study showed that sulfuric acid and dimethylamine produce new aerosols very efficiently and yield particle formation rates that are compatible with boundary layer observations. These previously published new particle formation (NPF) rates are reanalyzed in the present study with an advanced method. The results show that the NPF rates at 1.7 nm are more than a factor of 10 faster than previously published due to earlier approximations in correcting particle measurements made at a larger detection threshold. The revised NPF rates agree almost perfectly with calculated rates from a kinetic aerosol model at different sizes (1.7 and 4.3 nm mobility diameter). In addition, modeled and measured size distributions show good agreement over a wide range of sizes (up to ca. 30 nm). Furthermore, the aerosol model is modified such that evaporation rates for some clusters can be taken into account; these evaporation rates were previously published from a flow tube study. Using this model, the findings from the present study and the flow tube experiment can be brought into good agreement for the high base-to-acid ratios (∼ 100) relevant for this study. This confirms that nucleation proceeds at rates that are compatible with collision-controlled (a.k.a. kinetically controlled) NPF for the conditions during the CLOUD7 experiment (278 K, 38 % relative humidity, sulfuric acid concentration between 1 × 106 and 3 × 107 cm−3, and dimethylamine mixing ratio of ∼ 40 pptv, i.e., 1 × 109 cm−3).
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.
Convection-permitting models (CPMs) have proven their usefulness in representing precipitation on a sub-daily scale. However, investigations on sub-hourly scales are still lacking, even though these are the scales for which showers exhibit the most variability. A Lagrangian approach is implemented here to evaluate the representation of showers in a CPM, using the limited-area climate model COSMO-CLM. This approach consists of tracking 5‑min precipitation fields to retrieve different features of showers (e.g., temporal pattern, horizontal speed, lifetime). In total, 312 cases are simulated at a resolution of 0.01 ° over Central Germany, and among these cases, 78 are evaluated against a radar dataset. The model is able to represent most observed features for different types of convective cells. In addition, the CPM reproduced well the observed relationship between the precipitation characteristics and temperature indicating that the COSMO-CLM model is sophisticated enough to represent the climatological features of showers.
Chlorine monoxide (ClO) plays a key role in stratospheric ozone loss processes at midlatitudes. We present two balloon-borne in situ measurements of ClO conducted in northern hemisphere midlatitudes during the period of the maximum of total inorganic chlorine loading in the atmosphere. Both ClO measurements were conducted on board the TRIPLE balloon payload, launched in November 1996 in León, Spain, and in May 1999 in Aire sur l'Adour, France. For both flights a ClO daylight and night-time vertical profile was derived over an altitude range of approximately 15-35 km. ClO mixing ratios are compared to model simulations performed with the photochemical box model version of the Chemical Lagrangian Model of the Stratosphere (CLaMS). Simulations along 24-hour backward trajectories were performed to study the diurnal variation of ClO in the midlatitude lower stratosphere. Model simulations for the flight launched in Aire sur l'Adour 1999 show an excellent agreement with the ClO measurements. For the flight launched in León 1996, an overall good agreement is found, whereas the flight is characterized by a more complex dynamical situation due to a possible mixture of vortex and non-vortex air. We note that for both flights at solar zenith angles greater than 86°-87° simulated ClO mixing ratios are higher than observed ClO mixing ratios. However, the present findings indicate that no substantial uncertainties exist in midlatitude chlorine chemistry of the stratosphere.
Fractional release factors of long-lived halogenated organic compounds in the tropical stratosphere
(2009)
Fractional release factors (FRFs) of organic trace gases are time-independent quantities that influence the calculation of Global Warming Potentials and Ozone Depletion Potentials. We present the first set of vertically resolved FRFs for 15 long-lived halo carbons in the tropical stratosphere up to 34 km altitude. They were calculated from measurements on air samples collected on board balloons and a high altitude aircraft. We compare the derived dependencies of FRFs on the mean stratospheric transit times (the so-called mean ages of air) with similarly derived FRFs originating from measurements at higher latitudes and find significant differences. Moreover a comparison with averaged FRFs currently used by the World Meteorological Organisation revealed the latter to be imprecise measures due to their observed vertical and latitudinal variability. The presented data set could thus be used to improve future ozone level and climate projections.
Brachiopod shells are the most widely used geological archive for the reconstruction of the temperature and the oxygen isotope composition of Phanerozoic seawater. However, it is not conclusive whether brachiopods precipitate their shells in thermodynamic equilibrium. In this study, we investigated the potential impact of kinetic controls on the isotope composition of modern brachiopods by measuring the oxygen and clumped isotope compositions of their shells. Our results show that clumped and oxygen isotope compositions depart from thermodynamic equilibrium due to growth rate-induced kinetic effects. These departures are in line with incomplete hydration and hydroxylation of dissolved CO2. These findings imply that the determination of taxon-specific growth rates alongside clumped and bulk oxygen isotope analyses is essential to ensure accurate estimates of past ocean temperatures and seawater oxygen isotope compositions from brachiopods.
We generated transgenic tomato plants with altered expression of heat stress transcription factor HsfA1. Plants with 10-fold overexpression of HsfA1 (OE plants) were characterized by a singleHsfA1 transgene cassette, whereas plants harboring a tandem inverted repeat of the cassette showed cosuppression (CS plants) by posttranscriptional silencing of the HsfA1 gene connected with formation of small interfering RNAs. Under normal growth conditions, major developmental parameters were similar for wild-type (WT), OE, and CS plants. However, CS plants and fruits were extremely sensitive to elevated temperatures, because heat stress-induced synthesis of chaperones and Hsfs was strongly reduced or lacking. Despite the complexity of the plant Hsf family with at least 17 members in tomato, HsfA1 has a unique function as master regulator for induced thermotolerance. Using transient reporter assays with mesophyll protoplasts from WT tomato, we demonstrated that plasmid-encoded HsfA1 and HsfA2 were well expressed. However, in CS protoplasts the cosuppression phenomenon was faithfully reproduced. Only transformation with HsfA2 expression plasmid led to normal expression of the transcription factor and reporter gene activation, whereas even high amounts of HsfA1 expression plasmids were silenced. Thermotolerance in CS protoplasts was restored by plasmid-borne HsfA2, resulting in expression of chaperones, thermoprotection of firefly luciferase, and assembly of heat stress granules.
Size-resolved measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a full seasonal cycle (Mar 2014–Feb 2015). In a companion part 1 paper, we presented an in-depth CCN characterization based on annually as well as seasonally averaged time intervals and discuss different parametrization strategies to represent the Amazonian CCN cycling in modelling studies (M. Pöhlker et al., 2016b). The present part 2 study analyzes the aerosol and CCN variability in original time resolution and, thus, resolves aerosol advection and transformation for the following case studies, which represent the most characteristic states of the Amazonian atmosphere:
1. Near-pristine (NP) conditions, defined as the absence of detectable black carbon (< 0.01 µg m−3), showed their highest occurrence (up to 30 %) in the wet season (i.e., Mar–May). On average, the NP episodes are characterized by a bimodal aerosol size distribution (strong Aitken mode: DAit = 70 nm, NAit = ~ 200 cm−3 vs. weaker accumulation mode: Dacc = 170 nm, Nacc = ~ 60 cm−3), a mostly organic particle composition, and relatively low hygroscopicity levels (κAit = 0.12 vs. κacc = 0.18). The NP CCN efficiency spectrum shows that the CCN population is sensitive to changes in supersaturation (S) over a wide S range.
2. Long-range transport (LRT) conditions frequently mix Saharan dust, African combustion smoke, and sea spray aerosols into the Amazonian wet season atmosphere. The LRT episodes (i.e., Feb–Apr) are characterized by an accumulation mode dominated size distribution (DAit = 80 nm, NAit = 120 cm−3 vs. Dacc = 180 nm, Nacc = 300 cm−3), a clearly increased abundance of dust and salt compounds, and relatively high hygroscopicity levels (κAit = 0.18, κacc = 0.34). The LRT CCN efficiency spectrum shows that the CCN population is highly sensitive to changes in S in the low S regime.
3. Biomass burning (BB) conditions dominate the Amazonian dry season. A selected characteristic BB episode shows a very strong accumulation mode (DAit = 70 nm, NAit = ~ 140 cm−3 vs. Dacc = 170 nm, Nacc = ~ 3400 cm−3), particles with very high organic fractions (> 90 %), and correspondingly low hygroscopicity levels (κAit = 0.14, κacc = 0.17). The BB CCN efficiency spectrum shows that the CCN population is highly sensitive to changes in S in the low S regime.
4. Mixed pollution conditions show the superposition of African (i.e., volcanic) and Amazonian (i.e., biomass burning) aerosol emissions during the dry season. The African aerosols showed a broad monomodal distribution (D = 130 nm, N = ~ 1300 cm−3), with very high sulfate fractions (20 %), and correspondingly high hygroscopicity (κAit = 0.14, κacc = 0.22). This was superimposed by fresh smoke from nearby fires with one strong mode (D = 113 nm, Nacc = ~ 2800 cm−3), an organic-dominated aerosol, and sharply decreased hygroscopicity (κAit = 0.10, κacc = 0.20). These conditions underline the rapidly changing pollution regimes with clear impacts on the aerosol and CCN properties.
Overall, this study provides detailed insights into the CCN cycling in relation to aerosol-cloud interaction in the vulnerable and climate-relevant Amazon region. The detailed analysis of aerosol and CCN key properties and particularly the extracted CCN efficiency spectra with the associated fit parameters provide a basis for an in-depth analysis of aerosol-cloud interaction in the Amazon and beyond.
We developed a coupled regional climate system model based on the CCLM regional climate model. Within this model system, using OASIS3-MCT as a coupler, CCLM can be coupled to two land surface models (the Community Land Model (CLM) and VEG3D), the NEMO-MED12 regional ocean model for the Mediterranean Sea, two ocean models for the North and Baltic seas (NEMO-NORDIC and TRIMNP+CICE) and the MPI-ESM Earth system model.
We first present the different model components and the unified OASIS3-MCT interface which handles all couplings in a consistent way, minimising the model source code modifications and defining the physical and numerical aspects of the couplings. We also address specific coupling issues like the handling of different domains, multiple usage of the MCT library and exchange of 3-D fields.
We analyse and compare the computational performance of the different couplings based on real-case simulations over Europe. The usage of the LUCIA tool implemented in OASIS3-MCT enables the quantification of the contributions of the coupled components to the overall coupling cost. These individual contributions are (1) cost of the model(s) coupled, (2) direct cost of coupling including horizontal interpolation and communication between the components, (3) load imbalance, (4) cost of different usage of processors by CCLM in coupled and stand-alone mode and (5) residual cost including i.a. CCLM additional computations.
Finally a procedure for finding an optimum processor configuration for each of the couplings was developed considering the time to solution, computing cost and parallel efficiency of the simulation. The optimum configurations are presented for sequential, concurrent and mixed (sequential+concurrent) coupling layouts. The procedure applied can be regarded as independent of the specific coupling layout and coupling details.
We found that the direct cost of coupling, i.e. communications and horizontal interpolation, in OASIS3-MCT remains below 7 % of the CCLM stand-alone cost for all couplings investigated. This is in particular true for the exchange of 450 2-D fields between CCLM and MPI-ESM. We identified remaining limitations in the coupling strategies and discuss possible future improvements of the computational efficiency.
The assessment of water balance components using global hydrological models is subject to climate forcing uncertainty as well as to an increasing intensity of human water use within the 20th century. The uncertainty of five state-of-the-art climate forcings and the resulting range of cell runoff that is simulated by the global hydrological model WaterGAP is presented. On the global land surface, about 62 % of precipitation evapotranspires, whereas 38 % discharges into oceans and inland sinks. During 1971–2000, evapotranspiration due to human water use amounted to almost 1 % of precipitation, while this anthropogenic water flow increased by a factor of approximately 5 between 1901 and 2010. Deviation of estimated global discharge from the ensemble mean due to climate forcing uncertainty is approximately 4 %. Precipitation uncertainty is the most important reason for the uncertainty of discharge and evapotranspiration, followed by shortwave downward radiation. At continental levels, deviations of water balance components due to uncertain climate forcing are higher, with the highest discharge deviations occurring for river discharge in Africa (−6 to 11 % from the ensemble mean). Uncertain climate forcings also affect the estimation of irrigation water use and thus the estimated human impact of river discharge. The uncertainty range of global irrigation water consumption amounts to approximately 50 % of the global sum of water consumption in the other water use sector.