On solution of the inverse problem for confined aquifer flow via maximum likelihood

Joint estimation of transmissivity (T) and storativity (S) in a confined aquifer is done via maximum likelihood (ML). The differential equation of groundwater flow is discretized by the finite-element method, leading to equation _t+x _t=u _t. Elements of matrices and , as well as estimated covariance matrix of noise termu _t, are functions of T and S. By minimizing the negative loglikelihood function corresponding to discretized groundwater flow equation with respect to T and S, ML estimators are obtained. The ML approach is found to yield accurate estimates of T and S (within 9 and 10% of their actual values, respectively) and showed quadratic convergence in Newton's search technique. Prediction of aquifer response, using ML estimators, results in estimated piezometric heads accurate to ±0.5 m from their actual, exact values. Statistical properties of ML estimators are derived and some basic results for statistical inference are given.     

Three million years of monsoon variability over the northern Sahara

We present a 3 million year record of aeolian dust supply into the eastern Mediterranean Sea, based on hematite contents derived from magnetic properties of sediments from Ocean Drilling Program Site 967. Our record has an average temporal resolution of 400 years. Geochemical data validate this record of hematite content as a proxy for the supply of aeolian dust from the Sahara. We deduce that the aeolian hematite in eastern Mediterranean sediments derives from the eastern Algerian, Libyan, and western Egyptian lowlands located north of the central Saharan watershed (21°N). In corroboration of earlier work, we relate dust flux minima to penetration of the African summer monsoon front to the north of the central Saharan watershed. This would have enhanced soil humidity and vegetation cover in the source regions, in agreement with results from green Sahara climate models. Our results indicate that this northward monsoon penetration recurred during insolation maxima throughout the last 3 million years. As would be expected, this orbital precession-scale mechanism is modulated on both short (100-kyr) and long (400-kyr) eccentricity time scales. We also observe a strong expression of the 41-kyr (obliquity) cycle, which we discuss in terms of high- and low-latitude mechanisms that involve Southern Hemisphere meridional temperature contrasts and shifts in the latitudes of the tropics, respectively. We also observe a marked increase in sub-Milankovitch variability around the mid-Pleistocene transition (0.95 Ma), which suggests a link between millennial-scale climate variability, including monsoon dynamics, and the size of northern hemisphere ice sheets.     

Stable Eocene Magnetization Carried by Magnetite and Iron Sulphides in Marine Marls (Pamplona-Arguis Formation,Southern Pyrenees,Northern Spain)

In order to establish the magnetic carriers and assess the reliability of previous paleomagnetic results obtained for Eocene marine marls from the south Pyrenean basin, we carried out a combined paleo- and rock-magnetic study of the Pamplona-Arguis Formation, which crops out in the western sector of the southern Pyrenees (N Spain). The unblocking temperatures suggest that the characteristic remanent magnetization (ChRM) is carried by magnetite and iron sulphides. The ChRM has both normal and reversed polarities regardless of whether it resides in magnetite or iron sulphides, and represents a primary Eocene magnetization acquired before folding. Rock magnetic results confirm the presence of magnetite and smaller amounts of magnetic iron sulphides, most likely pyrrhotite, in all the studied samples. Framboidal pyrite is ubiquitous in the marls and suggests that iron sulphides formed during early diagenesis under sulphate-reducing conditions. ChRM directions carried by magnetic iron sulphides are consistent with those recorded by magnetite. These observations suggest that magnetic iron sulphides carry a chemical remanent magnetization that coexists with a remanence residing in detrital magnetite. We suggest that the south Pyrenean Eocene marls are suitable for magnetostratigraphic and tectonic purposes but not for studies of polarity transitions, secular variations and geomagnetic excursions, because it is difficult to test for short time differences in remanence lock-in time for the two minerals. The presence of iron sulphide minerals contributing to the primary magnetization in Eocene marine marls reinforces the idea that these minerals can persist over long periods of time in the geological record.     

A paleoenvironmental reconstruction of Laguna Babícora, Chihuahua, Mexico based on ostracode paleoecology and trace element shell chemistry

