Use of electrical resistivity methods for detecting subsurface fresh and saline water and delineating their interfacial configuration: a case study of the eastern Dead Sea coastal aquifers, Jordan

The alluvial aquifer is the primary source of groundwater along the eastern Dead Sea shoreline, Jordan. Over the last 20 years, salinity has risen in some existing wells and several new wells have encountered brackish water in areas thought to contain fresh water. A good linear correlation exists between the water resistivity and the chloride concentration of groundwater and shows that the salinity is the most important factor controlling resistivity. Two-dimensional electrical tomography (ET) integrated with geoelectrical soundings were employed to delineate different water-bearing formations and the configuration of the interface between them. The present hydrological system and the related brines and interfaces are controlled by the Dead Sea base level, presently at 410 m b.s.l. Resistivity measurements show a dominant trend of decreasing resistivity (thus increasing salinity) with depth and westward towards the Dead Sea. Accordingly, three zones with different resistivity values were detected, corresponding to three different water-bearing formations: (1) strata saturated with fresh to slightly brackish groundwater; (2) a transition zone of brine mixed with fresh to brackish groundwater; (3) a water-bearing formation containing Dead Sea brine. In addition, a low resistivity unit containing brine was detected above the 1955 Dead Sea base level, which was interpreted as having remained unflushed by infiltrating rain.     

Site Effects Study in the Serravalle Village (NE Italy) by Means of Spectral Ratios Analyses Relatively to Multiple Reference Sites

The area of Serravalle, sited in the northern part of the town of Vittorio Veneto (TV), NE Italy, has been the target of a seismic microzonation campaign. 10 seismic stations have been deployed for a 7 months period to record in continuous mode. Three stations were installed on bedrock outcrops and seven on sedimentary sites with variable cover thickness. Spectral analyses have been performed on the collected data-set using the Generalized Inversion Technique (GIT, e.g. Andrews, 1986). In particular, spectral ratios have been calculated for each station relatively to the average of the three reference, bedrock sites. The spectral ratios provide quantitative estimates of the seismic motion amplifications which occur in each of the monitored sites. Two sites show high values of amplification, 5 times larger than signal amplitude at the reference sites, in correspondence of well discernible peak frequencies of 5 Hz. Results for the other stations show smaller amounts of site amplification spreading over a broad range of frequencies. Sites where the highest amplifications were recorded all lie on the left bank of the Meschio River and in areas farther away from its outlet into the plain correlating with the presence of thick layers of Quaternary deposits.     

New data on Sr-and C-isotopic chemostratigraphy of the Upper Riphean type section (Southern Urals)

New data on Sr-and C-isotopic systematics of carbonate rocks from the Upper Riphean stratotype (Karatau Group of the southern Urals) are obtained for several southwestern sections of the Bashkirian meganticlinorium, which have not been studied before. The results obtained supplement the Sr-and C-isotopic information for the group upper horizons thus detailing chemostratigraphic characterization of the entire succession. Limestone and dolostone samples used to analyze the Sr isotope composition satisfy strict geochemical criteria of the isotopic system retentivity and have been subjected to preliminary treatment in ammonium acetate to remove secondary carbonate phases. Data on 255 samples of carbonate rocks (171 studied for the first time) show that δ¹³C value varies in the Karatau Group succession from ?2.8 to +5.9 ‰ V-PDB with several in-phase excursions from the general trend in all the sections studied in the area 90 × 130 km. The δ¹³C variation trend demarcates several levels in the carbonate succession of the Karatau Group suitable for objectives of regional stratigraphy and for C-isotope chemostratigraphic subdivision of the Upper Riphean. The results of Sr isotopic analysis of 121 samples (51 unstudied before) from the Karatau Group imply that rocks in its lower part (the Katav Formation and basal horizon of the Inzer Formation) experienced considerable secondary alterations, while limestones and dolostones of the overlying interval of the group are frequently unaltered. In the “best” samples satisfying geochemical criteria of the isotopic system retentivity, the ⁸⁷Sr/⁸⁶Sr initial ratio increases from 0.70521–0.70532 in the lower Inzer deposits to 0.70611 in the upper Min’yar carbonates, decreasing to <0.70600 near the top of the latter. Above the regional hiatus separating the Min’yar and Uk formation, this ratio grows from 0.70533 to 0.70605–0.70609 in the limestone succession of the last formation.     

New Generation of Probabilistic Seismic Hazard Assessment for the Area Cologne/Aachen Considering the Uncertainties of the Input Data

