Use of unit response functions for management of regional multilayered aquifers: application to the North Aquitaine Tertiary system (France)

Although there are many numerical applications used in aquifer management, few methodologies seem to be sufficiently adapted to provide solutions for regional-scale problems. Their applications for multilayered aquifers are also too limited. A new groundwater management tool adapted to such large physical systems has been developed, using the unit response function (URF) approach. This tool is based on the link between a numerical groundwater flow model and a mathematical optimisation tool. Pre- and post processing are based on a geographical information system (GIS). The developed tool has been applied to a real management case—the regional Aquitaine multilayered aquifer (France). A representative hydrodynamic model was built using MODFLOW 2000 code. First tests show that the linearity condition needed for the application of URF methodology is respected despite the high complexity of the system. Two management scenarios were tested. Optimal solutions thus obtained allow for viable exploitation alternatives to be proposed, which integrate complex environmental constraints.

---

Résumé Malgré un nombre important d’applications existant dans le domaine de la gestion des aquifères, peu de méthodologies semblent adaptées pour proposer des solutions à des problèmes réels de dimension régionale. Leurs application pour des systèmes multicouches sont souvent trop limitées. Un outil de gestion adapté à ces systèmes de large taille a été développé, basé sur la théorie des Fonctions de Réponses Unitaires (FRU). Cet outil est basé sur le couplage d’un modèle numérique d’écoulement et d’un code d’optimisation mathématique. Les entrées et sorties du modèle sont réalisées à l’aide d’un Système d’Information Géographique (SIG). L’outil développé est appliqué à un cas de gestion réel—le système aquifère régional multicouche aquitain (France). Un modèle hydrodynamique représentatif a été construit, basé sur le code MODFLOW 2000. Les premiers tests montrent que les conditions de linéarité nécessaires à l’application de la méthode FRU sont respectées, malgré l’importante complexité du système. Deux scénarii de gestion sont testés. Les solutions optimales ainsi obtenues permettent de proposer des solutions d’exploitation alternatives viables, intégrant des contraintes environnementales complexes.

---

Resumen Aunque hay muchas aplicaciones numéricas usadas en la gestión de acuíferos, unas pocas metodologías parecen estar adaptadas suficientemente para dar las soluciones a los problemas de escala regional. También son limitadas sus aplicaciones para los acuíferos multicapa. Se ha desarrollado una nueva herramienta para el manejo del agua subterránea, adaptada a tales sistemas físicos grandes, usando el avance de la Función de Respuesta Unitaria (URF). Esta herramienta se basa en la combinación entre un modelo numérico de flujo de agua subterránea y una herramienta de optimización matemática. El pre-proceso y el post-proceso se basan en un Sistema de Información Geográfico (GIS). La herramienta desarrollada se ha aplicado a un caso real de manejo en el acuífero regional multicapa de Aquitania (Francia). Se construyó un modelo hidrodinámico representativo, que usó el código MODFLOW 2000. Las primeras pruebas muestran que la condición de linearidad, necesaria para la aplicación de la metodología URF, se respeta a pesar de la complejidad alta del sistema. Se probaron dos escenarios de manejo. Las soluciones óptimas así obtenidas permiten la proyección de alternativas viables de explotación, las cuales integran restricciones medioambientales complejas.

     

Inferring Glacier Mass Balance Using Radarsat: Results From Peyto Glacier, Canada

