Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Oxygen isotopic constraints on the origin of Mg-rich olivines from chondritic meteorites

Chondrules are the major high temperature components of chondritic meteorites which accreted a few millions years after the oldest solids of the solar system, the calcium–aluminum-rich inclusions, were condensed from the nebula gas. Chondrules formed during brief heating events by incomplete melting of solid dust precursors in the protoplanetary disk. Petrographic, compositional and isotopic arguments allowed the identification of metal-bearing Mg-rich olivine aggregates among the precursors of magnesian type I chondrules. Two very different settings can be considered for the formation of these Mg-rich olivines: either a nebular setting corresponding mostly to condensation–evaporation processes in the nebular gas or a planetary setting corresponding mostly to differentiation processes in a planetesimal. An ion microprobe survey of Mg-rich olivines of a set of type I chondrules and isolated olivines from unequilibrated ordinary chondrites and carbonaceous chondrites revealed the existence of several modes in the distribution of the ?¹⁷O values and the presence of a large range of mass fractionation (several ‰) within each mode. The chemistry and the oxygen isotopic compositions indicate that Mg-rich olivines are unlikely to be of nebular origin (i.e., solar nebula condensates) but are more likely debris of broken differentiated planetesimals (each of them being characterized by a given ?¹⁷O). Mg-rich olivines could have crystallized from magma ocean-like environments on partially molten planetesimals undergoing metal–silicate differentiation processes. Considering the very old age of chondrules, Mg-rich olivine grains or aggregates might be considered as millimeter-sized fragments from disrupted first-generation differentiated planetesimals. Finally, the finding of only a small number of discrete ?¹⁷O modes for Mg-rich olivines grains or aggregates in a given chondrite suggests that these shattered fragments have not been efficiently mixed in the disk and/or that chondrite formation occurred in the first vicinity of the breakup of these planetary bodies.     

The lacustrine microbial carbonate factory of the successive Lake Bonneville and Great Salt Lake,Utah, USA

The Bonneville Basin is a continental lacustrine system accommodating extensive microbial carbonate deposits corresponding to two distinct phases: the deep Lake Bonneville (30 000 to 11 500 ¹⁴C bp ) and the shallow Great Salt Lake (since 11 500 ¹⁴C bp ). A characterization of these microbial deposits and their associated sediments provides insights into their spatio‐temporal distribution patterns. The Bonneville phase preferentially displays vertical distribution of the microbial deposits resulting from high‐amplitude lake level variations. Due to the basin physiography, the microbial deposits were restricted to a narrow shoreline belt following Bonneville lake level variations. Carbonate production was more efficient during intervals of relative lake level stability as recorded by the formation of successive terraces. In contrast, the Great Salt Lake microbial deposits showed a great lateral distribution, linked to the modern flat bottom configuration. A low vertical distribution of the microbial deposits was the result of the shallow water depth combined with a low amplitude of lake level fluctuations. These younger microbial deposits display a higher diversity of fabrics and sizes. They are distributed along an extensive ‘shore to lake’ transect on a flat platform in relation to local and progressive accommodation space changes. Microbial deposits are temporally discontinuous throughout the lake history showing longer hiatuses during the Bonneville phase. The main parameters controlling the rate of carbonate production are related to the interaction between physical (kinetics of the mineral precipitation, lake water temperature and runoff), chemical (Ca²⁺, Mg²⁺ and HCO₃^? concentrations, Mg/Ca ratio, dilution and depletion) and/or biological (trophic) factors. The contrast in evolution of Lake Bonneville and Great Salt Lake microbial deposits during their lacustrine history leads to discussions on major chemical and climatic changes during this interval as well as the role of physiography. Furthermore, it provides novel insights into the composition, structure and formation of microbialite‐rich carbonate deposits under freshwater and hypersaline conditions.     

Capture-zone design in an aquifer influenced by cyclic fluctuations in hydraulic gradients

Design of a groundwater pumping and treatment system for a wood-treatment facility adjacent to the tidally influenced Fraser River estuary required the development of methodologies to account for cyclic variations in hydraulic gradients. Design of such systems must consider the effects of these cyclic fluctuations on the capture of dissolved-phase contaminants. When the period of the cyclic fluctuation is much less than the travel time of the dissolved contaminant from the source to the discharge point, the hydraulic-gradient variations resulting from these cycles can be ignored. Capture zones are then designed based on the average hydraulic gradient determined using filter techniques on continuous groundwater-level measurements. When the period of cyclic fluctuation in hydraulic gradient is near to or greater than the contaminant travel time, the resulting hydraulic-gradient variations cannot be ignored. In these instances, procedures are developed to account for these fluctuations in the capture-zone design. These include proper characterization of the groundwater regime, assessment of the average travel time and period of the cyclic fluctuations, and numerical techniques which allow accounting for the cyclic fluctuations in the design of the capture zone. Electronic Publication     

Modelling interactions between fold–thrust belt deformation, foreland flexure and surface mass transport

