Water mass variation in the Mediterranean and Black Seas

The mass-induced sea level variability and the net mass transport between Mediterranean Sea and Black Sea are derived for the interval between August 2002 and July 2008 from satellite-based observations and from model data. We construct in each basin two time series representing the basin mean mass signal in terms of equivalent water height. The first series is obtained from steric-corrected altimetry while the other is deduced from GRACE data corrected for the contamination by continental hydrology. The series show a good agreement in terms of annual and inter-annual signals, which is in line with earlier works, although different model corrections influence the consistency in terms of seasonal signal and trend.In the Mediterranean Sea, we obtain the best agreement using a steric correction from the regional oceanographic model MFSTEP and a continental hydrological leakage correction derived from the global continental hydrological model WaterGAP2. The inter-annual time series show a correlation of 0.85 and a root mean square (RMS) difference of 15 mm. The two estimates have similar accuracy and their annual amplitude and phase agree within 3 mm and 23 days respectively. The GRACE-derived mass-induced sea level variability yields an annual amplitude of 27 ± 5 mm peaking in December and a trend of 5.3 ± 1.9 mm/yr, which deviates within 3 mm/yr from the altimetry-derived estimate.In the Black Sea, the series are less consistent, with lower accuracy of the GRACE-derived estimate, but still show a promising agreement considering the smaller size of the basin. The best agreement is realized choosing the corrections from WaterGAP2 and from the regional oceanographic model NEMO. The inter-annual time series have a correlation and RMS differences of 0.68 and 55 mm, their annual amplitude and phase agree within 4 mm and 6 days respectively. The GRACE-derived seawater mass signal has an annual amplitude of 32 ± 4 mm peaking in April. On inter-annual time scales, the mass-induced sea level variability is stronger than in the Mediterranean Sea, with an increase from 2003 to 2005 followed by a decrease from 2006 to 2008.Based on mass conservation, the mass-induced sea level variations, river runoff and precipitation minus evaporation are combined to derive the strait flows between the basins and with the Atlantic Ocean. At the Gibraltar strait, the net inflow varies annually with an amplitude of 52 ± 10 × 10⁻³ Sv peaking end of September (1 Sv = 10⁶ m³ s⁻¹). The inflow through the Bosphorus strait displays an annual amplitude of 13 ± 3 ×10⁻³ Sv peaking in the middle of March. Additionally, an increase of the Gibraltar net inflow (3.4 ± 0.8 × 10⁻³ Sv/yr) is detected.     

Os isotope heterogeneity of the upper mantle: Evidence from the Mayarí–Baracoa ophiolite belt in eastern Cuba

In an attempt to quantify the extent of geochemical heterogeneity within a restricted and well dated portion of the upper mantle, 27 chromite separates from the 90 My old chromite deposits in the Mayarí–Baracoa ophiolite belt in eastern Cuba have been investigated for platinum group element (PGE) concentrations and Re–Os isotopic systematics. The samples are characterized by systematically subchondritic initial ¹⁸⁷Os/¹⁸⁸Os ratios and substantial heterogeneity. The initial ¹⁸⁷Os/¹⁸⁸Os ratios vary with chromite chemistry and with geographical distribution, reflecting differences in the Os isotopic evolution for the different upper mantle sections represented by the ophiolite. Accordingly, the Os isotope data might be divided into three groups. In the Moa–Baracoa district, where the chromite bodies are located in the mantle–crust transition zone, the calculated initial γOs values average − 0.97 ± 0.69 (n = 13). In the Sagua de Tanamo district, where chromite chemistry is highly variable and their location in relation the mantle sequence is less clear, the initial γOs values are intermediate, with an average of − 1.77 ± 0.80 (n = 7). In the Mayarí district, where the chromite bodies are located in the lower part of the mantle sequence, initial γOs values average − 2.66 ± 0.29 (n = 7). These subchondritic (i.e. negative) initial γOs values are most simply explained by Re depletion during ancient partial melting and/or melt percolation events.The Os isotope heterogeneity documented here indicates a high degree of geochemical complexity on small to intermediate length scales in the upper mantle. Our results, in combination with data on chromites from the literature, show that an “average present-day Os isotopic composition” for the hypothetical depleted MORB mantle (DMM) reservoir cannot be precisely established beyond the statement that it is “broadly chondritic”. Indeed, the upper mantle cannot be considered a sufficiently homogeneous geochemical “reservoir” to serve meaningfully as a baseline against which geochemical “anomalies” are evaluated. On the other hand, our findings are consistent with the “Statistical Upper Mantle Assemblage” or “SUMA”-concept, according to which a high level of geochemical heterogeneity is maintained in the upper mantle at all relevant length scales, as a result of the plate-tectonic cycle and intra-mantle processes such as melt-migration and metasomatism.     

