Mineral dust mobilisation, transport, and deposition in different climate epochs

Mineral dust is an important component of the Earth's climate system and provides essential nutrientsrnto oceans and rain forests. During atmospheric transport, dust particles directly and indirectly influencernweather and climate. The strength of dust sources and characteristics of the transport, in turn, mightrnbe subject to climatic changes. Earth system models help for a better understanding of these complexrnmechanisms.rnrnThis thesis applies the global climate model ECHAM5/MESSy Atmospheric Chemistry (EMAC) for simulationsrnof the mineral dust cycle under different climatic conditions. The prerequisite for suitable modelrnresults is the determination of the model setup reproducing the most realistic dust cycle in the recentrnclimate. Simulations with this setup are used to gain new insights into properties of the transatlanticrndust transport from Africa to the Americas and adaptations of the model's climate forcing factors allowrnfor investigations of the impact of climatic changes on the dust cycle.rnrnIn the first part, the most appropriate model setup is determined through a number of sensitivity experiments.rnIt uses the dust emission parametrisation from Tegen et al. 2002 and a spectral resolutionrnof T85, corresponding to a horizontal grid spacing of about 155 km. Coarser resolutions are not able tornaccurately reproduce emissions from important source regions such as the Bodele Depression in Chad orrnthe Taklamakan Desert in Central Asia. Furthermore, the representation of ageing and wet deposition ofrndust particles in the model requires a basic sulphur chemical mechanism. This setup is recommended forrnfuture simulations with EMAC focusing on mineral dust.rnrnOne major branch of the global dust cycle is the long-range transport from the world's largest dustrnsource, the Sahara, across the Atlantic Ocean. Seasonal variations of the main transport pathways to thernAmazon Basin in boreal winter and to the Caribbean during summer are well known and understood,rnand corroborated in this thesis. Both Eulerian and Lagrangian methods give estimates on the typicalrntransport times from the source regions to the deposition on the order of nine to ten days. Previously, arnhuge proportion of the dust transported across the Atlantic Ocean has been attributed to emissions fromrnthe Bodele Depression. However, the contribution of this hot spot to the total transport is very low inrnthe present results, although the overall emissions from this region are comparable. Both model resultsrnand data sets analysed earlier, such as satellite products, involve uncertainties and this controversy aboutrndust transport from the Bodele Depression calls for future investigations and clarification.rnrnAforementioned characteristics of the transatlantic dust transport just slightly change in simulationsrnrepresenting climatic conditions of the Little Ice Age in the middle of the last millennium with meanrnnear-surface cooling of 0.5 to 1 K. However, intensification of the West African summer monsoon duringrnthe Little Ice Age is associated with higher dust emissions from North African source regions and wetterrnconditions in the Sahel. Furthermore, the Indian Monsoon and dust emissions from the Arabian Peninsula,rnwhich are affected by this circulation, are intensified during the Little Ice Age, whereas the annual globalrndust budget is similar in both climate epochs. Simulated dust emission fluxes are particularly influencedrnby the surface parameters. Modifications of the model do not affect those in this thesis, to be able tornascribe all differences in the results to changed forcing factors, such as greenhouse gas concentrations.rnDue to meagre comparison data sets, the verification of results presented here is problematic. Deeperrnknowledge about the dust cycle during the Little Ice Age can be obtained by future simulations, based onrnthis work, and additionally using improved reconstructions of surface parameters. Better evaluation ofrnsuch simulations would be possible by refining the temporal resolution of reconstructed dust depositionrnfluxes from existing ice and marine sediment cores.

Global modeling of pollutant transport and deposition from anthropogenic emission hot spots on global scale

