Sensitivity analysis of sediment flux derived by laser diffraction and acoustic backscatter within a reservoir

A combination of two indirect methods to measure sediment flux is presented in this study to evaluate suspended sediment transport in a hydropower reservoir. The acoustic backscatter signal (ABS) from an Acoustic Doppler Current Profiler (ADCP) is therefore applied in pre-defined transects within the reservoir in combination with a Laser In-Situ Scattering Transmissometry – stream lined device (LISST-SL). The stationary LISST-SL derived suspended sediment concentration (SSC) measurements are used to calibrate the ABS. From the LISST-SL measurements a time series of SSC is obtained. This enables, in addition, a comprehensive data analysis to evaluate the influence of natural fluctuations of the SSC on the calculated sediment flux, which should be taken into account when assessing sediment transport. Furthermore SSC measurements are done with the LISST-SL close to the reservoir bed. In areas close to the bed no information regarding the ABS is available from the ADCP measurements due to the side-lobe interference. In various studies the information from the last three valid cells is used for extrapolation. However, as result of a comparison of the LISST-SL measurements with extrapolated SSC values from the ADCP measurements it can be seen that, especially in deep reservoirs, this method has to be adapted to the in-situ conditions.     

Mesoscale and submesoscale mechanisms behind asymmetric cooling and phytoplankton blooms induced by hurricanes: a comparison between an open ocean case and a continental shelf sea case

Right-side bias in both sea surface cooling and phytoplankton blooms is often observed in the wake of hurricanes in the Northern Hemisphere. This idealized hurricane modeling study uses a coupled biological-physical model to understand the underlying mechanisms behind hurricane-induced cooling and phytoplankton bloom asymmetry. Both a deep ocean case and a continental shelf sea case are considered and contrasted. Model analyses show that while right-side asymmetric mixing due to inertial oscillations and restratification from strong right-side recirculation cells contributes to bloom asymmetry in the open ocean, the well-mixed condition in the continental shelf sea inhibits formation of recirculation cells, and the convergence of water onto the shelf is a more important process for bloom asymmetry.     

Water exchange between the inner and outer archipelago areas of the Finnish Archipelago Sea in the Baltic Sea

Ocean Dynamics - We studied the water age and transport of passive tracers in the Archipelago Sea, Baltic Sea, using the COHERENS 3D hydrodynamic model and the OpenDrift Lagrangian particle model....     

Surface-wave analysis for static corrections in mineral exploration: A case study from central Sweden

In mineral exploration, increased interest towards deeper mineralizations makes seismic methods attractive. One of the critical steps in seismic processing workflows is the static correction, which is applied to correct the effect of the shallow, highly heterogeneous subsurface layers, and improve the imaging of deeper targets. We showed an effective approach to estimate the statics, based on the analysis of surface waves (groundroll) contained in the seismic reflection data, and we applied it to a legacy seismic line acquired at the iron-oxide mining site of Ludvika in Sweden. We applied surface-wave methods that were originally developed for hydrocarbon exploration, modified as a step-by-step workflow to suit the different geologic context of hard-rock sites. The workflow starts with the detection of sharp lateral variations in the subsurface, the existence of which is common at hard-rock sites. Their location is subsequently used, to ensure that the dispersion curves extracted from the data are not affected by strong lateral variations of the subsurface properties. The dispersion curves are picked automatically, windowing the data and applying a wavefield transform. A pseudo-2D time-average S-wave velocity and time-average P-wave velocity profile are obtained directly from the dispersion curves, after inverting only a reference curve. The time-average P-wave velocity profile is then used for the direct estimation of the one-way traveltime, which provides the static corrections. The resulting P-wave statics from the field data were compared with statics computed through conventional P-wave tomography. Their difference was mostly negligible with more than 91% of the estimations being in agreement with the conventional statics, proving the effectiveness of the proposed workflow. The application of the statics obtained from surface waves provided a stacked section comparable with that obtained by applying tomostatics.     

Importance of Spatial Resolution in Global Groundwater Modeling