Nighttime E region plasma irregularities over an equatorial station Trivandrum

In situ measurements of electron density were made over Trivandrum (8.5°N, 76.9°E) during nighttime to study E-region plasma density irregularities. Irregularities, with vertical scale sizes from a few km to 15 cm, were detected during rocket ascent and descent. Electron density profiles during ascent and descent of an earlier nighttime rocket flight from Trivandrum are also presented. Some of the important results are as follows: (i) horizontal gradients in electron density exist in 110–120 km region with horizontal scale size of at least 40 km, (ii) based on the presence/absence of electron density structures during ascent and descent of both flights, the horizontal distance over which the gradient drift instability operates is found to be at least 80 km and 90 km, for both the flights, (iii) observed irregularities in regions of negative density gradient are suggested to be produced through the gradient drift instability (GDI) driven by vertical polarization electric field as well as by electric field produced through wind shears and those in positive gradient regions by wind driven GDI, (iv) largest irregularity amplitude (≈30%) was associated with steepest gradients and so was the presence of smallest vertical scale sizes (12 m to 15 cm), which were absent at other altitudes, (v) the spectral index of irregularities was in the range of ?2.2±0.2 for large scales (few kilometers>λ>50 m), ?3.25±0.25 for medium scales (50 m>λ>10 m) and ?2.6±0.1 for smaller scales (10 m>λ>1 m) and (vi) irregularities in large and medium scales are expected to be produced directly through GDI and the small and sum-meter scales through non-linear GDI.     

Temporal and spatial complexity in post-glacial sedimentation on the tectonically active,Poverty Bay continental margin of New Zealand

On the eastern Raukumara Ranges of the New Zealand East Coast, active tectonics, vigorous weather systems, and human colonisation have combined to cause widespread erosion of the mudstone- and sandstone-dominated hinterland. The Waipaoa River sedimentary dispersal system is an example that has responded to environmental change, and is now New Zealand's second largest river in terms of suspended sediment discharge. This paper presents new sediment accumulation rates for the continental shelf and slope that span century to post-glacial time scales. These data are derived from radiochemical tracer, palynological, tephrostratigraphic, and seismic methods. We hypothesise on the temporal and spatial complexity of post-glacial sedimentation across the margin and identify the broad extent of sediment dispersal from the Waipaoa system. The ∼15 km³ Poverty Bay mid-shelf basin lies adjacent to the mouth of the Waipaoa River, reaching a maximum thickness of ∼45 m. A post-glacial mud lobe of an additional ∼3 km³ extends through the Poverty Gap and out onto the uppermost slope, attaining 40 m thickness in a structurally controlled sub-basin. Here, an offset in the last-glacial erosion surface indicates that deposition was sympathetic with fault activity and the creation of accommodation space, implying that sedimentation was not supply limited. Contrary to classical shelf sedimentation models, the highest modern accumulation rate of 1 cm y⁻¹ occurs on the outer-shelf sediment lobe, approximately ∼2 times the rate recorded at the mid-shelf basin depocentre, and ∼10 times faster than the excess ²¹⁰Pb rates estimated from the slope. Pollen records from slope cores fingerprint Polynesian then European settlement, and broaden the spatial extent of post-settlement sedimentation initially documented from the Poverty Bay mid-shelf. Changes in sub-millennial sedimentation infer a 2–3-times increase in post-settlement accumulation on the shelf but a smaller 1–2 times increase on the slope. Over longer time scales, seismic evidence infers slower but steady sedimentation since the last transgression, and that significant cross-shelf sediment pathways pre-date the increase in sedimentation resulting from colonisation and deforestation. From a summation of coastal bedload, shelf and slope sediment mass accumulation, the total sediment budget for the Holocene is ∼1 Mt y⁻¹. Under modern conditions a larger proportion of the Waipaoa sediment dispersal system likely extends onto the slope and beyond.     

