Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

Application du réseau des neurones artificiels à la prévision des débits horaires: Cas du bassin versant de l’Eure,France

RÉSUMÉ

Les modèles pluie–débit sont fortement utilisés dans la gestion des risques hydrologiques et la prévision des crues. Dans cette étude, nous présentons un modèle pluie–débit pour la prévision des débits horaires basé sur la technique des réseaux de neurones artificiels (RNA). Ce modèle a été développé et appliqué sur le bassin versant de l’Eure au Nord-Ouest de la France afin de dépasser les problèmes dus à la non-linéarité de la relation pluie–débit et à l’imprécision des données collectées. La création de ce modèle a nécessité plusieurs étapes pendant lesquelles nous avons pu déterminer les paramètres du modèle permettant la compréhension de la complexité hydrologique et la production des informations nécessaires à la prévision. Elles ont abouti à un modèle de réseau de neurones artificiels capable d’effectuer, en quelques secondes, des prévisions des crues efficaces jusqu’à un horizon de prévision de 48 h. Ces résultats confirment que les modèles RNA peuvent jouer un rôle important dans le domaine de la prévision car capables de modéliser la non-linéarité des relations pluie–débit rencontrées sur certains bassins hydrologiques.

Editor Z.W. Kundzewicz Associate editor E. Gargouri     

Direct dating of thick‐ and thin‐skin thrusts in the Peruvian Subandean zone through apatite (U–Th)/He and fission track thermochronometry

Although the structure of the central Peruvian Subandean zone is well defined, the timing of thrust‐related exhumation and Cenozoic sedimentation remain poorly constrained. In this study, we report new apatite (U–Th)/He (AHe) and fission track (AFT) ages from thrust‐belt and foreland strata along three published balanced cross sections. AHe data from the northern, thick‐skinned domain (i.e. Shira Mountain, Otishi Cordillera and Ucayali Basin) show young AHe ages (ranging from 2.6 ± 0.2 to 13.1 ± 0.8 Ma) compared with AFT ages (ranging from 101 ± 5 to 133 ± 11 Ma). In the southern Camisea Basin, where deformation is mainly thin‐skinned, AHe and AFT ages have been both reset and show young cooling ages (3.7 ± 0.8 Ma and 8 ± 2 Ma respectively). Using low‐temperature thermochronology data and the latest fission track annealing and He diffusion codes, the thermal history of the study area has been reconstructed using inverse modelling. This history includes two steps of erosion: Early Cretaceous and late Neogene, but only Neogene sedimentation and exhumation varies in the different sectors of the study area. From a methodological point of view, large AHe data dispersion point to the need for refinement of AHe damage and annealing models. The influence of grain chemistry on damage annealing, multiple age components and the possibility of fission tracks as traps for He need further consideration. For the central Peruvian Subandes, AHe and AFT ages combined with balanced cross sections emphasize the dominant control of Paleozoic inheritance rather than climate on Cenozoic infilling and exhumation histories. Finally, our data provide the first field example of how thick‐skinned thrust‐related deformation and exhumation in the Subandes can be directly dated through AHe thermochronology.     

Hyperspectral remote sensing analysis of short rotation woody crops grown with controlled nutrient and irrigation treatments

Hyperspectral remote sensing research was conducted to document the biophysical and biochemical characteristics of controlled forest plots subjected to various nutrient and irrigation treatments. The experimental plots were located on the Savannah River Site near Aiken, SC. AISA hyperspectral imagery were analysed using three approaches, including: (1) normalized difference vegetation index based simple linear regression (NSLR), (2) partial least squares regression (PLSR) and (3) machine-learning regression trees (MLRT) to predict the biophysical and biochemical characteristics of the crops (leaf area index, stem biomass and five leaf nutrients concentrations). The calibration and cross-validation results were compared between the three techniques. The PLSR approach generally resulted in good predictive performance. The MLRT approach appeared to be a useful method to predict characteristics in a complex environment (i.e. many tree species and numerous fertilization and/or irrigation treatments) due to its powerful adaptability.     

Bertin’s matrix concepts reconsidered: transformations of semantics and semiotics to support geovisualization use

