550 Geowissenschaften
Refine
Year of publication
- 2016 (88) (remove)
Document Type
- Article (67)
- Doctoral Thesis (15)
- Contribution to a Periodical (3)
- Part of a Book (1)
- Master's Thesis (1)
- Working Paper (1)
Has Fulltext
- yes (88)
Is part of the Bibliography
- no (88)
Keywords
- Atmospheric chemistry (3)
- Anthropogene Klimaänderung (2)
- Modellierung (2)
- BSRN (1)
- Bayesian network (1)
- Biogeochemistry (1)
- Bochum (1)
- CCLM (1)
- CORDEX-EAS (1)
- COSMO-CLM (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.
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.
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 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
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.
A three-dimensional gridded climatology of carbon monoxide (CO) has been developed by trajectory mapping of global MOZAIC-IAGOS in situ measurements from commercial aircraft data. CO measurements made during aircraft ascent and descent, comprising nearly 41 200 profiles at 148 airports worldwide from December 2001 to December 2012, are used. Forward and backward trajectories are calculated from meteorological reanalysis data in order to map the CO measurements to other locations and so to fill in the spatial domain. This domain-filling technique employs 15 800 000 calculated trajectories to map otherwise sparse MOZAIC-IAGOS data into a quasi-global field. The resulting trajectory-mapped CO data set is archived monthly from 2001 to 2012 on a grid of 5° longitude × 5° latitude × 1 km altitude, from the surface to 14 km altitude.
The mapping product has been carefully evaluated, firstly by comparing maps constructed using only forward trajectories and using only backward trajectories. The two methods show similar global CO distribution patterns. The magnitude of their differences is most commonly 10 % or less and found to be less than 30 % for almost all cases. Secondly, the method has been validated by comparing profiles for individual airports with those produced by the mapping method when data from that site are excluded. While there are larger differences below 2 km, the two methods agree very well between 2 and 10 km with the magnitude of biases within 20 %. Finally, the mapping product is compared with global MOZAIC-IAGOS cruise-level data, which were not included in the trajectory-mapped data set, and with independent data from the NOAA aircraft flask sampling program. The trajectory-mapped MOZAIC-IAGOS CO values show generally good agreement with both independent data sets.
Maps are also compared with version 6 data from the Measurements Of Pollution In The Troposphere (MOPITT) satellite instrument. Both data sets clearly show major regional CO sources such as biomass burning in Central and southern Africa and anthropogenic emissions in eastern China. While the maps show similar features and patterns, and relative biases are small in the lowermost troposphere, we find differences of ∼ 20 % in CO volume mixing ratios between 500 and 300 hPa. These upper-tropospheric biases are not related to the mapping procedure, as almost identical differences are found with the original in situ MOZAIC-IAGOS data. The total CO trajectory-mapped MOZAIC-IAGOS column is also higher than the MOPITT CO total column by 12–16 %.
The data set shows the seasonal CO cycle over different latitude bands and altitude ranges as well as long-term trends over different latitude bands. We observe a decline in CO over the northern hemispheric extratropics and the tropics consistent with that reported by previous studies using other data sources.
We anticipate use of the trajectory-mapped MOZAIC-IAGOS CO data set as an a priori climatology for satellite retrieval and for air quality model validation and initialization.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Die vorliegende Arbeit beschäftigt sich mit der Entwicklung von regionalen Klimasimulationen für die Region Ostasien. Hierfür werden zwei verschiedene Modellierungsansätze verwendet. Der dynamische Regionalmodellierungsansatz, vertreten durch COSMO CLM (CCLM), und der statistische Modellierungsansatz, vertreten durch STARS. Die Simulationen erfolgten unter den Rahmenbedingungen des Coordinated Regional Climate Downscaling Experiment (CORDEX). Beide Regionalmodelle wurden im Rahmen dieser Arbeit umfassend für die Region CORDEX-Ostasien kalibriert und evaluiert. Das statistische Modell STARS wurde hierbei erstmals auf kontinentaler Ebene angewendet. Auf Basis der kalibrierten Modelle wurden Projektionen der zukünftigen klimatischen Entwicklung der Region durchgeführt.
Zur Auswertung der einzelnen Kalibrierungsläufe wurde ein komplexes Evaluierungsschema, mit einem Gütekennzahlensystem basierend auf einer linearisierten Form der relativen Modelldifferenz, entwickelt. Neben den etablierten univariaten statistischen Kennwerten (Mittelwert, Varianz, Trend) enthält das Gütekennzahlensystem auch ein bivariates statistisches Maß, welches die zweidimensionalen Stichprobenverteilungen zweier Variablen (beispielsweise Temperatur und Niederschlag) bewertet.
Im Rahmen der Kalibrierung konnte ein Großteil des Parameterraums des statistischen Modells STARS systematisch untersucht werden. Es zeigte sich, dass nur wenige Parameter einen Einfluss auf die Simulationen haben. Die meisten Parameter zeigten eine geringe und teilweise unsystematische Beeinflussung. Es konnte zudem eine Schwachstelle des Modells in Bezug auf die Variablenkorrelationen identifiziert werden. Bei der Kalibrierung des dynamischen Regionalmodells CCLM zeigte sich, dass aufgrund der groben horizontalen Auflösung des Modells eine signifikante Verbesserung der Simulationen durch eine Anpassung der physikalischen Parametrisierungen erfolgen kann.
Im Rahmen einer abschließenden Evaluierung wurden beide Modelle hinsichtlich ihres räumlichen Bias, des simulierten Jahresgangs und der Abbildung des asiatischen Monsunphänomens untersucht. Im ersten Punkt ergab sich kein qualitativer Unterschied zwischen CCLM und STARS. Beide Modelle zeigen eine deutliche Überschätzung der 2m-Temperatur im Winter über dem nördlichen Teil CORDEX-Ostasiens und eine Überschätzung des Luftdrucks über dem Hochland von Tibet im Sommer. Unterschiede zwischen beiden Modellen ergaben sich hingegen beim simulierten Jahresgang.
In Bezug auf die Modellierung des Monsunphänomens zeigt CCLM eine Unterschätzung der Intensität des indischen Sommermonsuns und eine Überschätzung des Sommermonsuns über dem westlichen Nordpazifik. Das statistische Modell STARS zeigte eine Auffälligkeit bei der Simulation des Jahresgangs sowie der räumlichen und zeitlichen Entwicklung des Sommermonsuns. Aufgrund der Konzeption des Modells ergab sich in einzelnen Regionen eine systematische Deformation des Jahresgangs. Trotz der identifizierten Schwachstellen von CCLM und STARS, bilden beide Modelle das Klima über der Region CORDEX-Ostasien qualitativ ähnlich gut ab wie aktuelle Reanalysen (ERA-Interim).
Auf Basis der kalibrierten und evaluierten Modelle wurden Klimaprojektionen für einen nahen (2020-2046), mittleren (2041-2070), und späten (2071-2100) Projektionszeitraum unter den Emissionsszenarien RCP2.6, RCP4.5 und RCP8.5 durchgeführt. Aufgrund von Modellbeschränkungen begrenzen sich die Rechnungen des Modells STARS auf den nahen Projektionszeitraum und die Emissionsszenarien RCP2.6 und RCP4.5. Die Projektionen beider Modelle zeigen eine deutliche und statistisch signifikante Erhöhung der 2m-Temperatur über der gesamten Region mit einer stärkeren Erwärmung über dem Kontinent gegenüber dem Meer. Aufgrund der relativ großen interannulären Variabilität des Niederschlags und des Luftdrucks werden statistisch nicht signifikante Änderungssignale und teils widersprüchliche Änderungen für den nahen Projektionszeitraum simuliert. Für den späten Projektionszeitraum ergeben sich jedoch deutliche Änderungssignale in den Simulationen des Modells CCLM. Insbesondere über dem Hochland von Tibet wird für den Zeitraum von 2071-2100 eine Temperaturerhöhung von über 7.0°C simuliert. Der Luftdruck und der Niederschlag zeigen räumlich heterogene Änderungssignale. Die spezifische Ausprägung der Luftdruckänderungen deutet auf eine Abschwächung der indischen Sommermonsunzirkulation und eine deutlichen Intensivierung des Sommermonsun über dem westlichen Nordpazifik hin. Die Niederschlagsänderungen über dem ostasiatischen Monsungebiet lassen auf eine Entkopplung der östlichen Monsunsysteme schließen. Trotz der heterogenen Änderungssignale im Niederschlag wird in den meisten Regionen eine Zunahme der Intensität von Extremniederschlägen simuliert. Dies gilt selbst für Regionen mit einer simulierten Abnahme der jährlichen Niederschlagssumme wie Westindonesien.