Marine quality assessment by use of benthic species-abundance distributions: a proposed new protocol within the European Union Water Framework Directive

The aim of this study is to develop a new method for classification of marine benthic quality according to the European Union Water Framework Directive. Tolerance values to environmental disturbance were determined in an objective analysis for benthic species along the Swedish west coast by using 4676 samples from 257 stations. Based on a combination of the species tolerance values, abundance and diversity, a benthic quality index (BQI) was calculated for the assessment of environmental status at a particular station. The qualification of BQI was evaluated in relation to known spatial and temporal gradients of disturbance.     

Spectra of poorly resolved ocean data

Standard temperature and pressure sensors on Aanderaa RCM8 current meters have a resolution of 0.024 °C and 0.6 bar, each equal to 1 digital number (value) over a range of 1024. It is shown that an 11-month deep-ocean temperature record using only four values can contain useful spectral information on internal wave motions. This is partially due to the modulation of high-frequency data by non-zero low-frequency (subinertial) variations. This result follows from the comparison of this record with artificial three- and four-value data constructed from temperature records observed in stronger stratified waters nearby. These artificial records show main features of the internal wave band similar to those observed in the original data spectra. Peaks at tidal harmonic frequencies and enhancements at sum-tidal-inertial interaction frequencies are preserved in the artificial data, but overall noise level (and thus the continuum spectral slope) is enhanced with respect to the properly resolved records (using 15 and 100 values). As a demonstration of the stable accuracy of the temperature sensors, the poorly resolved records provided an estimate of mean stratification to within 5% of the estimate using Seabird CTD data.Responsible Editor: Hans Burchard     

A multi-model ensemble of downscaled spatial climate change scenarios for the Dommel catchment, Western Europe

Regional or local scale hydrological impact studies require high resolution climate change scenarios which should incorporate some assessment of uncertainties in future climate projections. This paper describes a method used to produce a multi-model ensemble of multivariate weather simulations including spatial–temporal rainfall scenarios and single-site temperature and potential evapotranspiration scenarios for hydrological impact assessment in the Dommel catchment (1,350 km²) in The Netherlands and Belgium. A multi-site stochastic rainfall model combined with a rainfall conditioned weather generator have been used for the first time with the change factor approach to downscale projections of change derived from eight Regional Climate Model (RCM) experiments for the SRES A2 emission scenario for the period 2071–2100. For winter, all downscaled scenarios show an increase in mean daily precipitation (catchment average change of +9% to +40%) and typically an increase in the proportion of wet days, while for summer a decrease in mean daily precipitation (−16% to −57%) and proportion of wet days is projected. The range of projected mean temperature is 7.7°C to 9.1°C for winter and 19.9°C to 23.3°C for summer, relative to means for the control period (1961–1990) of 3.8°C and 16.8°C, respectively. Mean annual potential evapotranspiration is projected to increase by between +17% and +36%. The magnitude and seasonal distribution of changes in the downscaled climate change projections are strongly influenced by the General Circulation Model (GCM) providing boundary conditions for the RCM experiments. Therefore, a multi-model ensemble of climate change scenarios based on different RCMs and GCMs provides more robust estimates of precipitation, temperature and evapotranspiration for hydrological impact assessments, at both regional and local scale.     

Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

Consistency of observed winter precipitation trends in northern Europe with regional climate change projections

Often it is claimed that the recent changes in northern European climate are at least partly anthropogenic even though a human influence has not yet been successfully detected. Hence we investigate whether the recent changes are consistent with regional climate change projections. Therefore, trends in winter (DJF) mean precipitation in northern Europe are compared to human induced changes as predicted by a set of four regional climate model simulations. The patterns of recent trends and predicted changes match reasonably well as indicated by pattern correlation and the similarity is very likely not random. However, the model projections generally underestimate the recent change in winter precipitation. That is, the signal-to-noise ratio of the anthropogenic precipitation change is either rather low or the presently used simulations are significantly flawed in their ability to project changes into the future. European trends contain large signals related to the North Atlantic Oscillation (NAO), of which a major unknown part may be unrelated to the anthropogenic signal. Therefore, we also examine the consistency of recent and projected changes after subtracting the NAO signal in both the observations and in the projections. It turns out that even after the removal of the NAO signal, the pattern of trends in the observations is similar to those projected by the models. At the same time, the magnitude of the trends is considerably reduced and closer to the magnitude of the change in the projections.     

Modeled seasonality of glacial abrupt climate events