Interactions between fold and thrust belt deformation, foreland flexure and surface mass transport are investigated using a newly developed mathematical model incorporating fully dynamic coupling between mechanics and surface processes. The mechanical model is two dimensional (plane strain) and includes an elasto‐visco‐plastic rheology. The evolving model is flexurally compensated using an elastic beam formulation. Erosion and deposition at the surface are treated in a simple manner using a linear diffusion equation. The model is solved with the finite element method using a Lagrangian scheme with marker particles. Because the model is particle based, it enables straightforward tracking of stratigraphy and exhumation paths and it can sustain very large strain. It is thus ideally suited to study deformation, erosion and sedimentation in fold–thrust belts and foreland basins. The model is used to investigate how fold–thrust deformation and foreland basin development is influenced by the non‐dimensional parameter , which can be interpreted as the ratio of the deformation time scale to the time scale for surface processes. Large values of imply that the rate of surface mass transport is significantly greater than the rate of deformation. When , the rates of surface processes are so slow that one observes a classic propagating fold–thrust belt with well‐developed wedge top basins and a largely underfilled foreland flexural depression. Increasing causes (1) deposition to shift progressively from the wedge top into the foredeep, which deepens and may eventually become filled, (2) widespread exhumation of the fold–thrust belt, (3) reduced rates of frontal thrust propagation and possible attainment of a steady‐state orogen width and (4) change in the style and dynamics of deformation. Together, these effects indicate that erosion and sedimentation, rather than passively responding to tectonics, play an active and dynamic role in the development of fold–thrust belts and foreland basins. Results demonstrate that regional differences in the relative rates of surface processes (e.g. because of different climatic settings) may lead to fold–thrust belts and foreland basins with markedly different characteristics. Results also imply that variations in the efficiency of surface processes through time (e.g., because of climate change or the emergence of orogens above sea level) may cause major temporal changes in orogen and basin dynamics.     

Brest sea level record: a time series construction back to the early eighteenth century

The completeness and the accuracy of the Brest sea level time series dating from 1807 make it suitable for long-term sea level trend studies. New data sets were recently discovered in the form of handwritten tabulations, including several decades of the eighteenth century. Sea level observations have been made in Brest since 1679. This paper presents the historical data sets which have been assembled so far. These data sets span approximately 300 years and together constitute the longest, near-continuous set of sea level information in France. However, an important question arises: Can we relate the past and the present-day records? We partially provide an answer to this question by analysing the documents of several historical libraries with the tidal data using a ‘data archaeology’ approach advocated by Woodworth (Geophys Res Lett 26:1589–1592, 1999b). A second question arises concerning the accuracy of such records. Careful editing was undertaken by examining the residuals between tidal predictions and observations. It proved useful to remove the worst effects of timing errors, in particular the sundial correction to be applied prior to August 1, 1714. A refined correction based on sundial literature [Savoie, La gnomique, Editions Les Belles Lettres, Paris, 2001] is proposed, which eliminates the systematic offsets seen in the discrepancies in timing of the sea level measurements. The tidal analysis has also shown that shallow-water tidal harmonics at Brest causes a systematic difference of 0.023 m between mean sea level (MSL) and mean tide level (MTL). Thus, MTL should not be mixed with the time series of MSL because of this systematic offset. The study of the trends in MTL and MSL however indicates that MTL can be used as a proxy for MSL. Three linear trend periods are distinguished in the Brest MTL time series over the period 1807–2004. Our results support the recent findings of Holgate and Woodworth (Geophys Res Lett) of an enhanced coastal sea level rise during the last decade compared to the global estimations of about 1.8 mm/year over longer periods (Douglas, J Geophys Res 96:6981–6992, 1991). The onset of the relatively large global sea level trends observed in the twentieth century is an important question in the science of climate change. Our findings point out to an ‘inflexion point’ at around 1890, which is remarkably close to that in 1880 found in the Liverpool record by Woodworth (Geophys Res Lett 26:1589–1592, 1999b).     

Construction of three-phase data to model multiphase flow in porous media: Comparing an optimization approach to the finite element approach

Multiphase flow modelling is a major issue in the assessment of groundwater pollution. Three-phase flows are commonly governed by mathematical models that associate a pressure equation with two saturation equations. These equations involve a number of secondary variables that reflect the fluid behaviour in a porous medium. To improve the computational efficiency of multiphase flow simulators, several simplified reformulations of three-phase flow equations have been proposed. However, they require the construction of new secondary variables adapted to the reformulated flow equations. In this article, two different approaches are compared to quantify these variables. A numerical example is given for a typical fine sand.     

The retrieval of cloud microphysical properties using satellite measurements and an in situ database

By combining AVHRR data from the NOAA satellites with information from a database of in situ measurements, large-scale maps can be generated of the microphysical parameters most immediately significant for the modelling of global circulation and climate. From the satellite data, the clouds can be classified into cumuliform, stratiform and cirrus classes and then into further sub-classes by cloud top temperature. At the same time a database of in situ measurements made by research aircraft is classified into the same sub-classes and a statistical analysis is used to derive relationships between the sub-classes and the cloud microphysical properties. These two analyses are then linked to give estimates of the microphysical properties of the satellite observed clouds. Examples are given of the application of this technique to derive maps of the probability of occurrence of precipitating clouds and of precipitating water content derived from a case study within the International Cirrus Experiment (ICE) held in 1989 over the North Sea.