A similarity approach to determine response times to steady-state saturation in landscapes

Relationships describing response times for landscape saturation and subsurface flow for idealised hillslopes after a change in water balance are derived in terms of similarity parameters given by their topographic, soil and climatic attributes. The study was carried out under a set of specific assumptions.The work quantitatively describes effects of each of these attributes on travel time. For example, divergent hillslopes can exhibit travel times that are double those of convergent hillslopes, and concave slopes tend to have lower travel times than planar or convex slopes. It is shown that the hillslope travel time T is dependent on hillslope length L, hydraulic conductivity K and slope S such that T = f(L/KS, B, CR, Δq/smd) and if saturation does not occur in the lower reaches of a hillslope, this can be simplified to T = f(L/KS, B). Here, Δq is change in net external flux and smd is soil moisture deficit; CR and B are the convergence ratio and profile factor respectively. The relative importance of these attributes on response time is discussed.It is shown that subsurface flow hydrographs of simple hillslopes of different scale, derived numerically, can be collapsed into a single curve by normalising them by means of a form of response time, the time constant.Results from the analytical derivation were compared with numerical analysis with good agreement. The theory was applied to a natural catchment using observed data from a bushfire event. Comparison of analytically calculated response time with observed response time for the event showed a large discrepancy. The reason is discussed.     

Atmospheric polychlorinated biphenyl (pcb) inputs to a coastal city near the marmara sea

Recent studies have shown that polychlorinated biphenyl (PCB) levels are high in coastal sites. Various types of deposition events, including dry deposition, air–water exchange and wet deposition, were analyzed as part of the study to illustrate the pollution level. The atmospheric levels and deposition mechanisms of 82 PCBs were investigated in Mudanya, a coastal city in Turkey. A total of 175 samples and 112 blanks were collected. Air–water exchange and particle phase dry deposition fluxes were 120 ± 90 ng/m²-d and 125 ± 160 ng/m²-d, respectively. The wet-dry deposition sampler (WDDS) consisted of two reservoirs, wet and dry. Collected with the WDDS, the dry deposition flux was 21 ± 20 ng/m²-d, and the wet deposition fluxes (during the rainy period) were 145 ± 130 ng/m²-d and 195 ± 270 ng/m²-d for the dissolved and particulate phases, respectively. Mass transfer coefficients (MTCs), dry deposition velocities and washout ratios were calculated and compared with the literature values.     

Quantification of the gas mass emitted during single explosions on Stromboli with the SO2 imaging camera

We performed measurements using an SO₂ imaging camera of the SO₂ gas mass emitted during five discrete explosive events on Stromboli volcano on 3 October 2006. The SO₂ gas mass released during discrete explosions was 15–40 kg per explosion, producing 3–8% of the total daily SO₂ gas emission, demonstrating that in terms of gas flux Strombolian explosions are a second-order phenomenon compared with quiescent degassing. Using the typical gas composition measured with OP-FTIR allows us to determine the total gas mass released during an explosion as 360–960 kg with a volume of 1500–4100 m³ at 1 bar. At the probable source pressure of gas slug formation of 75 MPa this gas amount would occupy a volume equivalent to a sphere with a radius of 0.8–1 m, comparable with estimates of Stromboli's conduit geometry.     