Greenland ice cores, as well as many other paleo-archives from the northern hemisphere, recorded a series of 25 warm interstadial events, the so-called Dansgaard-Oeschger (D-O) events, during the last glacial period. We use the three-dimensional coupled global ocean–atmosphere–sea ice model ECBILT-CLIO and force it with freshwater input into the North Atlantic to simulate abrupt glacial climate events, which we use as analogues for D-O events. We focus our analysis on the Northern Hemisphere. The simulated events show large differences in the regional and seasonal distribution of the temperature and precipitation changes. While the temperature changes in high northern latitudes and in the North Atlantic region are dominated by winter changes, the largest temperature increases in most other land regions are seen in spring. Smallest changes over land are found during the summer months. Our model simulations also demonstrate that the temperature and precipitation change patterns for different intensifications of the Atlantic meridional overturning circulation are not linear. The extent of the transitions varies, and local non-linearities influence the amplitude of the annual mean response as well as the response in different seasons. Implications for the interpretation of paleo-records are discussed.     

Ion composition and chemistry in the coma of Comet 1P/Halley—A comparison between Giotto's Ion Mass Spectrometer and our ion-chemical network

In order to understand the cometary plasma environment it is important to track the closely linked chemical reactions that dominate ion evolution. We used a coupled MHD ion-chemistry model to analyze previously unpublished Giotto High Intensity Ion Mass Spectrometer (HIS-IMS) data. In this way we study the major species, but we also try to match some minor species like the CH_x and the NH_x groups. Crucial for this match is the model used for the electrons since they are important for ion-electron recombination. To further improve our results we included an enhanced density of supersonic electrons in the ion pile-up region which increases the local electron impact ionization. In this paper we discuss the results for the following important ions: C⁺, CH⁺, CH⁺₂, CH⁺₃, N⁺, NH⁺, NH⁺₂, NH⁺₃, NH⁺₄, O⁺, OH⁺, H₂O⁺, H₃O⁺, CO⁺, HCO⁺, H₃CO⁺, and CH₃OH⁺₂. We also address the inner shock which is very distinctive in our MHD model as well as in the IMS data. It is located just inside the contact surface at approximately 4550 km. Comparisons of the ion bulk flow directions and velocities from our MHD model with the data measured by the HIS-IMS give indication for a solar wind magnetic field direction different from the standard Parker angle at Halley's position. Our ion-chemical network model results are in a good agreement with the experimental data. In order to achieve the presented results we included an additional short lived inner source for the C⁺, CH⁺, and CH⁺₂ ions. Furthermore we performed our simulations with two different production rates to better match the measurements which is an indication for a change and/or an asymmetric pattern (e.g. jets) in the production rate during Giotto's fly-by at Halley's comet.     

Polarization Models for Rayleigh Scattering Planetary Atmospheres

We present Monte Carlo simulations for the polarization of light reflected from planetary atmospheres. We investigate dependencies of intensity and polarization on three main parameters: single scattering albedo, optical depth of a scattering layer, and albedo of a Lambert surface underneath. The main scattering process considered is Rayleigh scattering, but isotropic scattering and enhanced forward scattering on haze particles are also investigated. We discuss disk integrated results for all phase angles and radial profiles of the limb polarization at opposition. These results are useful to interpret available limb polarization measurements of solar system planets and to predict the polarization of extra-solar planets as a preparation for VLT/SPHERE. Most favorable for a detection are planets with an optically thick Rayleigh-scattering layer. The limb polarization of Uranus and Neptune is especially sensitive to the vertically stratified methane mixing ratio. From limb polarization measurements constraints on the polarization at large phase angles can be set.     

Late Holocene flood magnitudes in the Lower Rhine river valley and upper delta resolved by a two-dimensional hydraulic modelling approach

Palaeoflood hydraulic modelling is essential for quantifying ‘millennial flood’ events not covered in the instrumental record. Palaeoflood modelling research has largely focused on one-dimensional analysis for geomorphologically stable fluvial settings because two-dimensional analysis for dynamic alluvial settings is time consuming and requires a detailed representation of the past landscape. In this study, we make the step to spatially continuous palaeoflood modelling for a large and dynamic lowland area. We applied advanced hydraulic model simulations (1D–2D coupled set-up in HEC-RAS with 950 channel sections and 108 × 10³ floodplain grid cells) to quantify the extent and magnitude of past floods in the Lower Rhine river valley and upper delta. As input, we used a high-resolution terrain reconstruction (palaeo-DEM) of the area in early mediaeval times, complemented with hydraulic roughness values. After conducting a series of model runs with increasing discharge magnitudes at the upstream boundary, we compared the simulated flood water levels with an inventory of exceeded and non-exceeded elevations extracted from various geological, archaeological and historical sources. This comparison demonstrated a Lower Rhine millennial flood magnitude of approximately 14,000 m³/s for the Late Holocene period before late mediaeval times. This value exceeds the largest measured discharges in the instrumental record, but not the design discharges currently accounted for in flood risk management.