Distribution functions of spring discharges according to their lithologies and the influence of lower limit to flow: an example from Spain

Specific functions of distribution of the number of springs and of their contributions according to flow for each of nine lithological groups are established. These functions confirm the suitability of a general type of distribution (Sanz 1996), which is a borderline case of the lognormal distribution. The statistical methods proposed for the lithological types presented can be applied to any region if data on flow and lithology of the aquifers drained by springs are available. These methods have been applied to Spain, a representative region with varied geology, climate and topography; 71.2% of spring flow is supplied by limestones, 19.17% by alluvial sediments and marls, 6.7% by conglomerates and sandstones, and 3% by slates, plutonic rocks, quartzites, and other groups. Springs with discharge rates exceeding 2000 L/s exist only in limestones. The majority of springs with low flow occur in marls. If we consider springs with flow greater than 0.5 L/s, alluvial sediments, sandstones, marls, and limestones have the greatest density of springs per surface area, although the average flow varies greatly from one lithology to the next. The cumulative estimated discharge and the estimated number of springs depend upon the definition of the lower limit of flow. For example, total discharge is reduced by 30% if we eliminate all springs with flow lower than 0.5 L/min.     

Transport of three herbicides in ground water at Twin Lake test site,Chalk River,Ontario, Canada

A field tracer test performed under natural flow conditions at the Twin Lake test site, Chalk River Laboratories of the Atomic Energy of Canada Ltd. in Chalk River, Ontario, Canada, using tritium and three herbicides (Chlortoluron, Terbuthylazine, and Pendimethalin) was interpreted using the dispersion equation with a combined reaction model. The reaction model couples an instantaneous equilibrium reaction governed by a linear adsorption isotherm with a reversible or irreversible kinetic reaction of the first order, and decay. An improved interpretation method consists of a simultaneous fitting of theoretical concentration and mass-recovery curves to the experimental data, which leads to a more reliable determining of reaction models and improves the accuracy of fitting. Tritium served as the reference tracer to determine the flow velocity, dispersivity, and the recovery of the herbicides. Chlortoluron was slightly delayed by equilibrium exchange with strongly reduced concentration due to an irreversible kinetic reaction and/or decay. Terbuthilazine was slightly delayed by equilibrium exchange, with strongly reduced concentration due to a reversible kinetic reaction with some influence of decay. A strong equilibrium reaction and a strong reversible kinetic reaction without degradation governed the transport of Pendimethalin, reducing considerably its concentration. The results obtained show that simulations based only on Kd and decay constant, especially if these parameters are found in the laboratory, may considerably differ from those performed with reaction parameters determined in properly performed field tests. The dominant reaction types, and the values of parameters found in the study, supply useful information on the transport of the investigated herbicides in sandy aquifers under natural flow conditions.     

A data‐based mechanistic modelling (DBM) approach to understanding dynamic sediment transmission through Wyresdale Park Reservoir,Lancashire, UK

This paper outlines the application of a new data‐based mechanistic (DBM) modelling methodology to the characterization of the sediment transmission dynamics in a small upland reservoir, Wyresdale Park, Lancashire. The DBM modelling strategy exploits advanced statistical procedures to infer the dynamic model structure and its associated parameters directly from the instrumented data, producing a parametrically efficient, continuous time, transfer function model which relates suspended sediment load at the reservoir inflow to the outflow at the event scale. The associated differential equation model parameters have physical attributes which can be interpreted in terms of sediment transmission processes and associated reservoir trap efficiency. Sedigraph analysis suggests that wind‐induced resuspension episodically supplies an additional load to the reservoir outlet. The stochastic nature of the DBM model makes it ideal for evaluating the effects of uncertainty through Monte Carlo simulations (MCS) for discharge and sediment transmission. Copyright © 2000 John Wiley & Sons, Ltd.     

Flow prediction in high altitude ungauged catchments: A case study in the Italian Alps (Pantano Basin,Adamello Group)