The GGOS as the backbone for global observing and local monitoring: A user driven perspective

A key geodetic contribution to both the three Global Observing Systems and initiatives like the European Global Monitoring for Environment and Security is an accurate, long-term stable, and easily accessible reference frame as the backbone. Many emerging scientific as well as non-scientific high-accuracy applications require access to an unique, technique-independent reference frame decontaminated for short-term fluctuations due to global Earth system processes. Such a reference frame can only be maintained and made available through an observing system such as the Global Geodetic Observing System (GGOS), which is currently implemented and expected to provide sufficient information on changes in the Earth figure, its rotation and its gravity field. Based on a number of examples from monitoring of infrastructure, point positioning, maintenance of national references frames to global changes studies, likely future accuracy requirements for a global terrestrial reference frame are set up as function of time scales. Expected accuracy requirements for a large range of high-accuracy applications are less than 5 mm for diurnal and sub-diurnal time scales, 2–3 mm on monthly to seasonal time scales, better than 1 mm/year on decadal to 50 years time scales. Based on these requirements, specifications for a geodetic observing system meeting the accuracy requirements can be derived.     

On the visibility of airborne volcanic ash and mineral dust from the pilot’s perspective in flight

In April 2010, volcanic ash from the Eyjafjalla volcano in Iceland strongly impacted aviation in Europe. In order to prevent a similar scenario in the future, a threshold value for safe aviation based on actual mass concentrations was introduced (2 mg m⁻³ in Germany). This study contrasts microphysical and optical properties of volcanic ash and mineral dust and assesses the detectability of potentially dangerous ash layers (mass concentration larger than 2 mg m⁻³) from a pilot’s perspective during a flight. Also the possibility to distinguish between volcanic ash and other aerosols is investigated. The visual detectability of airborne volcanic ash is addressed based on idealized radiative transfer simulations and on airborne observations with the DLR Falcon gathered during the Eyjafjalla volcanic ash research flights in 2010 and during the Saharan Mineral Dust Experiments in 2006 and 2008. Mineral dust and volcanic ash aerosol both show an enhanced coarse mode (>1 μm) aerosol concentration, but volcanic ash aerosol additionally contains a significant number of Aitken mode particles (<150 nm) not present in mineral dust. Under daylight clear-sky conditions and depending on the viewing geometry, volcanic ash is visible already at mass concentrations far below what is currently considered dangerous for aircraft engines. However, it is not possible to visually distinguish volcanic ash from other aerosol layers or to determine whether a volcanic ash layer is potentially dangerous (mass concentration larger or smaller than 2 mg m⁻³). Different appearances due to microphysical differences of both aerosol types are not detectable by the human eye. Nonetheless, as ash concentrations can vary significantly over distances travelled by an airplane within seconds, this visual threat evaluation may contribute greatly to the short-term response of pilots in ash-contaminated air space.     

Is the distribution of the lancelet Branchiostoma caribaeum affected by sewage discharges? An analysis at multiple scales of variability

Spatial variation in the density and biomass of Branchiostoma caribaeum was analyzed along a sewage contamination gradient identified by fecal steroids in a subtropical estuary, southern Brazil. Sampling, repeated in the austral winter and summer, followed a hierarchical design nested at four spatial scales (sector > 1 km; area > 100 m; site > 10 m; replicate < 1 m). Density and biomass were significantly lower at sites characterized by high concentrations of fecal steroids. The best combinations of variables that explained the biological similarities among sites involved contamination indicators. Most of the variation of biological data was found at the smallest scales and could be related with the sediment texture. Our study highlighted the usefulness of a multi-scale perspective to evaluate distribution patterns of benthic invertebrates as a biological indication of environmental pollution. Gradient analyses at larger spatial scales may be invalidated by the patchy distribution of benthic fauna if they do not account for such small scale variability.     

