PM(10) remote sensing from geostationary SEVIRI and polar-orbiting MODIS sensors over the complex terrain of the European Alpine region

The subject of this study is to investigate the capability of spaceborne remote sensing data to predict ground concentrations of PM10 over the European Alpine region using satellite derived Aerosol Optical Depth (AOD) from the geostationary Spinning Enhanced Visible and InfraRed Imager (SEVIRI) and the polar-orbiting MODerate resolution Imaging Spectroradiometer (MODIS). The spatial and temporal resolutions of these aerosol products (10 km and 2 measurements per day for MODIS, ∼ 25 km and observation intervals of 15 min for SEVIRI) permit an evaluation of PM estimation from space at different spatial and temporal scales. Different empirical linear relationships between coincident AOD and PM10 observations are evaluated at 13 ground-based PM measurement sites, with the assumption that aerosols are vertically homogeneously distributed below the planetary Boundary Layer Height (BLH). The BLH and Relative Humidity (RH) variability are assessed, as well as their impact on the parameterization. The BLH has a strong influence on the correlation of daily and hourly time series, whilst RH effects are less clear and smaller in magnitude. Despite its lower spatial resolution and AOD accuracy, SEVIRI shows higher correlations than MODIS (rSEV∼ 0.7, rMOD∼ 0.6) with regard to daily averaged PM10. Advantages from MODIS arise only at hourly time scales in mountainous locations but lower correlations were found for both sensors at this time scale (r∼ 0.45). Moreover, the fraction of days in 2008 with at least one satellite observation was 27% for SEVIRI and 17% for MODIS. These results suggest that the frequency of observations plays an important role in PM monitoring, while higher spatial resolution does not generally improve the PM estimation. Ground-based Sun Photometer (SP) measurements are used to validate the satellite-based AOD in the study region and to discuss the impact of aerosols' micro-physical properties in the empirical models. A lower error limit of 30 to 60% in the PM10 assessment from space is estimated in the study area as a result of AOD uncertainties, variability of aerosols properties and the heterogeneity of ground measurement sites. It is concluded that SEVIRI has a similar capacity to map PM as sensors on board polar-orbiting platforms, with the advantage of a higher number of observations. However, the accuracy represents a serious limitation to the applicability of satellites for ground PM mapping, especially in mountainous areas.

Ground-based and airborne in-situ measurements of the Eyjafjallajökull volcanic aerosol plume in Switzerland in spring 2010

The volcanic aerosol plume resulting from the Eyjafjallajökull eruption in Iceland in April and May 2010 was detected in clear layers above Switzerland during two periods (17–19 April 2010 and 16–19 May 2010). In-situ measurements of the airborne volcanic plume were performed both within ground-based monitoring networks and with a research aircraft up to an altitude of 6000 m a.s.l. The wide range of aerosol and gas phase parameters studied at the high altitude research station Jungfraujoch (3580 m a.s.l.) allowed for an in-depth characterization of the detected volcanic aerosol. Both the data from the Jungfraujoch and the aircraft vertical profiles showed a consistent volcanic ash mode in the aerosol volume size distribution with a mean optical diameter around 3 ± 0.3 μm. These particles were found to have an average chemical composition very similar to the trachyandesite-like composition of rock samples collected near the volcano. Furthermore, chemical processing of volcanic sulfur dioxide into sulfate clearly contributed to the accumulation mode of the aerosol at the Jungfraujoch. The combination of these in-situ data and plume dispersion modeling results showed that a significant portion of the first volcanic aerosol plume reaching Switzerland on 17 April 2010 did not reach the Jungfraujoch directly, but was first dispersed and diluted in the planetary boundary layer. The maximum PM10 mass concentrations at the Jungfraujoch reached 30 μgm−3 and 70 μgm−3 (for 10-min mean values) duri ng the April and May episode, respectively. Even low-altitude monitoring stations registered up to 45 μgm−3 of volcanic ash related PM10 (Basel, Northwestern Switzerland, 18/19 April 2010). The flights with the research aircraft on 17 April 2010 showed one order of magnitude higher number concentrations over the northern Swiss plateau compared to the Jungfraujoch, and a mass concentration of 320 (200–520) μgm−3 on 18 May 2010 over the northwestern Swiss plateau. The presented data significantly contributed to the time-critical assessment of the local ash layer properties during the initial eruption phase. Furthermore, dispersion models benefited from the detailed information on the volcanic aerosol size distribution and its chemical composition.