We develop a simple model to evaluate the daily flow discharges in the ablation season for the 11 km² Pantano basin in the Retiche Italian Alps, based upon the data gathered during a three years field campaign. The Pantano basin embeds the Venerocolo debris covered and the Avio debris free glaciers, covering 2.14 km² in the Adamello Group, where the widest Italian glacier Adamello is located. First, degree-day models based upon air temperature are tuned to calculate snow and ice melt at daily scale. Glaciers’ meteorological data are collected from an automatic weather station (AWS), operating on the glacier during summer 2007. The melt factors in the debris covered areas of the glacier are estimated against debris thickness, using a data driven parameterization. The flow discharge from the catchment is estimated using semi distributed flow routing for the ablation seasons of four years, from 2006 to 2009. The predicted discharges are compared to those derived from inverse reservoir's routing at the Benedetto lake, catching the basin outflow. The proposed approach is valuable as a tool to investigate the hydrology of poorly gauged glacierized areas, including those with debris covered ice, widely diffused and yet poorly understood. Pending accurate parameterization the approach is usable for water resources evaluation and for long term assessment of the climate change impact on the glacierized areas within the Alps.     

Lake Winnipeg: geological setting and sediment seismostratigraphy

Lake Winnipeg, the seventh largest lake in North America, is located at the boundary between the Interior Plains and the Canadian Shield in Manitoba, Canada. Seismic profiles were obtained in Lake Winnipeg on two geoscientific cruises in 1994 and 1996. These data indicate the morphology of the bedrock surface. In most cases, a clear distinction between low relief Paleozoic carbonate rock and high relief Precambrian rock can be made. In northern Lake Winnipeg, the eastern limit of Paleozoic rock is clearly demarcated 30 km west of the previous estimate of its position. In southern Lake Winnipeg, all or most of the Paleozoic sequence terminates at a prominent buried escarpment in the centre of the lake. This indicates that Paleozoic rock on the eastern shore, known from drilling and outcrops, is an outlier. Major moraines are apparent as abrupt, large ridges having a chaotic internal reflection pattern. These include the Pearson Reef Moraine, the George Island Moraine and the offshore extension of The Pas Moraine. Little evidence for extensive or thick till was observed. Instead, fine-grained sediments deposited in glacial Lake Agassiz rest directly on bedrock over most of the lake basin. Hence an episode of erosion to bedrock was associated with glaciation and/or deglaciation. The Agassiz Sequence sediments are well-stratified, drape underlying relief and in some areas are over 100 m thick. In places, stratification in these sediments is disrupted, perhaps by dewatering. Evidence of erosion of Agassiz Sequence sediments by recent currents was observed. The contact between the Agassiz Sequence and the overlying Winnipeg Sequence sediments is a marked angular unconformity. The Agassiz Unconformity indicates up to 10 m of erosion in places. The low-relief character of this unconformity precludes subaerial erosion and the lack of till, moraines, or extensive deformation precludes glacial erosion. Waves appear to be the most likely erosional agent, either in waning Lake Agassiz or early Lake Winnipeg time. Winnipeg Sequence sediments, in places very thin, mantle most of the lakefloor. These sediments were deposited in the present Lake Winnipeg and are faintly stratified to massive and reach about 10 m in thickness in deep water. On the surface of the Winnipeg Sequence, vigorous, episodic currents are thought to contribute to the construction of flow-transverse sand waves as much as 6 m high in a deep, narrow constriction in the lake.     

Optimized and robust experimental design: a non-linear application to EM sounding

Pragmatic experimental design requires objective consideration of several classes of information including the survey goals, the range of expected Earth responses, acquisition costs, instrumental capabilities, experimental conditions and logistics. In this study we consider the ramifications of maximizing model parameter resolution through non-linear experimental design. Global optimization theory is employed to examine and rank different EM sounding survey designs in terms of model resolution as defined by linearized inverse theory. By studying both theoretically optimal and heuristic experimental survey configurations for various quantities of data, it is shown that design optimization is critical for minimizing model variance estimates, and is particularly important when the inverse problem becomes nearly underdetermined. We introduce the concept of robustness so that survey designs are relatively immune to the presence of potential bias errors in important data. Bias may arise during practical measurement, or from designing a survey using an appropriate model.     

The Taconite Inlet Lakes Project: a systems approach to paleoclimatic reconstruction

