The precession phase of hydrological variability in the Western Pacific Warm Pool during the past 400 ka

The low-latitude hydrological cycle is a key climate parameter on different timescales, as it contributes to various feedback processes. Modelling studies suggest that the interhemispheric insolation contrast is the major factor controlling the cycle, although the influence of glacial conditions and the phase relationships relative to insolation forcing remain undetermined. In this work, we studied precipitation variability over Papua New Guinea (PNG, 3°S) for the past 400 ka using terrigenous fractions transported by the Sepik River to the Western Pacific Warm Pool (WPWP). A multi-decadal to centennial resolution of the elemental content was obtained using X-ray fluorescence scanning of a marine sediment core using an age model based on ¹⁴C dates and benthic foraminiferal δ¹⁸O. Indicators of the coarse river particulate fraction (bulk and CaCO₃-free basis Ti concentrations, the log intensity ratios of Ti/K and Ti/Ca) displayed a dominant 23 ka periodicity without a clear glacial–interglacial trend. Our precipitation records showed a tight relationship with local summer insolation (3°S, January) with time-dependent lag of 0 to 4 ka. They were generally in anti-phase for U–Th dated Chinese speleothem δ¹⁸O records. Based on an analogy to modern climate, we propose that precipitation over PNG was primarily determined by interhemispheric insolation contrast, and the contribution of austral fall/winter precipitation added second-order variability that formed the lags. For the last four climate cycles, the WPWP hydrological cycle was closely associated with the eastern Asian monsoon, and the influence of glacial conditions on the low-latitude hydrological cycle was estimated to be limited.     

From paleostresses to paleoburial in fold-thrust belts: Preliminary results from calcite twin analysis in the Outer Albanides

This paper presents a new approach to constrain paleoburial and subsequent uplift by folding in fold-thrust belts, combining differential stress estimates from mechanically-induced calcite twins with the assumption that stress in the upper crust is in frictional equilibrium. Calcite twin data were collected from pre-folding veins in late Cretaceous limestones from the Ionian zone in Albania in order to (1) determine Paleogene-Neogene stresses associated with the development of the major vein sets in the frontal anticlines of the Outer Albanides and (2) estimate paleoburial of the Cretaceous reservoir rocks during pre-folding flexural subsidence of the foreland. The first vein set (set I) trends N140 (± 20) and the second set (set II) is oriented N060 (± 20). Calcite twinning analysis from set I veins reveals a pre-folding N030° extension likely related to foreland flexure; a later pre-folding, NE-directed compression (LPS) is identified either from one or from both vein sets in the samples from the Saranda anticline; this NE compression is instead recorded by twinning in set II veins from the Kremenara anticline during late stage fold tightening. This NE compression well agrees with independent microtectonic data, regional transport direction and contemporary stress. The differential stress values related to this NE compression are combined with the hypothesis of crustal frictional stress equilibrium to derive first-order estimates of paleoburial of the Cretaceous limestones just before they were uplifted by folding. The ~ 4 km paleoburial of these limestones estimated in the Saranda anticline is consistent with independent paleoburial estimates from stratigraphy, maturity rank of organic matter, paleotemperature/paleogeothermal gradients from fluid inclusions and predictions of kinematic modelling of the Albanian foreland. Our results therefore place reliable constraints on the amount and rate of vertical uplift of these Cretaceous limestones and yield a promising methodology for better constraining paleoburial and therefore erosion and uplift in fold-thrust belts.     

Impact of vertical structure on water mass circulation in a tropical lagoon (Ebrié, Ivory Coast)

A one-dimensional vertical model has been developed to simulate the water mass circulation along the vertical structure in all deep coastal areas. The model has hydrodynamic and transport components solved using finite difference scheme. The one-dimensional vertical model results are coupled to the vertically averaged two-dimensional model results at each point of a horizontal grid. A theoretical salinity profile is introduced for each vertically integrated value obtained from the 2DH model results. A viscosity profile, simulating a viscosity value close to zero at the surface and with large viscosity gradients, is applied along the water column. The model is applied to the Vridi channel, connecting the Ebrié lagoon to the sea (Ivory Coast).The response of the Ebrié lagoon is studied in terms of inflow and outflow of water in the system through the Vridi channel. Due to the abrupt variation of the surface slope, vertical velocities along the water column show an anticlockwise spiral from bottom to surface during a tidal cycle. Due to the bottom friction and to the vertical viscosity profile, velocities decrease from surface to bottom. However, the freshwater inflow slows down the tidal propagation during the flood and causes the surface velocity to be smaller than the bottom velocity at mid-tide. Close to the bottom, velocities follow an anticlockwise movement due to the tidal propagation. At the water surface, velocities follow only an alternative movement of either ebb or flood, along the channel direction. No cross shore velocities can develop at the surface in the channel.     