Long-Range Pollution Transport: Trans-Atlantic Mechanisms and Lagrangian Modeling Methods

Over the past several decades, it has become apparent that anthropogenic activities have resulted in the large-scale enhancement of the levels of many trace gases throughout the troposphere. More recently, attention has been given to the transport pathway taken by these emissions as they are dispersed throughout the atmosphere. The transport pathway determines the physical characteristics of emissions plumes and therefore plays an important role in the chemical transformations that can occur downwind of source regions. For example, the production of ozone (O3) is strongly dependent upon the transport its precursors undergo. O3 can initially be formed within air masses while still over polluted source regions. These polluted air masses can experience continued O3 production or O3 destruction downwind, depending on the air mass's chemical and transport characteristics. At present, however, there are a number of uncertainties in the relationships between transport and O3 production in the North Atlantic lower free troposphere. The first phase of the study presented here used measurements made at the Pico Mountain observatory and model simulations to determine transport pathways for US emissions to the observatory. The Pico Mountain observatory was established in the summer of 2001 in order to address the need to understand the relationships between transport and O3 production. Measurements from the observatory were analyzed in conjunction with model simulations from the Lagrangian particle dispersion model (LPDM), FLEX-PART, in order to determine the transport pathway for events observed at the Pico Mountain observatory during July 2003. A total of 16 events were observed, 4 of which were analyzed in detail. The transport time for these 16 events varied from 4.5 to 7 days, while the transport altitudes over the ocean ranged from 2-8 km, but were typically less than 3 km. In three of the case studies, eastward advection and transport in a weak warm conveyor belt (WCB) airflow was responsible for the export of North American emissions into the FT, while transport in the FT was governed by easterly winds driven by the Azores/Bermuda High (ABH) and transient northerly lows. In the fourth case study, North American emissions were lofted to 6-8 km in a WCB before being entrained in the same cyclone's dry airstream and transported down to the observatory. The results of this study show that the lower marine FT may provide an important transport environment where O3 production may continue, in contrast to transport in the marine boundary layer, where O3 destruction is believed to dominate. The second phase of the study presented here focused on improving the analysis methods that are available with LPDMs. While LPDMs are popular and useful for the analysis of atmospheric trace gas measurements, identifying the transport pathway of emissions from their source to a receptor (the Pico Mountain observatory in our case) using the standard gridded model output, particularly during complex meteorological scenarios can be difficult can be difficult or impossible. The transport study in phase 1 was limited to only 1 month out of more than 3 years of available data and included only 4 case studies out of the 16 events specifically due to this confounding factor. The second phase of this study addressed this difficulty by presenting a method to clearly and easily identify the pathway taken by only those emissions that arrive at a receptor at a particular time, by combining the standard gridded output from forward (i.e., concentrations) and backward (i.e., residence time) LPDM simulations, greatly simplifying similar analyses. The ability of the method to successfully determine the source-to-receptor pathway, restoring this Lagrangian information that is lost when the data are gridded, is proven by comparing the pathway determined from this method with the particle trajectories from both the forward and backward models. A sample analysis is also presented, demonstrating that this method is more accurate and easier to use than existing methods using standard LPDM products. Finally, we discuss potential future work that would be possible by combining the backward LPDM simulation with gridded data from other sources (e.g., chemical transport models) to obtain a Lagrangian sampling of the air that will eventually arrive at a receptor.

Laminar and turbulent nozzle-jet flows and their acoustic near-field

We investigate numerically the effects of nozzle-exit flow conditions on the jet-flow development and the near-field sound at a diameter-based Reynolds number of Re D = 18 100 and Mach number Ma = 0.9. Our computational setup features the inclusion of a cylindrical nozzle which allows to establish a physical nozzle-exit flow and therefore well-defined initial jet-flow conditions. Within the nozzle, the flow is modeled by a potential flow core and a laminar, transitional, or developing turbulent boundary layer. The goal is to document and to compare the effects of the different jet inflows on the jet flow development and the sound radiation. For laminar and transitional boundary layers, transition to turbulence in the jet shear layer is governed by the development of Kelvin-Helmholtz instabilities. With the turbulent nozzle boundary layer, the jet flow development is characterized by a rapid changeover to a turbulent free shear layer within about one nozzle diameter. Sound pressure levels are strongly enhanced for laminar and transitional exit conditions compared to the turbulent case. However, a frequency and frequency-wavenumber analysis of the near-field pressure indicates that the dominant sound radiation characteristics remain largely unaffected. By applying a recently developed scaling procedure, we obtain a close match of the scaled near-field sound spectra for all nozzle-exit turbulence levels and also a reasonable agreement with experimental far-field data.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