Global-scale gradient-based groundwater models are a new endeavor for hydrologists who wish to improve global hydrological models (GHMs). In particular, the integration of such groundwater models into GHMs improves the simulation of water flows between surface water and groundwater and of capillary rise and thus evapotranspiration. Currently, these models are not able to simulate water table depth adequately over the entire globe. Unsatisfactory model performance compared to well observations suggests that a higher spatial resolution is required to better represent the high spatial variability of land surface and groundwater elevations. In this study, we use New Zealand as a testbed and analyze the impacts of spatial resolution on the results of global groundwater models. Steady-state hydraulic heads simulated by two versions of the global groundwater model G³M, at spatial resolutions of 5 arc-minutes (9 km) and 30 arc-seconds (900 m), are compared with observations from the Canterbury region. The output of three other groundwater models with different spatial resolutions is analyzed as well. Considering the spatial distribution of residuals, general patterns of unsatisfactory model performance remain at the higher resolutions, suggesting that an increase in model resolution alone does not fix problems such as the systematic overestimation of hydraulic head. We conclude that (1) a new understanding of how low-resolution global groundwater models can be evaluated is required, and (2) merely increasing the spatial resolution of global-scale groundwater models will not improve the simulation of the global freshwater system.     

Simulating Groundwater-Streamflow Connections in the Upper Colorado River Basin

In mountain, snow driven catchments, snowmelt is supposed to be the primary contribution to river streamflows during spring. In these catchments the contribution of groundwater is not well documented because of the difficulty to monitor groundwater in such complex environment with deep aquifers. In this study we use an integrated hydrologic model to conduct numerical experiments that help quantify the effect of lateral groundwater flow on total annual and peak streamflow in predevelopment conditions. Our simulations focus on the Upper Colorado River Basin (UCRB; 2.8 × 10⁵ km²) a well-documented mountain catchment for which both streamflow and water table measurements are available for several important sub-basins. For the simulated water year, our results suggest an increase in peak flow of up to 57% when lateral groundwater flow processes are included—an unexpected result for flood conditions generally assumed independent of groundwater. Additionally, inclusion of lateral groundwater flow moderately improved the model match to observations. The correlation coefficient for mean annual flows improved from 0.84 for the no lateral groundwater flow simulation to 0.98 for the lateral groundwater flow one. Spatially we see more pronounced differences between lateral and no lateral groundwater flow cases in areas of the domain with steeper topography. We also found distinct differences in the magnitude and spatial distribution of streamflow changes with and without lateral groundwater flow between Upper Colorado River Sub-basins. A sensitivity test that scaled hydraulic conductivity over two orders of magnitude was conducted for the lateral groundwater flow simulations. These results show that the impact of lateral groundwater flow is as large or larger than an order of magnitude change in hydraulic conductivity. While our results focus on the UCRB, we feel that these simulations have relevance to other headwaters systems worldwide.     

Pay the PIED Piper: Guidelines to Visualize Large Geochemical Datasets on Piper Diagrams

The Piper diagram has increased in popularity since its 1944 introduction and is now one of the most familiar and effective tools in the hydrogeologist's toolbox. Within the Piper diagram, three points on three related plots fully display the major ionic species of a water sample. Recently the size and availability of datasets have increased as additional field measurements and modeling results are shared more effectively in online databases. This growth presents opportunities and challenges for data analysis and conveyance—larger and longer datasets increase the potential to identify trends and patterns, but traditional Piper diagrams are quickly overwhelmed by large datasets as dense points overlap and become obscured. We present guidelines for effectively displaying large geochemical datasets on traditional Piper diagrams and new code that adds novel functionality for following these generic guidelines. This code, plotting interesting environmental data with Piper diagrams (PIED Piper), can be run within the Matlab environment or through a stand-alone graphical user interface, and is the first Matlab code to generate Piper diagrams. The illustrative examples herein demonstrate (1) how limitations in displays of large datasets may be overcome with translucent symbology, contours, and heatmaps to identify trends and patterns, (2) how clusters of similar points can be identified and differentiated with convex hulls, and (3) how temporal-and-spatial patterns may be visually diagnosed with image groups and movies. The guidelines discussed in these examples will aid PIED Piper users to achieve the two goals of effective big data visualization: analysis and communication.     

Application of amplitude spectrum and genetic inversion to determine shaly facies distribution: a case study: Murzuq Basin,Libya

