Refine
Year of publication
- 2016 (73) (remove)
Document Type
- Article (61)
- Doctoral Thesis (11)
- Working Paper (1)
Language
- English (73) (remove)
Has Fulltext
- yes (73)
Is part of the Bibliography
- no (73) (remove)
Keywords
- Atmospheric chemistry (3)
- climate change (2)
- Andes (1)
- Anthropogene Klimaänderung (1)
- BSRN (1)
- Bayesian network (1)
- Biogeochemistry (1)
- COSMO-CLM (1)
- Clausius–Clapeyron scaling (1)
- Climate-change impacts (1)
Institute
- Geowissenschaften (73) (remove)
We discuss applications of a recently developed method for model reduction based on linear response theory of weakly coupled dynamical systems. We apply the weak coupling method to simple stochastic differential equations with slow and fast degrees of freedom. The weak coupling model reduction method results in general in a non-Markovian system; we therefore discuss the Markovianization of the system to allow for straightforward numerical integration. We compare the applied method to the equations obtained through homogenization in the limit of large timescale separation between slow and fast degrees of freedom. We numerically compare the ensemble spread from a fixed initial condition, correlation functions and exit times from a domain. The weak coupling method gives more accurate results in all test cases, albeit with a higher numerical cost.
Water is scarce in semi-arid and arid regions. Using alternative water sources (i.e. non-conventional water sources), such as municipal reuse water and harvested rain, contributes to using existing water resources more efficiently and productively. The aim of this study is to evaluate the two alternative water sources reuse water and harvested rain for the irrigation of small-holder agriculture from a system perspective. This helps decision and policy makers to have proper information about which system and technology to adopt under local conditions. For this, the evaluation included ecologic, societal, economic, institutional and political as well as technical aspects. For the evaluation, the study area in central-northern Namibia was chosen in the frame of the research and development project CuveWaters. The main methods used include a mathematical material flow analysis, the computation and modelling of crop requirements, a multi-criteria decision analysis using the Analytical Hierarchy Process (AHP) method and a financial cost-benefit analysis. From a systemic perspective, the proposed novel systems were compared to the exciting conventional infrastructure. The results showed that both water reuse and rainwater harvesting systems for the irrigation of small-holder horticulture offer numerous technological, ecologic, economic, societal, institutional and political benefits. Rainwater harvesting based gardens have a positive benefit-cost ratio under favorable conditions. Government programs could fund the infrastructure investment costs, while the micro-entrepreneur can assume a micro-credit to finance operation and maintenance costs. Installing sanitation in informal settlements and reusing municipal water for irrigation reduces the overall water demand of households and agriculture by 39%, compared to improving sanitation facilities in informal settlements without reusing the water for agriculture. Given that water is the limiting factor for crop fertigation, the generated nutrient-rich reuse water is sufficient to annually irrigate about 10 m2 to 13 m2 per sanitation user. Compared to crop nutrient requirements, there are too many nutrients in the reuse water. Thus when using nutrient-rich reuse water, no use of fertilizers and a careful salt management is necessary. When comparing this novel system with improved sanitation, advanced wastewater treatment and nutrient-rich water reuse to the conventional and to two adapted systems, results showed that the novel CuveWaters system is the best option for the given context in a semi-arid developing country. Therefore, the results of this study suggest a further roll-out of the novel CuveWaters system. The methodology developed and the results of this study demonstrated that taking sanitation users into consideration plays a major role for the planning of an integrated water reuse infrastructure because they are the determinant factor for the amount of available nutrient-rich reuse water. In addition, it could be shown that water reuse and rainwater harvesting systems for the irrigation of small-scale gardens provide a wide range of benefits and can be key to using scarce water resources more efficiently and to contributing to the Sustainable Development Goals.
This PhD thesis has been carried out within an interdisciplinary cooperational project between the Deutsches Bergbau-Museum Bochum and the Goethe-Universität Frankfurt, which is dedicated to ancient Pb-Ag mining and metal production in the hinterland of the municipium Ulpiana in central Kosovo. Geochemical analysis (OM, XRD, EMP, MC-ICP-MS) of ores, metallurgical (by-) products and metal artefacts allowed to reconstruct the local chaîne opératoire and to decipher significant chronological differences between presumably Roman/late antique and medieval/early modern metallurgical processing. Pb isotope provenance studies documented the relevance of local metal production within the Roman Empire and confirmed the actual existence of a Metalla Dardanica district, which until now solely has been suspected on basis of epigraphy.
The predominant abundance of the by-products matte (Cu, Pb, Fe and Zn sulphides) and speiss (ferrous speiss: Fe-As compounds; base metal speiss: ~(Cu,Ni,Fe,Ag )x(Sb,Sn,As )y ) at smelting sites with a preliminary Roman/late antique dating points to treatment of complex polymetallic ore. Pb isotope analysis demonstrated that the mining district of Shashkoc-Janjevo (partially) supplied six of the ten investigated metallurgical sites. In this mineralisation, parageneses with elevated Cu, As and Sb abundances comprise significant proportions of particularly tennantite-tetrahedrite minerals, chalcopyrite, arsenopyrite and were generated during the early and main stages of ore formation. Later precipitated ore in contrast is marked by a significantly less versatile mineralogy and consists almost exclusively of galena, sphalerite and pyrite/marcasite. Besides increased Cu, As and Sb contents, ore from the main formation stage also exhibits generally higher Ag abundances, which are mainly hosted by fahlore and locally abundant secondary Cu sulphides (chalcocite, digenite and covellite) and oxidised phases (e.g. malachite, azurite). The higher precious metal grades of this ore type, whose geochemical signature (i.e. higher proportions of Cu, As and Sb) is mirrored by the abundance of the metallurgical by-products matte and speiss (almost exclusively found at potentially Roman/late antique smelting sites; see above), presumably were a pivotal factor leading to its preferential exploitation in earlier times. Matte and base metal-rich speiss contain notable amounts of Ag, which are mainly present in Cu-(Fe) sulphides and particularly antimonides ((Cu,Ni)2Sb, Ag3Sb), respectively. While the speiss compounds due to their close association with Pb bullion presumably were cupelled automatically, the metallurgical treatment of matte could not have been proven unambiguously, but overall certainly is highly likely.
The beneficiated ore (i.e. crushed and sorted, potentially also treated by more lavish techniques such as grinding, sieving or wet-mechanical methods) possibly was partially roasted and subsequently together with fluxes and charcoal submitted to the furnaces. The working temperatures approximately ranged between 1100 and 1400 °C. Slags from all presumably Roman/late antique dated and few of their potentially medieval/early modern analogues were produced from smelting of (partially roasted) ore with charcoal and added siliceous material, thus resulting in fayalite-dominant phase assemblages or rarely observed glassy parageneses. Even though several subtypes of fayalite slags have been established on basis of the abundance of Fe-rich oxide phases (i.e. spinel ss and wüstite), late clinopyroxene and the general solidification sequence of the slags, the process conditions (i.e. temperature, fO2, added fluxing agents) must have been widely similar; chemical variations could be explained by varying degrees of interaction of the slag melt with charcoal ash and furnace material. The other investigated metallurgical remains indicate employment of a calcareous flux, which led to formation of Ca-rich olivine-, olivine+clinopyroxene-, clinopyroxene- or melilite-type slags. These types as well as glassy slags were generated at more oxidising conditions outside the fayalite stability field (FMQ buffer equilibrium, cf. Lindsley, 1976) than their olivine-dominant analogues. Conclusions on the furnace construction could be drawn on basis of the typology of the slags, which mostly were tapped into a basin located outside the furnace, but partially (at two presumably medieval/early modern sites) also accumulated in a reservoir within the smelter.
Lead artefacts excavated in Ulpiana could be isotopically related to ores from mineralisations in its vicinity and demonstrate that the resources were at least utilised for local metal production. However, also ship wreck cargo from Israel - including several lead ingots with the inscription 'MET DARD' (Raban, 1999) - and late antique lead-glazed pottery from Serbia and Romania (Walton & Tite, 2010) could be related to a possible Kosovarian/Serbian provenance of the raw material and thus indicate flourishing trade of metal from the Metalla Dardanica district within the Roman Empire.
References:
Lindsley, D. H. (1976). Experimental studies of oxide minerals. In D. Rumble, III (Hrsg.), Oxide minerals (61-88). Reviews in Mineralogy, Volume 3. Washington, DC: Mineralogical Society of America.
Raban, A. (1999). The lead ingots from the wreck site (area K8). Journal of Roman Archaeology, Supplementary Series, 35, 179-188.
Walton, M. S., & Tite, M. S. (2010). Production technology of Roman lead-glazed pottery and its continuance into late antiquity. Archaeometry, 52(5), 733-759.
We present the prototype of a regional climate system model based on the COSMO-CLM regional climate model coupled with several model components, analyze the performance of the couplings and present a strategy to find an optimum configuration with respect to computational costs and time to solution.
The OASIS3-MCT coupler is used to couple COSMO-CLM with two land surface models (CLM and VEG3D), a regional ocean model for the Mediterranean Sea (NEMO-MED12), two ocean models for the North and Baltic Sea (NEMO-NORDIC and TRIMNP+CICE) and the atmospheric component of an earth system model (MPI-ESM). We present a unified OASIS3-MCT interface which handles all couplings in a similar way, minimizes the model source code modifications and describes the physics and numerics of the couplings. Furthermore, we discuss solutions for specific regional coupling problems like handling of different domains, multiple usage of MCT interpolation library and efficient exchange of 3D fields.
A series of real-case simulations over Europe has been conducted and the computational performance of the couplings has been analyzed. The usage of the LUCIA tool of the OASIS3-MCT coupler enabled separation of the direct costs of: coupling, load imbalance and additional computations. The resulting limits for time to solution and costs are shown and the potential of further improvement of the computational efficiency is summarized for each coupling.
It was found that the OASIS3-MCT coupler keeps the direct coupling costs of communication and horizontal interpolation small in comparison with the costs of the additional computations and load imbalance for all investigated couplings. For the first time this could be demonstrated for an exchange of approximately 450 2D fields per time step necessary for the atmosphere-atmosphere coupling between COSMO-CLM and MPI-ESM.
A procedure for finding an optimum configuration for each of the couplings was developed considering the time to solution and costs of the simulations. The optimum configurations are presented for sequential and concurrent coupling layouts. The procedure applied can be regarded as independent on the specific coupling layout and coupling details.
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.
Amines are potentially important for atmospheric new particle formation, but their concentrations are usually low with typical mixing ratios in the pptv range or even smaller. Therefore, the demand for highly sensitive gas-phase amine measurements has emerged in the last several years. Nitrate chemical ionization mass spectrometry (CIMS) is routinely used for the measurement of gasphase sulfuric acid in the sub-pptv range. Furthermore, extremely low volatile organic compounds (ELVOCs) can be detected with a nitrate CIMS. In this study we demonstrate that a nitrate CIMS can also be used for the sensitive measurement of dimethylamine (DMA, (CH3/2NH) using the NO3−•(HNO3)1 − 2• (DMA) cluster ion signal. Calibration measurements were made at the CLOUD chamber during two different measurement campaigns. Good linearity between 0 and ~120 pptv of DMA as well as a sub-pptv detection limit of 0.7 pptv for a 10 min integration time are demonstrated at 278K and 38% RH.
Coccolithophores are an abundant phytoplankton group that exhibit remarkable diversity in their biology, ecology, and calcitic exoskeletons (coccospheres). Their extensive fossil record is testament to their important biogeochemical role and is a valuable archive of biotic responses to environmental change stretching back over 200 million years. However, to realise the potential of this archive requires an understanding of the physiological processes that underpin coccosphere architecture. Using culturing experiments on four modern coccolithophore species (Calcidiscus leptoporus, Calcidiscus quadriperforatus, Helicosphaera carteri and Coccolithus braarudii) from three long-lived families, we investigate how coccosphere architecture responds to shifts from exponential (rapid cell division) to stationary (slowed cell division) growth phases as cell physiology reacts to nutrient depletion. These experiments reveal statistical differences in cell size and the number of coccoliths per cell between these two growth phases, specifically that cells in exponential-phase growth are typically smaller with fewer coccoliths, whereas cells experiencing growth-limiting nutrient depletion have larger coccosphere sizes and greater numbers of coccoliths per cell. Although the exact numbers are species-specific, these growthphase
shifts in coccosphere geometry demonstrate that the core physiological responses of cells to nutrient depletion results in increased cell sizes and coccoliths per cell across four different coccolithophore families (Calcidiscaceae, Coccolithaceae, Isochrysidaceae, Helicosphaeraceae), a representative diversity of this phytoplankton group. Building on this, direct comparison of coccosphere geometries in modern and fossil coccolithophores enables a proxy for growth phase to be developed that allows growth responses to environmental change to be investigated throughout their evolutionary history. Our data also shows that changes in growth rate and coccoliths per cell associated with growth-phase shifts can substantially alter cellular calcite production. Coccosphere geometry is therefore a valuable tool for accessing growth information in the fossil record that provides unprecedented insights into the biotic responses of species to environmental change and its potential biogeochemical consequences.
Ice nucleating particles over the Eastern Mediterranean measured by unmanned aircraft systems
(2016)
During an intensive field campaign on aerosol, clouds and ice nucleation in the Eastern Mediterranean in April 2016, we have measured the abundance of ice nucleating particles (INP) in the lower troposphere from unmanned aircraft systems (UAS). Aerosol samples were collected by miniaturized electrostatic precipitators onboard the UAS at altitudes up to 2.5 km. The number of INP in these samples, which are active in the deposition and condensation modes at temperatures from −20 to −30 ◦C, were analyzed immediately after collection on site using the ice nucleus counter FRIDGE. During the one month campaign we encountered a series of Saharan dust plumes that traveled at several kilometers altitude. Here we present INP data from 42 individual flights, together with aerosol number concentrations, observations of lidar backscattering, dust concentrations derived by the dust transport model DREAM (Dust Regional Atmospheric Model), and results from scanning electron microscopy. The effect of the dust plumes is reflected by the coincidence of INP with the particulate mass (PM), the lidar signal and with the predicted dust mass of the model. This suggests that mineral dust or a constituent related to dust was a major contributor to the ice nucleating properties of the aerosol. Peak concentrations of above 100 INP std.l -1 were measured at −30 ◦C. The INP concentration in elevated plumes was on average a factor of 10 higher than at ground level. Since desert dust is transported for long distances over wide areas of the globe predominantly at several km altitude we conclude that INP measurements at ground level may be of limited significance for the situation at the level of cloud formation.
Recently significant advances have been made in the collection, detection and characterization of ice nucleating particles (INPs). Ice nuclei are particles that facilitate the heterogeneous formation of ice within the atmospheric aerosol by lowering the free energy barrier to spontaneous nucleation and growth of ice from atmospheric water and/or vapor. The Frankfurt isostatic diffusion chamber (FRankfurt Ice nucleation Deposition freezinG Experiment: FRIDGE) is an INP collection and offline detection system that has become widely deployed and shows additional potential for ambient measurements. Since its initial development FRIDGE has gone through several iterations and improvements. Here we describe improvements that have been made in the collection and analysis techniques. We detail the uncertainties inherent in the measurement method and suggest a systematic method of error analysis for FRIDGE measurements. Thus what is presented herein should serve as a foundation for the dissemination of all current and future measurements using FRIDGE instrumentation.
In-situ single particle analysis of ice particle residuals (IPR) and out-of-cloud aerosol particles was conducted by means of laser ablation mass spectrometry during the intensive INUIT-JFJ/CLACE campaign at the high alpine research station Jungfraujoch (3580 m a.s.l.) in January/February 2013. During the four week campaign more than 70000 out-of-cloud aerosol particles and 595 IPR were analyzed covering a particle size diameter range from 100 nm to 3 μm. The IPR were sampled during 273 hours while the station was covered by mixed-phase clouds at ambient temperatures between -27 °C and -6 °C. The identification of particle types is based on laboratory studies of different types of biological, mineral and anthropogenic aerosol particles. As outcome instrument specific marker peaks for the different investigated particle types were obtained and applied to the field data. The results show that the sampled IPR contain a larger relative amount of natural, primary aerosol, like soil dust (13 %) and minerals (11 %), in comparison to out-of-cloud aerosol particles (2 % and <1 %, respectively). Additionally, anthropogenic aerosol particles, like particles from industrial emissions and lead-containing particles, were found to be more abundant in the IPR than in the out-of-cloud aerosol. The out of-cloud aerosol contained a large fraction of aged particles (30 %, including organic material and secondary inorganics), whereas this particle type was much less abundant (3 %) in the IPR. In a selected subset of the data where a direct comparison between out-of-cloud aerosol particles and IPR in air masses with similar origin was possible, a pronounced enhancement of biological particles was found in the IPR.
Despite mounting evidence of the anthropogenic influence on the Earth's climate, underlying mechanisms of climate change often remain elusive. The investigation of periods of rapid climate change from geological archives may provide crucial information about magnitude, duration, teleconnections of and regional responses to global and hemispheric scale climate perturbations. Thus, paleoclimate reconstructions may help in mitigating and adapting to the challenges of the coming decades. The '8.2 kyr B.P. climatic event' has previously been proposed as a possible analogue for the future climatic scenario of a reduced Atlantic Meridional Overturning Circulation (AMOC). The catastrophic drainage of the Laurentide meltwater lakes through the Hudson Bay and into the Labrador Sea, that occurred ca. 8.47 kyr B.P., caused the slowdown of the AMOC around 8.2 kyr B.P.. Subsequently, reduced heat transfer towards Europe triggered a substantial decline in (winter) temperature and pronounced changes in atmospheric circulation patterns in many regions of the northern hemisphere, especially the North Atlantic realm and Europe. Among the regions affected by the 8.2 kyr B.P. climatic event, the Eastern Mediterranean region is of particular interest for both past and future climate developments. Traditionally characterized as a region highly sensitive to variations in the climate systems of the high and low latitudes, abrupt climate changes have the potential to strongly alter atmospheric circulation patterns and thus precipitation distribution in the region that may have severe socioeconomical consequences. The analysis of stable hydrogen (δD) and oxygen isotopes (δ18O) in precipitation is an excellent tool to trace changes in atmospheric circulation. Here, we present a comparative study of δD and δ18O in precipitation from the Eastern Mediterranean region both in a present day scenario and during the 8.2 kyr B.P. climatic event. We analyze the influences of topography, air mass trajectory, climate and seasonality among others the stable isotopic compositions of meteoric waters from the Central Anatolian Plateau (CAP), Turkey, in order to create a first-order template which may serve as a reference against which paleoenvironmental proxy data may be more accurately interpreted and tested. Further, we employ a multiproxy approach on the early Holocene peat deposits of the classical site of Tenaghi Philippon (TP), NE Greece, to investigate paleoenvironmental responses to northern hemisphere cooling during the 8.2 kyr B.P. climatic event and aim to determine changes atmospheric circulation from δD of leaf wax n-alkanes (δDwax).
Based on δD and δ18O data from more than 480 surface water samples from the CAP, we characterize moisture sources affecting the net isotopic budget of precipitation, manifesting in a systematic north-south difference in near-sea level moisture compositions. Rainout, induced by the major orographic barriers of the plateau, the Pontic Mountains to the north and the Taurus Mountains to the south, strongly shape the modern patterns of δD and δ18O. Stable isotope data from the semi-arid plateau interior provide clear evidence for an evaporitic regime that drastically affects surface water compositions. Strong evaporative enrichment contrasts rainfall patterns along the plateau margins, in part obfuscating the effects of topography and air mass trajectory.
Consequently, in order to address possible influences of evaporation on δD and δ18O in paleoprecipitation from TP, we analyze n-alkane abundances and distributions along with stable carbon isotope compositions of total organic carbon (δ13CTOC) and palynological data to estimate surface moisture conditions during the early Holocene (ca. 8.7 - 7.5 kyr B.P.) and especially during the 8.2 kyr B.P. climatic event. A period of relatively dry surface conditions from ca. 8.7 to 8.2 kyr B.P., indicated by low values of the 'aquatic index' (Paq ) and by elevated Average Chain Length (ACL) values, in concert with elevated δ13CTOC values, precedes the 8.2 kyr B.P. climatic event. The event itself is characterized by slightly wetter, more humid conditions, as suggested by an increase in Paq values as well as reduced ACL and δ13CTOC values between ca. 8.2 and 7.9 kyr B.P.. In the upper section of the core, a distinct change in paleohydrology becomes. A steep increase in Paq and a decrease in ACL values as well as variations in δ13CTOC from 7.9 kyr B.P. onwards imply considerably elevated surface moisture levels, likely caused by the increased activity of the karstic system of the surrounding mountains. Collectively, the biomarker proxies presented here, reveal a concise picture of changing moisture conditions at TP that is consistent with palynological data and provide detailed paleoenvironmental information for the analysis of δDwax as a paleoprecipitation proxy. The long-term decline in δDwax values characterizes the lower section of the core until ca. 8.2 kyr B.P.. The 8.2 kyr B.P. climatic event itself is connected to two distinct positive hydrogen isotope excursions: a minor shift in δDwax around 8.2 kyr B.P. and a major shift in δDwax between ca. 8.1 and 8.0 kyr B.P.. The upper part of the section shows a progressive trend towards higher δDwax values. With no indication of increased evaporitic conditions at TP during the 8.2 kyr B.P. climatic event, as evident in biomarker proxies and pollen data, we link shifts in δDwax to changes in Mediterranean air mass trajectories supplying precipitation to northeastern Greece, with variations in the relative contributions of northerly derived, D-depleted moisture and southerly-derived, D-enriched moisture. Possible control mechanisms include changes in the influence of the Siberian High and differences in the influence of the African and Asian monsoon circulations on anticyclonic conditions in the Mediterranean region as well as regional inflow of moist air from the Aegean Sea.
Within the framework of the Transboundary Waters Assessment Programme (TWAP), initiated by the Global Environment Facility (GEF), we contributed to a comprehensive baseline assessment of transboundary aquifers (TBAs) by quantifying different groundwater indicators using the global water resources and water use model WaterGAP 2.2. All indicators were computed under current (2010) and projected conditions in 2030 and 2050 for 91 selected TBAs larger than 20,000 km2 and for each nation’s share of the TBAs (TBA-CU: country unit). TBA outlines were provided by the International Groundwater Resources Assessment Centre (IGRAC). The set of indicators comprises groundwater recharge, groundwater depletion, per-capita groundwater recharge, dependency on groundwater, population density, and groundwater development stress (groundwater withdrawals to groundwater recharge). Only the latter four indicators were projected to 2030 and 2050. Current-state indicators were quantified using the Watch Forcing Data climate dataset, while projections were based on five climate scenarios that were computed by five global climate models for the high-emissions scenario RCP 8.5. Water use projections were based on the Shared Socio-economic Pathway SSP2 developed within ISI-MIP. Furthermore, two scenarios of future irrigated areas were explored. For individual water use sectors, the fraction of groundwater abstraction was assumed to remain at the current level.
According to our assessment, aquifers with the highest current groundwater depletion rates worldwide are not transboundary. Exceptions are the Neogene Aquifer System (Syria) with 53 mm/yr between 2000 and 2009 and the Indus River Plain aquifer (India) with 28 mm/yr. For current conditions, we identified 20 out of 258 TBA-CUs suffering from medium to very high groundwater development stress, which are located in the Middle East and North Africa region, in South Asia, China, and the USA. Considering projections, ensemble means of per-cent changes or percent point changes to current conditions were determined. Per-capita groundwater recharge is projected to decrease in 80-90% of all TBA-CUs until 2030/2050. Due to the strongly varying projections of the global climate models, we applied a worst-case scenario approach to define future hotspots of groundwater development stress, taking into account the strongest computed increase until either 2030 or 2050 among all scenarios and individual GCMs. Based on this approach, the number of TBA-CUs under at least medium groundwater development stress increases from 20 to 58, comprising all hotspots under current conditions. New hotspots are projected to develop mainly in Sub-Saharan Africa, China, and Mexico.
The CARIBIC (Civil Aircraft for the Regular Investigation of the Atmosphere Based on an Instrument Container) passenger aircraft observatory performed in situ measurements at 10–12 km altitude in the South Asian summer monsoon anticyclone between June and September 2008. These measurements enable us to investigate this atmospheric region (which so far has mostly been observed from satellites) using the broad suite of trace gases and aerosol particles measured by CARIBIC. Elevated levels of a variety of atmospheric pollutants (e.g. carbon monoxide, total reactive nitrogen oxides, aerosol particles, and several volatile organic compounds) were recorded. The measurements provide detailed information about the chemical composition of air in different parts of the monsoon anticyclone, particularly of ozone precursors. While covering a range of 3500 km inside the monsoon anticyclone, CARIBIC observations show remarkable consistency, i.e. with distinct latitudinal patterns of trace gases during the entire monsoon period.
Using the CARIBIC trace gas and aerosol particle measurements in combination with the Lagrangian particle dispersion model FLEXPART, we investigated the characteristics of monsoon outflow and the chemical evolution of air masses during transport. The trajectory calculations indicate that these air masses originated mainly from South Asia and mainland Southeast Asia. Estimated photochemical ages of the air were found to agree well with transport times from a source region east of 90–95° E. The photochemical ages of the air in the southern part of the monsoon anticyclone were systematically younger (less than 7 days) and the air masses were mostly in an ozone-forming chemical mode. In its northern part the air masses were older (up to 13 days) and had unclear ozone formation or destruction potential. Based on analysis of forward trajectories, several receptor regions were identified. In addition to predominantly westward transport, we found evidence for efficient transport (within 10 days) to the Pacific and North America, particularly during June and September, and also of cross-tropopause exchange, which was strongest during June and July. Westward transport to Africa and further to the Mediterranean was the main pathway during July.
The Earth’s surface condition we find today is a result of long exposure to metabolism of life forms. Particularly, molecular oxygen in the atmosphere is a feature which developed over time. The first substantial and lasting rise of atmospheric oxygen level happened ≈ 2.5 Ga ago, but localities are reported where transiently elevated oxygen levels appeared before this time-point. To trace the timing and circumstances of the earliest availability of free oxygen in the atmosphere is important to understand the habitats of early microbial life forms on Earth.
This thesis focuses to obtain information of oxygen levels and the related atmospheric cycling of metals in sediments of the 3.5 to 3.2 Ga Barberton Greenstone Belt. First, as iron was a ubiquitous constituent of Archean seawater, I investigated its isotopic composition in minerals of chemical sediments. Hereby, I tried to resolve the changes within the water basin on small scale sedimentary sequence cycles. Second, I focused on the minor constituents of Archean seawater. The Re-Os geochronologic system and the abundance patterns of the platinum-group elements were chosen to integrate information of oxygen promoted weathering of a large source area. To integrate information of a large time interval, the isotopes of uranium were investigated over a large stratigraphic section.
The two key findings of this thesis are:
• Quantitative oxidation of ferrous iron in surface layers of Paleoarchean seawater occurred during the onset and termination of hydrothermal FeIIaq delivery into shallow waters.
• Paleoarchean sedimentary successions of the Barberton Greenstone Belt lack any evidence of transient basin-scale oxygenation.
The Manzimnyama Iron Formation (IF, Fig Tree Group, Barberton Greenstone Belt, South Africa) has been deciphered to exist of cyclic stacks of lithostratigraphic units with varying amounts of iron oxide and carbonate minerals. In-situ femtosecond-Laser-Ablation ICP-MS iron isotope measurements showed that the majority of siderite (γ56Fe ≈ −0.5 ‰) precipitated directly from seawater of γ56Fe ≈ 0 ‰. Ferric iron from the surface layers is preserved in ≤ 1μ m hematite and in magnetite that has been grown within the consolidated sediment. During FeIIaq events, fine-grained hematite (γ56Fe ≈ 2.2 ‰) and magnetite (γ56Fe 0.5 to 0.8 ‰) indicate oxygen levels in surface waters of lower than 0.0002 μM. Upon onset and termination of iron oxide abundance, magnetite with γ56Fe ≈ 0 ‰ indicates that low concentrations of FeIIaq in surface waters were oxidized quantitatively. These observations demonstrate the existence of iron oxidation in Paleoarchean surface waters independent of FeIIaq concentration. This is the first investigation of Paleoarchean IF showing that lithostratigraphic cyclicity can be traced in iron isotopic composition of oxide minerals.
ID-ICP-MS measurement of Re, Ir, Ru, Pt and Pd, trace element (SF-ICP-MS) and ID-MCICP- MS uranium isotope determination have been applied to carbonaceous shale of the Mapepe Fm. (Fig Tree Group) after inverse Aqua Regia leaching and bulk digestion. The sediments reveal a silicified fraction which exhibits a seawater REE signature and a mixture of detrital and meteoritic PGE. Neither enrichment of the redox-sensitive elements Re or Mo nor fractionated uranium isotopes have been found on a stratigraphic interval of several hundred meters. The non-silica fraction shows no depletion of Re which indicates that the detrital material had no contact to oxidizing fluids. ID-TIMS measurements of Re and Os after the CrO3-SO4 Carius Tube method of two sample intervals showed that the Re-Os isotopic systems of the non-silica fractions are identical to two komatiite occurrences. Weltevreden Fm. and Komati Fm. rocks were uplifted, eroded and transported to the deep part of the sedimentary basin without any change to the Re-Os system. Negative fractionated uranium isotopes (γ238U = −0.41 ± 0.01 ‰) associated with detrital Ba-Cr-U occurrences suggest the existence of distal redox-processes that involve uranium species. This study demonstrates that over the time of exposure and deposition of the Mapepe Fm. sedimentation, free oxygen was not available for weathering in the catchment area.
The Fire Modeling Intercomparison Project (FireMIP), phase 1: experimental and analytical protocols
(2016)
The important role of fire in regulating vegetation community composition and contributions to emissions of greenhouse gases and aerosols make it a critical component of dynamic global vegetation models and Earth system models. Over two decades of development, a wide variety of model structures and mechanisms have been designed and incorporated into global fire models, which have been linked to different vegetation models. However, there has not yet been a systematic examination of how these different strategies contribute to model performance. Here we describe the structure of the first phase of the Fire Model Intercomparison Project (FireMIP), which for the first time seeks to systematically compare a number of models. By combining a standardized set of input data and model experiments with a rigorous comparison of model outputs to each other and to observations, we will improve the understanding of what drives vegetation fire, how it can best be simulated, and what new or improved observational data could allow better constraints on model behavior. Here we introduce the fire models used in the first phase of FireMIP, the simulation protocols applied, and the benchmarking system used to evaluate the models. The works published in this journal are distributed under the Creative Commons Attribution 3.0 License. This license does not affect the Crown copy-right work, which is re-usable under the Open Government Licence (OGL). The Creative Commons Attribution 3.0 License and the OGL are interoperable and do not conflict with, reduce, or limit each other.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations as well as hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a one-year period and full seasonal cycle (March 2014 - February 2015). The presented measurements provide a climatology of CCN properties for a characteristic central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The observed mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), elevated values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter κ exhibits remarkably little temporal variability: no pronounced diurnal cycles, weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modelling purposes, we compare different approaches of predicting CCN number concentration and present a novel parameterization, which allows accurate CCN predictions based on a small set of input data.
Size-resolved long-term measurements of atmospheric aerosol and cloud condensation nuclei (CCN) concentrations and hygroscopicity were conducted at the remote Amazon Tall Tower Observatory (ATTO) in the central Amazon Basin over a 1-year period and full seasonal cycle (March 2014–February 2015). The measurements provide a climatology of CCN properties characteristic of a remote central Amazonian rain forest site.
The CCN measurements were continuously cycled through 10 levels of supersaturation (S = 0.11 to 1.10 %) and span the aerosol particle size range from 20 to 245 nm. The mean critical diameters of CCN activation range from 43 nm at S = 1.10 % to 172 nm at S = 0.11 %. The particle hygroscopicity exhibits a pronounced size dependence with lower values for the Aitken mode (κAit = 0.14 ± 0.03), higher values for the accumulation mode (κAcc = 0.22 ± 0.05), and an overall mean value of κmean = 0.17 ± 0.06, consistent with high fractions of organic aerosol.
The hygroscopicity parameter, κ, exhibits remarkably little temporal variability: no pronounced diurnal cycles, only weak seasonal trends, and few short-term variations during long-range transport events. In contrast, the CCN number concentrations exhibit a pronounced seasonal cycle, tracking the pollution-related seasonality in total aerosol concentration. We find that the variability in the CCN concentrations in the central Amazon is mostly driven by aerosol particle number concentration and size distribution, while variations in aerosol hygroscopicity and chemical composition matter only during a few episodes.
For modeling purposes, we compare different approaches of predicting CCN number concentration and present a novel parametrization, which allows accurate CCN predictions based on a small set of input data.
We present an estimation of the uptake coefficient (γ) and yield of nitryl chloride (ClNO2) (f) for the heterogeneous processing of dinitrogen pentoxide (N2O5) using simultaneous measurements of particle and trace gas composition at a semi-rural, non-coastal, mountain site in the summer of 2011. The yield of ClNO2 varied between (0.035 ± 0.027) and (1.38 ± 0.60) with a campaign average of (0.49 ± 0.35). The large variability in f reflects the highly variable chloride content of particles at the site. Uptake coefficients were also highly variable with minimum, maximum and average γ values of 0.004, 0.11 and 0.028 ± 0.029, respectively, with no significant correlation with particle composition, but a weak dependence on relative humidity. The uptake coefficients obtained are compared to existing parameterizations based on laboratory datasets and with other values obtained by analysis of field data.
We present an estimation of the uptake coefficient (γ) and yield of nitryl chloride (ClNO2) (f) for the heterogeneous processing of dinitrogen pentoxide (N2O5) using simultaneous measurements of particle and trace gas composition at a semi-rural, non-coastal, mountain site in the summer of 2011. The yield of ClNO2 varied between (0.035 ± 0.027) and (1.38 ± 0.60) with a campaign average of (0.49 ± 0.35). The large variability in f reflects the highly variable chloride content of particles at the site. Uptake coefficients were also highly variable with minimum, maximum and average γ values of 0.004, 0.11 and 0.028 ± 0.029, respectively, with no significant correlation with particle composition, but a weak dependence on relative humidity. The uptake coefficients obtained are compared to existing parameterisations based on laboratory datasets and with other values obtained by analysis of field data.
Clouds and precipitation are essential climate variables. Because of their high spatial and temporal variability, their observation and modeling is difficult. In this thesis multiple observational data sources, including satellite data and station data are globally analyzed to understand the distribution and variability of clouds and precipitation, while a special focus is on the diurnal cycle of both variables. Substantial diurnal cycles of clouds and precipitation are observed in the tropics, with different properties over land and ocean. But also in Europe cloud diurnal cycles are observed in the summer season. Overall the maximum cloud cover and also the maximum precipitation is observed in the afternoon over land, and in the morning over ocean. The analyzed climate model simulations and the model-based reanalysis fail to simulate the observed diurnal cycles. Owing to their limited resolution, models can not fully resolve the processes responsible for the existence of diurnal cycles of clouds and precipitation.
The fractional release factor (FRF) gives information on the amount of a halocarbon that is released at some point in the stratosphere from its source form to the inorganic form, which can harm the ozone layer through catalytic reactions. The quantity is of major importance because it directly affects the calculation of the Ozone Depletion Potential (ODP). To apply FRF in this context, steady-state values are needed, thus representing a molecular property for a given atmospheric situation. In particular, these values should be independent of the tropospheric trends of the respective halogenated trace gases.
We analyzed the temporal evolution of FRF from ECHAM/MESSy Atmospheric Chemistry (EMAC) model simulations for several halocarbons and nitrous oxide between 1965–2011 on different mean age levels and found that the current formulation of FRF yields highly time-dependent values. We show that this is caused by the way that the tropospheric trend is handled in the current calculation method of FRF.
Taking into account chemical loss in the calculation of stratospheric mixing ratios reduces the time-dependence in correlations of different tracers. Therefore we implemented a loss term in the formulation of FRF and applied the parameterization of a "mean arrival time" to our data set.
We find that the time-dependence in FRF can almost be compensated by applying a new trend correction in the calculation of FRF. We suggest that this new method should be used to calculate time-independent FRF, which can then be used e.g. for the calculation of ODP
We present a compact and versatile cryofocusing– thermodesorption unit, which we developed for quantitative analysis of halogenated trace gases in ambient air. Possible applications include aircraft-based in situ measurements, in situ monitoring and laboratory operation for the analysis of flask samples. Analytes are trapped on adsorptive material cooled by a Stirling cooler to low temperatures (e.g. -80°C) and subsequently desorbed by rapid heating of the adsorptive material (e.g. 200°C). The set-up involves neither the exchange of adsorption tubes nor any further condensation or refocusing steps. No moving parts are used that would require vacuum insulation. This allows for a simple and robust design. Reliable operation is ensured by the Stirling cooler, which neither contains a liquid refrigerant nor requires refilling a cryogen. At the same time, it allows for significantly lower adsorption temperatures compared to commonly used Peltier elements. We use gas chromatography – mass spectrometry (GC–MS) for separation and detection of the preconcentrated analytes after splitless injection. A substance boiling point range of approximately -80 to +150°C and a substance mixing ratio range of less than 1 ppt (pmol mol−1)to more than 500 ppt in preconcentrated sample volumes of 0.1 to 10 L of ambient air is covered, depending on the application and its analytical demands. We present the instrumental design of the preconcentration unit and demonstrate capabilities and performance through the examination of analyte breakthrough during adsorption, repeatability of desorption and analyte residues in blank tests. Examples of application are taken from the analysis of flask samples collected at Mace Head Atmospheric Research Station in Ireland using our laboratory GC–MS instruments and by data obtained during a research flight with our in situ aircraft instrument GhOSTMS (Gas chromatograph for the Observation of Tracers – coupled with a Mass Spectrometer).
We present a compact and versatile cryofocusing–thermodesorption unit, which we developed for quantitative analysis of halogenated trace gases in ambient air. Possible applications include aircraft-based in situ measurements, in situ monitoring and laboratory operation for the analysis of flask samples. Analytes are trapped on adsorptive material cooled by a Stirling cooler to low temperatures (e.g. −80 °C) and subsequently desorbed by rapid heating of the adsorptive material (e.g. +200 °C). The set-up involves neither the exchange of adsorption tubes nor any further condensation or refocusing steps. No moving parts are used that would require vacuum insulation. This allows for a simple and robust design. Reliable operation is ensured by the Stirling cooler, which neither contains a liquid refrigerant nor requires refilling a cryogen. At the same time, it allows for significantly lower adsorption temperatures compared to commonly used Peltier elements. We use gas chromatography – mass spectrometry (GC–MS) for separation and detection of the preconcentrated analytes after splitless injection. A substance boiling point range of approximately −80 to +150 °C and a substance mixing ratio range of less than 1 ppt (pmol mol−1) to more than 500 ppt in preconcentrated sample volumes of 0.1 to 10 L of ambient air is covered, depending on the application and its analytical demands. We present the instrumental design of the preconcentration unit and demonstrate capabilities and performance through the examination of analyte breakthrough during adsorption, repeatability of desorption and analyte residues in blank tests. Examples of application are taken from the analysis of flask samples collected at Mace Head Atmospheric Research Station in Ireland using our laboratory GC–MS instruments and by data obtained during a research flight with our in situ aircraft instrument GhOST-MS (Gas chromatograph for the Observation of Tracers – coupled with a Mass Spectrometer).
We present the characterization and application of a new gas chromatography time-of-flight mass spectrometry instrument (GC-TOFMS) for the quantitative analysis of halocarbons in air samples. The setup comprises three fundamental enhancements compared to our earlier work (Hoker et al., 2015): (1) full automation, (2) a mass resolving power R = m/Δm of the TOFMS (Tofwerk AG, Switzerland) increased up to 4000 and (3) a fully accessible data format of the mass spectrometric data. Automation in combination with the accessible data allowed an in-depth characterization of the instrument. Mass accuracy was found to be approximately 5 ppm in mean after automatic recalibration of the mass axis in each measurement. A TOFMS configuration giving R = 3500 was chosen to provide an R-to-sensitivity ratio suitable for our purpose. Calculated detection limits are as low as a few femtograms by means of the accurate mass information. The precision for substance quantification was 0.15 % at the best for an individual measurement and in general mainly determined by the signal-to-noise ratio of the chromatographic peak. Detector non-linearity was found to be insignificant up to a mixing ratio of roughly 150 ppt at 0.5 L sampled volume. At higher concentrations, non-linearities of a few percent were observed (precision level: 0.2 %) but could be attributed to a potential source within the detection system. A straightforward correction for those non-linearities was applied in data processing, again by exploiting the accurate mass information. Based on the overall characterization results, the GC-TOFMS instrument was found to be very well suited for the task of quantitative halocarbon trace gas observation and a big step forward compared to scanning, quadrupole MS with low mass resolving power and a TOFMS technique reported to be non-linear and restricted by a small dynamical range.
Ziel der Arbeit war es, die Flugzeitmassenspektrometrie als neue Analysemethode für die instrumentelle Analytik halogenierter Spurengase in der Luft zu etablieren. Die grundle-gende Motivation dafür ist, dass anthropogene Emissionen vieler Vertreter dieser Sub-stanzklasse einen negativen Einfluss auf die Umwelt zeigen: in der Atmosphäre agieren die Substanzen bzw. ihre Abbauprodukte als Katalysatoren für den stratosphärischen Ozonab-bau und verstärken den Strahlungsantrieb der Erde durch Absorption elektromagnetischer Strahlung im sogenannten atmosphärischen Fenster. Um diese Effekte und deren Auswir-kung quantifizieren zu können, ist es notwendig, Konzentrationen und Trends der Substan-zen in der Atmosphäre zu überwachen. Nur so können Gegenmaßnahmen wie Produktions-reglementierungen geplant und bewertet werden. In Kombination mit inverser Modellie-rung können zudem Rückschlüsse über tatsächlich emittierten Mengen gezogen werden. Dies stellt den Anspruch an die Analytik, sehr geringe Mengen dieser Gase sehr präzise quantifizieren zu können, um auch schwache Trends zu erkennen. Zudem muss die Analy-semethode die Möglichkeit zu bieten, mit der wachsenden Anzahl bekannter und zu über-wachender Substanzen Schritt zu halten. Besonders für letzteren Aspekt bietet die Flug-zeitmassenspektrometrie einen entscheidenden Vorteil gegenüber der „konventionellen“ Methode, der Quadrupolmassenspektrometrie: sie zeichnet das gesamten Massenspektrum auf ohne dadurch an Empfindlichkeit einzubüßen. Um das atmosphärische Mischungsver-hältnis von Substanzen im Bereich von pmol mol−1 bis fmol mol−1 bestimmen zu können, muss das Quadrupolmassenspektrometer im Single Ion Monitoring Modus betrieben wer-den – so wird zwar eine hohe Sensitivität erreicht, es wird aber auch nur die Intensität eines bestimmten Masse zu Ladungsverhältnisses (kurz: Masse) zu einem Zeitpunkt aufgezeich-net. Ein Flugzeitmassenspektrometer hingegen extrahiert Ionen mit einer Frequenz im Ki-loherzbereich und zeichnet für jede Extraktion das vollständige Flugzeitspektrum und da-mit Massenspektrum auf.
Aufgabe dieser Arbeit war es, ein Flugzeitmassenspektrometer mit vorgeschalteter Pro-benanreicherungseinheit sowie Gaschromatograph zur Trennung des Subtanzgemisches vor der Detektion aufzubauen und Werkzeuge zur Datenauswertung zu entwickeln. Um einen zukünftigen Feldeinsatz vorzubereiten, sollte der Aufbau möglichst kompakt, mobil und vollständig automatisiert sein. Anschließend sollte Empfindlichkeit, Präzision und dynami-scher Messbereich geprüft, optimiert und die Anwendbarkeit zur Analyse halogenierter Spurengase gezeigt werden. Die Ergebnisse aus der in der vorliegenden Arbeit präsentier-ten Geräteentwicklung finden sich in drei Publikationen wieder, welche in thematischer Reihenfolge die Probenanreicherung (Obersteiner et al., 2016b), den Vergleich von Quadrupol- und Flugzeitmassenspektrometrie (Hoker et al., 2015) sowie Eigenschaften und Anwendung des neuen Aufbaus (Obersteiner et al., 2016a) behandeln. Mit den genannten Aufsätzen ist die Arbeitsgruppe Engel weltweit die erste, welche hochpräzise Analytik ha-logenierter Spurengase routinemäßig mittels Flugzeitmassenspektrometrie durchführt. Der nächste Schritt ist der Übergang von der Laboranwendung zur Feldmessung, z.B. in Form von bodenbasierter in situ Analyse troposphärischer Luftmassen am Taunus Observatorium auf dem Kleinen Feldberg. Da es bisher keine Messstation für die hier beschriebene analy-tische Fragestellung in Deutschland gibt, könnte eine deutliche Verbesserung der Überwa-chung halogenierter Treibhausgase und ozonzerstörender Substanzen in Europa erzielt wer-den. Weiterhin wäre eine Flugzeugapplikation in Zukunft denkbar, welche neben der durch das Flugzeitmassenspektrometer abgedeckten Substanzbandbreite auch von dessen hoher möglicher Spektrenrate profitieren könnte. In Kombination mit Hochgeschwindigkeitsgas-chromatographie könnte eine bisher unerreichte Zeitauflösung der Beprobung der Atmo-sphäre mittels Gaschromatographie-Massenspektrometrie erzielt werden.
Influence of sea surface roughness length parameterization on Mistral and Tramontane simulations
(2016)
The Mistral and Tramontane are mesoscale winds in southern France and above the Western Mediterranean Sea. They are phenomena well suited for studying channeling effects as well as atmosphere–land/ocean processes. This sensitivity study deals with the influence of the sea surface roughness length parameterizations on simulated Mistral and Tramontane wind speed and wind direction. Several simulations with the regional climate model COSMO-CLM were performed for the year 2005 with varying values for the Charnock parameter α. Above the western Mediterranean area, the simulated wind speed and wind direction pattern on Mistral days changes depending on the parameterization used. Higher values of α lead to lower simulated wind speeds. In areas, where the simulated wind speed does not change much, a counterclockwise rotation of the simulated wind direction is observed.
Mistral and tramontane wind speed and wind direction patterns in regional climate simulations
(2016)
The Mistral and Tramontane are important wind phenomena that occur over southern France and the northwestern Mediterranean Sea. Both winds travel through constricting valleys before flowing out towards the Mediterranean Sea. The Mistral and Tramontane are thus interesting phenomena, and represent an opportunity to study channeling effects, as well as the interactions between the atmosphere and land/ocean surfaces. This study investigates Mistral and Tramontane simulations using five regional climate models with grid spacing of about 50 km and smaller. All simulations are driven by ERA-Interim reanalysis data. Spatial patterns of surface wind, as well as wind development and error propagation along the wind tracks from inland France to offshore during Mistral and Tramontane events, are presented and discussed. To disentangle the results from large-scale error sources in Mistral and Tramontane simulations, only days with well simulated large-scale sea level pressure field patterns are evaluated. Comparisons with the observations show that the large-scale pressure patterns are well simulated by the considered models, but the orographic modifications to the wind systems are not well simulated by the coarse-grid simulations (with a grid spacing of about 50 km), and are reproduced slightly better by the higher resolution simulations. On days with Mistral and/or Tramontane events, most simulations underestimate (by 13 % on average) the wind speed over the Mediterranean Sea. This effect is strongest at the lateral borders of the main flow—the flow width is underestimated. All simulations of this study show a clockwise wind direction bias over the sea during Mistral and Tramontane events. Simulations with smaller grid spacing show smaller biases than their coarse-grid counterparts.
Cloud microphysical processes involving the ice phase in tropospheric clouds are among the major uncertainties in cloud formation, weather, and general circulation models. The detection of aerosol particles, liquid droplets, and ice crystals, especially in the small cloud particle-size range below 50 μm, remains challenging in mixed phase, often unstable environments. The Cloud Aerosol Spectrometer with Polarization (CASPOL) is an airborne instrument that has the ability to detect such small cloud particles and measure the variability in polarization state of their backscattered light. Here we operate the versatile Cosmics Leaving OUtdoor Droplets (CLOUD) chamber facility at the European Organization for Nuclear Research (CERN) to produce controlled mixed phase and other clouds by adiabatic expansions in an ultraclean environment, and use the CASPOL to discriminate between different aerosols, water, and ice particles. In this paper, optical property measurements of mixed-phase clouds and viscous secondary organic aerosol (SOA) are presented. We report observations of significant liquid–viscous SOA particle polarization transitions under dry conditions using CASPOL. Cluster analysis techniques were subsequently used to classify different types of particles according to their polarization ratios during phase transition. A classification map is presented for water droplets, organic aerosol (e.g., SOA and oxalic acid), crystalline substances such as ammonium sulfate, and volcanic ash. Finally, we discuss the benefits and limitations of this classification approach for atmospherically relevant concentrations and mixtures with respect to the CLOUD 8–9 campaigns and its potential contribution to tropical troposphere layer analysis.
Recent analysis of long-term balloon-borne measurements of Antarctic stratospheric condensation nuclei (CN) and temperature combined with global model calculations showed the wide extent of a mid stratospheric layer of new particles. Here the nucleation model SAWNUC is used to derive Antarctic stratospheric gaseous sulfuric acid profiles and to investigate the nucleation process of this CN layer. The sulfuric acid profiles were derived for an altitude range of 18-32 km between July and October by simulating air parcel trajectories that descend inside the polar vortex and calculating the sulfuric acid amount that reproduces the observations. The derived sulfuric acid concentrations (volume mixing ratios) are of the order of magnitude of 104 cm-3 (10-14) in July. In the following months the concentrations increase to about 107 cm-3 (10 exp -11) in October. They depend strongly on the temperature because a given temperature leaves only a small sulfuric acid range to reproduce the observed magnitude of CN. Ion-induced nucleation occurs. However, while it dominates nucleation at higher temperatures it has no significant influence on the nucleation rates at lower temperatures. Preexisting particles significantly reduce nucleation at sulfuric acid mixing ratios below 1 ppt. First estimates of sulfuric acid production rates range from 0.5 to about 500 molecules cm-3 s-1. A production mechanism for gaseous sulfuric acid during the Antarctic winter seems to be necessary to fully explain the observations. The derived sulfuric acid profiles compare well with midlatitude and Arctic sulfuric acid concentrations.
The transport of air masses originating from the Asian monsoon anticyclone into the extratropical upper troposphere and lower stratosphere (Ex-UTLS) above potential temperatures Θ = 380 K was identified during the HALO aircraft mission TACTS in August and September 2012. In situ measurements of CO, O3 and N2O during TACTS flight 2 on 30 August 2012 show the irreversible mixing of aged stratospheric air masses with younger (recently transported from the troposphere) ones within the Ex-UTLS. Backward trajectories calculated with the trajectory module of CLaMS indicate that these tropospherically affected air masses originate from the Asian monsoon anticyclone. These air masses are subsequently transported above potential temperatures Θ = 380 K from the monsoon circulation region into the Ex-UTLS, where they subsequently mix with stratospheric air masses. The overall trace gas distribution measured during TACTS shows that this transport pathway had affected the chemical composition of the Ex-UTLS during boreal summer and autumn 2012. This leads to an intensification of the tropospheric influence on the extratropical lower stratosphere with PV > 8 pvu within 3 weeks during the TACTS mission. During the same time period a weakening of the tropospheric influence on the lowermost stratosphere (LMS) is determined. The study shows that the transport of air masses originating from the Asian summer monsoon region within the lower stratosphere affects the change in the chemical composition of the Ex-UTLS over Europe and thus contributes to the flushing of the LMS during summer 2012.
Evaluation of radiation components in a global freshwater model with station-based observations
(2016)
In many hydrological models, the amount of evapotranspired water is calculated using the potential evapotranspiration (PET) approach. The main driver of several PET approaches is net radiation, whose downward components are usually obtained from meteorological input data, whereas the upward components are calculated by the model itself. Thus, uncertainties can be large due to both the input data and model assumptions. In this study, we compare the radiation components of the WaterGAP Global Hydrology Model, driven by two meteorological input datasets and two radiation setups from ERA-Interim reanalysis. We assess the performance with respect to monthly observations provided by the Baseline Surface Radiation Network (BSRN) and the Global Energy Balance Archive (GEBA). The assessment is done for the global land area and specifically for energy/water limited regions. The results indicate that there is no optimal radiation input throughout the model variants, but standard meteorological input datasets perform better than those directly obtained by ERA-Interim reanalysis for the key variable net radiation. The low number of observations for some radiation components, as well as the scale mismatch between station observations and 0.5° × 0.5° grid cell size, limits the assessment.
When assessing global water resources with hydrological models, it is essential to know about methodological uncertainties. The values of simulated water balance components may vary due to different spatial and temporal aggregations, reference periods, and applied climate forcings, as well as due to the consideration of human water use, or the lack thereof. We analyzed these variations over the period 1901–2010 by forcing the global hydrological model WaterGAP 2.2 (ISIMIP2a) with five state-of-the-art climate data sets, including a homogenized version of the concatenated WFD/WFDEI data set. Absolute values and temporal variations of global water balance components are strongly affected by the uncertainty in the climate forcing, and no temporal trends of the global water balance components are detected for the four homogeneous climate forcings considered (except for human water abstractions). The calibration of WaterGAP against observed long-term average river discharge Q significantly reduces the impact of climate forcing uncertainty on estimated Q and renewable water resources. For the homogeneous forcings, Q of the calibrated and non-calibrated regions of the globe varies by 1.6 and 18.5 %, respectively, for 1971–2000. On the continental scale, most differences for long-term average precipitation P and Q estimates occur in Africa and, due to snow undercatch of rain gauges, also in the data-rich continents Europe and North America. Variations of Q at the grid-cell scale are large, except in a few grid cells upstream and downstream of calibration stations, with an average variation of 37 and 74 % among the four homogeneous forcings in calibrated and non-calibrated regions, respectively. Considering only the forcings GSWP3 and WFDEI_hom, i.e., excluding the forcing without undercatch correction (PGFv2.1) and the one with a much lower shortwave downward radiation SWD than the others (WFD), Q variations are reduced to 16 and 31 % in calibrated and non-calibrated regions, respectively. These simulation results support the need for extended Q measurements and data sharing for better constraining global water balance assessments. Over the 20th century, the human footprint on natural water resources has become larger. For 11–18% of the global land area, the change of Q between 1941–1970 and 1971–2000 was driven more strongly by change of human water use including dam construction than by change in precipitation, while this was true for only 9–13 % of the land area from 1911–1940 to 1941–1970.
When assessing global water resources with hydrological models, it is essential to know the methodological uncertainties in the water resources estimates. The study presented here quantifies effects of the uncertainty in the spatial and temporal patterns of meteorological variables on water balance components at the global, continental and grid cell scale by forcing the global hydrological model WaterGAP 2.2 (ISI-MIP 2.1) with five state-of-the-art climate forcing input data-sets. While global precipitation over land during 1971–2000 varies between 103 500 and 111 000 km3 yr−1, global river discharge varies between 39 200 and 42 200 km3 yr−1. Temporal trends of global wa- ter balance components are strongly affected by the uncertainty in the climate forcing (except human water abstractions), and there is a need for temporal homogenization of climate forcings (in particular WFD/WFDEI). On about 10–20 % of the global land area, change of river discharge between two consecutive 30 year periods was driven more strongly by changes of human water use including dam construction than by changes in precipitation. This number increases towards the end of the 20th century due to intensified human water use and dam construction. The calibration approach of WaterGAP against observed long-term average river discharge reduces the impact of climate forcing uncertainty on estimated river discharge significantly. Different homgeneous climate forcings lead to a variation of Q of only 1.6 % for the 54 % of global land area that are constrained by discharge observations, while estimated renewable water resources in the remaining uncalibrated regions vary by 18.5 %. Uncertainties are especially high in Southeast Asia where Global Runoff Data Centre (GRDC) data availability is very sparse. By sharing already available discharge data, or installing new streamflow gauging stations in such regions, water balance uncertainties could be reduced which would lead to an improved assessment of the world’s water resources.
This thesis addresses the reconstruction of the topographic evolution and the climate dynamics of the Early Cenozoic North American Cordillera through integrated geochronology, sedimentology, stable isotope, and clumped isotope thermometry studies. It encompasses the scientific disciplines of geochemistry, tectonics, and Earth surface processes.
AirCore-HR: a high resolution column sampling to enhance the vertical description of CH₄ and CO₂
(2016)
An original and innovative sampling system called AirCore was presented by NOAA in 2010 (Karion et al., 2010). It consists of a long (> 100 m) and narrow (< 1 cm) stainless steel tube that can retain a profile of atmospheric air. The captured air sample has then to be analyzed with a gas analyzer for trace mole fraction. In this study, we introduce a new AirCore aiming at improved resolution along the vertical with the objectives to: (i) better capture the vertical distribution of CO2 and CH4, (ii) provide a tool to compare AirCores and validate the estimated vertical resolution achieved by AirCores. This AirCore-HR (high resolution) consists of a 300 m tube, combining 200 m of 1/8 in. (3.175 mm) tube and a 100 m of 1/4 in. (6.35 mm) tube. This new configuration allows to achieve a vertical resolution of 300 m up to 15 km and better than 500 m up to 22 km (if analysis of the retained sample is performed within 3 hours). The AirCore-HR was flown for the first time during the annual StratoScience campaign from CNES in August 2014 from Timmins (Ontario, Canada). High-resolution vertical profiles of CO2 and CH4 up to 25 km were successfully retrieved. These profiles revealed well defined transport structures in the troposphere (also seen in CAMS-ECMWF high resolution forecasts of CO2 and CH4 profiles) and captured the decrease of CO2 and CH4 in the stratosphere. The multi-instruments gondola from the flight carried two other low-resolution AirCore-GUF that allowed to perform direct comparisons and study the underlying processing method used to convert the sample of air to greenhouse gases vertical profiles. In particular, degrading the AirCore-HR derived profiles to the low resolution of AirCore-GUF yields an excellent match between both sets of CH4 profiles, and shows a good consistency between vertical structures of CO2 and CH4. These results fully validate the theoretical vertical resolution achievable by AirCores. Finally, the uncertainties associated with the measurements are assessed, yielding an average uncertainty below 3 ppb for CH4 and 0.25 ppm for CO2 with the major source of uncertainty coming from the potential loss of air sample on the ground and the choice of the starting and ending point of the collected air sample inside the tube. In an ideal case where the sample would be fully retained, it would be possible to know precisely the pressure at which air was sampled last and thus to improve the overall uncertainty to about 0.1 ppm for CO2 and 2 ppb for CH4.
Reconstructions of biomass burning from sediment-charcoal records to improve data–model comparisons
(2016)
The location, timing, spatial extent, and frequency of wildfires are changing rapidly in many parts of the world, producing substantial impacts on ecosystems, people, and potentially climate. Paleofire records based on charcoal accumulation in sediments enable modern changes in biomass burning to be considered in their long-term context. Paleofire records also provide insights into the causes and impacts of past wildfires and emissions when analyzed in conjunction with other paleoenvironmental data and with fire models. Here we present new 1000-year and 22 000-year trends and gridded biomass burning reconstructions based on the Global Charcoal Database version 3 (GCDv3), which includes 736 charcoal records (57 more than in version 2). The new gridded reconstructions reveal the spatial patterns underlying the temporal trends in the data, allowing insights into likely controls on biomass burning at regional to global scales. In the most recent few decades, biomass burning has sharply increased in both hemispheres but especially in the north, where charcoal fluxes are now higher than at any other time during the past 22 000 years. We also discuss methodological issues relevant to data–model comparisons and identify areas for future research. Spatially gridded versions of the global data set from GCDv3 are provided to facilitate comparison with and validation of global fire simulations.
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.
New particle formation driven by acid–base chemistry was initiated in the CLOUD chamber at CERN by introducing atmospherically relevant levels of gas-phase sulfuric acid and dimethylamine (DMA). Ammonia was also present in the chamber as a gas-phase contaminant from earlier experiments. The composition of particles with volume median diameters (VMDs) as small as 10 nm was measured by the Thermal Desorption Chemical Ionization Mass Spectrometer (TDCIMS). Particulate ammonium-to-dimethylaminium ratios were higher than the gas-phase ammonia-to-DMA ratios, suggesting preferential uptake of ammonia over DMA for the collected 10–30 nm VMD particles. This behavior is not consistent with present nanoparticle physicochemical models, which predict a higher dimethylaminium fraction when NH3 and DMA are present at similar gas-phase concentrations. Despite the presence in the gas phase of at least 100 times higher base concentrations than sulfuric acid, the recently formed particles always had measured base : acid ratios lower than 1 : 1. The lowest base fractions were found in particles below 15 nm VMD, with a strong size-dependent composition gradient. The reasons for the very acidic composition remain uncertain, but a plausible explanation is that the particles did not reach thermodynamic equilibrium with respect to the bases due to rapid heterogeneous conversion of SO2 to sulfate. These results indicate that sulfuric acid does not require stabilization by ammonium or dimethylaminium as acid–base pairs in particles as small as 10 nm.
New particle formation driven by acid-base chemistry was initiated in the CLOUD chamber at CERN by introducing atmospherically relevant levels of gas phase sulfuric acid and dimethylamine (DMA). Ammonia was also present in the chamber as a gas-phase contaminant from earlier experiments. The composition of particles with volume median diameters (VMDs) as small as 10 nm was measured by the Thermal Desorption Chemical Ionization Mass Spectrometer (TDCIMS). Particulate ammonium-to-dimethylaminium ratios were higher than the gas phase ammonia-to-DMA ratios, suggesting preferential uptake of ammonia over DMA for the collected 10-30 nm VMD particles. This behavior is not consistent with present nanoparticle physico-chemical models, which predict a higher dimethylaminium fraction when NH3 and DMA are present at similar gas phase concentrations. Despite the presence in the gas phase of at least 100 times higher base concentrations than sulfuric acid, the recently formed particles always had measured base:acid ratios lower than 1:1. The lowest base fractions were found in particles below 15 nm VMD, with a strong size-dependent composition gradient that suggests a change to a mixed-phase state as the particles grew beyond this size. The reasons for the very acidic composition remain uncertain, but a possible explanation is that the particles did not reach thermodynamic equilibrium with respect to the bases due to rapid heterogeneous conversion of SO2 to sulfate. These results indicate that sulfuric acid does not require stabilization by ammonium or dimethylaminium as acid-base pairs in particles as small as 10 nm.
The abiotic and biotic drivers of rapid diversification in Andean bellflowers (Campanulaceae)
(2016)
The tropical Andes of South America, the world's richest biodiversity hotspot, are home to many rapid radiations. While geological, climatic, and ecological processes collectively explain such radiations, their relative contributions are seldom examined within a single clade. We explore the contribution of these factors by applying a series of diversification models that incorporate mountain building, climate change, and trait evolution to the first dated phylogeny of Andean bellflowers (Campanulaceae: Lobelioideae). Our framework is novel for its direct incorporation of geological data on Andean uplift into a macroevolutionary model. We show that speciation and extinction are differentially influenced by abiotic factors: speciation rates rose concurrently with Andean elevation, while extinction rates decreased during global cooling. Pollination syndrome and fruit type, both biotic traits known to facilitate mutualisms, played an additional role in driving diversification. These abiotic and biotic factors resulted in one of the fastest radiations reported to date: the centropogonids, whose 550 species arose in the last 5 million yr. Our study represents a significant advance in our understanding of plant evolution in Andean cloud forests. It further highlights the power of combining phylogenetic and Earth science models to explore the interplay of geology, climate, and ecology in generating the world's biodiversity.
The exact mechanisms for new particle formation (NPF) under different boundary layer conditions are not known yet. One important question is whether amines and sulfuric acid lead to efficient NPF in the atmosphere. Furthermore, it is not clear to what extent highly oxidized organic molecules (HOMs) are involved in NPF. We conducted field measurements at a rural site in central Germany in the proximity of three larger dairy farms to investigate whether there is a connection between NPF and the presence of amines and/or ammonia due to the local emissions from the farms. Comprehensive measurements using a nitrate chemical ionization–atmospheric pressure interface time-of-flight (CI-APi-TOF) mass spectrometer, a proton-transfer-reaction mass spectrometer (PTR-MS), particle counters and differential mobility analyzers (DMAs), as well as measurements of trace gases and meteorological parameters, were performed. We demonstrate here that the nitrate CI-APi-TOF is suitable for sensitive measurements of sulfuric acid, amines, a nitrosamine, ammonia, iodic acid and HOMs. NPF was found to correlate with sulfuric acid, while an anti-correlation with RH, amines and ammonia is observed. The anti-correlation between NPF and amines could be due to the efficient uptake of these compounds by nucleating clusters and small particles. Much higher HOM dimer (C19/C20 compounds) concentrations during the night than during the day indicate that these HOMs do not efficiently self-nucleate as no nighttime NPF is observed. Observed iodic acid probably originates from an iodine-containing reservoir substance, but the iodine signals are very likely too low to have a significant effect on NPF.
The exact mechanisms for new particle formation (NPF) under different boundary layer conditions are not known yet. One important question is if amines and sulfuric acid lead to efficient NPF in the atmosphere. Furthermore, it is not clear to what extent highly oxidized organic molecules (HOM) are involved in NPF. We conducted field measurements at a rural site in central Germany in the proximity of three larger dairy farms to investigate if there is a connection between NPF and the presence of amines and/or ammonia due to the local emissions from the farms. Comprehensive measurements using a nitrate Chemical Ionization-Atmospheric Pressure interface-Time Of Flight (CI-APi-TOF) mass spectrometer, a Proton Transfer Reaction-Mass Spectrometer (PTR-MS), particle counters and Differential Mobility Analyzers (DMAs) as well as measurements of trace gases and meteorological parameters were performed. It is shown that the nitrate CI-APi-TOF is suitable for sensitive measurements of sulfuric acid, amines, a nitrosamine, ammonia, iodic acid and HOM. NPF was found to correlate with sulfuric acid, while an anti-correlation with RH, amines and ammonia is observed. The anti-correlation between NPF and amines could be due to the efficient uptake of these compounds by nucleating clusters and small particles. Much higher HOM dimer (C19/C20 compounds) concentrations during the night than during the day indicate that these HOM do not efficiently self-nucleate as no night-time NPF is observed. Observed iodic acid probably originates from an iodine-containing reservoir substance but the iodine signals are very likely too low to have a significant effect on NPF.
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.
Hygroscopicity of nanoparticles produced from homogeneous
nucleation in the CLOUD experiments
(2016)
Sulfuric acid, amines and oxidized organics have been found to be important compounds in the nucleation and initial growth of atmospheric particles. Because of the challenges involved in determining the chemical composition of objects with very small mass, however, the properties of the freshly nucleated particles and the detailed pathways of their formation processes are still not clear. In this study, we focus on a challenging size range, i.e., particles that have grown to diameters of 10 and 15 nm following nucleation, and measure their water uptake. Water uptake is useful information for indirectly obtaining chemical composition of aerosol particles. We use a nanometer-hygroscopicity tandem differential mobility analyzer (nano-HTDMA) at subsaturated conditions (ca. 90 % relative humidity at 293 K) to measure the hygroscopicity of particles during the seventh Cosmics Leaving OUtdoor Droplets (CLOUD7) campaign performed at CERN in 2012. In CLOUD7, the hygroscopicity of nucleated nanoparticles was measured in the presence of sulfuric acid, sulfuric acid–dimethylamine, and sulfuric acid–organics derived from α-pinene oxidation. The hygroscopicity parameter κ decreased with increasing particle size, indicating decreasing acidity of particles. No clear effect of the sulfuric acid concentration on the hygroscopicity of 10 nm particles produced from sulfuric acid and dimethylamine was observed, whereas the hygroscopicity of 15 nm particles sharply decreased with decreasing sulfuric acid concentrations. In particular, when the concentration of sulfuric acid was 5.1 × 106 molecules cm−3 in the gas phase, and the dimethylamine mixing ratio was 11.8 ppt, the measured κ of 15 nm particles was 0.31 ± 0.01: close to the value reported for dimethylaminium sulfate (DMAS) (κDMAS ∼ 0.28). Furthermore, the difference in κ between sulfuric acid and sulfuric acid–imethylamine experiments increased with increasing particle size. The κ values of particles in the presence of sulfuric acid and organics were much smaller than those of particles in the presence of sulfuric acid and dimethylamine. This suggests that the organics produced from α-pinene ozonolysis play a significant role in particle growth even at 10 nm sizes.
Development of the flash-heating method for measuring melting temperatures in the diamond anvil cell
(2016)
A new ‘laser flash-heating’ method has been developed for measuring melting temperatures above 2000 K in a diamond anvil cell at gigapascals of pressure. It overcomes the general difficulties in detecting an onset of melting in a diamond anvil cell. It also circumvents the notorious experimental difficulties associated with the long heating durations of the CW laser-heating and the short timescales in the pulsed laser-heating and shock-compression experiments.
In this method, the duration of heating a sample is tuned to avoid chemical reactions of the sample with the diamond anvils and the surrounding pressure medium, while maintaining the accuracy of the temperature measurements. The absence of chemical reactions is confirmed by the EDS technique. Melt detection is now unambiguous from the analysis of textures on the surface and in depth of the recovered samples using the SEM and FIBM techniques, respectively. Using this method, the following has been achieved.
1. The melting curve of hcp-Re has been measured to 48 GPa, 4200 K for the first time. It has a significantly steeper slope than those observed for other transition metals like W and Mo with bcc structures. Above 20 GPa, Re becomes the most refractory metal surpassing W.
2. The melting curve of bcc-Mo has been measured to 45 GPa, 3100 K. It agrees with previous melt-slopes approaching zero value with pressure as reported in the LHDAC experiments using ADXRD and visual observation techniques for inferring the onset of melting. Flash-heating experiments at pressures higher than 50 GPa are required to further corroborate the flat melt-slope and resolve the long standing controversy about melting of Mo.
3. The melting curve of bcc-Ta has been measured to 85 GPa, 4300 K. Unlike in previous experiments using ADXRD and visual observation as probes, it has been tightly bracketed with an unambiguous detection of the onset of melting, without any chemical reaction. The present melting curve cannot be reconciled with shock measurements and theoretical predictions, and the precision of measurements calls for a reevaluation of theoretical, shock compression, and other DAC approaches to determine melting at high pressures. A further analysis with TEM technique for investigating the structure of the heated portion below and above melting temperatures of Ta may benefit in resolving various phase transitions predicted to explain the vast discrepancies in the reported melt-slopes.
When extrapolated to one atmosphere pressure, all the measured flashmelting curves agree with the known melting points.
Melt segregation inside the earth consists of two different processes: 1) Generation of partially molten rock and 2) separation of melt, produced from partially molten rock, from the solid residual matrix. This thesis focuses on the later process. The 2 phase flow dynamics combines the study of flow dynamics of melt and matrix. Several studies have given the background theoretical frameworks for the flow dynamics of melt inside the earth. [McKenzie, 1984] summarizes the studies of [Ahern and Turcotte, 1979; Frank, 1968; Sleep, 1975] and gives a complete set of governing equations for the 2-phase flow problem.
[Bercovici et al., 2001] gives a general formulation considering the univariate system of equations related to matrix and melt flow which includes the interfacial surface force. The assumption of melt having negligible viscosity compare to the matrix has been abandoned. Therefore, based on these formulations, we have constructed our numerical model and thereafter a fortran code PERCOL2D to get an insight of melt percolation process through porous media. Additionally, we have used the Helmhotz decomposition, which splits a smooth and rapidly decaying vector field into an irrotational vector field and an incompressible vector field [Srámek, 2007], for matrix and fluid viscosity in order to lower the number of linearly independent variables to minimize the computational complications. The melt residing at inter-granular areas of lithosphere, forms an interconnected network even at low porosity. Therefore, being less dense than the matrix, melt moves up through porous media due to its buoyancy. Compaction of matrix, which occurs to compensate the melt separation, is considered in this thesis, where the effective bulk and shear viscosity of matrix are function of melt fraction. We have effective bulk viscosity of matrix as inversely proportional to melt fraction. Porosity dependence of effective bulk and shear viscosity leads to stronger melt focusing in highly porous region like mid ocean ridges [Katz, 2008] since the ratio of bulk and shear viscosity is smaller (< 10) than the constant viscosity case for the porous waves having non dimensional amplitude 5% or higher. Moreover, it is observed in [Richard et al., 2012] that the solitary wave formed in porosity dependent viscous matrix settings are steeper than the one formed in the constant matrix viscosity setting.
Firstly some 1D numerical experiments with PERCOL2D have been carried out using fixed and periodic boundary conditions for zero source term (i.e. no melting or no freezing) and negligible surface tension.
3 series of model setups with different initial conditions have been carried out varying the width, non-dimensional amplitude and the background porosity value of the initial input of porous wave.
A mathematical derivation for 1D solitary wave solution for the two phase flow through porosity dependent compacting media, is obtained in this thesis which is different than the study of [Barcilon and Lovera, 1989; Barcilon and Richter, 1986; Scott and Stevenson, 1984; Spiegelman, 1993a,b] as the effective viscosity of matrix is constant there.
Although [Simpson and Spiegelman, 2011] gives the solitary wave solutions in 1D, 2D and 3D considering the porosity dependent effective viscosity of the matrix, but using the small background porosity approximation, they neglect the background porosity (i.e φ0) and therefore the effect of variation of compaction lengths, which causes variation in the shape and dynamics of the solitary wave. Therefore, the study [This thesis, Richard et al., 2012] can be used for more general purpose. Solitary waves in varying viscous medium, are steeper (cf fig.5.1) compared to the one in constant viscous medium and their speed decreases as an inverse function of the background porosity. Additionally, this analytical solution is used in our code PERCOL2D and also in FDCON for numerical benchmarking (1D) of PERCOL2D.
The role of melt grain contiguity is considered in the revised viscosity formulation [Schmeling et al., 2012] based on elastic moduli theory of a fluid filled poro-elastic medium. This formulation is used in this thesis to produce a comparative dispersion relationship between speed of the wave and the non dimensional amplitude of porous wave, based on both the viscosity formulations (fig. 6.20) where one can see that the model based on [Bercovici et al., 2001] formulation, converges to the same dispersion relationship obtained from [Simpson and Spiegelman, 2011]. Whereas, the dispersion relationship using [Schmeling et al., 2012] formulations, shows time-dependent decrease of phase velocity with increasing amplitude and it is not yet clear that whether these solutions converge to steady state porosity waves before the porosity becomes 1.
Under certain conditions, secondary organic aerosol (SOA) particles can exist in the atmosphere in an amorphous solid or semi-solid state. To determine their relevance to processes such as ice nucleation or chemistry occurring within particles requires knowledge of the temperature and relative humidity (RH) range for SOA to exist in these states. In the Cosmics Leaving Outdoor Droplets (CLOUD) experiment at The European Organisation for Nuclear Research (CERN), we deployed a new in situ optical method to detect the viscous state of α-pinene SOA particles and measured their transition from the amorphous highly viscous state to states of lower viscosity. The method is based on the depolarising properties of laboratory-produced non-spherical SOA particles and their transformation to non-depolarising spherical particles at relative humidities near the deliquescence point. We found that particles formed and grown in the chamber developed an asymmetric shape through coagulation. A transition to a spherical shape was observed as the RH was increased to between 35 % at −10 °C and 80 % at −38 °C, confirming previous calculations of the viscosity-transition conditions. Consequently, α-pinene SOA particles exist in a viscous state over a wide range of ambient conditions, including the cirrus region of the free troposphere. This has implications for the physical, chemical, and ice-nucleation properties of SOA and SOA-coated particles in the atmosphere.
There are strong indications that particles containing secondary organic aerosol (SOA) exhibit amorphous solid or semi-solid phase states in the atmosphere. This may facilitate heterogeneous ice nucleation and thus influence cloud properties. However, experimental ice nucleation studies of biogenic SOA are scarce. Here, we investigated the ice nucleation ability of viscous SOA particles.
The SOA particles were produced from the ozone initiated oxidation of α-pinene in an aerosol chamber at temperatures in the range from −38 to −10 °C at 5–15 % relative humidity with respect to water to ensure their formation in a highly viscous phase state, i.e. semi-solid or glassy. The ice nucleation ability of SOA particles with different sizes was investigated with a new continuous flow diffusion chamber. For the first time, we observed heterogeneous ice nucleation of viscous α-pinene SOA for ice saturation ratios between 1.3 and 1.4 significantly below the homogeneous freezing limit. The maximum frozen fractions found at temperatures between −39.0 and −37.2 °C ranged from 6 to 20 % and did not depend on the particle surface area. Global modelling of monoterpene SOA particles suggests that viscous biogenic SOA particles are indeed present in regions where cirrus cloud formation takes place. Hence, they could make up an important contribution to the global ice nucleating particle budget.
The growth of aerosol due to the aqueous phase oxidation of sulfur dioxide by ozone was measured in laboratory-generated clouds created in the Cosmics Leaving OUtdoor Droplets (CLOUD) chamber at the European Organization for Nuclear Research (CERN). Experiments were performed at 10 and −10 °C, on acidic (sulfuric acid) and on partially to fully neutralised (ammonium sulfate) seed aerosol. Clouds were generated by performing an adiabatic expansion – pressurising the chamber to 220 hPa above atmospheric pressure, and then rapidly releasing the excess pressure, resulting in a cooling, condensation of water on the aerosol and a cloud lifetime of approximately 6 min. A model was developed to compare the observed aerosol growth with that predicted using oxidation rate constants previously measured in bulk solutions. The model captured the measured aerosol growth very well for experiments performed at 10 and −10 °C, indicating that, in contrast to some previous studies, the oxidation rates of SO2 in a dispersed aqueous system can be well represented by using accepted rate constants, based on bulk measurements. To the best of our knowledge, these are the first laboratory-based measurements of aqueous phase oxidation in a dispersed, super-cooled population of droplets. The measurements are therefore important in confirming that the extrapolation of currently accepted reaction rate constants to temperatures below 0 °C is correct.
The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated (nine chemical transport models and two chemistry–climate models) by simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993–2012). Except for three model simulations, all others were driven offline by (or nudged to) reanalysed meteorology. The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences due to emissions and transport processes. Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA's long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements – including high-altitude observations from the NASA Global Hawk platform.
The models generally capture the observed seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model–measurement correlation (r ≥ 0.7) at most sites. In a given model, the absolute model–measurement agreement at the surface is highly sensitive to the choice of emissions. Large inter-model differences are apparent when using the same emission inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve best agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2Br2). In general, the models reproduce observations of CHBr3 and CH2Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific well. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2Br2) most elevated over the tropical western Pacific during boreal winter. The models also indicate the Asian monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models.
We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2Br2 of 2.0 (1.2–2.5) ppt, ∼ 57 % larger than the best estimate from the most recent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. The transport-driven interannual variability in the annual mean bromine SGI is of the order of ±5 %, with SGI exhibiting a strong positive correlation with the El Niño–Southern Oscillation (ENSO) in the eastern Pacific. Overall, our results do not show systematic differences between models specific to the choice of reanalysis meteorology, rather clear differences are seen related to differences in the implementation of transport processes in the models.
The first concerted multi-model intercomparison of halogenated very short-lived substances (VSLS) has been performed, within the framework of the ongoing Atmospheric Tracer Transport Model Intercomparison Project (TransCom). Eleven global models or model variants participated, simulating the major natural bromine VSLS, bromoform (CHBr3) and dibromomethane (CH2Br2), over a 20-year period (1993-2012). The overarching goal of TransCom-VSLS was to provide a reconciled model estimate of the stratospheric source gas injection (SGI) of bromine from these gases, to constrain the current measurement-derived range, and to investigate inter-model differences
due to emissions and transport processes. Models ran with standardised idealised chemistry, to isolate differences due to transport, and we investigated the sensitivity of results to a range of VSLS emission inventories. Models were tested in their ability to reproduce the observed seasonal and spatial distribution of VSLS at the surface, using measurements from NOAA’s long-term global monitoring network, and in the tropical troposphere, using recent aircraft measurements - including high altitude observations from the NASA Global Hawk platform.
The models generally capture the seasonal cycle of surface CHBr3 and CH2Br2 well, with a strong model measurement correlation (r ≥0.7) and a low sensitivity to the choice of emission inventory, at most sites. In a given model, the absolute model-measurement agreement is highly sensitive to the choice of emissions and inter-model differences are also apparent, even when using the same inventory, highlighting the challenges faced in evaluating such inventories at the global scale. Across the ensemble, most consistency is found within the tropics where most of the models (8 out of 11) achieve optimal agreement to surface CHBr3 observations using the lowest of the three CHBr3 emission inventories tested (similarly, 8 out of 11 models for CH2Br2). In general, the models are able to reproduce well observations of CHBr3 and CH2Br2 obtained in the tropical tropopause layer (TTL) at various locations throughout the Pacific. Zonal variability in VSLS loading in the TTL is generally consistent among models, with CHBr3 (and to a lesser extent CH2Br2) most elevated over the tropical West Pacific during boreal winter. The models also indicate the Asian Monsoon during boreal summer to be an important pathway for VSLS reaching the stratosphere, though the strength of this signal varies considerably among models.
We derive an ensemble climatological mean estimate of the stratospheric bromine SGI from CHBr3 and CH2Br2 of 2.0 (1.2-2.5) ppt, ∼57% larger than the best estimate from the most re- cent World Meteorological Organization (WMO) Ozone Assessment Report. We find no evidence for a long-term, transport-driven trend in the stratospheric SGI of bromine over the simulation period. However, transport-driven inter-annual variability in the annual mean bromine SGI is of the order of a ±5%, with SGI exhibiting a strong positive correlation with ENSO in the East Pacific
A new global synthesis and biomization of long (> 40 kyr) pollen-data records is presented and used with simulations from the HadCM3 and FAMOUS climate models and the BIOME4 vegetation model to analyse the dynamics of the global terrestrial biosphere and carbon storage over the last glacial–interglacial cycle. Simulated biome distributions using BIOME4 driven by HadCM3 and FAMOUS at the global scale over time generally agree well with those inferred from pollen data. Global average areas of grassland and dry shrubland, desert, and tundra biomes show large-scale increases during the Last Glacial Maximum, between ca. 64 and 74 ka BP and cool substages of Marine Isotope Stage 5, at the expense of the tropical forest, warm-temperate forest, and temperate forest biomes. These changes are reflected in BIOME4 simulations of global net primary productivity, showing good agreement between the two models. Such changes are likely to affect terrestrial carbon storage, which in turn influences the stable carbon isotopic composition of seawater as terrestrial carbon is depleted in 13C.
Knowledge about mass discrimination effects in a chemical ionization mass spectrometer (CIMS) is crucial for quantifying, e.g., the recently discovered extremely low volatile organic compounds (ELVOCs) and other compounds for which no calibration standard exists so far. Here, we present a simple way of estimating mass discrimination effects of a nitrate-based chemical ionization atmospheric pressure interface time-of-flight (CI-APi-TOF) mass spectrometer. Characterization of the mass discrimination is achieved by adding different perfluorinated acids to the mass spectrometer in amounts sufficient to deplete the primary ions significantly. The relative transmission efficiency can then be determined by comparing the decrease of signals from the primary ions and the increase of signals from the perfluorinated acids at higher masses. This method is in use already for PTR-MS; however, its application to a CI-APi-TOF brings additional difficulties, namely clustering and fragmentation of the measured compounds, which can be treated with statistical analysis of the measured data, leading to selfconsistent results. We also compare this method to a transmission estimation obtained with a setup using an electrospray ion source, a high-resolution differential mobility analyzer and an electrometer, which estimates the transmission of the instrument without the CI source. Both methods give different transmission curves, indicating non-negligible mass discrimination effects of the CI source. The absolute transmission of the instrument without the CI source was estimated with the HR-DMA method to plateau between the m=z range of 127 and 568 Th at around 1.5 %; however, for the CI source included, the depletion method showed a steady increase in relative transmission efficiency from the m=z range of the primary ion (mainly at 62 Th) to around 550 Th by a factor of around 5. The main advantages of the depletion method are that the instrument is used in the same operation mode as during standard measurements and no knowledge of the absolute amount of the measured substance is necessary, which results in a simple setup.
Many natural minerals exist in the form of a solid solution. The systematic changes in structural and physical properties of oxide solid solutions are of geological importance and allow for wide applications. In order to understand the composition-structureproperty relations, substitutional solid solutions of CuxZn2−xTiO4, ZnxMg1−xTi2O5 and CuxMg1−xTi2O5 have been synthesised by mechanochemical activation assisted solid state synthesis. Self-propagating high-temperature synthesis has been employed to achieve the interstitial solid solutions of Ti5Si3Zx (Z refers to the element boron or oxygen).
The changes in the crystal structure and physical properties due to the formation of solid solutions are investigated by employing X-ray diffraction, neutron diffraction, Raman spectroscopy, low-temperature heat capacity, thermal expansion, scanning electron microscopy, UV-vis spectroscopy, plane-wave ultrasound spectroscopy and density functional theory calculations.
Biomass burning impacts vegetation dynamics, biogeochemical cycling, atmospheric chemistry, and climate, with sometimes deleterious socio-economic impacts. Under future climate projections it is often expected that the risk of wildfires will increase. Our ability to predict the magnitude and geographic pattern of future fire impacts rests on our ability to model fire regimes, either using well-founded empirical relationships or process-based models with good predictive skill. A large variety of models exist today and it is still unclear which type of model or degree of complexity is required to model fire adequately at regional to global scales. This is the central question underpinning the creation of the Fire Model Intercomparison Project - FireMIP, an international project to compare and evaluate existing global fire models against benchmark data sets for present-day and historical conditions. In this paper we summarise the current state-of-the-art in fire regime modelling and model evaluation, and outline what lessons may be learned from FireMIP.
Biomass burning impacts vegetation dynamics, biogeochemical cycling, atmospheric chemistry, and climate, with sometimes deleterious socio-economic impacts. Under future climate projections it is often expected that the risk of wildfires will increase. Our ability to predict the magnitude and geographic pattern of future fire impacts rests on our ability to model fire regimes, using either well-founded empirical relationships or process-based models with good predictive skill. While a large variety of models exist today, it is still unclear which type of model or degree of complexity is required to model fire adequately at regional to global scales. This is the central question underpinning the creation of the Fire Model Intercomparison Project (FireMIP), an international initiative to compare and evaluate existing global fire models against benchmark data sets for present-day and historical conditions. In this paper we review how fires have been represented in fire-enabled dynamic global vegetation models (DGVMs) and give an overview of the current state of the art in fire-regime modelling. We indicate which challenges still remain in global fire modelling and stress the need for a comprehensive model evaluation and outline what lessons may be learned from FireMIP.
This study aims to assess the skill of regional climate models (RCMs) at reproducing the climatology of Mediterranean cyclones. Seven RCMs are considered, five of which were also coupled with an oceanic model. All simulations were forced at the lateral boundaries by the ERA-Interim reanalysis for a common 20-year period (1989–2008). Six different cyclone tracking methods have been applied to all twelve RCM simulations and to the ERA-Interim reanalysis in order to assess the RCMs from the perspective of different cyclone definitions. All RCMs reproduce the main areas of high cyclone occurrence in the region south of the Alps, in the Adriatic, Ionian and Aegean Seas, as well as in the areas close to Cyprus and to Atlas mountains. The RCMs tend to underestimate intense cyclone occurrences over the Mediterranean Sea and reproduce 24–40 % of these systems, as identified in the reanalysis. The use of grid nudging in one of the RCMs is shown to be beneficial, reproducing about 60 % of the intense cyclones and keeping a better track of the seasonal cycle of intense cyclogenesis. Finally, the most intense cyclones tend to be similarly reproduced in coupled and uncoupled model simulations, suggesting that modeling atmosphere–ocean coupled processes has only a weak impact on the climatology and intensity of Mediterranean cyclones.
This paper is a contribution to the special issue on Med-CORDEX, an international coordinated initiative dedicated to the multi-component regional climate modelling (atmosphere, ocean, land surface, river) of the Mediterranean under the umbrella of HyMeX, CORDEX, and Med-CLIVAR and coordinated by Samuel Somot, Paolo Ruti, Erika Coppola, Gianmaria Sannino, Bodo Ahrens, and Gabriel Jordà.
MIPAS-Envisat is a satellite-borne sensor which measured vertical profiles of a wide range of trace gases from 2002 to 2012 using IR emission spectroscopy. We present geophysical validation of the MIPAS-Envisat operational retrieval (version 6.0) of N2O, CH4, CFC-12, and CFC-11 by the European Space Agency (ESA). The geophysical validation data are derived from measurements of samples collected by a cryogenic whole air sampler flown to altitudes of up to 34 km by means of large scientific balloons. In order to increase the number of coincidences between the satellite and the balloon observations, we applied a trajectory matching technique. The results are presented for different time periods due to a change in the spectroscopic resolution of MIPAS in early 2005. Retrieval results for N2O, CH4, and CFC-12 show partly good agreement for some altitude regions, which differs for the periods with different spectroscopic resolution. The more recent low spectroscopic resolution data above 20 km altitude show agreement with the combined uncertainties, while there is a tendency of the earlier high spectral resolution data set to underestimate these species above 25 km. The earlier high spectral resolution data show a significant overestimation of the mixing ratios for N2O, CH4, and CFC-12 below 20 km. These differences need to be considered when using these data. The CFC-11 results from the operation retrieval version 6.0 cannot be recommended for scientific studies due to a systematic overestimation of the CFC-11 mixing ratios at all altitudes.
We present the first 3-D model of seismic P and S velocities in the crust and uppermost mantle beneath the Gulf of Aqaba and surrounding areas based on the results of passive travel time tomography. The tomographic inversion was performed based on travel time data from ∼ 9000 regional earthquakes provided by the Egyptian National Seismological Network (ENSN), and this was complemented with data from the International Seismological Centre (ISC). The resulting P and S velocity patterns were generally consistent with each other at all depths. Beneath the northern part of the Red Sea, we observed a strong high-velocity anomaly with abrupt limits that coincide with the coastal lines. This finding may indicate the oceanic nature of the crust in the Red Sea, and it does not support the concept of gradual stretching of the continental crust. According to our results, in the middle and lower crust, the seismic anomalies beneath the Gulf of Aqaba seem to delineate a sinistral shift (∼ 100 km) in the opposite flanks of the fault zone, which is consistent with other estimates of the left-lateral displacement in the southern part of the Dead Sea Transform fault. However, no displacement structures were visible in the uppermost lithospheric mantle.
Binary nucleation of sulphuric acid-water particles is expected to be an important process in the free troposphere at low temperatures. SAWNUC (Sulphuric Acid Water Nucleation) is a model of binary nucleation that is based on laboratory measurements of the binding energies of sulphuric acid and water in charged and neutral clusters. Predictions of SAWNUC are compared for the first time comprehensively with experimental binary nucleation data from the CLOUD chamber at European Organization for Nuclear Research. The experimental measurements span a temperature range of 208–292 K, sulphuric acid concentrations from 1·106 to 1·109 cm−3, and distinguish between ion-induced and neutral nucleation. Good agreement, within a factor of 5, is found between the experimental and modeled formation rates for ion-induced nucleation at 278 K and below and for neutral nucleation at 208 and 223 K. Differences at warm temperatures are attributed to ammonia contamination which was indicated by the presence of ammonia-sulphuric acid clusters, detected by an Atmospheric Pressure Interface Time of Flight (APi-TOF) mass spectrometer. APi-TOF measurements of the sulphuric acid ion cluster distributions (math formula with i = 0, 1, ..., 10) show qualitative agreement with the SAWNUC ion cluster distributions. Remaining differences between the measured and modeled distributions are most likely due to fragmentation in the APi-TOF. The CLOUD results are in good agreement with previously measured cluster binding energies and show the SAWNUC model to be a good representation of ion-induced and neutral binary nucleation of sulphuric acid-water clusters in the middle and upper troposphere.
Profiles of CFC-11 (CCl3F) and CFC-12 (CCl2F2) of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) aboard the European satellite Envisat have been retrieved from versions MIPAS/4.61 to MIPAS/4.62 and MIPAS/5.02 to MIPAS/5.06 level-1b data using the scientific level-2 processor run by Karlsruhe Institute of Technology (KIT), Institute of Meteorology and Climate Research (IMK) and Consejo Superior de Investigaciones Científicas (CSIC), Instituto de Astrofísica de Andalucía (IAA). These profiles have been compared to measurements taken by the balloon-borne cryosampler, Mark IV (MkIV) and MIPAS-Balloon (MIPAS-B), the airborne MIPAS-STRatospheric aircraft (MIPAS-STR), the satellite-borne Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) and the High Resolution Dynamic Limb Sounder (HIRDLS), as well as the ground-based Halocarbon and other Atmospheric Trace Species (HATS) network for the reduced spectral resolution period (RR: January 2005–April 2012) of MIPAS. ACE-FTS, MkIV and HATS also provide measurements during the high spectral resolution period (full resolution, FR: July 2002–March 2004) and were used to validate MIPAS CFC-11 and CFC-12 products during that time, as well as profiles from the Improved Limb Atmospheric Spectrometer, ILAS-II. In general, we find that MIPAS shows slightly higher values for CFC-11 at the lower end of the profiles (below ∼ 15 km) and in a comparison of HATS ground-based data and MIPAS measurements at 3 km below the tropopause. Differences range from approximately 10 to 50 pptv ( ∼ 5–20 %) during the RR period. In general, differences are slightly smaller for the FR period. An indication of a slight high bias at the lower end of the profile exists for CFC-12 as well, but this bias is far less pronounced than for CFC-11 and is not as obvious in the relative differences between MIPAS and any of the comparison instruments. Differences at the lower end of the profile (below ∼ 15 km) and in the comparison of HATS and MIPAS measurements taken at 3 km below the tropopause mainly stay within 10–50 pptv (corresponding to ∼ 2–10 % for CFC-12) for the RR and the FR period. Between ∼ 15 and 30 km, most comparisons agree within 10–20 pptv (10–20 %), apart from ILAS-II, which shows large differences above ∼ 17 km. Overall, relative differences are usually smaller for CFC-12 than for CFC-11. For both species – CFC-11 and CFC-12 – we find that differences at the lower end of the profile tend to be larger at higher latitudes than in tropical and subtropical regions. In addition, MIPAS profiles have a maximum in their mixing ratio around the tropopause, which is most obvious in tropical mean profiles. Comparisons of the standard deviation in a quiescent atmosphere (polar summer) show that only the CFC-12 FR error budget can fully explain the observed variability, while for the other products (CFC-11 FR and RR and CFC-12 RR) only two-thirds to three-quarters can be explained. Investigations regarding the temporal stability show very small negative drifts in MIPAS CFC-11 measurements. These instrument drifts vary between ∼ 1 and 3 % decade−1. For CFC-12, the drifts are also negative and close to zero up to ∼ 30 km. Above that altitude, larger drifts of up to ∼ 50 % decade−1 appear which are negative up to ∼ 35 km and positive, but of a similar magnitude, above.
In this study we investigate the scaling of precipitation extremes with temperature in the Mediterranean region by assessing against observations the present day and future regional climate simulations performed in the frame of the HyMeX and MED-CORDEX programs. Over the 1979–2008 period, despite differences in quantitative precipitation simulation across the various models, the change in precipitation extremes with respect to temperature is robust and consistent. The spatial variability of the temperature–precipitation extremes relationship displays a hook shape across the Mediterranean, with negative slope at high temperatures and a slope following Clausius–Clapeyron (CC)-scaling at low temperatures. The temperature at which the slope of the temperature–precipitation extreme relation sharply changes (or temperature break), ranges from about 20 °C in the western Mediterranean to <10 °C in Greece. In addition, this slope is always negative in the arid regions of the Mediterranean. The scaling of the simulated precipitation extremes is insensitive to ocean–atmosphere coupling, while it depends very weakly on the resolution at high temperatures for short precipitation accumulation times. In future climate scenario simulations covering the 2070–2100 period, the temperature break shifts to higher temperatures by a value which is on average the mean regional temperature change due to global warming. The slope of the simulated future temperature–precipitation extremes relationship is close to CC-scaling at temperatures below the temperature break, while at high temperatures, the negative slope is close, but somewhat flatter or steeper, than in the current climate depending on the model. Overall, models predict more intense precipitation extremes in the future. Adjusting the temperature–precipitation extremes relationship in the present climate using the CC law and the temperature shift in the future allows the recovery of the temperature–precipitation extremes relationship in the future climate. This implies negligible regional changes of relative humidity in the future despite the large warming and drying over the Mediterranean. This suggests that the Mediterranean Sea is the primary source of moisture which counteracts the drying and warming impacts on relative humidity in parts of the Mediterranean region.
In this study, we aim to reconstruct a relevant and new database of monthly zonal mean distribution of carbon dioxide (CO2) at global scale extending from the upper-troposphere (UT) to stratosphere (S). This product can be used for model and satellite validation in the UT/S, as a prior for inversion modelling and mainly to analyse a plausible feature of the stratospherictropospheric exchange as well as the stratospheric circulation and its variability. To do so, we investigate the ability of a Lagrangian trajectory model guided by ERA-Interim reanalysis to construct the CO2 abundance in the UT/S. From 10 year backward trajectories and tropospheric observations of CO2, we reconstruct upper-tropospheric and stratospheric CO2 over the period 2000–2010. The inter-comparisons of the reconstructed CO2 with mid-latitude vertical profiles measured by balloon samples as well as quasi-horizontal air samples from ER-2 aircraft during SOLVE and CONTRAIL campaigns exhibit a remarkable agreement. That demonstrates the potential of Lagrangian model to reconstruct CO2 in the UT/S. The zonal mean distribution exhibits relatively large CO2 in the tropical stratosphere due to the seasonal variation of the tropical upwelling of Brewer-Dobson circulation. During winter and spring, the tropical pipe is relatively isolated but is less confined during summer and autumn so that high CO2 values are more readily transported out of the tropics to the mid- and high latitude stratosphere. The shape of the vertical profiles suggests that relatively high CO2 above 20 km altitude mainly enter the stratosphere through tropical upwelling. CO2 mixing ratio is relatively low in the polar and tropical regions above 25 km. On average the CO2 mixing ratio decreases with altitude by 6-8 ppmv from the UT to stratosphere (e.g. up to 35 km) and is nearly constant with altitude.
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.
Small-scale thermal upwellings under the northern East African Rift from S travel time tomography
(2016)
There is a long-standing debate over how many and what types of plumes underlie the East African Rift and whether they do or do not drive its extension and consequent magmatism and seismicity. Here we present a new tomographic study of relative teleseismic S and SKS residuals that expands the resolution from previous regional studies below the northern East African Rift to image structure from the surface to the base of the transition zone. The images reveal two low-velocity clusters, below Afar and west of the Main Ethiopian Rift, that extend throughout the upper mantle and comprise several smaller-scale (about 100 km diameter), low-velocity features. These structures support those of our recent P tomographic study below the region. The relative magnitude of S to P residuals is around 3.5, which is consistent with a predominantly thermal nature of the anomalies. The S and P velocity anomalies in the low-velocity clusters can be explained by similar excess temperatures in the range of 100–200°C, consistent with temperatures inferred from other seismic, geochemical, and petrological studies. Somewhat stronger VS anomalies below Afar than west of the Main Ethiopian Rift may include an expression of volatiles and/or melt in this region. These results, together with a comparison with previous larger-scale tomographic models, indicate that these structures are likely small-scale upwellings with mild excess temperatures, rising from a regional thermal boundary layer at the base of the upper mantle.
We report on HCFC-22 data acquired by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) in the reduced spectral resolution nominal observation mode. The data cover the period from January 2005 to April 2012 and the altitude range from the upper troposphere (above cloud top altitude) to about 50 km. The profile retrieval was performed by constrained nonlinear least squares fitting of modelled spectra to the measured limb spectral radiances. The spectral ν4-band at 816.5 ± 13 cm−1 was used for the retrieval. A Tikhonov-type smoothing constraint was applied to stabilise the retrieval. In the lower stratosphere, we find a global volume mixing ratio of HCFC-22 of about 185 pptv in January 2005. The rate of linear growth in the lower latitudes lower stratosphere was about 6 to 7 pptv year−1 in the period 2005–2012. The profiles obtained were compared with ACE-FTS satellite data v3.5, as well as with MkIV balloon profiles and cryosampler balloon measurements. Between 13 and 22 km, average agreement within −3 to +5 pptv (MIPAS – ACE) with ACE-FTS v3.5 profiles is demonstrated. Agreement with MkIV solar occultation balloon-borne measurements is within 10–20 pptv below 30 km and worse above, while in situ cryosampler balloon measurements are systematically lower over their full altitude range by 15–50 pptv below 24 km and less than 10 pptv above 28 km. MIPAS HCFC-22 time series below 10 km altitude are shown to agree mostly well to corresponding time series of near-surface abundances from the NOAA/ESRL and AGAGE networks, although a more pronounced seasonal cycle is obvious in the satellite data. This is attributed to tropopause altitude fluctuations and subsidence of polar winter stratospheric air into the troposphere. A parametric model consisting of constant, linear, quasi-biennial oscillation (QBO) and several sine and cosine terms with different periods has been fitted to the temporal variation of stratospheric HCFC-22 for all 10°-latitude/1-to-2-km-altitude bins. The relative linear variation was always positive, with relative increases of 40–70 % decade−1 in the tropics and global lower stratosphere, and up to 120 % decade−1 in the upper stratosphere of the northern polar region and the southern extratropical hemisphere. Asian HCFC-22 emissions have become the major source of global upper tropospheric HCFC-22. In the upper troposphere, monsoon air, rich in HCFC-22, is instantaneously mixed into the tropics. In the middle stratosphere, between 20 and 30 km, the observed trend is inconsistent with the trend at the surface (corrected for the age of stratospheric air), hinting at circulation changes. There exists a stronger positive trend in HCFC-22 in the Southern Hemisphere and a more muted positive trend in the Northern Hemisphere, implying a potential change in the stratospheric circulation over the observation period.
Convection-permitting climate model are promising tools for improved representation of extremes, but the number of regions for which these models have been evaluated are still rather limited to make robust conclusions. In addition, an integrated interpretation of near-surface characteristics (typically temperature and precipitation) together with cloud properties is limited. The objective of this paper is to comprehensively evaluate the performance of a ‘state-of-the-art’ regional convection-permitting climate model for a mid-latitude coastal region with little orographic forcing. For this purpose, an 11-year integration with the COSMO-CLM model at Convection-Permitting Scale (CPS) using a grid spacing of 2.8 km was compared with in-situ and satellite-based observations of precipitation, temperature, cloud properties and radiation (both at the surface and the top of the atmosphere). CPS clearly improves the representation of precipitation, in especially the diurnal cycle, intensity and spatial distribution of hourly precipitation. Improvements in the representation of temperature are less obvious. In fact the CPS integration overestimates both low and high temperature extremes. The underlying cause for the overestimation of high temperature extremes was attributed to deficiencies in the cloud properties: The modelled cloud fraction is only 46 % whereas a cloud fraction of 65 % was observed. Surprisingly, the effect of this deficiency was less pronounced at the radiation balance at the top of the atmosphere due to a compensating error, in particular an overestimation of the reflectivity of clouds when they are present. Overall, a better representation of convective precipitation and a very good representation of the daily cycle in different cloud types were demonstrated. However, to overcome remaining deficiencies, additional efforts are necessary to improve cloud characteristics in CPS. This will be a challenging task due to compensating deficiencies that currently exist in ‘state-of-the-art’ models, yielding a good representation of average climate conditions. In the light of using the CPS models to study climate change it is necessary that these deficiencies are addressed in future research.
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Traditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies over several decades. Recently some feldspar species were identified to be ice active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands, and the Sinai Peninsula). Additionally, 11 dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore, we investigated how representative surface-collected dust is for the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC setup to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 and 250 K. Dust particles were collected in parallel on filters for offline cold-stage ice nucleation experiments at 253–263 K. To help the interpretation of the ice nucleation experiments the mineralogical composition of the dusts was investigated. We find that a higher ice nucleation activity in a given sample at 253 K can be attributed to the K-feldspar content present in this sample, whereas at temperatures between 238 and 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content, in contrast, is associated with lower ice nucleation activity. This confirms the importance of feldspar above 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower temperatures as found by earlier studies for monomineral dusts. The airborne samples show on average a lower ice nucleation activity than the surface-collected ones. Furthermore, we find that under certain conditions milling can lead to a decrease in the ice nucleation ability of polymineral samples due to the different hardness and cleavage of individual mineral phases causing an increase of minerals with low ice nucleation ability in the atmospherically relevant size fraction. Comparison of our data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modeling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key, while other minerals are only of minor importance.
Desert dust is one of the most abundant ice nucleating particle types in the atmosphere. Tra ditionally, clay minerals were assumed to determine the ice nucleation ability of desert dust and constituted the focus of ice nucleation studies. Only recently some feldspar species were identified to be ice-active at much higher temperatures than clay minerals, redirecting studies to investigate the contribution of feldspar to ice nucleation on desert dust. However, so far no study has shown the atmospheric relevance of this mineral phase.
For this study four dust samples were collected after airborne transport in the troposphere from the Sahara to different locations (Crete, the Peloponnese, Canary Islands and the Sinai Peninsula). Additionally, eleven dust samples were collected from the surface from nine of the biggest deserts worldwide. The samples were used to study the ice nucleation behavior specific to different desert dusts. Furthermore we investigated how representative ice nucleation on surface-collected dust is for that in the atmosphere by comparing to the ice nucleation activity of the airborne samples. We used the IMCA-ZINC set-up to form droplets on single aerosol particles which were subsequently exposed to temperatures between 233 - 250 K. Dust particles were collected in parallel on filters for offline cold stage ice nucleation experiments at 253 - 263 K. To help the interpretation of the results from the ice nucleation experiments the mineralogical composition of the dusts was investigated.We found that a higher ice nucleation activity in a given sample can be attributed at 253 K to the K-feldspar content present in this sample whereas at temperatures between 238 - 245 K it is attributed to the sum of feldspar and quartz content present. A high clay content on the other hand is associated with a lower ice nucleation activity of a sample. This confirms the importance of feldspar at T > 250 K and the role of quartz and feldspars determining the ice nucleation activities at lower T as found by earlier studies for monomineral dust surrogates. Furthermore, we find that milling may lead to a decrease in the ice nucleation ability of polymineral samples due to a change in mineralogical composition in the atmospherically relevant size fraction arising from the different hardness and cleavage of individual mineral phases. Comparison of our comprehensive data set to an existing desert dust parameterization confirms its applicability for climate models. Our results suggest that for an improved prediction of the ice nucleation ability of desert dust in the atmosphere, the modelling of emission and atmospheric transport of the feldspar and quartz mineral phases would be key while other minerals are only of minor importance.
The oceans at the time of the Cenomanian–Turonian transition were abruptly perturbed by a period of bottom-water anoxia. This led to the brief but widespread deposition of black organic-rich shales, such as the Livello Bonarelli in the Umbria–Marche Basin (Italy). Despite intensive studies, the origin and exact timing of this event are still debated. In this study, we assess leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world by providing a 6 Myr long astronomically tuned timescale across the Cenomanian–Turonian boundary. We procure insights into the relationship between orbital forcing and the Late Cretaceous carbon cycle by deciphering the imprint of astronomical cycles on lithologic, physical properties, and stable isotope records, obtained from the Bottaccione, Contessa and Furlo sections in the Umbria–Marche Basin. The deposition of black shales and cherts, as well as the onset of oceanic anoxia, is related to maxima in the 405 kyr cycle of eccentricity-modulated precession. Correlation to radioisotopic ages from the Western Interior (USA) provides unprecedented age control for the studied Italian successions. The most likely tuned age for the base of the Livello Bonarelli is 94.17 ± 0.15 Ma (tuning 1); however, a 405 kyr older age cannot be excluded (tuning 2) due to uncertainties in stratigraphic correlation, radioisotopic dating, and orbital configuration. Our cyclostratigraphic framework suggests that the exact timing of major carbon cycle perturbations during the Cretaceous may be linked to increased variability in seasonality (i.e. a 405 kyr eccentricity maximum) after the prolonged avoidance of seasonal extremes (i.e. a 2.4 Myr eccentricity minimum). Volcanism is probably the ultimate driver of oceanic anoxia, but orbital periodicities determine the exact timing of carbon cycle perturbations in the Late Cretaceous. This unites two leading hypotheses about the inception of oceanic anoxia in the Late Cretaceous greenhouse world.
The organic rich Livello Bonarelli formed as a result of oxygen deficiency and carbonate dissolution in the oceans during the Cenomanian/Turonian (C/T) transition. During this Ocean Anoxic Event 2 (OAE2), a combination of factors caused increased productivity, incomplete decomposition of organic matter and widespread deposition of black shales. Although these sediments are extensively studied, the exact extent, cause, timing and duration of oceanic anoxia are debated (Sinton and Duncan, 1997; Mitchell et al., 2008). Contrasting causal mechanisms have been suggested, including stratification of the water column (Lanci et al., 2010) versus intensification of the hydrological cycle driving a dynamic ocean circulation (Trabucho-Alexandre et al., 2010). Studies on trace-elemental and (radiogenic) isotope compositions of Cenomanian marine successions have suggested a volcanic origin of OAE2, by delivering nutrients to the semi-enclosed proto-North Atlantic (Zheng et al., 2013, and references therein; Du Vivier et al., 2014). Deciphering the importance of volcanic and oceanographic processes requires tight constraints on their relative timing. Regularly occurring black cherts and shales below the Livello Bonarelli demonstrate that oceanic conditions in the Umbria-Marche Basin were punctuated by episodes of regional anoxia from the mid-Cenomanian onwards. Their hierarchical stacking pattern suggests an orbital control on the deposition of organic rich horizons (Mitchell et al., 2008; Lanci et al., 2010). Stable carbon isotope data reveal that long-term 15 variations in eccentricity paced the carbon cycle (Sprovieri et al., 2013) and sea level changes (Voigt et al., 2006) of the Late Cretaceous. Here we investigate the role of orbital forcing on climate and the carbon cycle, and, specifically, on organic-rich sedimentation prior, during, and after OAE2.
We also explore the potential for establishing an anchored astrochronology for the C/T interval in Europe. Recent improvements in the astronomical solution (La2011; Laskar et al., 2011b) and in the intercalibration of radiometric and astronomical dating techniques (Kuiper et al., 2008; Renne et al., 2013) allow the extension of the astronomical time scale into the Cretaceous. The C/T boundary in the Western Interior (USA) has been dated at 93.90 ± 0.15 Ma by intercalibration of radio-isotopic and astrochronologic time scales (Meyers et al., 2012b). Also, reinterpretation of proxy records spanning the C/T interval seems to resolve discrepancies in reported durations of the OAE2 (Sageman et al., 2006; Meyers et al., 2012a). The well-documented Italian rhythmic successions, reference sections for climatic processes in the Tethyan realm, need to be tied in with the absolute time scale. Biostratigraphic correlation to radioisotopically-dated ash beds in the Western Interior is complicated by the provinciality of faunas and floras. However, δ13C stratigraphy provides a reliable correlation tool (Gale et al., 2005) and we present a new 40Ar/39Ar age for the Thatcher bentonite from the Western Interior occurring within the mid-Cenomanian δ13C event (MCE). This study integrates the well-developed cyclostratigraphy from the Umbria-Marche Basin with radioisotopic ages from the Western Interior and derives a numerical timescale for this critical interval in Earth’s history.
Chemical reduction-oxidation mechanisms within mantle rocks link to the terrestrial carbon cycle by influencing the depth at which magmas can form, their composition, and ultimately the chemistry of gases released into the atmosphere. The oxidation state of the uppermost mantle has been widely accepted to be unchanged over the past 3800 m.y., based on the abundance of redox-sensitive elements in greenstone belt–associated samples of different ages. However, the redox signal in those rocks may have been obscured by their complex origins and emplacement on continental margins. In contrast, the source and processes occurring during decompression melting at spreading ridges are relatively well constrained. We retrieve primary redox conditions from metamorphosed mid-oceanic ridge basalts (MORBs) and picrites of various ages (ca. 3000–550 Ma), using V/Sc as a broad redox proxy. Average V/Sc values for Proterozoic suites (7.0 ± 1.4, 2σ, n = 6) are similar to those of modern MORB (6.8 ± 1.6), whereas Archean suites have lower V/Sc (5.2 ± 0.4, n = 5). The lower Archean V/Sc is interpreted to reflect both deeper melt extraction from the uppermost mantle, which becomes more reduced with depth, and an intrinsically lower redox state. The pressure-corrected oxygen fugacity (expressed relative to the fayalite-magnetite-quartz buffer, ΔFMQ, at 1 GPa) of Archean sample suites (ΔFMQ –1.19 ± 0.33, 2σ) is significantly lower than that of post-Archean sample suites, including MORB (ΔFMQ –0.26 ± 0.44). Our results imply that the reducing Archean atmosphere was in equilibrium with Earth’s mantle, and further suggest that magmatic gases crossed the threshold that allowed a build-up in atmospheric O2 levels ca. 3000 Ma, accompanied by the first “whiffs” of oxygen in sediments of that age.
In Belize, which is well known for the Belize Barrier Reef and its offshore atolls, coastal lagoons are frequent morphological features along the coast. They represent transitional environments between siliciclastic and carbonate settings. In order to shed light into the Holocene evolution of coastal lagoon environments, five localities along the central coast of Belize were selected as coring sites. These include four coastal lagoons and one marsh area, namely Mantatee Lagoon, Mullins River Beach, Colson Point Lagoon, Commerce Bight Lagoon and Sapodilla Lagoon. A total of 26 sediment cores with core lengths ranging from 109 cm to 500 cm, were drilled using a Lanesky-vibracorer. Overall, 73 m of Holocene sediments and Pleistocene soil were recovered. Together with 58 radiocarbon dates the sediments reveal details on the sediment architecture and depositional features of the localities.
Marine inundation of the mainland and coastal lagoon formation started around 6 kyrs cal BP.
As a response to sea-level rise during the Holocene transgression, facies retrograded towards the coast, as seen in marginal marine overlying brackish mollusk faunas. Evidence for late Holocene progradation of facies due to sea-level stagnation is largely lacking. The occurrence of landward thinning sand beds, hiatuses and marine fauna in lagoonal successions are indications of event (overwash) sedimentation. Sediments recovered are largely of Holocene age (<7980 cal BP), overlying Pleistocene sections. Analyses of sediment composition and texture, radiocarbon dating and mollusk shell identification were used to describe and correlate sedimentary facies.
XRD analyses have identified quartz as the dominant mineral, with the Maya Mountains as main source of coastal lagoon sediments. The most common sedimentary facies include peat and peaty sediment, mud, sand, and poorly sorted sediments. Pleistocene soil forms the basement of Holocene sediments. Holocene mud represents lagoon background permanent sedimentation.
Peats and peat-rich sequences were deposited in mangrove swamp environments, whereas sandy facies mainly occur in the shoreface, beach, barriers, bars, barrier spits and overwash deposits. Facies successions could be identified for each locality, but it has proven difficult to correlate the stratigraphic sequences, especially among localities. These differences among the five locations studied suggest that apart from regional influence such as sea-level rise, local environmental factors such as small-scale variation in geomorphology and resulting facies heterogeneity, connectivity of the lagoon with the sea, antecedent topography and river discharge, were responsible for coastal sedimentation and lagoon development in the Holocene of Belize.
Faunal composition and distribution patterns of mollusk assemblages from 20 shell concentrations in cores collected in coastal lagoons, a mangrove-fringed tidal inlet and the marginal marine area (shallow subtidal) show considerable variation due to environmental heterogeneity and the interplay of several environmental factors in the course of the mid-late Holocene (ca. 6000 cal BP to modern). The investigated fauna ≥2 mm comprises 2246 bivalve, gastropod and 11 scaphopod specimens. Fifty-three mollusk species, belonging to 42 families, were identified. The bivalve Anomalocardia cuneimeris and cerithid gastropods are the dominant species and account for 78% of the total fauna. Diversity indices are low in concentrations from lagoons and relatively high in the marginal marine and tidal inlet areas.
Based on cluster analysis and nonmetric multidimensional scaling (NMDS), seven lagoonal assemblages and three marginal marine/tidal inlet assemblages were defined. A separation between lagoonal and marginal marine/tidal inlet assemblages seen in ordination indicates a lagoon-onshore gradient. The statistical separation among lagoonal assemblages demonstrates environmental changes during the Holocene evolution of the coastal lagoons, which is probably related to the formation of barriers and spits. The controlling factors of species distribution patterns are difficult to figure out, probably due to the heterogeneity of the barrier-lagoon systems and the interaction of paleoecological and paleoenvironmental factors. In addition to the taxonomic analysis, a taphonomic analysis of 1827 valves of A. cuneimeris from coastal lagoons was carried out. There is no relationship between depth and age of shells and their taphonomic condition. Size-frequency distributions and right-left valve ratios of A. cuneimeris suggest that valves were not transported over long distances but were deposited parautochthonously in their original habitat. Shells from tidal inlet and marginal marine environments were also predominantly deposited in their original habitats.
Since the Belize coast was repeatedly affected by hurricanes and the paleohurricane record for this region is poor, the sediment cores have been examined in order to identify storm deposits.
The paleohurricane record presented in this study spans the past 8000 years and exhibits three periods with increased evidences of hurricane strikes occurring at 6000-4900 cal yr BP, 4200-3600 cal yr BP and 2200-1500 cal yr BP. Two earlier events around 7100 and 7900 cal yr BP and more recent events around 180 cal yr BP and during modern times have been detected. Sand layers, redeposited corals and lagoon shell concentrations have been used as proxies for storm deposition. Additionally, hiatuses and reversed ages may indicate storm influence. While sand layers and corals represent overwash deposits, the lagoon shell concentrations, which mainly comprise the bivalve Anomalocardia cuneimeris and cerithid gastropods, have been deposited due to changes in lagoon salinity during and after storm landfalls. Comparison with other studies reveals similarities with one record from Belize, but hardly any matches with other published records. The potential for paleotempestology reconstructions of the barrier-lagoon complexes along the central Belize coast differs depending on geomorphology, and deposition of washovers in the lagoon basins is limited, probably due to the interplay of biological, geological and geomorphological processes.