Paleoecology of Laguna Babícora, Chihuahua, Mexico was reconstructed using ostracode faunal assemblages and shell chemistry. The paleolimnological record is used to show the magnitude of paleoclimatic changes in the area from 25,000 years to the present.Faunal assemblages consist of four species of the genus Limnocythere: L. sappaensis, L. ceriotuberosa, L. bradburyi and L. platyforma, all associated with Candona caudata, Candona patzcuaro and Cypridopsis vidua. A paleosalinity index developed from these assemblages indicates that the lake's salinity fluctuated frequently from oligo- to meso-haline conditions during the last 25,000 years. This pattern and low salinity range are in good agreement with modern TDS (here used as an indicator of salinity) values recorded from 26 wells and one spring from the area (258–975 mg l^–1). To estimate paleotemperature we examined the trace element content (Mg/Ca ratios) from individual valves of L. ceriotuberosa and L. platyforma, the two species most commonly recorded in Laguna Babícora.Shell Mg/Ca ratios of 204 specimens of these two species were used to estimate water temperature (Mg/Ca) by means of experimental standard coefficients. Our data show that paleowater temperature ranged from 5.6–21.3 °C (with 2 values ranging from 0.2–4.8 °C), which suggest a close correlation with atmospheric temperatures around the lake. These results are in good agreement with a modern mean winter temperature (3.5 °C) and mean summer temperature (20 °C) recorded in the area between 1970 and 1980.     

Radio Emission from Three-dimensional Relativistic Hydrodynamic Jets: Observational Evidence of Jet Stratification

High-resolution, high signal-to-noise ratio, blue-violet spectra of three red giant branch tip stars in M15 have been obtained with the Keck I High-Resolution Echelle Spectrograph. These spectra have been analyzed to determine the abundances of several neutron-capture elements, including the radioactive chronometer element thorium. There are two principal results of this study. First, the abundances of the heavier (Z>/=56) elements for each of the three stars is well matched by a scaled solar system r-process abundance distribution. Second, a weighted mean-observed Th/Eu ratio for the stars implies an age for the neutron-capture material in M15 stars of 14+/-3 Gyr, in reasonable agreement with other recent age estimates for Galactic globular clusters.     

An audio-magnetotelluric investigation in Terceira Island (Azores)

Ten audio-magnetotelluric soundings have been carried out along a profile crossing the Serra do Cume caldera in the eastern part of the Terceira Island (Azores). The main objectives of this investigation were to detect geoelectrical features related with tectonic structures and to characterize regional hydrological and hydrothermal aspects mainly those related to geothermal fluid dynamics.Three-dimensional numerical investigation showed that the data acquired at periods shorter than 1 s are not significantly affected by ocean effect. The data was analysed using the Smith's decomposition method in order to investigate possible distortions caused by superficial structures and to estimate a global regional strike. The results suggest that in general the soundings were not distorted. A regional N55°W strike was chosen for the two-dimensional data inversion.The low-resistivity zones (10–30 ohm-m) displayed in the central part of the 2-D geoelectrical model have been interpreted as caused by hydrothermal circulation. The low-resistivity anomalies at the ends of the profile might be attributed to alteration zones with interaction of seawater intrusion. High-resistivity (> 300 ohm-m) values have been related with less permeable zones in the SW of Cinco Picos and Guilherme Moniz caldera walls.     

Trophic Assessment in Chinese coastal systems-review of methods and application to the Changjiang (Yangtze) Estuary and Jiaozhou Bay

Coastal eutrophication has become one of the main threats to Chinese coastal areas during the last two decades. High nutrient loads from human activities have modified the natural background water quality in coastal water bodies, resulting in a range of undesirable effects. There is a need to assess the eutrophic level in coastal systems and to identify the extent of this impact to guide development of appropriate management efforts. Traditional Chinese assessment methods are discussed and compared with other currently-used methods, such as the Oslo-Paris Convention for the Protection of the North Sea (OSPAR) Comprehensive Procedure and Assessment of Estuarine Trophic Status (ASSETS). The ASSETS method and two Chinese methods were tested on two Chinese systems: the Changjiang (Yangtze) Estuary and Jiaozhou Bay. ASSETS is process based, and uses a pressure-state-response model based on three main indices: Influencing Factors, Overall Eutrophic Condition, and Future Outlook. The traditional methods are based on a nutrient index. ASSETS was successfully applied to both systems, classifying the Changjiang Estuary as Bad (high eutrophication) and Jiaozhou Bay as High (low eutrophication). The traditional methods led to ambiguous results, particularly for Jiaozhou Bay, due to the high spatial variability of data and a failure to assess the role of shellfish aquaculture in nutrient control. An overview of the Chinese coastal zone identifies 50 estuaries and bays that should form part of a national assessment. A comparison of methods and results suggests that ASSETS is a promising tool for evaluating the eutrophication status of these systems.     