Physical oceanography from the TRACTOR project

Primary Objectives - Describe and quantify the present strength and variability of the circulation and oceanic processes of the Nordic Seas regions using primarily observations of the long term spread of a tracer purposefully released into the Greenland Sea Gyre in 1996. - Improve our understanding of ocean processes critical to the thermaholine circulation in the Nordic Seas regions so as to be able to predict how this region may respond to climate change. - Assess the role of mixing and ageing of water masses on the carbon transport and the role of the thermohaline circulation in carbon storage using water transports and mixing coefficients derived from the tracer distribution. Specific Objectives Perform annual hydrographic, chemical and SF6 tracer surveys into the Nordic regions in order to: - Measure lateral and diapycnal mixing rates in the Greenland Sea Gyre and in the surrounding regions. - Document the depth and rates of convective mixing in the Greenland Sea using the SF6 and the water masses characteristics. - Measure the transit time and transport of water from the Greenland Sea to surrounding seas and outflows. Document processes of water mass transformation and entrainment occurring to water emanating from the central Greenland Sea. - Measure diapycnal mixing rates in the bottom and margins of the Greenland Sea basin using the SF6 signal observed there. Quantify the potential role of bottom boundary-layer mixing in the ventilation of the Greenland Sea Deep Water in absence of deep convection. Monitor the variability of the entrainment of water from the Greenland Sea using time series auto-sampler moorings at strategic positions i.e., sill of the Denmark Strait, Labrador Sea, Jan Mayen fracture zone and Fram Strait. Relate the observed variability of the tracer signal in the outflows to convection events in the Greenland Sea and local wind stress events. Obtain a better description of deepwater overflow and entrainment processes in the Denmark Strait and Faeroe Bank Channel overflows and use these to improve modelling of deepwater overflows. Monitor the tracer invasion into the North Atlantic using opportunistic SF6 measurements from other cruises: we anticipate that a number of oceanographic cruises will take place in the north-east Atlantic and the Labrador Sea. It should be possible to get samples from some cruises for SF6 measurements. Use process models to describe the spread of the tracer to achieve better parameterisation for three-dimensional models. One reason that these are so resistant to prediction is that our best ocean models are as yet some distance from being good enough, to predict climate and climate change.

Description of sampling sites in Northern Russia

The youngest ice marginal zone between the White Sea and the Ural mountains is the W-E trending belt of moraines called the Varsh-Indiga-Markhida-Harbei-Halmer-Sopkay, here called the Markhida line. Glacial elements show that it was deposited by the Kara Ice Sheet, and in the west, by the Barents Ice Sheet. The Markhida moraine overlies Eemian marine sediments, and is therefore of Weichselian age. Distal to the moraine are Eemian marine sediments and three Palaeolithic sites with many C-14 dates in the range 16-37 ka not covered by till, proving that it represents the maximum ice sheet extension during the Weichselian. The Late Weichselian ice limit of M. G. Grosswald is about 400 km (near the Urals more than 700 km) too far south. Shorelines of ice dammed Lake Komi, probably dammed by the ice sheet ending at the Markhida line, predate 37 ka. We conclude that the Markhida line is of Middle/Early Weichselian age, implying that no ice sheet reached this part of Northern Russia during the Late Weichselian. This age is supported by a series of C-14 and OSL dates inside the Markhida line all of >45 ka. Two moraine loops protrude south of the Markhida line; the Laya-Adzva and Rogavaya moraines. These moraines are covered by Lake Komi sediments, and many C-14 dates on mammoth bones inside the moraines are 26-37 ka. The morphology indicates that the moraines are of Weichselian age, but a Saalian age cannot be excluded. No post-glacial emerged marine shorelines are found along the Barents Sea coast north of the Markhida line.

Mineralogy and diagenesis of slope sediments at DSDP Leg 84 holes

