The graininess of dark matter haloes

We use the recently completed one billion particle Via Lactea II Λ cold dark matter simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as the structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions, the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo cannot be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.     

Multimass spherical structure models for N-body simulations

We present a simple and efficient method to set up spherical structure models for N -body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N -body simulations for a given computer run time. Here, we present several models with an effective resolution of up to  1.68 × 10⁹  particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen that in mergers of collisionless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10⁸ particles. We also find that in a core–core merger the central density approximately doubles whereas in the cusp–cusp case the central density only increases by approximately 50 per cent. This may suggest that the central regions of flat structures are better protected and get less energy input through the merger process.     

Can agricultural groundwater economies collapse? An inquiry into the pathways of four groundwater economies under threat

The aim of this paper is to investigate the notion of collapse of agricultural groundwater economies using the adaptive-cycle analytical framework. This framework was applied to four case studies in southern Europe and North Africa to question and discuss the dynamics of agricultural groundwater economies. In two case studies (Saiss in Morocco and Clain basin in France), the imminent physical or socio-economic collapse was a major concern for stakeholders and the early signs of collapse led to re-organization of the groundwater economy. In the other two cases (Biskra in Algeria and Almeria in Spain), collapse was either not yet a concern or had been temporarily resolved through increased efficiency and access to additional water resources. This comparative analysis shows the importance of taking the early signs of collapse into account. These signs can be either related to resource depletion or to environmental and socio-economic impacts. Beyond these four case studies, the large number of groundwater economies under threat in (semi-)arid areas should present a warning regarding their possible collapse. Collapse can have severe and irreversible consequences in some cases, but it can also mean new opportunities and changes.     

3D structural model and kinematic interpretation of the Panixer Pass Transverse Zone (Infrahelvetic Complex,eastern Switzerland)

The Panixer Pass Transverse Zone in the eastern Swiss Alps is oriented perpendicular to most alpine structures in the area. Its main element is the SSE-trending Crena-Martin Fold, a downward facing fold with Permian Verrucano in its core, which is cut by the Glarus Thrust. Hence Verrucano can be found below the Glarus Thrust in the Infrahelvetic Complex. Across the Panixer Pass Transverse Zone the structural buildup of the Infrahelvetic Complex changes considerably. Multiple published theories of the structural evolution are not satisfying particularly because traditional 2D geological cross-sections are insufficient to understand the 3D complexity. The main result and product of our study is a 3D structural model of the Panixer Pass Transverse Zone providing insight into its geometry. As modeling input, we produced a lithostratigraphic map and collected structural orientation data. The 3D structural model honors the observed surface geology and the expected 3D subsurface geometry. Our field data indicates that the shearing and transport direction was continuously NNW-directed, except for a phase of north-directed shearing during the early movement along the Glarus Thrust (late Calanda Phase) and related foliation development in the Helvetic Nappes. The Panixer Pass Transverse Zone developed prior to the penetrative foliation during a thrust-dominated deformation phase (Cavistrau Phase), for which we created a kinematic block model. According to this model, the Crena-Martin Fold is the result of multiple lateral ramps and related lateral fault-bend folds that all developed in a similar positon. In particular, we do not propose ENE-WSW-directed shortening to form the Crena-Martin Fold. The latter was finally cut at low angle by a dextral strike-slip fault to create the final geometry of the Panixer Pass Transverse Zone. Our kinematic model reproduces the main features of the 3D structural model and embeds well into previously proposed sequences of deformation phases.     

Seasonal climate information preserved in West Antarctic ice core water isotopes: relationships to temperature, large-scale circulation, and sea ice

As part of the United States’ contribution to the International Trans-Antarctic Scientific Expedition (ITASE), a network of precisely dated and highly resolved ice cores was retrieved from West Antarctica. The ITASE dataset provides a unique record of spatial and temporal variations of stable water isotopes (δ¹⁸O and δD) across West Antarctica. We demonstrate that, after accounting for water vapor diffusion, seasonal information can be successfully extracted from the ITASE cores. We use meteorological reanalysis, weather station, and sea ice data to assess the role of temperature, sea ice, and the state of the large-scale atmospheric circulation in controlling seasonal average water isotope variations in West Antarctica. The strongest relationships for all variables are found in the cores on and west of the West Antarctic Ice Sheet Divide and during austral fall. During this season positive isotope anomalies in the westernmost ITASE cores are strongly related to a positive pressure anomaly over West Antarctica, low sea ice concentrations in the Ross and Amundsen Seas, and above normal temperatures. Analyses suggest that this seasonally distinct climate signal is due to the pronounced meridional oriented circulation and its linkage to enhanced sea ice variations in the adjacent Southern Ocean during fall, both of which also influence local to regional temperatures.     

