The DRASTIC-Sg model: an extension to the DRASTIC approach for mapping groundwater vulnerability in aquifers subject to differential land subsidence,with application to Mexico City

Mexico City relies on groundwater for most of its domestic supply. Over the years, intensive pumping has led to significant drawdowns and, subsequently, to severe land subsidence. Tensile cracks have also developed or reactivated as a result. All such processes cause damage to urban infrastructure, increasing the risk of spills and favoring contaminant propagation into the aquifer. The effects of ground deformation are frequently ignored in groundwater vulnerability studies, but can be relevant in subsiding cities. This report presents an extension to the DRASTIC methodology, named DRASTIC-Sg, which focuses on evaluating groundwater vulnerability in urban aquifers affected by differential subsidence. A subsidence parameter is developed to represent the ground deformation gradient (Sg), and then used to depict areas where damage risk to urban infrastructure is higher due to fracture propagation. Space-geodetic SqueeSAR data and global positioning system (GPS) validation were used to evaluate subsidence rates and gradients, integrating hydrogeological and geomechanical variables into a GIS environment. Results show that classic DRASTIC approaches may underestimate groundwater vulnerability in settings such as the one at hand. Hence, it is concluded that the Sg parameter is a welcome contribution to develop reliable vulnerability assessments in subsiding basins.     

Hard X-ray emission from the SNR G337.2+0.1

We report hard X-ray emission of the non-thermal supernova remnant G337.2+0.1. The source presents centrally filled and diffuse X-ray emission. A spectral study confirms that the column density of the central part of the object is about N _H∼5.9(±1.5)×10²² cm⁻² and its X-ray spectrum is well represented by a single power-law with a photon index Γ=0.96±0.56. Detailed spectral analysis indicates that the outer region is highly absorbed and quite softer than the inner region. Characteristics already observed in other well-known X-ray plerions. Based on the gathered information, we confirm the SNR nature of G337.2+0.1, and suggest that the central region of the source is a pulsar wind nebula (PWN), originated by an energetic though yet undetected pulsar.     

Optimal low-thrust trajectories to asteroids through an algorithm based on differential dynamic programming

In this paper an optimisation algorithm based on Differential dynamic programming is applied to the design of rendezvous and fly-by trajectories to near Earth objects. Differential dynamic programming is a successive approximation technique that computes a feedback control law in correspondence of a fixed number of decision times. In this way the high dimensional problem characteristic of low-thrust optimisation is reduced into a series of small dimensional problems. The proposed method exploits the stage-wise approach to incorporate an adaptive refinement of the discretisation mesh within the optimisation process. A particular interpolation technique was used to preserve the feedback nature of the control law, thus improving robustness against some approximation errors introduced during the adaptation process. The algorithm implements global variations of the control law, which ensure a further increase in robustness. The results presented show how the proposed approach is capable of fully exploiting the multi-body dynamics of the problem; in fact, in one of the study cases, a fly-by of the Earth is scheduled, which was not included in the first guess solution.     

The growth and contamination mechanism of the Cana Brava layered mafic-ultramafic complex: new field and geochemical evidences

Mineralogy and Petrology - The Cana Brava complex is the northernmost of three layered complexes outcropping in the Goiás state (central Brasil). New field and geochemical evidences suggest...     

From Alluvial Fan to Playa: An Upper Jurassic Ephemeral Fluvial System, Neuquén Basin, Argentina

The Upper Jurassic fill of the back-arc Neuquén Basin includes a lowstand wedge known as the Tordillo Formation. The studied deposits crop out along a N-S oriented belt that runs parallel to the Andean magmatic arc. They are limited to the south by the east-west oriented positive structure of the Huincul arch formed as a result of Upper Jurassic tectonic inversion. The Tordillo deposits were formed in an arid fluvial-dominated system characterised by systematic downstream changes in architectural style. A gravelly and sandy bedload fluvial system is recognised in the southern upstream sector. The reduced thickness and the coarse grain size suggest steep gradients, excess of bedload supply and a low subsidence rate. Thicker and finer-grained deposits prevail to the north and northeast. They were formed under arid conditions in a wadi-sand flat-playa fluvial system. This distal facies association indicates increased accommodation owing to high rates of subsidence relative to coarse siliciclastic sedimentation rates. These low-gradient deposits are characterised by cyclic alternations of mud-dominated and sand-dominated packages interpreted as high- and low-accommodation systems tracts. The overall fining upward stacking pattern of the Tordillo Formation suggests a change towards higher accommodation rates. This is accompanied by frequent development of soil horizons and darker primary and reworked pyroclastic deposits. These attributes indicate a stronger explosive volcanic activity associated with increased precipitation and high water table emplacement towards the end of the Tordillo lowstand wedge.     

First lead isotopic data for cinnabar in the Almadén district (Spain): implications for the genesis of the mercury deposits

The Almadén district constitutes the largest and probably the most intriguing mercury concentration in the world. Two types of mineralization are recognized: 1) stratabound, of Lower Silurian age, well represented by the large Almadén deposit; and 2) fully discordant mercury deposits of minor importance in terms of size, and exemplified by the deposit of Las Cuevas. The latter ones can be found at different positions along the Almadén stratigraphic column. Both types of deposits are always associated with the so-called frailesca rocks (diatremes of alkali basaltic composition). This paper reports the first lead isotope compositions of cinnabar in the district. Whole samples and stepwise leaching cinnabar aliquots display relatively homogeneous isotopic compositions (²⁰⁶Pb/²⁰⁴Pb = 18.112–18.460; ²⁰⁷Pb/²⁰⁴Pb = 15.635–15.705; ²⁰⁸Pb/²⁰⁴Pb = 38.531–38.826). Taken together with Jébrak et al.s (2002) pyrite lead isotope results, the new cinnabar isotopic data define a steep array trend on the ²⁰⁷Pb/²⁰⁴Pb– ²⁰⁶Pb/²⁰⁴Pb diagram, indicating a mixed contribution of lead and probably mercury from different sedimentary sources in the Almadén basin. The Almadén Hg deposits are related to a contemporaneous mafic magmatism that might have provided part of the mercury. Hydrothermal leaching of organic matter from sedimentary rocks and formation of Hg organic complexes enhanced metal solubility, promoting transport from and within the volcanic units.Editorial handling: M. Chiaradia     

Theory of motion of saturn's coorbiting satellites

A simple analytic theory describing the 1:1 orbital resonance is presented and applied to Saturn's coorbiting pair, 1980S1 and 1980S3. These satellites are very small and can approach to within 15,000 km, but are prevented from passing each other by their mutual gravitational interaction. The long-term stability of the S1–S3 orbital configuration is discussed in this paper, and a tie between the 1966 and 1980 observations is establised.     

Numerical calculation of covariances between area means

When area means are given on a regular grid, their covariance functions, needed in least squares collocation to estimate quantities statistically related to the data, can be calculated as accurately as desired and very efficiently in the form of a Fourier series truncated at a suitably high frequency.     

Gravity field determination from satellite gradiometry

An orbiting gradiometer measures simultaneously several gravity quantities, ideally all six second-order derivatives of the gravitational potential. These contain information on the orbit, on the structure of the gravity field, and on the attitude of the space-craft. Due to the availability of several components simultaneously it is possible to separate orbit determination from attitude or gravity field recovery. This facilitates the analysis of the gradiometer measurements and allows the use of the principles of fast spherical harmonic analysis. The separation of gravity field recovery and orbit determination is tested numerically with a simplified gravity field (with a purely zonal spherical harmonic expansion) up to degree 300. For both the potential coefficients and for the orbit an almost exact recovery is attained after two iteration steps.