Textural Characterization Using Gy’s Heterogeneity Functions

A solid can be regarded as a set of contiguous elementary units. The distribution within the solid of any properties, measurable within each elementary unit, can be characterized using two parameters. These parameters are built using the constitution and distribution heterogeneities of Gy. The former account for the granularity of the elementary units, whereas the latter assess the spatial distribution of the property. A texture which definition involves several properties can be described using a diagram where both parameters work as variables. Potential applications encompass (i) the textural classification of soils, ore, breccia, and concrete and (ii) the monitoring of textural transformation during process like dolomitization, metamorphism, weathering, deformation, or annealing.     

A grassed waterway and earthen dams to control muddy floods from a cultivated catchment of the Belgian loess belt

Muddy floods, i.e. runoff from cultivated areas carrying large quantities of soil, are frequent and widespread in the European loess belt. They are mainly generated in dry zero-order valleys and are nowadays considered as the most likely process transferring material eroded from cultivated hillslopes during the Holocene to the flood plain. The huge costs of muddy flood damages justify the urgent installation of control measures. In the framework of the ‘Soil Erosion Decree’ of the Belgian Flemish region, a 12 ha-grassed waterway and three earthen dams have been installed between 2002–2004 in the thalweg of a 300-ha cultivated dry valley in the Belgian loess belt. The measures served their purpose by preventing any muddy flood in the downstream village, despite the occurrence of several extreme rainfall events (with a maximum return period of 150 years). The catchment has been intensively monitored from 2005–2007 and 39 runoff events were recorded in that period. Peak discharge (per ha) was reduced by 69% between the upstream and the downstream extremities of the grassed waterway (GWW). Furthermore, runoff was buffered for 5–12 h behind the dams, and the lag time at the outlet of the catchment was thereby increased by 75%. Reinfiltration was also observed within the waterway, runoff coefficients decreasing by a mean of 50% between both extremities of the GWW. Sediment discharge was also reduced by 93% between the GWW's inflow and the outlet. Before the installation of the control measures, specific sediment yield (SSY) of the catchment reached 3.5 t ha^− 1 yr^− 1 and an ephemeral gully was observed nearly each year in the catchment. Since the control measures have been installed, no (ephemeral) gully has developed and the SSY of the catchment dropped to a mean of 0.5 t ha^− 1 yr^− 1. Hence, sediment transfer from the cultivated dry valley to the alluvial plain should dramatically decrease. Total cost of the control measures that are built for a 20 year-period is very low (126 € ha^− 1) compared to the mean damage cost associated with muddy floods in the study area (54 € ha^− 1 yr^− 1). Similar measures should therefore be installed to protect other flooded villages of the Belgian loess belt and comparable environments.     

Operational performance of current synthetic aperture radar sensors in mapping soil surface characteristics in agricultural environments: application to hydrological and erosion modelling

Synthetic aperture radar (SAR) sensors are often used to characterize the surface of bare soils in agricultural environments. They enable the soil moisture and roughness to be estimated with constraints linked to the configurations of the sensors (polarization, incidence angle and radar wavelength). These key soil characteristics are necessary for different applications, such as hydrology and risk prediction. This article reviews the potential of currently operational SAR sensors and those planned for the near future to characterize soil surface as a function of users' needs. It details what it is possible to achieve in terms of mapping soil moisture and roughness by specifying optimal radar configurations and the precision associated with the estimation of soil surface characteristics. The summary carried out for the present article shows that mapping soil moisture is optimal with SAR sensors at low incidence angles (<35 ). This configuration, which enables an estimated moisture accuracy greater than 6% is possible several times a month taking into account all the current and future sensors. Concerning soil roughness, it is best mapped using three classes (smooth, moderately rough, and rough). Such mapping requires high‐incidence data, which is possible with certain current sensors (RADARSAT‐1 and ASAR both in band C). When L‐band sensors (ALOS) become available, this mapping accuracy should improve because the sensitivity of the radar signal to Soil Surface Characteristics (SSC) increases with wavelength. Finally, the polarimetric mode of certain imminent sensors (ALOS, RADARSAT‐2, TerraSAR‐X, etc.), and the possibility of acquiring data at very high spatial resolution (metre scale), offer great potential in terms of improving the quality of SSC mapping. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental modelling of shallow magma emplacement: Application to saucer-shaped intrusions