Observations of narrow bipolar events in East China

Narrow bipolar events (NBEs) are a distinct class of intra-cloud lightning discharge. In this paper we present observations of 10 negative and 67 positive such events in East China. Positive NBEs occurred at 7–12 km altitude above mean sea level (MSL) with a mean altitude of 9.5 km, and negative NBEs occurred at 14–16 km altitude. Electrical/channel characteristics of these events were derived from NBE pulse waveforms based on the transmission-line model. On average, the peak current moment and the charge moment change of a NBE event is 15 kA km, and 0.12 C km, respectively. The mean time for the propagation of current front along the channel is 2.2 μs. The upper limit on channel length for NBEs in this study is 510–1060 m, the lower limit on discharge current amplitude is 12.5–43.2 kA, and the minimum charge transfer is 0.1–0.3 C.     

40Ar/39Ar geochronology using a quadrupole mass spectrometer

Magnetic sector mass spectrometers dominate the field of ⁴⁰Ar/³⁹Ar geochronology. Recent advances in quadrupole mass spectrometer technology, especially improvements in resolution, have increased the performance of these instruments to the extent that they can be used for isotopic determinations. We describe a triple filter quadrupole mass spectrometer (Hiden HAL 3F Series Pulse Ion Counting Triple Filter QMS) linked to an automated furnace extraction and cleaning system dedicated to ⁴⁰Ar/³⁹Ar incremental heating experiments.The instrument produces peaks with broad flat tops and a width of 0.9 amu at 10 cps height and 0.84 amu at 500,000 cps height on a 1 million cps high peak (peak width at 0.01‰ and 50% peak height respectively). This allows measurement of ratios of the main Ar peaks in the 1‰ range. Measurements of 1.6 × 10^?12 mole of air reference gas over two years yields ⁴⁰Ar/³⁶Ar = 257.9 ± 1.3 (1σ, n = 34). The ability of the instrument to produce ⁴⁰Ar/³⁹Ar ages from rocks/minerals of a wide age range, reaching into the late Quaternary, are demonstrated by a series of tests and comparison with geochronological data from other studies and an in-house MAP 215-50 magnetic sector mass spectrometer. We demonstrate that high-end quadrupole systems can be used for routine ⁴⁰Ar/³⁹Ar dating purposes.     

Earthquake scaling,fault segmentation,and structural maturity

Slip and length measurements on earthquakes suggest large stress drop variability. We analyze an extended set of slip-length measurements for large earthquakes (M ≥ 6) to seek for the possible origin(s) of this apparent variability. We propose that such variability arises from earthquakes breaking a variable number of major fault segments. That number depends on the strength of the inter-segment zones, which itself depends on the structural maturity of the faults. We propose new D_max–L parameterizations based on that idea of multiple segment-ruptures. In such parameterizations, each broken segment roughly scales as a crack, while the total multi-segment rupture does not. Stress drop on individual segments is roughly constant, only varying between 3.5 to 9 MPa. The slight variation that is still observed depends on fault structural maturity; more mature faults have lower stress drops than immature ones. The new D_max–L functions that we propose reduce uncertainties with respect to available relationships. They thus provide a more solid basis to estimate seismic hazard by integrating fault properties revealed by geological studies.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Temporal and spatial variations of seismicity scaling behavior in Southern México