A system for managing and visualizing archaeological finds and their spatio-temporal properties is presented. The system is an attempt to combine archaeologists' interest to keep track of time with interests related to the representation and visualization of the spatio-temporal component of archaeological excavations, with the help of cartographic and GIS methods and techniques. Several parameters of a prehistoric excavation such as finds, their thematic characteristics, and their spatio-temporal distributions are examined. A three-dimensional GIS/cartographic environment and its user interface to the excavation database have been developed as a tool for enhancing archaeological interpretation in an exploratory environment. The development of the system is an ongoing task; the article describes its current use as a prototype.     

On the Importance of High-Resolution in Large-Scale Ocean Models

Eddying global ocean models are now routinely used for ocean prediction, and the value-added of a better representation of the observed ocean variability and western boundary currents at that resolution is currently being evaluated in climate models. This overview article begins with a brief summary of the impact on ocean model biases of resolving eddies in several global ocean–sea ice numerical simulations. Then, a series of North and Equatorial Atlantic configurations are used to show that an increase of the horizontal resolution from eddy-resolving to submesoscale-enabled together with the inclusion of high-resolution bathymetry and tides significantly improve the models' abilities to represent the observed ocean variability and western boundary currents. However, the computational cost of these simulations is extremely large, and for these simulations to become routine, close collaborations with computer scientists are essential to ensure that numerical codes can take full advantage of the latest computing architecture.     

Projections of severe droughts in future climate in Southeast Brazil: a case study in Southern Minas Gerais State,Brazil

Theoretical and Applied Climatology - South of Minas Gerais state, in Southeast Brazil, is known for the coffee crop production (more than 30% of country’s production) and hydroelectricity...     

Geology of the epithermal Ag–Au Huevos Verdes vein system and San José district, Deseado massif, Patagonia, Argentina

The San José district is located in the northwest part of the Deseado massif and hosts a number of epithermal Ag–Au quartz veins of intermediate sulfidation style, including the Huevos Verdes vein system. Veins are hosted by andesitic rocks of the Bajo Pobre Formation and locally by rhyodacitic pyroclastic rocks of the Chon Aike Formation. New ⁴⁰Ar/³⁹Ar constraints on the age of host rocks and mineralization define Late Jurassic ages of 151.3 ± 0.7 Ma to 144.7 ± 0.1 Ma for volcanic rocks of the Bajo Pobre Formation and of 147.6 ± 1.1 Ma for the Chon Aike Formation. Illite ages of the Huevos Verdes vein system of 140.8 ± 0.2 and 140.5 ± 0.3 Ma are 4 m.y. younger than the volcanic host rock unit. These age dates are among the youngest reported for Jurassic volcanism in the Deseado massif and correlate well with the regional context of magmatic and hydrothermal activity. The Huevos Verdes vein system has a strike length of 2,000 m, with several ore shoots along strike. The vein consists of a pre-ore stage and three main ore stages. Early barren quartz and chalcedony are followed by a mottled quartz stage of coarse saccharoidal quartz with irregular streaks and discontinuous bands of sulfide-rich material. The banded quartz–sulfide stage consists of sulfide-rich bands alternating with bands of quartz and bands of chlorite ± illite. Late-stage sulfide-rich veinlets are associated with kaolinite gangue. Ore minerals are argentite and electrum, together with pyrite, sphalerite, galena, chalcopyrite, minor bornite, covellite, and ruby silver. Wall rock alteration is characterized by narrow (< 3 m) halos of illite and illite/smectite next to veins, grading outward into propylitic alteration. Gangue minerals are dominantly massive quartz intergrown with minor to accessory adularia. Epidote, illite, illite/smectite, and, preferentially at deeper levels, Fe-chlorite gangue indicate near-neutral pH hydrothermal fluids at temperatures of >220°C. Kaolinite occurring with the late sulfide-rich veinlet stage indicates pH < 4 and a temperature of <200°C. The Huevos Verdes system has an overall strike of 325°, dipping on average 65° NE. The orientations of individual ore shoots are controlled by vein strike and intersecting north-northwest-striking faults. We propose a structural model for the time of mineralization of the San José district, consisting of a conjugate shear pair of sinistral north-northwest- and dextral west-northwest-striking faults that correspond to R and R′ in the Riedel shear model and that are related to master faults (M) of north-northeast-strike. Veins of 315° strike can be interpreted as nearly pure extensional fractures (T). Variations in vein strike predict an induced sinistral shear component for strike directions of >315°, whereas strike directions of <315° are predicted with an induced dextral strike–slip movement. The components of the structural model appear to be present on a regional scale and are not restricted to the San José district.