Verification of a coupled climate-hydrological model against Holocene palaeohydrological records

We have coupled a climate model (ECBilt-CLIO-VECODE) and a hydrological model (STREAM) offline to simulate palaeodischarge of nineteen rivers (Amazon, Congo, Danube, Ganges, Krishna, Lena, Mackenzie, Mekong, Meuse, Mississippi, Murray–Darling, Nile, Oder, Rhine, Sacramento–San Joaquin, Syr Darya, Volga, Volta, Zambezi) for three time-slices: Early Holocene (9000–8650 BP), Mid-Holocene (6200–5850 BP) and Recent (1750–2000 AD). To evaluate the model's skill in retrodicting broad changes in mean palaeodischarge we have compared the model results with palaeodischarge estimates from multi-proxy records. We have compared the general trends inferred from the proxy data with statistical differences in modelled discharge between the three periods, thereby developing a technique to assess the level of agreement between the model and proxy data. The quality of the proxy data for each basin has been classed as good, reasonable or low. Of the model runs for which the proxy data were good or reasonable, 72% were in good agreement with the proxy data, and 92% were in at least reasonable agreement. We conclude that the coupled climate-hydrological model performs well in simulating mean discharge in the time-slices studied. The discharge trends inferred from the proxy and model data closely follow latitudinal and seasonal variations in insolation over the Holocene. For a number of basins for which agreement was not good we have identified specific mechanisms which could be responsible for the discrepancy, primarily the absence of the Laurentide ice sheet in our model. In order to use the model in an operational sense within water management studies it would be useful to use a higher spatial resolution and a daily time-step.     

Holocene glacier and deep water dynamics,Adélie Land region,East Antarctica

This study presents a high-resolution multi-proxy investigation of sediment core MD03-2601 and documents major glacier oscillations and deep water activity during the Holocene in the Adélie Land region, East Antarctica. A comparison with surface ocean conditions reveals synchronous changes of glaciers, sea ice and deep water formation at Milankovitch and sub-Milankovitch time scales. We report (1) a deglaciation of the Adélie Land continental shelf from 11 to 8.5 cal ka BP, which occurred in two phases of effective glacier grounding-line retreat at 10.6 and 9 cal ka BP, associated with active deep water formation; (2) a rapid glacier and sea ice readvance centred around 7.7 cal ka BP; and (3) five rapid expansions of the glacier–sea ice systems, during the Mid to Late Holocene, associated to a long-term increase of deep water formation. At Milankovich time scales, we show that the precessionnal component of insolation at high and low latitudes explains the major trend of the glacier–sea ice–ocean system throughout the Holocene, in the Adélie Land region. In addition, the orbitally-forced seasonality seems to control the coastal deep water formation via the sea ice–ocean coupling, which could lead to opposite patterns between north and south high latitudes during the Mid to Late Holocene. At sub-Milankovitch time scales, there are eight events of glacier–sea ice retreat and expansion that occurred during atmospheric cooling events over East Antarctica. Comparisons of our results with other peri-Antarctic records and model simulations from high southern latitudes may suggest that our interpretation on glacier–sea ice–ocean interactions and their Holocene evolutions reflect a more global Antarctic Holocene pattern.     

Isotopic evidence for the origin of an acid sulphate alteration,Styrian basin,Austria

The Gossendorf volcanic body is the only one in the Styrian basin that shows extensive hydrothermal alteration. K‐Ar dating of primary volcanic biotite and alteration products (alunite) suggests that the emplacement of the volcanic body and hydrothermal alteration took place synchronously, 15 Ma ago. The stable isotope compositions of the alteration products such as opal, barite, pyrite and alunite combined with crystallographia investigations indicate temperatures between 150 and 200 °C for the formation of the alteration zones. The calculated stable isotopic compositions of the parent fluid, responsible for the alteration, show an exogene marine component, which interacted with the host rock. Sulphur isotopic compositions of sulphur, sulphides and sulphates indicate disequilibrium, and progressive oxidation. This fact, combined with the mineral zonation of the alteration zone, reflects not only change in the pH but also change in the fO₂ of the ascending fluids.     

The crustal evolution of South America from a zircon Hf‐isotope perspective

Hf‐isotope data of >1100 detrital zircon grains from the Palaeozoic, south‐central Andean Gondwana margin record the complete crustal evolution of South America, which was the predominant source. The oldest grains, with crustal residence ages of 3.8–4.0 Ga, are consistent with complete recycling of existing continental crust around 4 Ga. We confirm three major Archaean, Palaeoproterozoic (Transamazonian) and late Mesoproterozoic to early Neoproterozoic crust‐addition phases as well as six igneous phases during Proterozoic to Palaeozoic time involving mixing of juvenile and crustally reworked material. A late Mesoproterozoic to early Neoproterozoic, Grenville‐age igneous belt can be postulated along the palaeo‐margin of South America. This belt was the basement for later magmatic arcs and accreted allochthonous microcontinents as recorded by similar crustal residence ages. Crustal reworking likely dominated over juvenile addition during the Palaeozoic era, and Proterozoic and Archaean zircon was mainly crustally reworked from the eroding, thickened Ordovician Famatinian arc.     

On the causes of pressure oscillations in low‐permeable and low‐compressible porous media

