Modeling of transient groundwater flow,pollutant transport,and biodegradation in an aquifer with large hydraulic head variations

Industrially sourced dense non-aqueous phase liquids (DNAPLs) contaminated an alluvial aquifer in France decades ago. The location(s) and nature of the pollution source zone(s) were unknown, and the dissolved concentrations of volatile organic compounds in the monitoring wells varied greatly with time. The aquifer was in hydraulic equilibrium with an artificial canal whose water level was highly variable (up to 5 m). These variations propagated into the aquifer, causing changes in the groundwater flow direction; a transient numerical model of flow and solute transport showed that they correlate with the concentration variations because the changes in the flow direction resulted in the contaminant plume shifting. The transient hydrogeological numerical model was built, taking into account solvent biodegradation with first-order chain, since biodegradation has a significant influence on the pollutant concentration evolution. The model parameterization confirms the position of the source zones among the potential troughs in the bedrock where DNAPLs could have accumulated. The groundwater model was successfully calibrated to reproduce the observed concentration variations over several years and allowed a rapid validation of the hypotheses on the functioning of the polluted system.     

Investigating the impact of surface wave breaking on modeling the trajectories of drifters in the northern Adriatic Sea during a wind-storm event

An accurate numerical prediction of the oceanic upper layer velocity is a demanding requirement for many applications at sea and is a function of several near-surface processes that need to be incorporated in a numerical model. Among them, we assess the effects of vertical resolution, different vertical mixing parameterization (the so-called Generic Length Scale –GLS– set of k–ε, k–ω, gen, and the Mellor–Yamada), and surface roughness values on turbulent kinetic energy (k) injection from breaking waves.First, we modified the GLS turbulence closure formulation in the Regional Ocean Modeling System (ROMS) to incorporate the surface flux of turbulent kinetic energy due to wave breaking. Then, we applied the model to idealized test cases, exploring the sensitivity to the above mentioned factors. Last, the model was applied to a realistic situation in the Adriatic Sea driven by numerical meteorological forcings and river discharges. In this case, numerical drifters were released during an intense episode of Bora winds that occurred in mid-February 2003, and their trajectories compared to the displacement of satellite-tracked drifters deployed during the ADRIA02-03 sea-truth campaign.Results indicted that the inclusion of the wave breaking process helps improve the accuracy of the numerical simulations, subject to an increase in the typical value of the surface roughness z₀. Specifically, the best performance was obtained using α_CH = 56,000 in the Charnok formula, the wave breaking parameterization activated, k–ε as the turbulence closure model. With these options, the relative error with respect to the average distance of the drifter was about 25% (5.5 km/day). The most sensitive factors in the model were found to be the value of α_CH enhanced with respect to a standard value, followed by the adoption of wave breaking parameterization and the particular turbulence closure model selected.     

Validation of turbulence closure parameterisations for stably stratified flows using the PROVESS turbulence measurements in the North Sea

A number of parameterisations for the simulation of mixing processes in the thermocline are compared and tested against the microstructure data of the PROVESS campaigns, conducted in the northern part of the North Sea during the autumn of 1998. The transport term in the turbulent kinetic energy equation is parameterised via the introduction of a third stability function S_k for turbulent energy diffusion. The formulations are compared with a simpler scheme based upon limiting conditions for turbulence variables. Improved results are obtained with a new form of S_k. The best agreement is, however, found with the simpler limiting scheme. This is explained in terms of a turbulence length scale theory for stably stratified turbulence. In agreement with previous laboratory and ocean data it is found that the ratios of the Thorpe and Kolmogorov scales to the Ozmidov length scale approach critical limiting values in the thermocline. The first of these conditions is satisfied when limiting conditions are implemented into the scheme, providing the necessary minimum value for the dissipation rate, whereas the schemes without limiting conditions fail to produce this critical ratio. The basic reason for this failure is that the Thorpe scale is overestimated, which is shown to be connected to an even larger overprediction of the dissipation rate of temperature variance. To investigate the impact of non-resolved advective processes and salinity stratification on the turbulence predictions, additional numerical experiments were conducted using a simple scheme for data assimilation. The best agreement is found again with the limiting scheme, which is able to make reasonable predictions for the dissipation rate without knowing the detailed shape of the mean stratification profile. It is shown that advective transport due to tidally and wind-driven motions has a non-negligible impact on vertical mixing. This is seen in the data and the models by periodic enhancements of turbulent mixing inside the thermocline.     