Saucer-shaped dolerite and sandstone intrusions are common in sedimentary basins world-wide. We have conducted a series of scaled experiments simulating the process of magma emplacement in sedimentary basins, with particular attention on the formation of saucer-shaped sills. The model materials were (1) cohesive fine-grained silica flour, representing brittle crust; and (2) molten low-viscosity oil, representing magma. The experiments were performed in both homogeneous and layered models. In all the experiments, oil injection resulted in doming of the surface. In the homogeneous models, the injected oil formed cone sheets and sub-vertical dykes. Cone sheets formed for shallow injection (1–3 cm), and vertical dykes formed for deeper injection (4–5 cm). In layered models, the injected oil always formed saucer-shaped intrusions. Our experimental results show that (1) sill intrusion results in the formation of a dome, with melt erupting at the rim; (2) layering controls the formation of sills and saucer-shaped sills; (3) saucer-shaped sills are fed from the bottom and the fluid flows upward and outward; and (4) the diameter of saucer-shaped sills increase with increasing emplacement depth. The systematic relation between domes and sills and the depth-dependence of sill diameters show that saucer-shaped intrusions result from the interaction between a growing flat-lying shallow sill and doming of the free surface. We conclude that saucer-shaped intrusions represent fundamental geometries formed by shallow magma intrusion in stratified basins.     

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.     

The NAFE’06 data set: Towards soil moisture retrieval at intermediate resolution

The National Airborne Field Experiment 2006 (NAFE’06) was conducted during a three week period of November 2006 in the Murrumbidgee River catchment, located in southeastern Australia. One objective of NAFE’06 was to explore the suitability of the area for SMOS (Soil Moisture and Ocean Salinity) calibration/validation and develop downscaling and assimilation techniques for when SMOS does come on line. Airborne L-band brightness temperature was mapped at 1 km resolution 11 times (every 1–3 days) over a 40 by 55 km area in the Yanco region and 3 times over a 40 by 50 km area that includes Kyeamba Creek catchment. Moreover, multi-resolution, multi-angle and multi-spectral airborne data including surface temperature, surface reflectance (green, read and near infrared), lidar data and aerial photos were acquired over selected areas to develop downscaling algorithms and test multi-angle and multi-spectral retrieval approaches. The near-surface soil moisture was measured extensively on the ground in eight sampling areas concurrently with aircraft flights, and the soil moisture profile was continuously monitored at 41 sites. Preliminary analyses indicate that (i) the uncertainty of a single ground measurement was typically less than 5% vol. (ii) the spatial variability of ground measurements at 1 km resolution was up to 10% vol. and (iii) the validation of 1 km resolution L-band data is facilitated by selecting pixels with a spatial soil moisture variability lower than the point-scale uncertainty. The sensitivity of passive microwave and thermal data is also compared at 1 km resolution to illustrate the multi-spectral synergy for soil moisture monitoring at improved accuracy and resolution. The data described in this paper are available at www.nafe.unimelb.edu.au.     

Zinc stable isotopes in seafloor hydrothermal vent fluids and chimneys

Many of the heaviest and lightest natural zinc (Zn) isotope ratios have been discovered in hydrothermal ore deposits. However, the processes responsible for fractionating Zn isotopes in hydrothermal systems are poorly understood. In order to better assess the total range of Zn isotopes in hydrothermal systems and to understand the factors which are responsible for this isotopic fractionation, we have measured Zn isotopes in seafloor hydrothermal fluids from numerous vents at 9–10°N and 21°N on the East Pacific Rise (EPR), the TAG hydrothermal field on the Mid-Atlantic Ridge, and in the Guaymas Basin. Fluid δ⁶⁶Zn values measured at these sites range from + 0.00‰ to + 1.04‰. Of the many physical and chemical parameters examined, only temperature was found to correlate with fluid δ⁶⁶Zn values. Lower temperature fluids (< 250 °C) had both heavier and more variable δ⁶⁶Zn values compared to higher temperature fluids from the same hydrothermal fields. We suggest that subsurface cooling of hydrothermal fluids leads to precipitation of isotopically light sphalerite (Zn sulfide), and that this process is a primary cause of Zn isotope variation in hydrothermal fluids. Thermodynamic calculations carried out to determine saturation state of sphalerite in the vent fluids support this hypothesis with isotopically heaviest Zn found in fluids that were calculated to be saturated with respect to sphalerite. We have also measured Zn isotopes in chimney sulfides recovered from a high-temperature (383 °C) and a low-temperature (203 °C) vent at 9–10°N on the EPR and, in both cases, found that the δ⁶⁶Zn of chimney minerals was lighter or similar to the fluid δ⁶⁶Zn. The first measurements of Zn isotopes in hydrothermal fluids have revealed large variations in hydrothermal fluid δ⁶⁶Zn, and suggest that subsurface Zn sulfide precipitation is a primary factor in causing variations in fluid δ⁶⁶Zn. By understanding how chemical processes that occur beneath the seafloor affect hydrothermal fluid δ⁶⁶Zn, Zn isotopes may be used as a tracer for studying hydrothermal processes.