Effects of geology and human activity on the dynamics of salt-water intrusion in three coastal aquifers in southern Spain

 Mathematical modelling of salt-water intrusion processes in three aquifers on the southern coast of Spain (Río Verde, Río Vélez and Castell de Ferro) reveals that, although all three systems are subject to the same climate and seasonal over-exploitation, geological and human factors have very different effects on the dynamics of contamination. In the Río Verde aquifer, the most important influence is the high volume of extractions occurring during the dry season; in Río Vélez, the intrusion is strongly controlled by infiltration of water from the river to the aquifer, and, in the Castell de Ferro system, an intensely karstified carbonate massif lying in contact with both the sea water and the detrital aquifer represents the main entrance point for influx of sea water and subsequent washing of the aquifer. We have undertaken a mathematical simulation of various possible measures to counteract intrusion, according to the specific characteristics of the process in each aquifer. These measures include artificial recharge, use of natural recharge from the river as a hydraulic barrier, and the construction of a low-permeability barrier. Received: 5 December 1995 · Accepted: 12 April 1996     

Frequency distributions of major elements (Pb,Zn, Cu) in the main polymetallic deposits of the Baia Mare district (Romania)

The frequency distributions of Pb, Zn, Cu in the hydrothermal sulphide deposits of the Baia Mare district are established. It is shown that the distribution of Pb is similar to that of Zn, unlike that of Cu. This law appears to be valid whatever the composition and grade of the ore, the type of the host rock, and the type of hydrothermal alterations. This law can provide arguments to establish the genesis of mineralizations of unknown origin.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Seismic gaps and plate tectonics: Seismic potential for major boundaries

The theory of plate tectonics provides a basic framework for evaluating the potential for future great earthquakes to occur along major plate boundaries. Along most of the transform and convergent plate boundaries considered in this paper, the majority of seismic slip occurs during large earthquakes, i.e., those of magnitude 7 or greater. The concepts that rupture zones, as delineated by aftershocks, tend to abut rather than overlap, and large events occur in regions with histories of both long- and short-term seismic quiescence are used in this paper to delineate major seismic gaps.In detail, however, the distribution of large shallow earthquakes along convergent plate margins is not always consistent with a simple model derived from plate tectonics. Certain plate boundaries, for example, appear in the long term to be nearly aseismic with respect to large earthquakes. The identification of specific tectonic regimes, as defined by dip of the inclined seismic zone, the presence or absence of aseismic ridges and seamounts on the downgoing lithospheric plate, the age contrast between the overthrust and underthrust plates, and the presence or absence of back-arc spreading, have led to a refinement in the application of plate tectonic theory to the evaluation of seismic potential.The term seismic gap is taken to refer to any region along an active plate boundary that has not experienced a large thrust or strike-slip earthquake for more than 30 years. A region of high seismic potential is a seismic gap that, for historic or tectonic reasons, is considered likely to produce a large shock during the next few decades. The seismic gap technique provides estimates of the location, size of future events and origin time to within a few tens of years at best.The accompanying map summarizes six categories of seismic potential for major plate boundaries in and around the margins of the Pacific Ocean and the Caribbean, South Sandwich and Sunda (Indonesia) regions for the next few decades. These categories range from what we consider high to low potential for being the site of large earthquakes during that period of time. Categories 1, 2 and 6 define a time-dependent potential based on the amount of time elapsed since the last large earthquake. The remaining categories, 3, 4, and 5, are used for areas that have ambiguous histories for large earthquakes; their seismic potential is inferred from various tectonic criteria. These six categories are meant to be interpreted as forecasts of the location and size of future large shocks and should not be considered to be predictions in which a precise estimate of the time of occurrence is specified.Several of the segments of major plate boundaries that are assigned the highest potential, i.e., category 1, are located along continental margins, adjacent to centers of population. Some of them are hundreds of kilometers long. High priority should be given to instrumenting and studying several of these major seismic gaps since many are now poorly instrumented. The categories of potential assigned here provide a rationale for assigning prorities for instrumentation, for future studies aimed at predicting large earthquakes and for making estimates of tsunami potential.Lamont-Doherty Geological Observatory Contribution No. 2906.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

Interpretations of new features of time domain electric-field structures in the auroral acceleration region

Possible theoretical interpretations of the various nonlinear electric-field structures in the auroral acceleration region are provided.     

A storm and tidally-influenced prograding shoreline—Upper Cretaceous Milk River Formation of Southern Alberta, Canada