Nonphysical pressure oscillations are observed in finite element calculations of Biot's poroelastic equations in low‐permeable media. These pressure oscillations may be understood as a failure of compatibility between the finite element spaces, rather than elastic locking. We present evidence to support this view by comparing and contrasting the pressure oscillations in low‐permeable porous media with those in low‐compressible porous media. As a consequence, it is possible to use established families of stable mixed elements as candidates for choosing finite element spaces for Biot's equations. Copyright © 2011 John Wiley & Sons, Ltd.     

Preliminary evaluation of the runoff processes in a remote montane cloud forest basin using Mixing Model Analysis and Mean Transit Time

In this study, the Mean Transit Time and Mixing Model Analysis methods are combined to unravel the runoff generation process of the San Francisco River basin (73.5 km²) situated on the Amazonian side of the Cordillera Real in the southernmost Andes of Ecuador. The montane basin is covered with cloud forest, sub‐páramo, pasture and ferns. Nested sampling was applied for the collection of streamwater samples and discharge measurements in the main tributaries and outlet of the basin, and for the collection of soil and rock water samples. Weekly to biweekly water grab samples were taken at all stations in the period April 2007–November 2008. Hydrometric data, Mean Transit Time and Mixing Model Analysis allowed preliminary evaluation of the processes controlling the runoff in the San Francisco River basin. Results suggest that flow during dry conditions mainly consists of lateral flow through the C‐horizon and cracks in the top weathered bedrock layer, and that all subcatchments have an important contribution of this deep water to runoff, no matter whether pristine or deforested. During normal to low precipitation intensities, when antecedent soil moisture conditions favour water infiltration, vertical flow paths to deeper soil horizons with subsequent lateral subsurface flow contribute most to streamflow. Under wet conditions in forested catchments, streamflow is controlled by near surface lateral flow through the organic horizon. Exceptionally, saturation excess overland flow occurs. By absence of the litter layer in pasture, streamflow under wet conditions originates from the A horizon, and overland flow. Copyright © 2011 John Wiley & Sons, Ltd.     

Quaternary denudation of southern Fennoscandia – evidence from the marine realm

Throughout the last 1.1 million years repeated glaciations have modified the southern Fennoscandian landscape and the neighbouring continental shelf into their present form. The glacigenic erosion products derived from the Fennoscandian landmasses were transported to the northern North Sea and the SE Nordic Seas continental margin. The prominent sub‐marine Norwegian Channel trough, along the south coast of Norway, was the main transport route for the erosion products between 1.1 and 0.0 Ma. Most of these erosion products were deposited in the North Sea Fan, which reaches a maximum thickness of 1500 m and has nearly 40 000 km³ of sediments. About 90% of the North Sea Fan sediments have been deposited during the last 500 000 years, in a time period when fast‐moving ice streams occupied the Norwegian Channel during each glacial stage. Back‐stripping the sediment volumes in the northern North Sea and SE Nordic Seas sink areas, including the North Sea Fan, to their assumed Fennoscandian source area gives an average vertical erosion of 164 m for the 1.1–0.0 Ma time period. The average 1.1–0.0 Ma erosion rate in the Fennoscandian source area is estimated to be 0.15 mm a^?1. We suggest, however, that large variations in erosion rates have existed through time and that the most intense Fennoscandian landscape denudation occurred during the time period of repeated shelf edge ice advances, namely from Marine Isotope Stage 12 (c. 0.5 Ma) onwards.     

Observations of Comet 9P/Tempel 1 around the Deep Impact event by the OSIRIS cameras onboard Rosetta

The OSIRIS cameras on the Rosetta spacecraft observed Comet 9P/Tempel 1 from 5 days before to 10 days after it was hit by the Deep Impact projectile. The Narrow Angle Camera (NAC) monitored the cometary dust in 5 different filters. The Wide Angle Camera (WAC) observed through filters sensitive to emissions from OH, CN, Na, and OI together with the associated continuum. Before and after the impact the comet showed regular variations in intensity. The period of the brightness changes is consistent with the rotation period of Tempel 1. The overall brightness of Tempel 1 decreased by about 10% during the OSIRIS observations. The analysis of the impact ejecta shows that no new permanent coma structures were created by the impact. Most of the material moved with . Much of it left the comet in the form of icy grains which sublimated and fragmented within the first hour after the impact. The light curve of the comet after the impact and the amount of material leaving the comet ( of water ice and a presumably larger amount of dust) suggest that the impact ejecta were quickly accelerated by collisions with gas molecules. Therefore, the motion of the bulk of the ejecta cannot be described by ballistic trajectories, and the validity of determinations of the density and tensile strength of the nucleus of Tempel 1 with models using ballistic ejection of particles is uncertain.     

Dynamics of planets in retrograde mean motion resonance

In a previous paper (Gayon and Bois 2008a), we have shown the general efficiency of retrograde resonances for stabilizing compact planetary systems. Such retrograde resonances can be found when two-planets of a three-body planetary system are both in mean motion resonance and revolve in opposite directions. For a particular two-planet system, we have also obtained a new orbital fit involving such a counter-revolving configuration and consistent with the observational data. In the present paper, we analytically investigate the three-body problem in this particular case of retrograde resonances. We therefore define a new set of canonical variables allowing to express correctly the resonance angles and obtain the Hamiltonian of a system harboring planets revolving in opposite directions. The acquiring of an analytical “rail” may notably contribute to a deeper understanding of our numerical investigation and provides the major structures related to the stability properties. A comparison between our analytical and numerical results is also carried out.