A data assimilation approach to discharge estimation from space

River discharge is currently monitored by a diminishing network of gauges, which provide a spatially incomplete picture of global discharges. This study assimilated water level information derived from a fused satellite Synthetic Aperture Radar (SAR) image and digital terrain model (DTM) with simulations from a coupled hydrological and hydrodynamic model to estimate discharge in an un‐gauged basin scenario. Assimilating water level measurements led to a 79% reduction in ensemble discharge uncertainty over the coupled hydrological hydrodynamic model alone. Measurement bias was evident, but the method still provided a means of improving estimates of discharge for high flows. The study demonstrates the potential of currently available synthetic aperture radar imagery to reduce discharge uncertainty in un‐gauged basins when combined with model simulations in a data assimilation framework, where sufficient topographic data are available. The work is timely because in the near future the launch of satellite radar missions will lead to a significant increase in the volume of data available for space‐borne discharge estimation. Copyright © 2009 John Wiley & Sons, Ltd.     

Influence of cluster environments on X-ray cluster morphology

We investigated the influence of environment on cluster morphology by examining the alignment of the cluster position angle with respect to the nearest neighbour cluster. The cluster position angle and ellipticity were measured using high spatial resolution X-ray data taken from the Chandra archive, while the nearest neighbour sample was extracted from the Abell cluster catalogue. We found high confidence for an alignment when neighbour distance (Dn) was less than  70  h ⁻¹  Mpc, while we found no significant variations of the confidence level at various values of Dn less than  70  h ⁻¹  Mpc.     

Orbital and thermal evolutions of four potential targets for a sample return space mission to a primitive near-Earth asteroid

In this paper, we present our study of the orbital and thermal evolutions, due to solar radiative heating, of four near-Earth asteroids (NEAs) considered as potential target candidates for sample return space missions to primitive asteroids. We used a dynamical model of the NEA population to estimate the most likely source region and orbital history of these objects. Then, for each asteroid, we integrated numerically over their entire lifetime a set of 14 initially indistinguishable orbit (clones), obtained by small variations of the nominal initial conditions. Using a thermal model, we then computed surface and sub-surface temperatures of these bodies during their dynamical history. Our aim is to determine whether these bodies are likely to have experienced high temperature level, and whether great temperature changes can be expected due to the orbital changes as well as their maximum and minimum values. Such information is important in the framework of sample return space missions whose goal is to bring back pristine materials. The knowledge of the temperature range of materials at different depth over the orbital evolution of potential targets can help defining sampling strategies that ensure the likelihood that unaltered material will be brought back. Our results suggest that for all the considered potential targets, the surface has experienced for some time temperatures greater than 400 K and at most 500 K with 50% probability. This probability drops rapidly with increasing temperature. Sub-surface materials at a depth of only 3 cm are much more protected from high temperature and generally do not reach temperatures exceeding 450 K (with 50% probability). They should thus be unaltered at this depth at least from a Sun-driven heating point of view. On the other hand, surface material for some of the considered objects can have a range of temperature which can make them less reliable as pristine materials. However, it is assumed here that the same material is constantly exposed to solar heat, while regolith turnover may occur. The latter can be caused by different processes such as seismic shaking and/or impact cratering. This would reduce the total time that materials are exposed to a certain temperature. Thus, it is very likely that a sample collected from any of the four considered targets, or any primitive NEA with similar dynamical properties, will have components that will be thermally unaltered as long as some of it comes from only 3 to 5 cm depth. Such a depth is not considered difficult to reach with some of the current designs of sampling devices.     

Models of high velocity impacts into dust-covered ice: Application to Martian northern lowlands

The detection of fresh impact craters with bright floors and ejecta (arising from fresh clean water ice) in the northern lowlands of Mars (Byrne et al., 2009b, Science 325, 1674), together with observations of polygonal structures and evidence from the Phoenix probe, suggests that there are substantial water ice deposits just below the surface over large areas. Specifically in cases of the larger craters observed, the impacts themselves may have raised the temperature and the pressure of the water ice deposits locally to values which allow phase changes. In this paper, we use smoothed particle hydrodynamics to model hyper-velocity impacts. We estimate peak shock pressures in a solid water ice target covered by a layer of loose material, modeled by pre-damaged dunite. In addition, we account for the possibility of a thin layer of sub-surface water ice by using a three-layer model where the ice is surrounded by dunite. We find that the peak shock pressures reached in the simulated events are high enough to produce several 100-1000 kg of liquid water depending upon the impact parameters and the exact shock pressure needed for the phase change. A difficulty remains however in determining whether liquid is generated or whether a type of fluidized ice is produced (or indeed some combination of the two). We also note that the process can become rather complex as the number of layers increases because of reflections of the shock at sub-surface boundaries—a process which should lead to increased fluidization.     

Alpine tectonics in the Calabrian–Peloritan belt (southern Italy): New 40Ar/39Ar data in the Aspromonte Massif area

This study provides new ⁴⁰Ar/³⁹Ar geochronological constraints on the age of the Alpine tectonics in the Aspromonte Massif (southern part of the Calabrian–Peloritan belt). This massif exposes the upper units of the Calabride Complex which originated from the European continental margin. The Calabride Complex was incorporated in the Alpine orogenic wedge and then integrated into the Apennines and Maghrebides fold-and-thrust belts. Throughout the Calabride Complex there is evidence for a two stage tectonic history, which remains however rather poorly dated: Alpine nappe stacking is followed by extensional reworking along the former thrust contacts or along new detachment surfaces. Our new ages suggest that exhumation of the uppermost units, which accompanied nappe stacking, probably started at 45 Ma and that the deepest units were almost completely exhumed at 33 Ma. This kinematics probably corresponds to syn-orogenic extension while the end of exhumation is clearly related to the extensional tectonics dated at 28.6 Ma along detachment structures.Our geochronological data reveal a very short lag time between accretional and extensional processes in this part of the Mediterranean Alpine orogenic belt. The direction of extension, when the units are restored to their initial position (i.e. before the opening of the Western Mediterranean basins and the bending of the arc) is NNE–SSW. Such a direction does not fit with the eastward slab-retreat model generally put forward to explain extension in the Western Mediterranean. In contrast, we provide evidence for roughly N–S middle Oligocene extension in the accretionary prism, not previously described in this part of the Mediterranean domain.     