R/S analysis is used in this work to investigate the fractal correlations in terms of the Hurst exponent for the 1998–2011 seismicity data in Southern Mexico. This region is the most seismically active area in Mexico, where epicenters for severe earthquakes (e.g., September 19, 1985, M_w = 8.1) causing extensive damage in highly populated areas have been located. By only considering the seismic events that meet the Gutenberg–Ritcher law completeness requirement (b = 0.97, M_GR = 3.6), we found time clustering for scales of about 100 and 135 events. In both cases, a cyclic behavior with dominant spectral components at about one cycle per year is revealed. It is argued that such a one-year cycle could be related to tidal effects in the Pacific coast. Interestingly, it is also found that high-magnitude events (M_w ≥ 6.0) are more likely to occur under increased interevent correlations with Hurst exponent values H > 0.65. This suggests that major earthquakes can occur when the tectonic stress accumulates in preferential directions. In contrast, the high-magnitude seismic risk is reduced when stresses are uniformly distributed in the tectonic shell. Such cointegration between correlations (i.e., Hurst exponent) and macroseismicity is confirmed for spatial variations of the Hurst exponent. In this way, we found that, using the Hurst exponent standpoint, the former presumed Michoacan and the Guerrero seismic gaps are the riskiest seismic zones. To test this empirical finding, two Southern Mexico local regions with large earthquakes were considered. These are the Atoyac de Alvarez, Guerrero (M_w = 6.3), and Union Hidalgo, Oaxaca (M_w = 6.6), events. In addition, we used the Loma Prieta, California, earthquake (October 17, 1989, M_w = 6.9) to show that the high-magnitude earthquakes in the San Andreas Fault region can also be linked to the increments of determinism (quantified in terms of the Hurst exponent) displayed by the stochastic dynamics of the interevent period time series. The results revealed that the analysis of seismic activity by means of R/S analysis could provide further insights in the advent of major earthquakes.     

Scattering coefficients in the mantle revealed from the seismogram envelope analysis based on the multiple isotropic scattering model

At longer periods, scattered ScS waves sometimes dominate over coda waves at large lapse times. Examining recordings of seismic envelopes at 9 IRIS seismic stations of regional earthquakes with focal depths deeper than 150 km in periods from 1 to 20 s for a wide lapse time range up to 2000 s, we found significant frequency dependence. The coda decay gradient at short periods is steeper than that at longer periods; however, the change of coda gradient associated with the ScS arrival becomes distinct as the period becomes longer. In particular, a clear offset of coda amplitude appears in central Asia for 10 and 15 s period bands. The multiple isotropic scattering process of S-waves in the heterogeneous mantle can be simply simulated by using the Monte Carlo simulation method based on the radiative transfer theory in scattering media. Assuming a two-plane-layer attenuation structure and smoothed velocity model of the PREM, we estimated the average total scattering coefficients of S-waves such as 7.52 × 10^− 4∼1.32 × 10^− 3 km^− 1 and 2.08 × 10^− 4∼6.23 × 10^− 4 km^− 1 at 4 s, and 4.51 × 10^− 4∼7.37 × 10^− 4 km^− 1 and 2.80 × 10^− 5∼2.71 × 10^− 4 km^− 1 at 10 s, for the lithosphere and the upper mantle and for the lower mantle, respectively. Our results indicate that scattering occurs mostly in the lithosphere and the upper mantle and support that medium heterogeneity spreads over the whole mantle though its scattering power is small. Strong scattering occurs beneath central Asia and Papua New Guinea, whereas the scattering beneath Italy and regions of east Russia is much weaker. The numerical calculation enables us to confirm that much stronger scattering than intrinsic attenuation causes the offset behavior with coda decay gradient change after the ScS arrival for 4 and 10 s period bands in some regions.     

Effects of starvation on bacterial transport through porous media