Homogeneity of the Geochemical Reference Material BRP-1 (Paraná Basin Basalt) and Assessment of Minimum Mass

Reference materials (RM) are required for quantitative analyses and their successful use is associated with the degree of homogeneity, and the traceability and confidence limits of the values established by characterisation. During the production of a RM, the chemical characterisation can only commence after it has been demonstrated that the material has the required level of homogeneity. Here we describe the preparation of BRP-1, a proposed geochemical reference material, and the results of the tests to evaluate its degree of homogeneity between and within bottles. BRP-1 is the first of two geochemical RM being produced by Brazilian institutions in collaboration with the United States Geological Survey (USGS) and the International Association of Geoanalysts (IAG). Two test portions of twenty bottles of BRP-1 were analysed by wavelength dispersive-XRF spectrometry and major, minor and eighteen trace elements were determined. The results show that for most of the investigated elements, the units of BRP-1 were homogeneous at conditions approximately three times more rigorous than those strived for by the test of "sufficient homogeneity". Furthermore, the within bottle homogeneity of BRP-1 was evaluated using small beam (1 mm²) synchrotron radiation XRF spectrometry and, for comparison, the USGS reference materials BCR-2 and GSP-2 were also evaluated. From our data, it has been possible to assign representative minimum masses for some major constituents (1 mg) and for some trace elements (1-13 mg), except Zr in GSP-2, for which test portions of 74 mg are recommended.     

Analysis of the precipitation and cloudiness associated with COLs occurrence in the Iberian Peninsula

Summary The Iberian Peninsula is one of the regions in the world with higher occurrence of cut-off low systems (COL). The aim of this paper is to analyse the weather events (rainfall and cloudiness layer) associated to COLs in the Iberian Peninsula with tools not previously used: (a) the use of the new multidecadal COLs database developed by Nieto et al (2005) that permit us to study a 41 years period (1958–1998), (b) the checking of the expected weather events (rainfall and cloudiness layer) associated with COLs in a conceptual model (Winkler et al, 2005) and (c) the extensive use of radiosoundings to analyse convective instability in areas inside and close to the COL. Two points of view are used to make the analysis: (1) a source oriented method, when a particular COL is followed and its associated precipitation and cloudiness is analysed over four quadrants in which Iberia was divided and (2) a receptor oriented method, when the precipitation associated to COLs is analysed in given areas, defined by patterns of precipitation. Results reveal that the precipitation and cloudiness patterns associated to COLs in the conceptual model reproduce quite well the main characteristics found over the Iberian Peninsula. The generalized idea that most of the COLs produce intense convective rainfall is show to be misleading. Convective phenomena are important usually when the centre of the COL is located on the Mediterranean region. Most of the rainfall associated with COLs comes from the baroclinic shield; specially in cases located over the west half of the Iberian Peninsula. It is shown that nearly 30% of COLs do not induce any rainfall; most of them located in the southern half of the Peninsula, and mainly during autumn. Only 30% of COLs produce generalized rainfall over the whole analysed territory, being most of them (about 90%) located over the western half of the Iberian Peninsula.     

EJECTION–COLLISION ORBITS AND INVARIANT PUNCTURED TORI IN A RESTRICTED FOUR‐BODY PROBLEM

We study the motion of an infinitesimal mass point under the gravitational action of three mass points of masses μ, 1–2μ and μ moving under Newton's gravitational law in circular periodic orbits around their center of masses. The three point masses form at any time a collinear central configuration. The body of mass 1–2μ is located at the center of mass. The paper has two main goals. First, to prove the existence of four transversal ejection–collision orbits, and second to show the existence of an uncountable number of invariant punctured tori. Both results are for a given large value of the Jacobi constant and for an arbitrary value of the mass parameter 0<μ≤1/2. This revised version was published online in July 2006 with corrections to the Cover Date.