Martian gravity field model and its time variations from MGS and Odyssey data

We present results of several years of research and data processing aimed at modelling the Mars gravity field and its longest wavelength time variations. The new solution includes tracking data from Mars Global Surveyor (MGS) from 1998 to 2006 (end of mission) and from Mars Odyssey from 2002 to the spring of 2008; this is the longest analyzed data set from these two orbiter missions as compared to previous works. The new model has been obtained by a team working in Europe, independently from the works of groups at NASA Jet Propulsion Laboratory (JPL) and Goddard Space Flight Center (GSFC), also with totally independent software. Observations consist in two and three-way Doppler measurements (also one way for MGS), and range tracking data collected by the Deep Space Network and have been processed in 4 day arcs, taking into account all disturbing forces of gravitational and non-gravitational origins; for each arc the state vector, drag and solar pressure model multiplying factors, and angular momentum dump parameters are adjusted. The static field (MGGM08A) is represented in spherical harmonics up to degree and order 95 and is very close to previously published models (in terms of spectral components and also over specific features); correlations with the global Mars topography are established and apparent depths of compensation by degree are derived. Lumped zonal harmonics of degree two and three are solved for every 10 days, exhibiting variations in line with previous results (including authors’ ones); the work also shows the difficulty of finding clean signatures (annual and semi-annual) for the zonal coefficient of second degree. The k₂ Love number is also derived from the ensemble of data, as well as from subsets of them; values between 0.110 and 0.130 are found, which are consistent with the existence of a Martian fluid core of significant radius.     

Time-dependent simulations of steady C-type shocks

Using a time-dependent multifluid, magnetohydrodynamic code, we calculated the structure of steady perpendicular and oblique C-type shocks in dusty plasmas. We included relevant processes to describe mass transfer between the different fluids, radiative cooling by emission lines and grain charging, and studied the effect of single- and multiple-sized grains on the shock structure. Our models are the first of oblique fast-mode molecular shocks in which such a rigorous treatment of the dust grain dynamics has been combined with a self-consistent calculation of the thermal and ionization structures including appropriate microphysics. At low densities, the grains do not play any significant rôle in the shock dynamics. At high densities, the ionization fraction is sufficiently low that dust grains are important charge and current carriers and, thus, determine the shock structure. We find that the magnetic field in the shock front has a significant rotation out of the initial upstream plane. This is most pronounced for single-sized grains and small angles of the shock normal with the magnetic field. Our results are similar to previous studies of steady C-type shocks showing that our method is efficient, rigorous and robust. Unlike the method employed in the previous most detailed treatment of dust in steady oblique fast-mode shocks, ours allow a reliable calculation even when chemical or other conditions deviate from local statistical equilibrium. We are also able to model transient phenomena.     

Assessment of Raman spectroscopy as a tool for the non-destructive identification of organic minerals and biomolecules for Mars studies

Several characteristic geological features found on the surface of Mars by planetary rovers suggest that a possible extinct biosphere could exist based on similar sources of energy as occurred on Earth. For this reason, analytical instrumental protocols for the detection of biomarkers in suitable geological matrices unequivocally have to be elaborated for future unmanned explorations including the forthcoming ESA ExoMars mission. As part of the Pasteur suite of analytical instrumentation on ExoMars, the Raman/LIBS instrument will seek elemental and molecular information about geological, biological and biogeological markers in the Martian record. A key series of experiments on terrestrial Mars analogues, of which this paper addresses a particularly important series of compounds, is required to obtain the Raman spectra of key molecules and crystals, which are characteristic for each biomarker. Here, we present Raman spectra of several examples of organic compounds which have been recorded non-destructively—higher n-alkanes, polycyclic aromatic hydrocarbons, carotenoids, salts of organic acids, pure crystalline terpenes as well as oxygen-containing organic compounds. In addition, the lower limit of β-carotene detection in sulphate matrices using Raman microspectroscopy was estimated.     

Do extrusion ages reflect magma generation processes at depth? An example from the Neogene Volcanic Province of SE Spain

The high-K calc-alkaline volcanic rocks along the Neogene Volcanic Province of SE Spain represent crustal anatectic melts mixed with mantle components during the opening of the Alborán Sea. Partially melted metapelitic enclaves, along with the geochemical signature, provide evidence of their crustal source. U–Pb SHRIMP geochronology on monazite and zircon from enclaves and their hosting lavas in the localities of El Hoyazo, Mazarrón and Mar Menor reveals variable delays between the melting at depth and the eruption of the volcanics. These data indicate that: (1) the most important event of anatexis in the Neogene spanned at least the 3 m.y. interval between 12 and 9 Ma; (2) there is no trend in age of crustal melting; and (3) the delay between magma generation and extrusion varies from more than 3 m.y. at El Hoyazo to ~0.5 m.y. and possibly 2.5 m.y. at Mar Menor, with no significant delay measurable at Mazarrón. The variable time delay between anatexis and lava extrusion indicates that radiometric ages of volcanics may provide misleading information on the timing of magma genesis occurring at depth. This highlights the pitfall of basing detailed geodynamic models on volcanic extrusion ages alone. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The Mammoth Peak sheeted complex,Tuolumne batholith,Sierra Nevada,California: a record of initial growth or late thermal contraction in a magma chamber?