The Santonian-Campanian Milk River Formation of Southern Alberta represents the transition from an open shelf, through a storm-dominated shoreface into a non-marine sequence of shales and sandstones, with coal. The open shelf deposits consist of interbedded bioturbated mudstones with sharp-based hummocky cross-stratified sandstones. There are no indications of fairweather reworking of the sandstones, which are therefore interpreted as having been deposited below fairweather wavebase. The shoreface sequence consists of a 28 m thick sandstone. It has a very sharp, loaded base, and is dominated by swaley cross-stratification, a close relative of hummocky cross-stratification. Angle of repose cross-bedding is preserved in scattered patches only in the top 5 m of the sand body. Channels up to 180 m wide and 7 m deep are cut into this sand body, with channel margins characterized by lateral accretion surfaces. Regional dispersal trends, as well as local palaeocurrent readings suggest flow toward the NW. Within the channels there is some herringbone cross-bedding and at least two examples of neap-spring bundle cycles, suggesting that the channels are tidally-influenced. Above the channels there is a sequence of carbonaceous shales with in situ root casts and lignitic coal seams. No marine, brackish or lagoonal fauna was identified, and the sequence appears to represent a distal floodplain. The sequence from interbedded hummocky cross-stratified sandstones and bioturbated mudstones into a 10–20 m thick, sharp-based shoreface sandstone characterized by swaley cross-stratification is uncommon. The scarcity or absence of angle of repose cross-bedding in the shoreface, and the dominance of swaley cross-stratification suggests that the shoreface was so storm-dominated that almost no fairweather record was preserved. Other examples of swaley cross-stratified shorefaces are reviewed in the paper.     

青藏高原边缘积雪辐射特征

利用中国陆地生态系统通量观测研究网络的玛曲站观测的一次降雪过程的资料,对青藏高原东部边缘冬季的降雪、积雪过程的辐射特征进行了分析.研究结果表明;积雪期晴天和降雪过程的向上短波辐射的峰值分别约为降雪前晴天的3和2倍.无积雪晴天地表反射率主要分布在0.175～0.36,新雪地表反射率主要分布在0.8～0.9.大气逆辐射变化较小,降雪过程的最大,积雪时的最小.地表长波辐射则为降雪前最大,降雪时最小.积雪覆盖的晴天比无积雪时的净辐射变化幅度减小,且早上由负转正的时间推迟.     

Effects of chemical environment on the solubility and incorporation mechanism for hydrogen in olivine

To investigate the solubility and the sites of incorporation of hydrogen in olivine as a function of point defect concentration, two-stage high-temperature annealing experiments have been carried out. The first annealing stage (the dry preannealing stage) was conducted at a total pressure of 0.1 MPa, a temperature of 1300° C and various oxygen fugacities in the range 10^?11–10^?4 MPa for times > 12 h. In these heat treatments, the samples were buffered against either orthopyroxene or magnesiowustite, or they remained unbuffered. The second annealing stage (the hydrothermal annealing stage) was performed at 300 MPa and 900–1050 ° C under a hydrogen fugacity of ～ 158 MPa for 1–5 h. Infrared spectra from the annealed samples revealed two distinct groups of bands. Group I bands occurred at wavenumbers in the range 3450–3650 cm^?1, while Group II bands occurred in the range 3200–3450 cm^?1. The hydrogen solubility associated with Group I bands is proportional to f _{O 2} to the 1/6 power for samples preannealed in contact with orthopyroxene, to the 1/3 power for samples preannealed in contact with magnesiowustite, and to the 1/13 power for samples preannealed in the absence of a solid-state buffer. The hydrogen concentration for Group II bands varies with f _{o 2} to the 1/3 power for opxbuffered samples, to the 1/2 power for mw-buffered samples, and to the 1/3 power for unbuffered samples. The dependence of hydrogen solubility on oxygen fugacity and orthopyroxene activity suggests that hydrogen is incorporated into the olivine structure via association with point defects. The presence of two distinct groups of absorption bands indicates that hydrogen is associated with two distinct lattice defects. The following point defect model for the mechanism of incorporation of hydrogen in olivine is consistent with these results: Hydrogen ions responsible for the Group I bands are associated with doubly charged oxygen interstitials, while hydrogen ions responsible for the Group II bands are associated with singly charged oxygen interstitials. Furthermore, the infrared bands observed in naturally derived olivines are present in spectra from our hydrothermally annealed crystals. Thus, the mechanisms of incorporation of hydrogen in olivine under geological conditions are the same as those operative under laboratory conditions. The maximum solubility reached in these experiments was ～ 360H/10⁶Si, which corresponds to ～ 0.002 wt% of H₂O. This value is a lower bound for the solubility of hydrogen in olivine under upper mantle conditions.