Late Ordovician glacial deposits of the Mamuniyat Formation are the main oil reservoir in the Murzuq Basin in Libya. Autopicking the strong reflection at the base of the Silurian shales can be used to map the top of the Mamuniyat reservoir in the area where it is in direct contact with the Silurian shales. However, in areas where the Bir Tlacsin Formation, a mud‐prone unit, is between the Silurian shales and the Mamuniyat reservoir, the top of the Mamuniyat is difficult to pick because the units juxtaposed across the boundary are too similar to produce a strong reflection. Defining the Bir Tlacsin facies is important because it impacts hydrocarbon accumulation and migration. To predict the distribution of the shaly facies of Bir Tlacsin and enhance mapping of the top Mamuniyat reservoir, we utilized a continuous wavelet transform to identify the distinctive thickness of the Hot Shale and Bir Tlacsin units. We also used genetic inversion to distinguish the bulk density of the Bir Tlacsin facies. A 64 Hz frequency gave good time resolution to the amplitude spectrum and was used to predict the facies distribution of the Bir Tlacsin. In contrast, the 24 Hz frequency showed good frequency resolution of the amplitude spectrum and was used to estimate the temporal thickness of the non‐reservoir unit of Bir Tlacsin and Hot Shale. That estimate was then used to modify the autopick horizon for the base of the Silurian reflector to approximate the top of the Mamuniyat reservoir. Because of the large density contrast between the shaly facies of the Bir Tlacsin and the underlying and overlying units, inverted density also provides a way to predict the distribution of the Bir Tlacsin through estimated temporal thickness and to enhance mapping of the top Mamuniyat reservoir through mapping the base of the inverted density of the Bir Tlacsin. A comparison between mapping of the top reservoir using spectral decomposition and inverted density with respect to autopick shows that both methods improved the top of the Mamuniyat reservoir mapping. Prediction of the presence of Bir Tlacsin and improved accuracy of the top of the Mamuniyat reservoir mapping reduce the risk of drilling the shaly facies of Bir Tlacsin and provide a better estimate of the reservoir reserve.     

Storage and Preservation of Artificial Sweeteners in Groundwater Samples

Compound‐specific standardized sampling and storage methods are not available for artificial sweeteners found in groundwater. This study aimed to understand: (1) the appropriate length of storage time for samples containing acesulfame (ACE), sucralose (SUC), saccharin (SAC), and cyclamate (CYC) in simulated groundwater (SGW); (2) conditions of their stability; and (3) which sampling materials are appropriate for sample collection. The evaluated storage conditions included acidification, headspace, exposure to light, and refrigeration; the evaluated sampling materials included steel, stainless steel, aluminum, polyvinyl chloride, polyamide (nylon), polypropylene (PharMed BPT?) tubing, styrene‐ethylene‐butylene co‐polymer (MasterFlex?) tubing, and polytetrafluoroethylene (Teflon?) tubing. All compounds evaluated were stable in storage at 4 °C for 241 d (8 months). Concentrations of artificial sweeteners were consistently within 60% to 120% of original concentrations, except ACE and SAC that were substantially lower under acidified conditions at 25 °C after 241 d. Artificial sweetener concentrations remained nearly constant while in contact with all sampling materials except steel. SEM and TEM images showed oxidation of steel occurred; moreover, removal of all artificial sweetener compounds from aqueous solution had occurred after 289 d. These results suggest artificial sweetener analyses conducted within 14 d of sample collection produce optimal results; however, longer storage times may be acceptable under certain conditions. The results also suggest concentrations of artificial sweeteners in SGW are not affected by contact with typical well casing, sampling, and storage materials, with the exception of steel. The findings from this study will improve the use of artificial sweeteners as tracers in environmental studies.     

The 7 and 11 May 1984 earthquakes in Abruzzo-Latium (Central Italy): reappraisal of the existing macroseismic datasets according to the EMS98

The aim of this paper is to provide a complete and reliable macroseismic knowledge of the events that stroke a large area in Central Italy on 7 and 11 May 1984. Previous studies, together with original accounts integrated with new and unpublished information, have been gathered and examined in order to re-evaluate macroseismic intensities in terms of the European Macroseismic Scale (EMS98). New intensity maps have been compiled; the total number of localities with available information for both the shocks increases from 1254 of the previous study to 1576. On the basis of the new dataset, the macroseismic magnitude of the first shock is M_W 5.6 which is lower than the previous macroseismic computation (M_W 5.7). Moreover, the topic of assessing macroseismic intensity in the presence of multiple shocks has been also investigated, proposing an unconventional approach to presenting the macroseismic data: an overall picture of the cumulative effects produced by all the seismic sequence is given to support a partial but faithful reconstruction of the second shock. This approach is inspired by the common experience in interpreting historical seismic sequences and gives a picture of the impact of the 1984 events on the territory.