Comparison of climate field reconstruction techniques: application to Europe

This paper presents a comparison of principal component (PC) regression and regularized expectation maximization (RegEM) to reconstruct European summer and winter surface air temperature over the past millennium. Reconstruction is performed within a surrogate climate using the National Center for Atmospheric Research (NCAR) Climate System Model (CSM) 1.4 and the climate model ECHO-G 4, assuming different white and red noise scenarios to define the distortion of pseudoproxy series. We show how sensitivity tests lead to valuable “a priori” information that provides a basis for improving real world proxy reconstructions. Our results emphasize the need to carefully test and evaluate reconstruction techniques with respect to the temporal resolution and the spatial scale they are applied to. Furthermore, we demonstrate that uncertainties inherent to the predictand and predictor data have to be more rigorously taken into account. The comparison of the two statistical techniques, in the specific experimental setting presented here, indicates that more skilful results are achieved with RegEM as low frequency variability is better preserved. We further detect seasonal differences in reconstruction skill for the continental scale, as e.g. the target temperature average is more adequately reconstructed for summer than for winter. For the specific predictor network given in this paper, both techniques underestimate the target temperature variations to an increasing extent as more noise is added to the signal, albeit RegEM less than with PC regression. We conclude that climate field reconstruction techniques can be improved and need to be further optimized in future applications.     

Detailed debris flow hazard assessment in Andorra: A multidisciplinary approach

In many mountainous areas, the rapid development of urbanisation and the limited space in the valley floors have created a need to construct buildings in zones potentially exposed to debris flow hazard. In these zones, a detailed and coherent hazard assessment is necessary to provide an adequate urban planning. This article presents a multidisciplinary procedure to evaluate the debris flow hazard at a local scale. Our four-step approach was successfully applied to five torrent catchments in the Principality of Andorra, located in the Pyrenees. The first step consisted of a comprehensive geomorphologic and geologic analysis providing an inventory map of the past debris flows, a magnitude–frequency relationship, and a geomorphologic–geologic map. These data were necessary to determine the potential initiation zones and volumes of future debris flows for each catchment. A susceptibility map and different scenarios were the principal outcome of the first step, as well as essential input data for the second step, the runout analysis. A one-dimensional numerical code was applied to analyse the scenarios previously defined. First, the critical channel sections in the fan area were evaluated, then the maximum runout of the debris flows on the fan was studied, and finally simplified intensity maps for each defined scenario were established. The third step of our hazard assessment was the hazard zonation and the compilation of all the results from the two previous steps in a final hazard map. The base of this hazard map was the hazard matrix, which combined the intensity of the debris flow with its probability of occurrence and determined a certain hazard degree. The fourth step referred to the hazard mitigation and included some recommendations for hazard reduction. In Andorra, this four-step approach is actually being applied to assess the debris flow hazard. The final hazard maps, at 1 : 2000 scale, provide an obligatory tool for local land use planning. Experience achieved during the study showed that the collaboration between geologists, geomorphologists, engineers, and decision makers is essential and that only a multidisciplinary approach allows for solving all the problems of such a complex process as debris flows. Finally, we propose that our approach may be applied to other mountainous areas, adapting the hazard matrix to new local conditions.     

Partial melting processes above subducting plates: Constraints from Pa-U disequilibria