Urban centers significantly contribute to anthropogenic air pollution, although they cover only a minor fraction of the Earth's land surface. Since the worldwide degree of urbanization is steadily increasing, the anthropogenic contribution to air pollution from urban centers is expected to become more substantial in future air quality assessments. The main objective of this thesis was to obtain a more profound insight in the dispersion and the deposition of aerosol particles from 46 individual major population centers (MPCs) as well as the regional and global influence on the atmospheric distribution of several aerosol types. For the first time, this was assessed in one model framework, for which the global model EMAC was applied with different representations of aerosol particles. First, in an approach with passive tracers and a setup in which the results depend only on the source location and the size and the solubility of the tracers, several metrics and a regional climate classification were used to quantify the major outflow pathways, both vertically and horizontally, and to compare the balance between pollution export away from and pollution build-up around the source points. Then in a more comprehensive approach, the anthropogenic emissions of key trace species were changed at the MPC locations to determine the cumulative impact of the MPC emissions on the atmospheric aerosol burdens of black carbon, particulate organic matter, sulfate, and nitrate. Ten different mono-modal passive aerosol tracers were continuously released at the same constant rate at each emission point. The results clearly showed that on average about five times more mass is advected quasi-horizontally at low levels than exported into the upper troposphere. The strength of the low-level export is mainly determined by the location of the source, while the vertical transport is mainly governed by the lifting potential and the solubility of the tracers. Similar to insoluble gas phase tracers, the low-level export of aerosol tracers is strongest at middle and high latitudes, while the regions of strongest vertical export differ between aerosol (temperate winter dry) and gas phase (tropics) tracers. The emitted mass fraction that is kept around MPCs is largest in regions where aerosol tracers have short lifetimes; this mass is also critical for assessing the impact on humans. However, the number of people who live in a strongly polluted region around urban centers depends more on the population density than on the size of the area which is affected by strong air pollution. Another major result was that fine aerosol particles (diameters smaller than 2.5 micrometer) from MPCs undergo substantial long-range transport, with about half of the emitted mass being deposited beyond 1000 km away from the source. In contrast to this diluted remote deposition, there are areas around the MPCs which experience high deposition rates, especially in regions which are frequently affected by heavy precipitation or are situated in poorly ventilated locations. Moreover, most MPC aerosol emissions are removed over land surfaces. In particular, forests experience more deposition from MPC pollutants than other land ecosystems. In addition, it was found that the generic treatment of aerosols has no substantial influence on the major conclusions drawn in this thesis. Moreover, in the more comprehensive approach, it was found that emissions of black carbon, particulate organic matter, sulfur dioxide, and nitrogen oxides from MPCs influence the atmospheric burden of various aerosol types very differently, with impacts generally being larger for secondary species, sulfate and nitrate, than for primary species, black carbon and particulate organic matter. While the changes in the burdens of sulfate, black carbon, and particulate organic matter show an almost linear response for changes in the emission strength, the formation of nitrate was found to be contingent upon many more factors, e.g., the abundance of sulfuric acid, than only upon the strength of the nitrogen oxide emissions. The generic tracer experiments were further extended to conduct the first risk assessment to obtain the cumulative risk of contamination from multiple nuclear reactor accidents on the global scale. For this, many factors had to be taken into account: the probability of major accidents, the cumulative deposition field of the radionuclide cesium-137, and a threshold value that defines contamination. By collecting the necessary data and after accounting for uncertainties, it was found that the risk is highest in western Europe, the eastern US, and in Japan, where on average contamination by major accidents is expected about every 50 years.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global three-dimensional ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during the early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global 3-D ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

Development of room temperature operating single electron transistor using FIB etching and deposition technology

The single-electron transistor (SET) is one of the best candidates for future nano electronic circuits because of its ultralow power consumption, small size and unique functionality. SET devices operate on the principle of Coulomb blockade, which is more prominent at dimensions of a few nano meters. Typically, the SET device consists of two capacitively coupled ultra-small tunnel junctions with a nano island between them. In order to observe the Coulomb blockade effects in a SET device the charging energy of the device has to be greater that the thermal energy. This condition limits the operation of most of the existing SET devices to cryogenic temperatures. Room temperature operation of SET devices requires sub-10nm nano-islands due to the inverse dependence of charging energy on the radius of the conducting nano-island. Fabrication of sub-10nm structures using lithography processes is still a technological challenge. In the present investigation, Focused Ion Beam based etch and deposition technology is used to fabricate single electron transistors devices operating at room temperature. The SET device incorporates an array of tungsten nano-islands with an average diameter of 8nm. The fabricated devices are characterized at room temperature and clear Coulomb blockade and Coulomb oscillations are observed. An improvement in the resolution limitation of the FIB etching process is demonstrated by optimizing the thickness of the active layer. SET devices with structural and topological variation are developed to explore their impact on the behavior of the device. The threshold voltage of the device was minimized to ~500mV by minimizing the source-drain gap of the device to 17nm. Vertical source and drain terminals are fabricated to realize single-dot based SET device. A unique process flow is developed to fabricate Si dot based SET devices for better gate controllability in the device characteristic. The device vi parameters of the fabricated devices are extracted by using a conductance model. Finally, characteristic of these devices are validated with the simulated data from theoretical modeling.

LATTICE BOLTZMANN METHOD AND CELLULAR AUTOMATA SIMULATION OF PARTICLE MOTION AND DEPOSITION IN 2-D CASE