The distribution and composition of minerals in the silt and clay fraction of the fine-grained slope sediments were examined. Special interest was focused on diagenesis. The results are listed as follows. (1) Smectite, andesitic Plagioclase, quartz, and low-Mg calcite are the main mineral components of the sediment. Authigenic dolomite was observed in the weathering zones of serpentinites, together with aragonite, as well as in clayey silt. (2) The mineralogy and geochemistry of the sediments is analogous to that of the andesitic rocks of Costa Rica and Guatemala. (3) Unstable components like volcanic glass, amphiboles, and pyroxenes show increasing etching with depth. (4) The diagenetic alteration of opal-A skeletons from etching pits and replacement by opal-CT to replacement by chalcedony as a final stage corresponds to the typical opal diagenesis. (5) Clinoptilolite is the stable zeolite mineral according to mineral stability fields; its neoformation is well documented. (6) The early diagenesis of smectites is shown by an increase of crystallinity with depth. Only the smectites in the oldest sediments (Oligocene and early Eocene) contain nonexpanding illite layers.

Geological mapping of the Heber, Rhüdener Sattel, Northern Germany

This master thesis describes the geological mapping of an 8 km**2 area of the 300 m high elevation HEBER in Northern Germany which is part of the Rhüdener Sattel (Harzvorland). The geology consists of mesozoic rocks from Buntsandstein to Keuper; major parts are Muschelkalk, partly covered with a thin layer of Quaternary loess. The map includes measurements of strike and dip in various outcrops, drafting a scetch of the probable tectonic structure of Rühdener Sattel.

Whole-rock and mineral geochemistry of ODP Hole 176-735B gabbros

We report mineral chemistry, whole-rock major element compositions, and trace element analyses on Hole 735B samples drilled and selected during Leg 176. We discuss these data, together with Leg 176 shipboard data and Leg 118 sample data from the literature, in terms of primary igneous petrogenesis. Despite mineral compositional variation in a given sample, major constituent minerals in Hole 735B gabbroic rocks display good chemical equilibrium as shown by significant correlations among Mg# (= Mg/[Mg + Fe2+]) of olivine, clinopyroxene, and orthopyroxene and An (=Ca/[Ca + Na]) of plagioclase. This indicates that the mineral assemblages olivine + plagioclase in troctolite, plagioclase + clinopyroxene in gabbro, plagioclases + clinopyroxene + olivine in olivine gabbro, and plagioclase + clinopyroxene + olivine + orthopyroxene in gabbronorite, and so on, have all coprecipitated from their respective parental melts. Fe-Ti oxides (ilmenite and titanomagnetite), which are ubiquitous in most of these rocks, are not in chemical equilibrium with olivine, clinopyroxene, and plagioclase, but precipitated later at lower temperatures. Disseminated oxides in some samples may have precipitated from trapped Fe-Ti-rich melts. Oxides that concentrate along shear bands/zones may mark zones of melt coalescence/transport expelled from the cumulate sequence as a result of compaction or filter pressing. Bulk Hole 735B is of cumulate composition. The most primitive olivine, with Fo = 0.842, in Hole 735B suggests that the most primitive melt parental to Hole 735B lithologies must have Mg# 0.637, which is significantly less than Mg# = 0.714 of bulk Hole 735B. This suggests that a significant mass fraction of more evolved products is needed to balance the high Mg# of the bulk hole. Calculations show that 25%-45% of average Eastern Atlantis II Fracture Zone basalt is needed to combine with 55%-75% of bulk Hole 735B rocks to give a melt of Mg# 0.637, parental to the most primitive Hole 735B cumulate. On the other hand, the parental melt with Mg# 0.637 is far too evolved to be in equilibrium with residual mantle olivine of Fo > 0.89. Therefore, a significant mass fraction of more primitive cumulate (e.g., high Mg# dunite and troctolite) is yet to be sampled. This hidden cumulate could well be deep in the lower crust or simply in the mantle section. We favor the latter because of the thickened cold thermal boundary layer atop the mantle beneath slow-spreading ridges, where cooling and crystallization of ascending mantle melts is inevitable. These observations and data interpretation require reconsideration of the popular concept of primary mantle melts and relationships among the extent of mantle melting, melt production, and the composition and thickness of igneous crust.

Chemistry of altered basaltic glass shards and zeolites of volcaniclastic sediments from ODP Hole 157-953C

