Increasing levels and biomagnification of persistent organic pollutants (POPs) in Antarctic biota

Representatives of the Antarctic food web (krill, cephalopod, fish, penguin, seal) of the area around Elephant Island and from the Weddell Sea were analysed for the most recalcitrant organochlorine compounds. Due to sorption of the compounds to sinking particles and accumulation in sediments, two benthic fish species (Gobionotothen gibberifrons, Chaenocephalus aceratus) feeding on benthos invertebrates and fish reflected significantly increasing concentrations within a decade (1987-1996), while a benthopelagic species (Champsocephalus gunnari) feeding on krill did not. In the pelagic food chain, lipid normalised concentrations of all compounds increased from Antarctic krill to fish proving that biomagnification of highly lipophilic pollutants (log octanol-water partition coefficient>5) occurs in water-breathing animals. As top predators Weddell and southern elephant seals (Leptonychotes weddellii, Mirounga leonina) biomagnified the persistent organic pollutants relative to krill 30-160 fold with the exception of hexachlorobenzene, the levels of which were lower than in fish indicating its intense specific elimination.     

Broadband frequency-dependent amplification of seismic waves across Bucharest, Romania

The determination of seismic amplitude amplification is a fundamental contribution to seismic hazard assessment. While often only high-frequency amplitude variations (>1 Hz) are taken into account, we analyse broadband waveforms from 0.14 to 8.6 Hz using a temporary network of 32 stations in and around the earthquake-prone city of Bucharest. Spectral amplitudes are calculated with an adaptive multiple-taper approach. Across our network (aperture 25 km × 25 km), we find a systematic northwest/southeast-oriented structural influence on teleseismic P-wave amplitudes from 0.14 to 0.86 Hz that can be explained by constructive interference in the dipping Cenozoic sedimentary layers. For higher frequencies (1.4–8.75 Hz), more local site effects prevail and can be correlated partly among neighbouring stations. The transition between systematic and localised amplitude variations occurs at about 1 Hz.     

Flux Ropes Embedded In A Radial Magnetic Field: Analytic Solutions For The External Magnetic Field

Three sympathetic flares were observed with the Solar Magnetic Field Telescope (SMFT) at the Huairou Solar Observing Station of Beijing Astronomical Observatory on 29 August, and 1 September 1990. Each set of sympathetic flares had three ribbons. Two ribbons appeared in active region NOAA 6233 and one ribbon occurred in NOAA 6240 embedded in a single polarity area. Photospheric vector magnetograms were simultaneously obtained from both regions as well. We use a new numerical technique to reconstruct the chromospheric and coronal magnetic fields by making use of the observed vector magnetic fields in the photosphere as boundary conditions. Magnetic field loops linking both regions were identified from the reconstructed 3-D fields. The analysis of chromospheric filtergrams and reconstructed 3-D magnetic fields indicates that interaction between a sheared lower loop in the active region NOAA 6233 and a higher loop linking the two regions resulted in sympathetic flares. The analysis of the time delay between flare ribbons in NOAA 6233 and 6240 indicates that heat conduction along the higher loop from the primary energy release site is responsible for the sympathetic flaring in NOAA 6240. The events reported in this paper represent only one alternative as the cause of sympathetic flaring in which energy transport along coronal interconnecting loops plays the major role, and no in-situ energy release is required.     

Higher-order compositional modeling with Fickian diffusion in unstructured and anisotropic media

We present advances in compositional modeling of two-phase multi-component flow through highly complex porous media. Higher-order methods are used to approximate both mass transport and the velocity and pressure fields. We employ the Mixed Hybrid Finite Element (MHFE) method to simultaneously solve, to the same order, the pressure equation and Darcy's law for the velocity. The species balance equation is approximated by the discontinuous Galerkin (DG) approach, combined with a slope limiter. In this work we present an improved DG scheme where phase splitting is analyzed at all element vertices in the two-phase regions, rather than only as element averages. This approximation is higher-order than the commonly employed finite volume method and earlier DG approximations. The method reduces numerical dispersion, allowing for an accurate capture of shock fronts and lower dependence on mesh quality and orientation. Further new features are the extension to unstructured grids and support for arbitrary permeability tensors (allowing for both scalar heterogeneity, and shear anisotropy). The most important advancement in this work is the self-consistent modeling of two-phase multi-component Fickian diffusion. We present several numerical examples to illustrate the powerful features of our combined MHFE–dg method with respect to lower-order calculations, ranging from simple two component fluids to more challenging real problems regarding CO₂ injection into a vertical domain saturated with a multi-component petroleum fluid.     

Modelling and forecasting long-term dynamics of Western Baltic macrobenthic fauna in relation to climate signals and environmental change

Long-term macrobenthos data from Kiel Bight in the Western Baltic collected between 1968 and 2000 have been correlated with the winter NAO index (North Atlantic Oscillation Index) and other environmental data such as temperature, salinity and oxygen content in the bottom water in order to detect systematic patterns related to so far unexplained abiotic signals in the dynamics of zoobenthic species assemblages. The benthos data come from a cluster of five stations (Süderfahrt/ Millionenviertel) in Kiel Bay. Our investigations concentrated on the macrobenthic dynamics with a focus on the number of species m^− 2 (species richness). Using logarithms and the time series analysis approach of Box/Jenkins (ARIMA modelling, transfer function modelling) it was shown that species richness was strongly influenced by the winter NAO (adjusted for a linear time trend within the 1968-2000 period) and salinity (with a shift/lag of four years). Bootstrapping experiments (i.e. sampling from the error process) and analysis of prediction power (by means of the one- or more-years leaving-out method) showed that the parameter estimates behaved in a stable way, leading to a relatively robust model.     