This technical report discusses the application of the Lattice Boltzmann Method (LBM) and Cellular Automata (CA) simulation in fluid flow and particle deposition. The current work focuses on incompressible flow simulation passing cylinders, in which we incorporate the LBM D2Q9 and CA techniques to simulate the fluid flow and particle loading respectively. For the LBM part, the theories of boundary conditions are studied and verified using the Poiseuille flow test. For the CA part, several models regarding simulation of particles are explained. And a new Digital Differential Analyzer (DDA) algorithm is introduced to simulate particle motion in the Boolean model. The numerical results are compared with a previous probability velocity model by Masselot [Masselot 2000], which shows a satisfactory result.

IMPACT OF 2000-2050 CLIMATE CHANGE ON GLOBAL ATMOSPHERIC TRANSPORT AND DEPOSITION OF MERCURY

Since it is very toxic and accumulates in organisms, particularly in fish, mercury is a very important pollutant and one of the most studies. And this concern over the toxicity and human health risks of mercury has prompted efforts to regulate anthropogenic emissions. As mercury pollution problem is getting increasingly serious, we are curious about how serious this problem will be in the future. What is more, how the climate change in the future will affect the mercury concentration in the atmosphere. So we investigate the impact of climate change on mercury concentration in the atmosphere. We focus on the comparison between the mercury data for year 2000 and for year 2050. The GEOS-Chem model shows that the mercury concentrations for all tracers (1 to 3), elemental mercury (Hg(0)), divalent mercury (Hg(II)) and primary particulate mercury (Hg(P)) have differences between 2000 and 2050 in most regions over the world. From the model results, we can see the climate change from 2000 to 2050 would decrease Hg(0) surface concentration in most of the world. The driving factors of Hg(0) surface concentration changes are natural emissions(ocean and vegetation) and the transformation reactions between Hg(0) and Hg(II). The climate change from 2000 to 2050 would increase Hg(II) surface concentration in most of mid-latitude continental parts of the world while decreasing Hg(II) surface concentration in most of high-latitude part of the world. The driving factors of Hg(II) surface concentration changes is deposition amount change (majorly wet deposition) from 2000 to 2050 and the transformation reactions between Hg(0) and Hg(II). Climate change would increase Hg(P) concentration in most of mid-latitude area of the world and meanwhile decrease Hg(P) concentration in most of high-latitude regions of the world. For the Hg(P) concentration changes, the major driving factor is the deposition amount change (mainly wet deposition) from 2000 to 2050.

Modelling woody material transport and deposition in alpine rivers

Recent flood events in Switzerland and Western Austria in 2005 were characterised by an increase in impacts and associated losses due to the transport of woody material. As a consequence, protection measures and bridges suffered considerable damages. Furthermore, cross-sectional obstructions due to woody material entrapment caused unexpected flood plain inundations resulting in severe damage to elements at risk. Until now, the transport of woody material is neither sufficiently taken into account nor systematically considered, leading to prediction inaccuracies during the procedure of hazard mapping. To close this gap, we propose a modelling approach that (1) allows the estimation of woody material recruitment from wood-covered banks and flood plains; (2) allows the evaluation of the disposition for woody material entrainment and transport to selected critical configurations along the stream and that (3) enables the delineation of hazard process patterns at these critical configurations. Results from a case study suggest the general applicability of the concept. This contribution to woody material transport analysis refines flood hazard assessments due to the consideration of woody material transport scenarios.

CaCO3, density, sedimentation and accumulation rates at DSDP Hole 68-503B

Pelagic sediments from DSDP Hole 5O3B contain, in their carbonate abundance data, a clear record of glacial-interglacial cycles. The eolian component of those sediments was analyzed over the past four carbonate cycles, and the mass accumulation rate (MAR) and grain size of the eolian component was determined. Eolian MARs range from 24 to 169 mg/cm**2/10**3y. and commonly are higher by a factor of three to five during times of glacial retreat. Reduced contribution during periods of glaciation most likely reflects glacial-age humidity in the American source. Grain-size values (phi50) range from 8.25 to a minimum of 8.79phi-a variation in grain mass by a factor of 3.1. Larger grains reflect more vigorous atmospheric circulation, but sizes do not covary with the carbonate or eolian accumulation curves. These data suggest that the intensity of atmospheric circulation in the tropics may reflect the 42,000 y.-tilt cycle rather than the 100,000 y.-cycle of glacial advance.

Late Neogene carbonate sedimentation and stable isotope record of ODP Site 121-758

