Groundwater chemistry patterns in the phreatic aquifer of the central Belgian coastal plain

The Holocene geological evolution of the Belgian coastal plain is dominated by a transgression of the North Sea, silting up of the coastal plain and human intervention (impoldering). This has led to a typical pattern in groundwater quality which is discussed here for the central part of the coastal plain. Therefore, a database with available groundwater samples is composed. Water type according to the Stuyfzand classification is determined and different hydrosomes and their hydrochemical facies are identified. Based on this, the origin and evolution of the water types is explained using Piper plots and geochemical calculations with PHREEQC. Before the impoldering, salinising and freshening conditions alternated with a general salinisation of the aquifer after about 3400 BP. This results in a dominance of brackish and salt NaCl subtypes which are still found in the deeper part of the aquifer. The subsequent impoldering resulted in an major freshening of the aquifer leading to NaHCO₃, MgHCO₃ and CaHCO₃ subtypes. Overall, mixing, cation exchange, carbonate mineral dissolution and oxidation of organic matter are identified as the major processes determining the general water quality. The close link between geological evolution, water quality and what is still observable today is illustrated with this example of the Belgian coastal plain.     

Spatiotemporal dynamics: the need for an innovative approach in mountain hazard risk management

Starting with an overview on losses due to mountain hazards in the Russian Federation and the European Alps, the question is raised why a substantial number of events still are recorded—despite considerable efforts in hazard mitigation and risk reduction. The main reason for this paradox lies in a missing dynamic risk-based approach, and it is shown that these dynamics have different roots: firstly, neglecting climate change and systems dynamics, the development of hazard scenarios is based on the static approach of design events. Secondly, due to economic development and population dynamics, the elements at risk exposed are subject to spatial and temporal changes. These issues are discussed with respect to temporal and spatial demands. As a result, it is shown how risk is dynamic on a long-term and short-term scale, which has to be acknowledged in the risk concept if this concept is targeted at a sustainable development of mountain regions. A conceptual model is presented that can be used for dynamical risk assessment, and it is shown by different management strategies how this model may be converted into practice. Furthermore, the interconnectedness and interaction between hazard and risk are addressed in order to enhance prevention, the level of protection and the degree of preparedness.     

Horizontal rotation signals detected by “G-Pisa” ring laser for the M w = 9.0, March 2011, Japan earthquake

We report the observation of the ground rotation induced by the M _w = 9.0, 11th of March 2011, Japan earthquake. The rotation measurements have been conducted with a ring laser gyroscope operating in a vertical plane, thus detecting rotations around the horizontal axis. Comparison of ground rotations with vertical accelerations from a co-located force balance accelerometer shows excellent ring laser coupling at periods longer than 100?s. Under the plane wave assumption, we derive a theoretical relationship between horizontal rotation and vertical acceleration for Rayleigh waves. Due to the oblique mounting of the gyroscope with respect to the wave direction of arrival, apparent velocities derived from the acceleration/rotation rate ratio are expected to be always larger than or equal to the true wave propagation velocity. This hypothesis is confirmed through comparison with fundamental mode, Rayleigh-wave phase velocities predicted for a standard Earth model.     

Local ensemble assimilation scheme with global constraints and conservation

Ensemble assimilation schemes applied in their original, global formulation respect linear conservation properties if the ensemble perturbations are set up accordingly. For realistic ocean systems, only a relatively small number of ensemble members can be calculated. A localization of the ensemble increment is therefore necessary to filter out spurious long-range correlations. The conservation of the global properties will be lost if the assimilation is performed locally, since the conservation requires a coupling between all model grid points which is removed by the localization. The distribution of ocean observations is often highly inhomogeneous. Systematic errors of the observed parts of the ocean state can lead to spurious adjustment of the non-observed parts via data assimilation and thus to a spurious increase or decrease in long-term simulations of global properties which should be conserved. In this paper, we propose a local assimilation scheme (with different variants and assumptions) which can satisfy global conservation properties. The proposed scheme can also be used for non-local observation operators. Different variants of the proposed scheme are tested in an idealized model and compared to the traditional covariance localization with an ad-hoc step enforcing conservation. It is shown that the inclusion of the conservation property reduces the total RMS error and that the presented stochastic and deterministic schemes avoiding error space rotation provide better results than the traditional covariance localization.     

The generalized kinetic equation as a model for the nonlinear transfer in third-generation wave models

An alternative model for the nonlinear interaction term S_nl in spectral wave models, the so called generalized kinetic equation (Janssen J Phys Oceanogr 33(4):863–884, 2003; Annenkov and Shrira J Fluid Mech 561:181–207, 2006b; Gramstad and Stiassnie J Fluid Mech 718:280–303, 2013), is discussed and implemented in the third generation wave model WAVEWATCH-III. The generalized kinetic equation includes the effects of near-resonant nonlinear interactions, and is therefore able, in theory, to describe faster nonlinear evolution than the existing forms of S_nl which are based on the standard Hasselmann kinetic equation (Hasselmann J Fluid Mech 12:481–500, 1962). Numerical simulations with WAVEWATCH have been carried out to thoroughly test the performance of the new form of S_nl, and to compare it to the existing models for S_nl in WAVEWATCH; the DIA and WRT. Some differences between the different models for S_nl are observed. As expected, the DIA is shown to perform less well compared to the exact terms in certain situations, in particular for narrow wave spectra. Also for the case of turning wind significant differences between the different models are observed. Nevertheless, different from the case of unidirectional waves where the generalized kinetic equation represents a obvious improvement to the standard forms of S_nl (Gramstad and Stiassnie 2013), the differences seems to be less pronounced for the more realistic cases considered in this paper.     