The capability of RADARSAT synthetic aperture radar (SAR) for the purpose of snow-line/accumulation area mapping for a temperate alpine glacier is examined. In agreement with other orbital C-band SAR studies, RADARSAT can discriminate between firn and bare ice facies. Limited observations are reported with respect to the electromagnetic variability of the ice facies in the ablation area, but they are inconclusive. Operational considerations are discussed with respect to reconciling the uncertainties of late-summer weather and their possible impact on the dielectric and scattering properties of the glacier surface. Vagaries associated with other glacier settings, mass balance states and their associated facies configurations are discussed including the difficulty of using the transient snow-line to define the equilibrium line and the lower extent of the accumulation area for glaciers where superimposed ice may form.
The radar remote-sensing reconnaissance of equilibrium line altitude (ELA) and accumulation area ratio (AAR) for estimating glacier mass balance requires serious consideration in those instances where traditional ground measurements used in the direct glaciological method are absent. However, with respect to the ELA, such estimates can vary depending on the accuracy of the reference digital elevation information. Moreover, for many glacier configurations, where mass balance variations due to altitude are influenced or in some cases completely masked by local balance variations, defining the ELA may be an irreconcilable problem. Using the AAR may be more robust in this regard. It is further determined that the total error inherent in the reconnaissance method would have serious implications for the confident estimation of mass balance normals and climate-related trends if the method were to be utilized over the longer term.     

Conflicting structural and geochronological data from the Ibituruna quartz-syenite (SE Brazil): Effect of protracted “hot” orogeny and slow cooling rate?

The Ibituruna quartz-syenite was emplaced as a sill in the Ribeira-Araçuaí Neoproterozoic belt (Southeastern Brazil) during the last stages of the Gondwana supercontinent amalgamation. We have measured the Anisotropy of Magnetic Susceptibility (AMS) in samples from the Ibituruna sill to unravel its magnetic fabric that is regarded as a proxy for its magmatic fabric. A large magnetic anisotropy, dominantly due to magnetite, and a consistent magnetic fabric have been determined over the entire Ibituruna massif. The magmatic foliation and lineation are strikingly parallel to the solid-state mylonitic foliation and lineation measured in the country-rock. Altogether, these observations suggest that the Ibituruna sill was emplaced during the high temperature (~ 750 °C) regional deformation and was deformed before full solidification coherently with its country-rock. Unexpectedly, geochronological data suggest a rather different conclusion. LA-ICP-MS and SHRIMP ages of zircons from the Ibituruna quartz-syenite are in the range 530–535 Ma and LA-ICP-MS ages of zircons and monazites from synkinematic leucocratic veins in the country-rocks suggest a crystallization at ~ 570–580 Ma, i.e., an HT deformation > 35My older than the emplacement of the Ibituruna quartz-syenite. Conclusions from the structural and the geochronological studies are therefore conflicting. A possible explanation arises from ⁴⁰Ar–³⁹Ar thermochronology. We have dated amphiboles from the quartz-syenite, and amphiboles and biotites from the country-rock. Together with the ages of monazites and zircons in the country-rock, ⁴⁰Ar–³⁹Ar mineral ages suggest a very low cooling rate: < 3 °C/My between 570 and ~ 500 Ma and ~ 5 °C/My between 500 and 460 Ma. Assuming a protracted regional deformation consistent over tens of My, under such stable thermal conditions the fabric and microstructure of deformed rocks may remain almost unchanged even if they underwent and recorded strain pulses separated by long periods of time. This may be a characteristic of slow cooling “hot orogens” that rocks deformed at significantly different periods during the orogeny, but under roughly unchanged temperature conditions, may display almost indiscernible microstructure and fabric.     

Twin structure in coesite studied by high resolution electron microscopy

A synthetic crystal of coesite, SiO₂, has been found to contain microtwins. The (100) reflection twin has been made evident on incoherent (010) planes. The twin structure is directly imaged by high resolution electron microscopy at 2,5 Å resolution, allowing an atomic model to be proposed. This twin could be incoherent on large surfaces and the dislocation density is predicted. The X-ray crystallographic data are analyzed critically in view of the presence of such defects in single crystals. The mutual orientation between a silicon matrix and coesite precipitates in Czochralski-silicon annealed specimens is proposed to be non unique as a consequence of the easy twinning of coesite.     

Data Assimilation (4D-VAR) for Shallow-Water Flow: The Case of the Chicoutimi River