Trace element concentrations of altered basaltic glass shards (layer silicates) and zeolites in volcaniclastic sediments drilled in the volcanic apron northeast of Gran Canaria during Ocean Drilling Program (ODP) leg 157 document variable element mobilities during low-temperature alteration processes in a marine environment. Clay minerals (saponite, montmorillonite, smectite) replacing volcanic glass particles are enriched in transition metals and rare earth elements (REE). The degree of retention of REE within the alteration products of the basaltic glass is correlated with the field strength of the cations. The high field-strength elements are preferentially retained or enriched in the alteration products by sorption through clay minerals. Most trace elements are enriched in a boundary layer close to the interface mineral-altered glass. This boundary layer has a key function for the physico-chemical conditions of the subsequent alteration process by providing a large reactive surface and by lowering the fluid permeability. The release of most elements is buffered by incorporation into secondary precipitates (sodium-rich zeolites, phillipsite, Fe- and Mn-oxides) as shown by calculated distribution coefficients between altered glasses and authigenic minerals. Chemical fluxes change from an open to a closed system behavior during prograde low-temperature alteration of volcaniclastic sediments with no significant trace metal flux from the sediment to the water column.

Seawater carbonate chemistry and concentration boundary layers around complex assemblages of macroalgae in a laboratory experiment

Metabolic processes have the potential to modulate the effects of ocean acidification (OA) in nearshore macroalgal beds. We investigated whether natural mixed assemblages of the articulate coralline macroalgae Arthrocardia corymbosa and understory crustose coralline algae (CCA) altered pH and O2 concentrations within and immediately above their canopies. In a unidirectional flume, we tested the effect of water velocity (0-0.1 m/s), bulk seawater pH (ambient pH 8.05, and pH 7.65), and irradiance (photosynthetically saturating light and darkness) on pH and O2 concentration gradients, and the derived concentration boundary layer (CBL) thickness. At bulk seawater pH 7.65 and slow velocities (0 and 0.015 m/s), pH at the CCA surface increased to 7.90-8.00 in the light. Although these manipulations were short term, this indicates a potential daytime buffering capacity that could alleviate the effects of OA. Photosynthetic activity also increased O2 concentrations at the surface of the CCA. However, this moderating capacity was flow dependent; the CBL thickness decreased from an average of 26.8 mm from the CCA surface at 0.015 m/s to 4.1 mm at 0.04 m/s. The reverse trends occurred in the dark, with respiration causing pH and O2 concentrations to decrease at the CCA surface. At all flow velocities the CBL thicknesses (up to 68 mm) were much greater than those previously published, indicating that the presence of canopies can alter the CBL substantially. In situ, the height of macroalgal canopies can be an order of magnitude larger than those used here, indicating that the degree of buffering to OA will be context dependent.

Diffusion boundary layers ameliorate the negative effects of ocean acidification on the temperate coralline macroalga Arthrocardia corymbosa

Anthropogenically-modulated reductions in pH, termed ocean acidification, could pose a major threat to the physiological performance, stocks, and biodiversity of calcifiers and may devalue their ecosystem services. Recent debate has focussed on the need to develop approaches to arrest the potential negative impacts of ocean acidification on ecosystems dominated by calcareous organisms. In this study, we demonstrate the role of a discrete (i.e. diffusion) boundary layer (DBL), formed at the surface of some calcifying species under slow flows, in buffering them from the corrosive effects of low pH seawater. The coralline macroalga Arthrocardia corymbosa was grown in a multifactorial experiment with two mean pH levels (8.05 'ambient' and 7.65 a worst case 'ocean acidification' scenario projected for 2100), each with two levels of seawater flow (fast and slow, i.e. DBL thin or thick). Coralline algae grown under slow flows with thick DBLs (i.e., unstirred with regular replenishment of seawater to their surface) maintained net growth and calcification at pH 7.65 whereas those in higher flows with thin DBLs had net dissolution. Growth under ambient seawater pH (8.05) was not significantly different in thin and thick DBL treatments. No other measured diagnostic (recruit sizes and numbers, photosynthetic metrics, %C, %N, %MgCO3) responded to the effects of reduced seawater pH. Thus, flow conditions that promote the formation of thick DBLs, may enhance the subsistence of calcifiers by creating localised hydrodynamic conditions where metabolic activity ameliorates the negative impacts of ocean acidification.