The Mammoth Peak sheeted intrusive complex formed in the interior of a ~7–10 km deep magma chamber, specifically in the Half Dome granodiorite of the Tuolumne batholith, central Sierra Nevada, CA (USA). The sheets consist of fractionated melts with accumulated hornblende, biotite, magnetite, titanite, apatite, and zircon. The accumulation, especially of titanite, had a profound effect on minor and trace elements (Nb, Ta, Ti, REE, U, Th, P, Zr, Hf, etc.), increasing their contents up to five to six times. Our thermal–mechanical modeling using the finite element method shows that cooling-generated tensile stresses resulted in the inward propagation of two perpendicular sets of dilational cracks in the host granodiorite. We interpret the sheeted complex to have formed by a crack-seal mechanism in a high strength, crystal-rich mush, whereby outward younging pulses of fractionated magma were injected into these syn-magmatic cracks at the margin of an active magma chamber. Thermal–mechanical instabilities developed after the assembly of the sheeted complex, which was then overprinted by late ~NW–SE magmatic foliation. This case example provides a cautionary note regarding the interpretation that sheeted zones in large granitoid plutons imply a diking mechanism of growth because the sheeted/dike complexes in plutons (1) may display inverse growth directions from the growth of the overall intrusive sequence; (2) need not record initial chamber construction and instead may reflect late pulsing of magma within an already constructed magma chamber; (3) have an overprinting magmatic fabric indicating the continued presence of melt after construction of sheeted complexes and thus a prolonged thermal history as compared to dikes; and (4) because the scale of the observed sheeted complexes may be small (<1%) in comparison to large homogenous parts of plutons, in which there is no evidence for sheeting or diking. Thus, where extensive dike complexes in plutons are absent, such as in much of the Tuolumne batholith, the application of an incremental diking model to explain chamber construction is at best speculative.     

Fluidization of buried mass-wasting deposits in lake sediments and its relevance for paleoseismology: Results from a reflection seismic study of lakes Villarrica and Calafquén (South-Central Chile)

A dense grid of very-high resolution seismic profiles on Lake Villarrica provides a quasi-3D view on intercalated lenses of low-amplitude reflections, which are connected by acoustic wipe-out patches and fractures to an underlying voluminous mass-wasting deposit. The lenses are interpreted as being created by earthquake-triggered liquefaction in this buried mass-wasting deposit and subsequent sediment fluidization and extrusion at the paleo-lake bottom. These sediment volcanoes are mapped in detail. They have a rather uniform circular geometry and show a linear relationship between apparent width and maximum thickness on a seismic section. The largest sediment volcanoes are up to 80 m wide and 1.9 m thick. Their slope angles designate a syn- to post-depositional sagging of most sediment volcanoes. Sediment volcano detection and mapping from nearby Lake Calafquén further strengthen the revealed geometrical relationships. Locally, some of the sediment/fluid escape structures extend to a higher position in the stratigraphy, which points to a polyphase escape process associated with multiple multi-century spaced strong earthquakes. Thickness and morphology of the source layer seem to exert a dominant control in the production of sediment/fluid extrusions. This study shows that reflection seismic profiling allowed recognizing 4 different seismic events in the studied stratigraphic interval, which are evidenced by mass-wasting deposits and/or fluidization features.     

Regional stratigraphy and distribution of epigenetic stratabound celestine, fluorite, barite and Pb–Zn deposits in the MVT province of northeastern Mexico

Northeastern Mexico hosts numerous epigenetic stratabound carbonate-hosted low-temperature hydrothermal deposits of celestine, fluorite, barite and zinc-lead, which formed by replacement of Mesozoic evaporites or carbonate rocks. Such deposits can be permissively catalogued as Mississippi Valley-type (MVT) deposits. The deposits studied in the state of Coahuila are associated with granitic and metasedimentary basement highs (horsts) marginal or central to the Mesozoic Sabinas Basin. These horsts controlled the stratigraphy of the Mesozoic basins and subsequently influenced the Laramide structural pattern. The Sabinas Basin consists of ~6,000-m-thick Jurassic to Cretaceous siliciclastic, carbonate and evaporitic series. The MVT deposits are mostly in Barremian and in Aptian-Albian to Cenomanian formations and likely formed from basinal brines that were mobilized during the Laramide orogeny, although earlier diagenetic replacement of evaporite layers (barite and celestine deposits) and lining of paleokarstic cavities in reef carbonates (Zn–Pb deposits) is observed. Fluid inclusion microthermometry and isotopic studies suggest ore formation due to mixing of basinal brines and meteoric water. Homogenization temperatures of fluid inclusions range from 45°C to 210°C; salinities range from 0 to 26 wt.% NaCl equiv., and some inclusions contain hydrocarbons or bitumen. Sulfur isotope data suggest that most of the sulfur in barite and celestine is derived from Barremian to Cenomanian evaporites. Regional geology and a compilation of metallogenic features define the new MVT province of northeastern Mexico, which comprises most of the state of Coahuila and portions of the neighboring states of Nuevo León, Durango and, perhaps extends into Zacatecas and southern Texas. This province exhibits a regional metal zonation, with celestine deposits to the south, fluorite deposits to the north and barite and Zn–Pb deposits mostly in the central part.     