This is a web presentation of the work presented at the 10th Annual Conference of the CFD Society of Canada, “CFD 2002”, at the University of Windsor on June 9-11, 2002. This discussion paper presents the four-dimensional variational data assimilation (4D-VAR) technique as a tool to forecast floods. This discussion will be limited to hydrological forecast. We assume that the weather, here a large rainstorm, had already been forecasted by the meteorological services. In the 4D-VAR technique, we need to minimize, in the sense of Lagrange, a cost function which measures the difference between the forecast and the observations. The physical equations acts as a set of constraints. Here, the model is the shallow-water equations modified to include sediment transport. The minimum was found by using the steepest descent algorithm. This is made possible because the gradient of the cost function can be calculated analytically by using the adjoint equations of the model. To illustrate the 4D-VAR technique, the bypass of a simple theoretical dam as well as the more complex overflowing of the Chicoutimi River at the Chute-Garneau dam (during the 1996 flood) are investigated.     

Differential distribution of the two cryptic species, Pomatoschistus microps and P. marmoratus, in the lagoons of southern France, with an emphasis on the genetic organisation of P. microps

The presence and the distribution of two cryptic sedentary gobies, Pomatoschistus microps and Pomatoschistus marmoratus, were investigated in several lagoons stretching along the Golfe du Lion, southern France Mediterranean coast, and Corsica. Pomatoschistus microps is the only sedentary Pomatoschistus species in two shallow lagoons with large variations of salinity and temperature. In contrast, P. marmoratus, another sedentary species, can be found preferentially in a deep salty lagoon with small salinity variations. Both species occur sympatrically in the Rhone Delta and hybridise, producing fertile offspring. Using allozymes, competition, linked to fitness, is invoked to explain the distribution pattern of these sibling species.     

A new method for sampling and analyzing volcanic sublimates — application to Merapi volcano, Java

A new method for the sampling of sublimates from high-temperature volcanic gases has been used at Merapi volcano, Java, in 1978. The sublimates were collected on the inner walls of silica tubes introduced into fumarolic vents. Volcanic gases were allowed to move freely through the tubes and as they cooled, a fraction of the volatile components condensed on the inner walls of the tubes along the temperature gradient. The sublimates were then analyzed by a combination of light microscopy, scanning electron microscopy, electron microprobe and X-ray diffraction.Six successive zones of different compositions and mineralogical associations have been identified along the covered range of temperatures (900° to around 400°C). From the high to the low temperatures, these zones are composed of: (1) cristobalite, magnetite, hercynite; (2) molybdenite; (3) acmite; (4) halite, sylvite; (5) sphalerite, pyrite; and (6) galena. Equilibrium calculations show that these crystalline phases are stable for pS₂, pC₁, and pO₂, values typical of magma-buffered gases that have not been contaminated by atmospheric oxygen.The deposits observed in the tubes may be useful in aiding the understanding of the mechanisms acting during the cooling of the gaseous phase on its way to the surface and before its emission into the atmosphere.     

On different types of adjustment usable to calculate the parameters of the stream power law

Model parameterization through adjustment to field data is a crucial step in the modeling and the understanding of the drainage network response to tectonic or climatic perturbations. Using as a test case a data set of 18 knickpoints that materialize the migration of a 0.7-Ma-old erosion wave in the Ourthe catchment of northern Ardennes (western Europe), we explore the impact of various data fitting on the calibration of the stream power model of river incision, from which a simple knickpoint celerity equation is derived. Our results show that statistical least squares adjustments (or misfit functions) based either on the stream-wise distances between observed and modeled knickpoint positions at time t or on differences between observed and modeled time at the actual knickpoint locations yield significantly different values for the m and K parameters of the model. As there is no physical reason to prefer one of these approaches, an intermediate least-rectangles adjustment might at first glance appear as the best compromise. However, the statistics of the analysis of 200 sets of synthetic knickpoints generated in the Ourthe catchment indicate that the time-based adjustment is the most capable of getting close to the true parameter values. Moreover, this fitting method leads in all cases to an m value lower than that obtained from the classical distance adjustment (for example, 0.75 against 0.86 for the real case of the Ourthe catchment), corresponding to an increase in the non-linear character of the dependence of knickpoint celerity on discharge.