A major problem preventing widespread implementation of microbial injection strategies for bioremediation and/or microbially enhanced oil recovery is the tendency of bacteria to strongly adhere to surfaces in the immediate vicinity of the injection point. Long term (weeks to months) nutrient starvation of bacteria prior to injection can decrease attachment and enhance transport through porous media. This paper summarizes results of starvation-enhanced transport experiments in sand columns of 30 cm, 3 m, and 16 m in length. The 16 m column experiments compared transport, breakthrough and distribution of adhered cells for starved and vegetative cultures of Klebsiella oxytoca, a copious biofilm producer. Results from these experiments were subsequently used to design and construct a field-scale biofilm barrier using starved Pseudomonas fluorescens. The 30 cm and 3 m sand columns experiments investigated starvation-enhanced transport of Shewanella algae BrY, a dissimilatory metal-reducing bacterium. In both cases the vegetative cells adsorbed onto the sand in higher numbers than the starved cells, especially near the entrance of the column. These results, taken together with studies cited in the literature, indicate that starved cells penetrate farther (i.e. higher breakthrough concentration) and adsorb more uniformly along the flow path than vegetative cells.     

Some aspects of vanadium and chromium chemistry in the English Channel

Dissolved and particulate vanadium and chromium concentrations along a transect between Cherbourg and the Isle of Wight were investigated in the English Channel. Seawater samples were collected at two different depths (surface and bottom) at five sampling stations during five cruises carried out between September 1994 and September 1995. Calculated mass flows through this Channel section were 7600 T yr⁻¹ for vanadium (about 66% was in the dissolved phase) and 1650 T yr⁻¹ for chromium (about 50% was in the dissolved phase). Dissolved chromium concentrations do not vary significantly along the transect. Seasonal variations in chromium distribution linked to biotic parameters were noted during the September 1994 cruise, when a significant relationship between particulate chromium and algal organic carbon was found. In addition, dissolved Cr (III) and Cr (VI) were well correlated with both algal and terrestrial organic matter. In the May cruise, during phytoplanktonic decay, particulate chromium was correlated with the detritic and bacterial organic fractions. These observations suggest interactions between chromium and biotic material. During the winter period, no relationship between chromium and POC was found. Dissolved and particulate vanadium concentrations varied, respectively, between 15 and 28 nmol l⁻¹ and 2–32 nmol l⁻¹. Values of dissolved vanadium showed depletion in the Channel with respect to oceanic waters. This loss in dissolved vanadium was higher near the English coast, but was compensated for by the increase in the particulate vanadium. There was no clear relationship between dissolved vanadium and algal organic carbon and it is inferred that the vanadium transfer cannot be explained by trapping with biotic material. On the other hand, the dissolved vanadium depletion could be attributed to the presence of ferric oxyhydroxide phases in particles, which have strong adsorbing properties for a range of dissolved metal ions.     

Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

Experimental and numerical evaluation of the effectiveness of a stiff wave barrier in the soil

This paper discusses the design, the installation, and the experimental and numerical evaluation of the effectiveness of a stiff wave barrier in the soil as a mitigation measure for railway induced vibrations. A full scale in situ experiment has been conducted at a site in El Realengo (Spain), where a barrier consisting of overlapping jet grout columns has been installed along a railway track. This barrier is stiff compared to the soil and has a depth of 7.5 m, a width of 1 m, and a length of 55 m. Geophysical tests have been performed prior to the installation of the barrier for the determination of the dynamic soil characteristics. Extensive measurements have been carried out before and after installation of the barrier, including free field vibrations during train passages, transfer functions between the track and the free field, and the track receptance. Measurements have also been performed at a reference section adjacent to the test section in order to verify the effect of changing train, track, and soil conditions over time. The in situ measurements show that the barrier is very effective: during train passages, a reduction of vibration levels by 5 dB is already obtained from 8 Hz upwards, while a peak reduction of about 12 dB is observed near 30 Hz immediately behind the barrier. The performance decreases further away from the jet grouting wall, but remains significant. The experimental results are also compared to numerical simulations based on a coupled finite element–boundary element methodology. A reasonable agreement between experiments and predictions is found, largely confirming the initially predicted reduction. This in situ test hence serves as a ‘proof of concept׳, demonstrating that stiff wave barriers are capable of significantly reducing vibration levels, provided that they are properly designed.