The evolution of oceanic and climatic conditions the northeast Indian Ocean during the last 7 m.y. is revealed in the sediments from Site 758. We present detailed and continuous records of d18O and d13C from planktonic foraminifers, weight percent calcium carbonate, weight percent coarse fraction, magnetic susceptibility, and geomagnetic reversals. Sample spacing of the records ranges from 3 to 10 cm and is equivalent to an average time interval of 2000 to 6000 yr. Despite the fact that core recovery ranged between 100% and 105%, recovery gaps as large as 2.7 m occurred at nearly every break between advanced hydraulic piston cores. Approximately 12% of the late Neogene sequence was not recovered in each of the two holes drilled at Site 758. To circumvent the discontinuity introduced by the gaps, a composite depth section was constructed from multiple cores taken from offset holes at Site 758. The resulting composite depth section extends continuously from 0 to 116 mbsf, from the Holocene to the upper Miocene. A detailed chronostratigraphy is based on geomagnetic reversals which extend from the Brunhes Chron to Chron 6, and on d18O stages 1 through 105, which span from 0 to 2.5 Ma. The d18O record is dominated by a ~40-k.y. cycle in the late Pliocene and early Pleistocene, and is followed by a change to a ~100-k.y. cycle in the late Pleistocene. The mid-Pleistocene transition between these two modes of variability occurs between d18O stages 25 and 22 (between 860 and 800 Ka). Thirteen major volcanic ash horizons from the Indonesian arc are observed throughout the sedimentary section and are dated by their relative position within the geomagnetic reversals and the d18O chronostratigraphy. Since 5 Ma, there has been a long-term decline in weight percent CaCO3 and CaCO3 mass accumulation rates, and an associated rise in non-CaCO3 mass accumulation rates. We attribute these changes to a decrease in CaCO3 productivity and an increase in terrigenous sedimentation through enhanced riverine input. Such input may be linked to rapid tectonic uplift of the Himalayas and the Tibetan Plateau via mechanisms such as the intensification of the monsoonal rains, increased fluvial erosion, and regional glaciation. The long-term increase in percent coarse fraction since 5 Ma suggests a gradual increase in CaCO3 preservation. Higher frequency fluctuations in CaCO3 preservation are superimposed on the long-term trend and are related to climate fluctuations. The abrupt drop (-50%) in CaCO3 accumulation at 3.4 Ma signals a dramatic decrease in CaCO3 production that occurred over much of the Indian Ocean.

Sedimentation and sediment redistribution on the Cocos Ridge from two seismic acquisition cruises, NEMO-3 and MV1014

We use digital seismic reflection profiles within a 1° * 1° survey area on the Cocos Ridge (COCOS6N) to study the extent and timing of sedimentation and sediment redistribution on the Cocos Ridge. The survey was performed to understand how sediment focusing might affect paleoceanographic flux measurements in a region known for significant downslope transport. COCOS6N contains ODP Site 1241 to ground truth the seismic stratigraphy, and there is a seamount ridge along the base of the ridge that forms a basin (North Flank Basin) to trap sediments transported downslope. Using the Site 1241 seismic stratigraphy and densities extrapolated from wireline logging, we document mass accumulation rates (MARs) since 11.2 Ma. The average sediment thickness at COCOS6N is 196 m, ranging from outcropping basalt at the ridge crest to ~ 400 m at North Flank Basin depocenters. Despite significant sediment transport, the average sedimentation over the entire area is well correlated to sediment fluxes at Site 1241. A low mass accumulation rate (MAR) interval is associated with the 'Miocene carbonate crash' interval even though COCOS6N was at the equator at that time and relatively shallow. Highest MAR occurs within the late Miocene-early Pliocene biogenic bloom interval. Lowest average MAR is in the Pleistocene, as plate tectonic motions caused COCOS6N to leave the equatorial productivity zone. The Pliocene and Pleistocene also exhibit higher loss of sediment from the ridge crest and transport to North Flank Basin. Higher tidal energy on the ridge caused by tectonic movement toward the margin increased sediment focusing in the younger section.

(Table 2) Sedimentation and accumulation rates of DSDP Holes 38-337 and 38-338

The modern depositional environment of the deep Norwegian-Greenland Sea is highly asymmetric in an E-W direction because of the hydrography of the surface water masses and because of the more or less permanent pack ice cover of the East Greenland Current regime along the Greenland continental margin. By means of sedimentation rates we have tried to investigate whether this hydrographic asymmetry influenced the sediment input to the Norwegian-Greenland Sea over the past 60 m.y. Sediment input can be quantified if thicknesses of sediment sections accumulated over known time intervals can be measured and if some of their physical properties have been determined. Sedimentation rates have been estimated for Tertiary and Quaternary times, and their temporal as well as their spatial changes are discussed. Basin structure and morphology exerted an important influence on sediment distribution. During the Early Tertiary major sediment source regions in the southern Barents Sea and to the north and west of Iceland could be identified; these source regions supplied the bulk of the sediment fill of the Norwegian-Greenland Sea. Since inception of a "glacial" type sedimentation major elements of the sea surface circulation seem to have controlled the sediment input into this polar and subpolar deep-sea basin.