The interior structure of Mercury and its core sulfur content

Depth-dependent interior structure models of Mercury are calculated for several plausible chemical compositions of the core and of the mantle. For those models, we compute the associated libration amplitude, obliquity, tidal deformation, and tidal changes in the external potential. In particular we study the relation between the interior structure parameters for five different mantle mineralogies and two different temperature profiles together with two extreme crust density values. We investigate the influence of the core light element concentration, temperature, and melting law on core state and inner and outer core size. We show that a sulfur concentration above 10 wt% is unlikely if the temperature at the core-mantle boundary is above 1850 K and the silicate shell at least 240 km thick. The interior models can only have an inner core if the sulfur weight fraction is below 5 wt% for core-mantle boundary temperature in the 1850-2200 K range. Within our modeling hypotheses, we show that with the expected precision on the moment of inertia the core size can be estimated to a precision of about 50 km and the core sulfur concentration with an error of about 2 wt%. This uncertainty can only be reduced when more information on the mantle mineralogy of Mercury becomes available. However, we show that the uncertainty on the core size estimation can be greatly reduced, to about 25 km, if tidal surface displacements and tidal variations in the external potential are considered.     

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.     

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.     

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.     

Chondritic micrometeorites from the Transantarctic Mountains

Abstract– On the basis of morphological and petrographic characteristics, eight “giant” unmelted micrometeorites in the 300–1100 μm size range were selected from the Transantarctic Mountain micrometeorite collection, Victoria Land, Antarctica. Mineralogical and geochemical data obtained by means of scanning electron microscopy, electron probe microanalyses, and synchrotron X‐ray diffraction allow their classification as chondritic micrometeorites. The large size of the micrometeorites increases considerably the amount of mineralogical and geochemical information compared to micrometeorites in smaller size fractions, therefore allowing a better definition of their parent material. A large variety of material is observed: five micrometeorites are related to unequilibrated and equilibrated ordinary chondrite, one to CV chondrite, one to CM chondrite, and one to CI chondrite parent materials. Besides reporting the first occurrence of a CV‐like micrometeorite, our study shows that the abundance of chondritic material supports observations from recent studies on cosmic spherules that a large part of the micrometeorite flux in this size range is of asteroidal origin.     

Effects of riverbed incision on the hydrology of the Vietnamese Mekong Delta

The hydrogeomorphology of the Vietnamese Mekong Delta (VMD) has been significantly altered by natural and anthropogenic drivers. In this study, the spatiotemporal changes of the flow regime were examined by analysing the long-term daily, monthly, annual and extreme discharges and water levels from 1980 to 2018, supported by further investigation of the long-term annual sediment load (from the 1960s to 2015), river bathymetric data (in 1998, 2014 and 2017) and daily salinity concentration (from the 1990s to 2015) using various statistical methods and a coupled numerical model. Then, the effects of riverbed incision on the hydrology were investigated. The results show that the dry season discharge (i.e., in March–June) of the Tien River increased by up to 23% from the predam period (1980–1992) to the postdam period (1993–2018) but that the dry season water level at My Thuan decreased by up to −46%. The annual mean and monthly water levels in June at Tan Chau and in January and June–October at My Thuan in the Tien River decreased statistically, even though the respective discharges increased significantly. These decreased water levels instead of the increased discharges were attributed to the accelerated riverbed incision upstream from My Thuan, which increased by more than three times, from a mean rate of −0.16 m/year (−16.7 Mm³/year) in 1998–2014 to −0.5 m/year (−52.5 Mm³/year) in 2014–2017. This accelerated riverbed incision was likely caused by the reduction in the sediment load of the VMD (from 166.7 Mt/year in the predam period to 57.6 Mt/year in the postdam period) and increase in sand mining (from 3.9 Mm³ in 2012 to 13.43 Mm³ in 2018). Collectively, the decreased dry season water level in the Tien River is likely one of the main causes of the enhanced salinity intrusion.     

The dark matter environment of the Abell 901/902 supercluster: a weak lensing analysis of the HST STAGES survey

We present a high-resolution dark matter reconstruction of the   z = 0.165  Abell 901/902 supercluster from a weak lensing analysis of the Hubble Space Telescope STAGES survey. We detect the four main structures of the supercluster at high significance, resolving substructure within and between the clusters. We find that the distribution of dark matter is well traced by the cluster galaxies, with the brightest cluster galaxies marking out the strongest peaks in the dark matter distribution. We also find a significant extension of the dark matter distribution of Abell 901a in the direction of an infalling X-ray group Abell 901α. We present mass, mass-to-light and mass-to-stellar mass ratio measurements of the structures and substructures that we detect. We find no evidence for variation of the mass-to-light and mass-to-stellar mass ratio between the different clusters. We compare our space-based lensing analysis with an earlier ground-based lensing analysis of the supercluster to demonstrate the importance of space-based imaging for future weak lensing dark matter 'observations'.     

Spatially distributed modelling of soil erosion and sediment yield at regional scales in Spain

Initiated by the need to quantify erosion rates and the impacts of global changes on erosion, several attempts have been made to apply erosion models at regional scales. However, these models have often been directed towards on-site soil erosion estimates, emphasising sheet and rill erosion processes, and disregarding gully erosion, channel erosion and sediment transport. These models are therefore of limited use for the assessment of sediment yield, off-site impacts of erosion, and for the development of environmental management to control these impacts at regional scale. This study analyses and compares three spatially distributed models for the prediction of soil erosion and/or sediment yield at regional scales: the WATEM-SEDEM model that is based on the empirical Revised Universal Soil Loss Equation (RUSLE) in combination with a sediment transport equation, the physics-based Pan European Soil Erosion Risk Assessment model (PESERA), and a newly developed Spatially Distributed Scoring model (SPADS). The three models were applied to 61 Spanish drainage basins and model predictions were evaluated against data on measured reservoir sedimentation rates. Global data sets on land use, climate, elevation and soil characteristics were used as model input for WATEM-SEDEM and SPADS, whereas published soil erosion estimates of PESERA at 1 km² resolution were used directly. SPADS and WATEM-SEDEM provided best results after separate calibration for basins with a Sediment Delivery Ratio (SDR) higher than 5% and those with an SDR lower than 5%. In this way, SPADS explained 67% of variation in sediment yield, while WATEM-SEDEM explained 48% of the variation. PESERA represents a promising alternative to the use of empirical models at the regional scale as it can be applied to very diverse environments with little calibration. However, PESERA provides soil erosion rates and not sediment yield estimates. For most basins PESERA soil erosion rates vary between fifty and close to zero percent of total sediment yield. Two major factors may explain this discrepancy between modelled soil erosion rates and measured sediment yield. First, it may be that PESERA underestimates soil erosion under Mediterranean conditions, although PESERA soil erosion rates are of the same order of magnitude as erosion rates measured in erosion plot studies. Second, gully-, river channel erosion and sediment transport processes may be much more important than sheet- and rill erosion for regional scale sediment yield in these environments. These issues therefore require further attention in future model development. Although spatially lumped models provide better predictions of sediment yield at the basin scale, and while validation of the predicted spatial patterns of sources and sinks of sediment requires further research, spatially distributed models are expected to be of value to support management decisions regarding the assessment of on-site and off-site impacts of erosion at the regional scale.