Small-scale heterogeneities below the Gräfenberg array, Germany from seismic wavefield fluctuations of Hindu Kush events

Small-scale elastic heterogeneities (<5  km) are found in the upper lithosphere underneath the Gräfenberg array, southeast Germany. The results are based on the analysis of broadband recordings of 17 intermediate-depth (201–272  km) events from the Hindu Kush region. The wavefront of the first P arrival and the following 40  s coda are separated into coherent and incoherent (scattered) parts in the frequency range from 0.05 to 5  Hz. The frequency-dependent intensities of the mean and fluctuation wavefields are used to describe the scattering characteristics of the lithosphere underneath the receivers. It is possible to discriminate a weak-fluctuation regime of the wavefield in the frequency range below approximately 1.5–2.5  Hz and a strong-fluctuation regime starting at 2.0–2.5  Hz and continuing to higher frequencies. In order to explain the observed wavefield fluctuations, an approach with seismic scattering at random media-type structures is proposed. The preferred model contains heterogeneities with 3–7 per cent perturbations in seismic velocity and correlation lengths of 0.6–4.8  km in the crust. This is compatible with models from active seismic experiments. Scattering in the lithospheric mantle is not required, but cannot be excluded at weak velocity contrasts (<3 per cent).     

Spatially distributed hydrotope‐based modelling of evapotranspiration and runoff in mountainous basins

River basins in mountainous regions are characterized by strong variations in topography, vegetation, soils, climatic conditions and snow cover conditions, and all are strongly related to altitude. The high spatial variation needs to be considered when modelling hydrological processes in such catchments. A complex hydrological model, with a great potential to account for spatial variability, was developed and applied for the hourly simulation of evapotranspiration, soil moisture, water balance and the runoff components for the period 1993 and 1994 in 12 subcatchments of the alpine/pre‐alpine basin of the River Thur (area 1703 km²). The basin is located in the north‐east of the Swiss part of the Rhine Basin and has an elevation range from 350 to 2500 m a.s.l. A considerable part of the Thur Basin is high mountain area, some of it above the tree‐line and a great part of the basin is snow covered during the winter season. In the distributed hydrological model, the 12 sub‐basins of the Thur catchment were spatially subdivided into sub‐areas (hydrologically similar response units—HRUs or hydrotopes) using a GIS. Within the HRUs a hydrologically similar behaviour was assumed. Spatial interpolations of the meteorological input variables wereemployed for each altitudinal zone. The structure of the model components for snow accumulation and melt, interception, soil water storage and uptake by evapotranspiration, runoff generation and flow routing are briefly outlined. The results of the simulated potential evapotranspiration reflect the dominant role of altitudinal change in radiation and albedo of exposure, followed by the influence of slope. The actual evapotranspiration shows, in comparison with the potential evapotranspiration, a greater variability in the lower and medium altitudinal zones and a smaller variability in the upper elevation zones, which was associated with limitations of available moisture in soil and surface depression storages as well as with the evaporative demand of the local vegetation. The higher altitudinal dependency and variability of runoff results from the strong increase in precipitation and the decrease in evaporation with increased altitude. An increasing influence of snow cover on runoff as well as evapotranspiration with altitude is obvious. The computed actual evapotranspiration and runoff were evaluated against the observed values of a weighting lysimeter and against runoff hydrographs. Copyright © 1999 John Wiley & Sons, Ltd.     

Oxidation of CO on surface hematite in high CO2 atmospheres

We propose a mechanism for the oxidation of gaseous CO into CO₂ occurring on the surface mineral hematite (Fe₂O₃(s)) in hot, CO₂-rich planetary atmospheres, such as Venus. This mechanism is likely to constitute an important source of tropospheric CO₂ on Venus and could at least partly address the CO₂ stability problem in Venus’ stratosphere, since our results suggest that atmospheric CO₂ is produced from CO oxidation via surface hematite at a rate of 0.4 petagrammes (Pg) CO₂ per (Earth) year on Venus which is about 45% of the mass loss of CO₂ via photolysis in the Venusian stratosphere. We also investigated CO oxidation via the hematite mechanism for a range of planetary scenarios and found that modern Earth and Mars are probably too cold for the mechanism to be important because the rate-limiting step, involving CO(g) reacting onto the hematite surface, proceeds much slower at lower temperatures. The mechanism may feature on extrasolar planets such as Gliese 581c or CoRoT-7b assuming they can maintain solid surface hematite which, e.g. starts to melt above about 1200 K. The mechanism may also be important for hot Hadean-type environments and for the emerging class of hot Super-Earths with planetary surface temperatures between about 600 and 900 K.     

A comparison of chemical and chemodynamical models

Representative results from a comparison of the chemical evolution of spherical collapse models without and with a intercloud medium are presented. The hot metal-rich gas distributes quickly the metals produced in supernovae throughout the galaxy, thus leading to a more homogeneous chemical evolution and flatter metallicity gradients in the gas and the stars. The stellar population is somewhat less concentrated towards the centre. The strong outflow results in a substantial loss of metals from the galaxy to its surroundings, and a lower effective yield in the galaxy. This revised version was published online in August 2006 with corrections to the Cover Date.