The processes involved in the formation and transport of partial melts above subducting plates remain poorly constrained relative to those at mid-ocean ridges. In particular, ²³⁸U-²³⁰Th-²²⁶Ra disequilibria, that might normally be used to constrain melting dynamics, tend to be swamped by the effects of fluid addition from the down-going plate. The ²³¹Pa-²³⁵U system provides an exciting exception to this because the highly incompatible nature of Pa means that fractionation and in-growth during partial melting overwrite the effects of fluid U addition. We present ²³¹Pa-²³⁵U data on 50 well-characterised lavas from seven subduction zones in order to examine partial melting processes. Measured (²³¹Pa/²³⁵U) ratios are all >1 and 15% are >2. Overall (²³¹Pa/²³⁵U) shows broad positive correlations with (²³⁰Th/²³⁸U) and La/Yb and negative trends against Ba/Th and (²²⁶Ra/²³⁰Th). These systematics can differ from arc to arc but suggest that (²³¹Pa/²³⁵U) tends to be higher in sediment-rich arc lavas where the effects of fluid addition are muted and there is less of a ²³¹Pa deficit for melting to overprint. We have explored the effects of decompression melting, frictional drag dynamic melting with and without ageing subsequent to fluid U addition to the wedge as well as flux melting models. Globally, average (²³¹Pa/²³⁵U) appears to correlate negatively with convergence rate and so in the numerical models we use the local subduction rate for the rate of matrix flow through the melting zone. Using this assumption and reasonable values for other parameters, the melting models can simulate the overall range of (²³¹Pa/²³⁵U) and some of the data trends. However, it is clear that local variations in some parameters, especially source composition and extent of melting, exert a major influence on ²³¹Pa-²³⁵U disequilibria. Some data, which lie at a high angle to the modelled trends, may be explained by mixing between small degree hydrous melts formed near the slab and larger degree, decompression melts produced at shallow depth.     

Genesis of the Dordogne (France) high-purity silica placers

This article describes the formation of the principal French high-purity silica placer (source for electrometallurgy of silicon and ferrosilicon), the Boudeau deposit in Dordogne, and of one of its subordinate deposits, the Pays Brûlé deposit. Boudeau has always been considered to be a colluvial deposit, formed in a trap produced by extensional tectonism. The present work suggests a new genetic interpretation, based on a combined study of the deposits and their underlying rocks. It shows that the quartz gravel deposits, laid down by stream flow within braided channels, were reworked by karstic progressive sinking, and that the mechanism responsible for the increase in the thickness of the economic formations is cryptokarstic corrosion. This undercovers weathering process, which involves dissolution of karstifiable rocks (oolitic and bioclastic limestones) at their interface with a permeable, non-karstifiable cover composed of chemically pure quartz gravels in a sandy-argillaceous matrix. A typological study of the gravels provides a basis for tackling the problem of their origin and their purity. They seem to be derived mainly from the synmetamorphic quartz mobilisates of the stripped metamorphic units (the lower and upper gneiss units forming the western edge of the Massif Central). Lastly, a dynamic model is suggested for the genesis of the high-purity silica deposits of Dordogne.     

Ceilometer observations of aerosol layer structure above the Petit Lubéron during ESCOMPTE's IOP 2

A modified ceilometer has been used during the second Intensive Observation Period (IOP) of the “Expérience sur Site pour COntraindre les Modèles de Pollution atmosphériques et de Transport d'Émission” (ESCOMPTE) to perform continuous remote observations of aerosol accumulations in the first 3 km of the atmosphere. These observations encompassed an episode of intense particulate and photochemical pollution. The submicronic particles density, measured at an altitude of 600 m, went from a very low point of a few tens of particles per cubic centimeter (at the end of a Mistral episode in the free atmosphere) to a high point of more than 4500 particles per cubic centimeter (when pollutants were trapped by thermal inversions).The main result is that this instrument enables a fine documentation of the mixing layer height and of aerosol particles stratifications and circulation. Airborne aerosol measurements have been made above the mountainous region of Mérindol in order to validate in situ the remote sensing measurements. Ozone measurements near the summit of the mountains as well as in the valley were performed in order to correlate aerosol accumulation and ozone concentration. As a notable example, the two-layer aerosol stratification seen in the first 2 days of IOP 2b in that part of the ESCOMPTE domain confirms the results of another team which used backtrajectories. The low-altitude pollution for this timeframe had a local origin (the Fos industrial area), whereas above 500 m, the air masses had undergone regional-scale transport (from north-eastern Spain).The second major result is the highlighting of a pattern, in sea breeze conditions and in this part of the ESCOMPTE experiment zone, of nocturnal aerosol accumulation at an altitude of between 500 and 2000 m, followed by high ozone concentration the next day.