A feather hydrogen isoscape for Mexico

Developing useful biological isoscapes for areas of the world is a priority. This is the case for Mexico that hosts a large percentage of North America's Neotropical migrant birds. Here we investigated the use of House Sparrow (Passer domesticus) feathers to create a spatially explicit feather deuterium isoscape for that country using samples (n = 461) that were collected across Mexico. Considerable and useful spatial hydrogen isotopic structure was observed, suggesting that isotopes may be a potential forensic tool for evaluating origins of Mexican derived fauna and flora. The most positive feather δD values occurred in the northeast and most negative in the south-central part of the country, roughly matching δD patterns observed in groundwater. A weak negative isotopic relationship was found with altitude in both the Pacific and Atlantic drainage systems. The most parsimonious model describing isotopic spatial variation in feathers between 300 and 3000 m a.s.l. included groundwater δD (δD_gw; precipitation proxy), sex, amount of precipitation, and the coefficient of variation in amount of precipitation. Overall, δD_gw was a poor predictor of sparrow δD_f values for all of Mexico. However, this relationship was considerably strengthened when we considered sex separately, removed the Baja peninsula from our sample, and considered the Atlantic and Pacific drainage basins separately. The strongest relationship between δD_gw and δD_f was found for female sparrows in the Atlantic drainage basin (r² = 0.464). We recommend that researchers interested in inferring origins of migratory birds and other animals in Mexico create species specific isotopic basemaps that may be guided by the isotopic patterns we have observed for House Sparrows and groundwater.     

A groundwater isoscape (δD, δO) for Mexico

Numerous studies have shown that precipitation isocapes drive δD and δ¹⁸O patterns in surficial waters and in terrestrial food webs. While the GNIP (Global Network for Isotopes in Precipitation) dataset provided a key foundation for linking precipitation-terrestrial isoscapes globally, it has insufficient spatial coverage in many countries like Mexico. To overcome this limitation, we hypothesized that shallow phreatic groundwaters in Mexico could be used as an isotopic integrator of long-term seasonally weighted precipitation inputs to the landscape to aid in calibrating spatial H and O isotope datasets for terrestrial, biological and hydrological research. Groundwater was sampled from 234 sites in Mexico at ~ 50 km latitudinal spacing to obtain high spatial resolution and country-wide coverage for the construction of a groundwater isoscape. Our data revealed that shallow groundwater infiltration in Mexico appears largely unaffected by evaporation and reflects seasonally weighted precipitation inputs. These precipitation inputs are primarily biased to summertime when highest rainfall occurs, but a small degree of post-precipitation evaporation revealed a lower d-excess zone that corresponded to the interior semi-arid ecozone. We developed a predictive general linear model (GLM) for hydrogen and oxygen isotopic spatial patterns in Mexican groundwater and then compared the results to a validation subset of our field data, as well external data reported in the literature. The GLM used elevation, latitude, drainage basin (Atlantic vs. Pacific), and rainfall as the most relevant predictive variables. The GLM explained 81% of the overall isotopic variance observed in groundwater, 68% of the variance within our validation subset, and 77% of the variance in the external data set. Our predictive GLM is sufficiently accurate to allow for future ecological, hydrological and forensic isoscape applications in Mexico, and may be an approach that is applicable to other countries and regions where GNIP stations are lacking.     

A site evaluation campaign for a ground based atmospheric Cherenkov telescope in Romania

Around the world, several scientific projects share the interest of a global network of small Cherenkov telescopes for monitoring observations of the brightest blazars??the DWARF network. A small, ground based, imaging atmospheric Cherenkov telescope of last generation is intended to be installed and operated in Romania as a component of the DWARF network. To prepare the construction of the observatory, two support projects have been initiated. Within the framework of these projects, we have assessed a number of possible sites where to settle the observatory. In this paper we submit a brief report on the general characteristics of the best four sites selected after the local infrastructure, the nearby facilities and the social impact criteria have been applied.