Extensive characterization of Italian accelerometric stations from single-station ambient-vibration measurements

This paper describes the analyses of the single-station ambient-vibration measurements performed on the Italian accelerometric network to detect site resonance phenomena potentially affecting earthquake recordings. The use of low cost, high quality microtremor measurement can be helpful to discriminate among soil classes, since several classification schemes based on resonance frequencies were proposed in the last decades. Operatively, in the framework of the Italian Strong Motion Database project (DPC-INGV 2007–2009 S4; ), soil resonance frequencies have been evaluated from more than 200 ambient vibration measurements in correspondence of accelerometric stations included in ITACA (). The noise recordings have been analyzed using the same numerical protocol in order to standardize the results. Particular attention has been paid to evaluate the quality of measurements and to develop an on-purpose mathematical tool to automatically estimate the peaks in the horizontal-to-vertical spectral ratio (HVSR) curve. The reliability of the resonance frequencies from HVSR has been tested by comparing estimates provided by independent methods (modeling or earthquake recordings). The test confirmed the reliability of the microtremor HVSR for assessing the resonance frequencies of the examined sites.     

Use of biomarkers indices in a sediment core to evaluate potential pollution sources in a subtropical reservoir in Brazil

The presence of PAHs, n-alkanes, pristane, and phytanes in core sediment from the Vossoroca reservoir (Parana, southern Brazil) was investigated. The total concentration of the 16 PAHs varied from 15.5 to 1646 μg kg⁻¹. Naphthalene was present in all layers (3.34–74.0 μg kg⁻¹). The most abundant and dominant n-alkanes were n-C₁₅ and n-C₃₆, with average concentrations of 198.1 ± 46.8 and 522.9 ± 167.7 μg kg⁻¹, respectively. Lighter n-alkanes were distributed more evenly through the layers and showed less variation, specially n-C₉, n-C₁₂, and n-C₁₈, with average concentrations of 14.6 ± 3.0, 31.6 ± 1.9, and 95.0 ± 5.2 μg kg⁻¹, respectively; heavier n-alkanes were more unevenly distributed.     

The Radio and X-ray Mode-Switching Pulsar PSR B0943+10

Observations obtained in the last years challenged the widespread notion that rotation-powered neutron stars are steady X-ray emitters. Besides a few allegedly rotation-powered neutron stars that showed ‘magnetar-like’ variability, a particularly interesting case is that of PSR B0943+10. Recent observations have shown that this pulsar, well studied in the radio band where it alternates between a bright and a quiescent mode, displays significant X-ray variations, anticorrelated in flux with the radio emission. The study of such synchronous radio/X-ray mode switching opens a new window to investigate the processes responsible for the pulsar radio and high-energy emission. Here we review the main X-ray properties of PSR B0943+10 derived from recent coordinated X-ray and radio observations.     

The Italian strong motion network

The Italian Strong Motion Network is a permanent monitoring system run by the Italian national emergency management department (Dipartimento della Protezione Civile, DPC). The network is known as RAN (Rete Accelerometrica Nazionale). An extensive project for updating and improving the technology of RAN instruments as well as the number of recording points was performed in the last 10 years. A wide site selection survey was carried out from eastern Sicily along the Italian peninsula, covering high seismic risk areas. The recording station density and the choice of high-quality digital strong motion instruments ensure reliability of the RAN network in the long-term. At the end of 2008, the free field sites selection and instruments installation, planned in the project, were quite completed. In planning and drawing the new RAN, special attention has been devoted to the robustness of the transmission systems, and to the distribution of new stations in order to ensure plenty of data during a seismic emergency. We spent special care both in the estimation of the RAN site responses and in the diffusion of the strong motion data. In order to better identify damaged earthquake areas, improved ground motion parameters need to be set. Such parameters will also assist future progress for engineering seismic design techniques as well as disaster mitigation.