Multi-stage reaction history in different eclogite types from the Pakistan Himalaya and implications for exhumation processes

Metabasites were sampled from rock series of the subducted margin of the Indian Plate, the so-called Higher Himalayan Crystalline, in the Upper Kaghan Valley, Pakistan. These vary from corona dolerites, cropping out around Saif-ul-Muluk in the south, to coesite–eclogite close to the suture zone against rocks of the Kohistan arc in the north. Bulk rock major- and trace- element chemistry reveals essentially a single protolith as the source for five different eclogite types, which differ in fabric, modal mineralogy as well as in mineral chemistry. The study of newly-collected samples reveals coesite (confirmed by in situ Raman spectroscopy) in both garnet and omphacite. All eclogites show growth of amphiboles during exhumation. Within some coesite-bearing eclogites the presence of glaucophane cores to barroisite is noted whereas in most samples porphyroblastic sodic–calcic amphiboles are rimmed by more aluminous calcic amphibole (pargasite, tschermakite, and edenite). Eclogite facies rutile is replaced by ilmenite which itself is commonly surrounded by titanite. In addition, some eclogite bodies show leucocratic segregations containing phengite, quartz, zoisite and/or kyanite. The important implication is that the complex exhumation path shows stages of initial cooling during decompression (formation of glaucophane) followed by reheating: a very similar situation to that reported for the coesite-bearing eclogite series of the Tso Morari massif, India, 450 km to the south-east.     

Development of a Darcy-Brinkman model to simulate water flow and tracer transport in a heterogeneous karstic aquifer (Val d’Orléans, France)

Darcy’s law is the equation of reference widely used to model aquifer flows. However, its use to model karstic aquifers functioning with large pores is problematic. The physics occurring within the karstic conduits requires the use of a more representative macroscopic equation. A hydrodynamic model is presented which is adapted to the karstic aquifer of the Val d’Orléans (France) using two flow equations: (1) Darcy’s law, used to describe water flow within the massive limestone, and (2) the Brinkman equation, used to model water flow within the conduits. The flow equations coupled with the transport equation allow the prediction of the karst transfer properties. The model was tested by using six dye tracer tests and compared to a model that uses Darcy’s law to describe the flow in karstic conduits. The simulations show that the conduit permeability ranges from 5?×?10^?6 to 5.5?×?10^?5?m² and the limestone permeability ranges from 8?×?10^?11 to 6?×?10^?10?m². The dispersivity coefficient ranges from 23 to 53 m in the conduits and from 1 to 5 m in the limestone. The results of the simulations carried out using Darcy’s law in the conduits show that the dispersion towards the fractures is underestimated.     

A field study (Massachusetts, USA) of the factors controlling the depth of groundwater flow systems in crystalline fractured-rock terrain

Groundwater movement and availability in crystalline and metamorphosed rocks is dominated by the secondary porosity generated through fracturing. The distributions of fractures and fracture zones determine permeable pathways and the productivity of these rocks. Controls on how these distributions vary with depth in the shallow subsurface (<300 m) and their resulting influence on groundwater flow is not well understood. The results of a subsurface study in the Nashoba and Avalon terranes of eastern Massachusetts (USA), which is a region experiencing expanded use of the fractured bedrock as a potable-supply aquifer, are presented. The study logged the distribution of fractures in 17 boreholes, identified flowing fractures, and hydraulically characterized the rock mass intersecting the boreholes. Of all fractures encountered, 2.5% are hydraulically active. Boreholes show decreasing fracture frequency up to 300 m depth, with hydraulically active fractures showing a similar trend; this restricts topographically driven flow. Borehole temperature profiles corroborate this, with minimal hydrologically altered flow observed in the profiles below 100 m. Results from this study suggest that active flow systems in these geologic settings are shallow and that fracture permeability outside of the influence of large-scale structures will follow a decreasing trend with depth.     

Application of drone-captured thermal imagery in aiding in the recovery of meteorites within a snow-covered strewn field

The rapid recovery of meteorites mitigates the exposure of astromaterials to the terrestrial environment and subsequent contamination. Modern fireball observatories have enabled the more accurate triangulation of fireball trajectories, which has aided in the location of strewn fields, in the case of meteorite-producing events. Despite this advancement, most meteorite searches still use manual searching to locate any meteorite falls, which is often labor-intensive and has a slow coverage rate (km² day⁻¹). Recent work has begun exploring the application of drone technology to the recovery of meteorites; however, most of this work has focused on falls in arid environments. Our study examines the utilization of drones with thermal imaging technology to aid in the recovery of meteorites that have fallen on a snow-covered field. We created a simulated strewn field that included meteorite specimens as well as Earth rocks with similar properties (“meteowrongs”). Thermal imagery was utilized to determine whether the thermal contrast between meteorites and snow could aid in the identification of meteorites. We found that the thermal contrast was significant enough that meteorites were readily identifiable within thermal images; however, it was not significant enough to distinguish between the meteorites and the meteowrongs. The utilization of thermal imagery in conjunction with visible imagery has the potential to aid in the rapid recovery of meteorites in snow-covered landscapes.