Tracing offshore low-salinity plumes in the Northeastern Gulf of Mexico during the summer season by use of multispectral remote-sensing data

To trace offshore surface low-salinity water (LSW) in the northeastern Gulf of Mexico, a proxy was developed using the surface water beam attenuation coefficient (c _p), and salinity matched with synchronous Sea-viewing Wide Field-of-view Sensor (SeaWiFS) satellite data from three annual summer cruises (July 1998–August 2000) using a two-step empirical approach. First, a relationship between in-situ salinity and c _p was obtained. Second, in-situ c _p was matched with SeaWiFS radiance ratios of all available blue-to-green wavelengths. Finally, satellite-derived surface salinity was determined directly by combining the two empirical relationships, providing a robust estimate over a range of salinities (26–36). This significantly improves the limited spatial and temporal resolution of surface salinity distribution obtained by shipboard sampling. The resulting correlation is best explained as mixing between low-salinity plume waters and normal salinity Gulf waters. The empirical relationships were used to map satellite-derived salinity using the average of SeaWiFS images during each summer cruise. As expected for summer, spatial patterns of LSW plumes with high c _p, particulate matter (PM), particulate organic carbon (POC), and chlorophyll-a (Chl-a) were connected to the mouth of the Mississippi River Delta and extended to the east-southeast. Normal salinity Gulf water with lower c _p, PM, POC, and Chl-a was confined to the shelf and upper slope in the eastern part of the study area, outside the plumes. This proxy approach can be applied throughout the region of shipboard sampling for more detailed coverage and analysis.     

Possible variation in methane flux caused by gas hydrate formation on the northeastern continental slope off Sakhalin Island, Russia

The Sakhalin Slope Gas Hydrate Project (SSGH) is an international collaborative effort by scientists from Japan, Korea, and Russia to investigate natural gas hydrates (GHs) that have accumulated on the continental slope off Sakhalin Island, Okhotsk Sea. From 2009 to 2011, field operations of the SSGH-09, -10, and -11 projects were conducted. GH-bearing and -free sediment cores were retrieved using steel hydro- and gravity corers. The concentrations of sulfate ions in sediment pore waters were measured to investigate sulfate concentration–depth profiles. Seventeen cores showed linear depth profiles of sulfate concentrations. In contrast, eight cores and two cores showed concave-up and -down profiles plausibly explained by sudden increase and decrease in methane flux from below, respectively, presumably caused by the formation of gas hydrate adjacent to the core sampling sites.     

A meltwater origin for Antarctic shelf bedforms with special attention to megalineations

The geomorphology of troughs crossing the Antarctic shelf is described and interpreted in terms of ice-stream hydrology. The scale of tunnel channels on the inner shelf and the absence of sediment at their mouths are taken to infer catastrophic drainage. Drumlins on the inner and outer shelves with pronounced crescentic and hairpin scours are also interpreted as products of catastrophic flow. Gullies and channels on the continental slope and turbidites on the rise and abyssal plain point to abundant meltwater discharge across the shelf. Attempts to explain this morphology and sedimentology in terms of release or discharge of meltwater by pressure melting, strain heating, Darcian flow, or advection in deforming till are shown to be unrealistic. We suggest that meltwater flow across the middle and outer shelves might have been in broad, turbulent floods, which raises the possibility that megascale glacial lineations (MSGL) on the shelf might originate by erosion in turbulent flow. This possibility is explored by use of analogs for MSGL from flood and eolian landscapes and marine environments. An extended discussion reflects on objections that stand in the way of the flood hypothesis.     

Predicting changes in alluvial channel patterns in North-European Russia under conditions of global warming

Global climate change may have a noticeable impact on the northern environment, leading to changes in permafrost, vegetation and fluvial morphology. In this paper we compare the results from three geomorphological models and study the potential effects of changing climatic factors on the river channel types in North-European Russia. Two of the selected models by Romashin [Romashin, V.V., 1968. Variations of the river channel types under governing factors, Annals of the Hydrological Institute, vol. 155. Hydrometeoizdat, Leningrad, pp. 56–63.] and Leopold and Wolman [Leopold, L.B., Wolman, M.G., 1957. River channel pattern: braided, meandering and straight, Physiographic and hydraulic studies of rivers. USA Geological Survey Professional Paper 252, pp. 85–98.] are conventional QS-type models, which predict the existence of either multi-thread or single-tread channel types using data on discharge and channel slope. The more advanced model by Van den Berg [Van den Berg, J.H., 1995. Prediction of alluvial channel pattern of perennial rivers. Geomorphology 12, 259–270.] takes into account the size of the sediment material.We used data from 16 runoff gauges to validate the models and predict the channel types at selected locations under modern and predicted for the future climatic conditions. Two of the three models successfully replicated the currently existing channel types in all but one of the studied sites. Predictive calculations under the hypothetical scenarios of 10%, 15%, 20% and 35% runoff increase gave different results. Van den Berg's model predicted potential transformation of the channel types, from single- to multi-thread, at 4 of 16 selected locations in the next few decades, and at 5 locations by the middle of the 21st century. Each of the QS-type models predicted such transformation at one site only.Results of the study indicate that climatic warming in combination with other environmental changes may lead to transformation of the river channel types at selected locations in north-western Russia. Further efforts are needed to improve the performance of the fluvial geomorphological models and their ability to predict such changes.