A new earthquake catalogue for central, northern and northwestern Europe with unified M_w magnitudes, in part derived from chi-square maximum likelihood regressions, forms the basis for seismic hazard calculations for the Lower Rhine Embayment. Uncertainties in the various input parameters are introduced, a detailed seismic zonation is performed and a recently developed technique for maximum expected magnitude estimation is adopted and quantified. Applying the logic tree algorithm, resulting hazard values with error estimates are obtained as fractile curves (median, 16% and 84% fractiles and mean) plotted for pga (peak ground acceleration; median values for Cologne 0.7 and 1.2 m/s² for probabilities of exceedence of 10% and 2%, respectively, in 50 years), 0.4 s (0.8 and 1.5 m/s²) and 1.0 s (0.3 and 0.5 m/s²) pseudoacclerations, and intensity (I₀ = 6.5 and 7.2). For the ground motion parameters, rock foundation is assumed. For the area near Cologne and Aachen, maps show the median and 84% fractile hazard for 2% probability of exceedence in 50 years based on pga (maximum median value about 1.5 m/s²), and 0.4 s (>2 m/s²) and 1.0 s (about 0.8 m/s²) pseudoaccelerations, all for rock. The pga 84% fractile map also has a maximum value above 2 m/s² and shows similarities with the median map for 0.4 s. In all maps, the maximum values fall within the area 6.2–6.3° E and 50.8–50.9° N, i.e., east of Aachen.     

Relations between climatic variability and hydrologic time series from four alluvial basins across the southwestern United States

Hydrologic time series of groundwater levels, streamflow, precipitation, and tree-ring indices from four alluvial basins in the southwestern United States were spectrally analyzed, and then frequency components were reconstructed to isolate variability due to climatic variations on four time scales. Reconstructed components (RCs), from each time series, were compared to climatic indices like the Pacific Decadal Oscillation (PDO), North American Monsoon (NAM), and El Niño-Southern Oscillation (ENSO), to reveal that as much as 80% of RC variation can be correlated with climate variations on corresponding time scales. In most cases, the hydrologic RCs lag behind the climate indices by 1–36 months. In all four basins, PDO-like components were the largest contributors to cyclic hydrologic variability. Generally, California time series have more variation associated with PDO and ENSO than the Arizona series, and Arizona basins have more variation associated with NAM. ENSO cycles were present in all four basins but were the largest relative contributors in southeastern Arizona. Groundwater levels show a wide range of climate responses that can be correlated from well to well in the various basins, with climate responses found in unconfined and confined aquifers from pumping centers to mountain fronts.     

A new approach for the development and management of brackish karst springs

A new approach to the method of artificial upraising of the water outlet point, for management and development of brackish karst springs, uses the MODKARST model. Brackish karst springs simulation can be used to estimate the necessary upraising of the spring elevation, so that sea-water intrusion is blocked. The consequent freshwater loss to the sea, due to this upraising, can also be estimated. The method has been applied to the periodically brackish karst Almiros spring at Heraklion of Crete, Greece. The spring simulation showed that the sea-water intrusion could be prevented through an artificial upraising of the water-outlet point, realized by the construction of a dam. The exact upraising has been estimated. Freshwater loss to the sea because of this upraising has also been estimated. The model could also be used as a tool for the management of the spring. For example, it was used to assess management options for the spring during the depletion period of the year 1994. The best scenario for the development of the spring during this period has been estimated.     

Optimization of SAGD Well Elevation

The application of steam-assisted gravity drainage (SAGD) to recover heavy oil sands is becoming increasingly important in the northern Alberta McMurray Formation because of the vast resources/reserves accessible with this mechanism. Selecting the stratigraphic elevations of SAGD well pairs is a vital decision for reservoir evaluation and planning. The inherent uncertainty in the distribution of geological variables significantly influences this decision. Geostatistical simulation is used to capture geological uncertainty, which is used can be transformed into a distribution of the best possible well pair elevations. A simple exhaustive calculation scheme is used to determine the optimum stratigraphic location of a SAGD well pair where the recovery R is maximized. There are three basic steps to the methodology: (1) model the uncertainty in the top continuous bitumen (TCB) and bottom continuous bitumen (BCB) surfaces, (2) calculate the recovery at all possible elevation increments within the TCB and BCB interval, and (3) identify the elevation that maximizes R. This is repeated for multiple TCB/BCB pairs of surfaces to assess uncertainty. The methodology is described and implemented on a subset of data from the Athabasca Oilsands in Fort McMurray, Alberta.     

Spherical Slepian functions and the polar gap in geodesy

The estimation of potential fields such as the gravitational or magnetic potential at the surface of a spherical planet from noisy observations taken at an altitude over an incomplete portion of the globe is a classic example of an ill-posed inverse problem. We show that this potential-field estimation problem has deep-seated connections to Slepian's spatiospectral localization problem which seeks bandlimited spherical functions whose energy is optimally concentrated in some closed portion of the unit sphere. This allows us to formulate an alternative solution to the traditional damped least-squares spherical harmonic approach in geodesy, whereby the source field is now expanded in a truncated Slepian function basis set. We discuss the relative performance of both methods with regard to standard statistical measures such as bias, variance and mean squared error, and pay special attention to the algorithmic efficiency of computing the Slepian functions on the region complementary to the axisymmetric polar gap characteristic of satellite surveys. The ease, speed, and accuracy of our method make the use of spherical Slepian functions in earth and planetary geodesy practical.