Physical and chemical properties of the waters of saline lakes and their importance for deep-water renewal: Lake Issyk-Kul, Kyrgyzstan

The relationships between electrical conductivity, temperature, salinity, and density are studied for brackish Lake Issyk-Kul. These studies are based on a newly determined major ion composition, which for the open lake shows a mean absolute salinity of 6.06 g kg⁻¹. The conductivity-temperature relationship of the lake water was determined experimentally showing that the lake water is about 1.25 times less conductive than seawater diluted to the same absolute salinity as that of the lake water. Based on these results, an algorithm is presented to calculate salinity from in-situ conductivity measurements. Applied to the field data, this shows small but important vertical salinity variations in the lake with a salinity maximum at 200 m and a freshening of the surface water with increasing proximity to the shores. The algorithm we adopt to calculate density agrees well with earlier measurements and shows that at 20°C and 1 atm Lake Issyk-Kul water is about 530 g m⁻³ denser than seawater at the same salinity. The temperature of maximum density at 1 atm is about 0.15°C lower than that for seawater diluted to the same salinity. Despite its small variations, salinity plays an important role, together with temperature changes, in the static stability and in the production of deep-water in this lake. Changes in salinity may have had important consequences on the mixing regime and the fate of inflowing river water over geological time. Uncharged silicic acid is negligible for the stability of the water column except near an ∼15 m thick nepheloid layer observed at the bottom of the deep basin.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.     

Integrated Two-Dimensional Modeling of Fluid Flow and Compaction in a Sedimentary Basin

In this paper we employ mixed finite elements and numerically study an integrated two-dimensional model of fluid flow and compaction in a sedimentary basin. This model describes a single phase incompressible flow in a two-dimensional section of a sedimentary basin with vertical compaction. At each time step, an iterative algorithm is used to solve this model. The determination of the grid movement is based on the mass conservation and movement of sediments in the basin, while the mixed method is utilized to solve the fluid flow over the moving grid. Numerical experiments are presented to verify this iterative algorithm and show representative solutions for the model under consideration.     

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.     

Horizontal and vertical components of earthquake ground motions at liquefiable sites

Field observations on ground motions from recent earthquakes imply that current knowledge is limited with regard to relating vertical and horizontal motions at liquefiable sites. This paper describes a study with the purpose of clarifying this emerging issue to some extent. A series of numerical analyses is carried out on a liquefiable soil deposit with a verified, fully coupled, nonlinear procedure. It is shown that the transformation of vertical motions in the deposit differs considerably from the transformation of horizontal motions. Both the amplitude and frequency content of the horizontal motions are strongly dependent on the shaking level or the associated nonlinear soil behavior. The transfer function for vertical motions is however likely to be independent of the intensity of input motions; no reduction in the amplitude occurs even in the case of strong shaking. The results are shown to be in consistence with the laboratory observations on shaking table tests and recent field observations that less nonlinearity exists for vertical motions. It is also shown that the possibility exists for using information on spectral ratios between the horizontal and vertical surface motions to quickly identify in situ soil behavior and liquefaction that are not readily covered by conventional field or laboratory experimentation procedures.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

Relationships between tissue concentrations of chlorinated hydrocarbons (polychlorinated biphenyls and chlorinated pesticides) and antioxidative responses of marine mussels, Perna viridis

Marine mussels, Perna viridis, were transplanted from a reference site to various polluted sites around Hong Kong. After 30 d of exposure, antioxidative responses in the gills and hepatopancreas and tissue concentrations of chlorinated hydrocarbons [polychlorinated biphenyls (PCBs) and chlorinated pesticides (CPs)] were determined for individual mussels. Glutathione S transferase (GST) and glutathione (GSH) were positively correlated with tissue PCB concentrations. Only one of the enzymatic antioxidants, glutathione peroxidase (GPx), showed significant response to tissue PCB. No significant correlation was found between tissue concentrations of chlorinated hydrocarbons and other enzymatic antioxidants (superoxide dismutase (SOD), catalase (CAT), glutathione reductase (GR) and NADPH DT-diaphorase (DT-d). Oxidative stress, measured as thiobarbituric acid reactive substances, was correlated with chlorinated pesticide concentrations in tissues. This study demonstrated a correlation between GST/GSH and chlorinated hydrocarbons. The apparent lack of correlation between trace organic pollutants and some of the enzymatic antioxidants may be due to the inhibitory effects caused by these chemicals. The above results suggest that more investigations are needed before these enzymes can be used as biomarkers.     

Determining Q of near-surface materials from Rayleigh waves

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.     

African monsoon variability during the previous interglacial maximum

Little is known about centennial- to millennial-scale climate variability during interglacial times, other than the Holocene. We here present high-resolution evidence from anoxic (unbioturbated) sediments in the eastern Mediterranean Sea that demonstrates a sustained ∼800-yr climate disturbance in the monsoonal latitudes during the Eemian interglacial maximum (∼125 ka BP). Results imply that before and after this event, the Intertropical Convergence Zone (ITCZ) penetrated sufficiently beyond the central Saharan watershed (∼21°N) during the summer monsoon to fuel flooding into the Mediterranean along the wider North African margin, through fossil river/wadi systems that to date have been considered only within a Holocene context. Relaxation in the ITCZ penetration during the intra-Eemian event curtailed this flux, but flow from the Nile - with its vast catchment area - was not affected. Previous work suggests a concomitant Eurasian cooling event, with intensified impact of the higher-latitude climate on the Mediterranean basin. The combined signals are very similar to those described for the Holocene cooling event around 8 ka BP. The apparent type of concurrent changes in the monsoon and higher-latitude climate may reflect a fundamental mechanism for variability in the transfer of energy (latent heat) between the tropics and higher latitudes.