A comprehensive study of meteorological, hydrological, limnological and sedimentological conditions in the watersheds of density-stratified (meromictic) lakes around Taconite Inlet, Northern Ellesmere Island, N.W.T., Canada was carried out from 1990–1992. Lakes C1 and C2 contain seawater trapped by isostatic uplift as the former embayments became isolated from the sea. These lakes, and Lake C3, contain varved sediments which provide an annually resolvable paleoclimatic record. By studing the major systems influencing sedimentation in one of these lakes (Lake C2) a better understanding of the climatic controls on varve formation, and hence on the paleoclimatic signal in the varved sediment record, was obtained. The varves of Lake C2 provide a proxy record of summer temperature for the region.This is the first in a series of papers published in this issue on the Taconite Inlet Lakes Project. These papers were collected by Dr R. S. Bradley.     

Ccte d'Ivoire-Ghana margin: seismic imaging of passive rifted crust adjacent to a transform continental margin

During May 1990 and January-February 1991, an extensive geophysical data set was collected over the Côte d'Ivoire-Ghana continental margin, located along the equatorial coast of West Africa. The Ghana margin is a transform continental margin running subparallel to the Romanche Fracture Zone and its associated marginal ridge—the Côte d'Ivoire-Ghana Ridge. From this data set, an explosive refraction line running ∼ 150 km, ENE-WSW between 3°55'N, 3°21'W and 4°23'N, 2°4'W, has been modelled together with wide-angle airgun profiles, and seismic reflection and gravity data. This study is centred on the Côte d'Ivoire Basin located just to the north of the Côte d'Ivoire-Ghana Ridge, where bathymetric data suggest that a component of normal rifting occurred, rather than the transform motion observed along the majority of the equatorial West African margin.
Traveltime and amplitude modelling of the ocean-bottom seismometer data shows that the continental Moho beneath the margin rises in an oceanward direction, from ∼ 24 km below sea level to ∼ 17 km. In the centre of the line where the crust thins most rapidly, there exists a region of anomalously high velocity at the base of the crust, reaching some 8 km in thickness. This higher-velocity region is thought to represent an area of localized underplating related to rifting. Modelling of marine gravity data, collected coincident with the seismic line, has been used to test the best-fitting seismic model. This modelling has shown that the observed free-air anomaly is dominated by the effects of crustal thickness, and that a region of higher density is required at the base of the crust to fit the observed data. This higher-density region is consistent in size and location with the high velocities required to fit the seismic data.     

Periodicity of magnetic intensities in magnetic anomaly profiles: the Cenozoic of the South Atlantic

Spectral analyses of several published magnetic anomaly profiles from Candé & Kent (1992a) were undertaken prior to analysing, in the same way, raw magnetic anomaly data from similar parts of the South Atlantic. It was found that similar and distinct medium and short wavelengths were present in both the published and raw data. When these are converted into the time domain using the average rate of spreading for each profile, these periodicities appear similar, possibly identical, to those expected from the long-term eccentricity orbital parameters (Fischer, DeBoer & Premoli Silva 1990). While such correlations are not necessarily causative, they suggest that magnetohydro-dynamical processes near the core-mantle boundary may be affected by gravitational changes due to planetary orbital perturbations.     

Calculating the two-dimensional magnetotelluric Jacobian in finite elements using reciprocity

To speed up the calculation of the field Jacobian for 2-D magnetoteliuric inversion using finite elements, the principle of electromagnetic reciprocity is applied. The governing relationship for the Jacobian of the field along strike is obtained by differentiating the Helmholtz equation with respect to the resistivity of each region in the finite-element mesh. The result is a similar Helmholtz equation for the Jacobian, with new sources distributed over all nodes within the parameter medium. However, according to the principle of electromagnetic reciprocity, the roles of sources and receivers are interchangeable. Utilizing reciprocity, the field values obtained from the original forward problem and for new unit sources imposed at the receivers are then utilized in the calculation of the Jacobian by simple multiplication and summation with finite-element terms at each rectangle in the mesh. For the auxiliary (across-strike) fields, the Jacobian terms are obtained by solving source vectors loaded with parabola coefficients used in the approximation to Maxwell's equations. Jacobian terms for the apparent resistivity ( p _a), the impedance phase (φ) and the vertical magnetic field ( K _zy) are then calculated utilizing the parallel- and auxiliary-field Jacobians. Comparison of Jacobian values obtained from reciprocity calculations and by differencing two forward solutions show that the reciprocity method is accurate and can be used to decrease the number of calculations required to obtain sensitivities by one to two orders of magnitude.