MODFLOW as a Configurable Multi-Model Hydrologic Simulator     

Christian D. Langevin  Joseph D. Hughes  Alden M. Provost  Martijn J. Russcher  Sorab Panday 《Ground water》2024,62(1):111-123

MODFLOW 6 is the latest in a line of six “core” versions of MODFLOW released by the U.S. Geological Survey. The MODFLOW 6 architecture supports incorporation of additional hydrologic processes, in addition to groundwater flow, and allows interaction between processes. The architecture supports multiple model instances and multiple types of models within a single simulation, a flexible approach to formulating and solving the equations that represent hydrologic processes, and recent advances in interoperability, which allow MODFLOW to be accessed and controlled by external programs. The present version of MODFLOW 6 consolidates popular capabilities available in MODFLOW variants, such as the unstructured grid support in MODFLOW-USG, the Newton-Raphson formulation in MODFLOW-NWT, and the support for partitioned stress boundaries in MODFLOW-CDSS. The flexible multi-model capability allows users to configure MODFLOW 6 simulations to represent the local-grid refinement (LGR) capabilities available in MODFLOW-LGR, the multi-species transport capabilities in MT3DMS, and the coupled variable-density capabilities available in SEAWAT. This paper provides a new, holistic and integrated overview of simulation capabilities made possible by the MODFLOW 6 architecture, and describes how ongoing and future development can take advantage of the program architecture to integrate new capabilities in a way that is minimally invasive and automatically compatible with the existing MODFLOW 6 code.  相似文献   

Dense collection areas and terrestrial alteration of meteorites in the Atacama Desert     

Gabriel A. Pinto  Adrien Tavernier  Jérôme Gattacceca  Alexandre Corgne  Millarca Valenzuela  Béatrice Luais  Laura Flores  Felipe Olivares  Yves Marrocchi 《Meteoritics & planetary science》2024,59(2):351-367

In the last 15 years, more than 2700 meteorites have been recovered and officially classified from the Atacama Desert. Although the number of meteorites collected in the Atacama has risen, the physical and climatic properties of the dense collection areas (DCAs) have not been fully characterized. In this article, we compiled the published data of all classified meteorites found in the Atacama Desert to (i) describe the distribution by meteorite groups, (ii) compare the weathering degree of chondrites among different Atacama DCAs and other hot and cold deserts, and (iii) determine the preservation conditions of chondrites in the main Atacama DCAs in relation with the local climatic conditions. The 35 DCAs so far identified in the Atacama Desert are located in three main morphotectonic units: The Coastal Range (CR), Central Depression (CD), and Pre-Andean Range/Basement. A comparison with reported weathering data from other cold and hot deserts indicates that the mean terrestrial weathering of Atacama chondrites (W1–2), displays less alteration than other hot deserts (W2–3) and resembles the weathering distribution of the Antarctic meteorites (W1–2). The highest abundance of Atacama chondrites with low weathering (≤W2) is localized in the CD (78.8%, N = 1435), which is protected from the coastal fog influence and seasonal rainfalls and displays the oldest surfaces in the Atacama Desert. The morphogenetic classification based on present-day temperatures and precipitations of the main Atacama DCAs reveals similar regional/subregional climatic conditions in the most productive areas and a truly productive surface for meteorite recovery between 5% and 58% of the quadrangles formally defined for each Atacama DCA. Our morphogenetic classification lacks consideration of some meteorological parameters such as the coastal fog, so it cannot fully explain the differences in weathering patterns among CR chondrites. Future studies of chondrite preservation in the Atacama DCAs should consider other meteorological variables such as relative humidity, specific humidity, or dew point, in combination with exposure ages of meteorites and its surfaces.  相似文献   

Hydraulic-Head Formulation for Density-Dependent Flow and Transport     

Christian D. Langevin  Sorab Panday  Alden M. Provost 《Ground water》2020,58(3):349-362

Density-dependent flow and transport solutions for coastal saltwater intrusion investigations, analyses of fluid injection into deep brines, and studies of convective fingering and instabilities of denser fluids moving through less dense fluids typically formulate the groundwater flow equation in terms of pressure or equivalent freshwater head. A formulation of the flow equation in terms of hydraulic head is presented here as an alternative. The hydraulic-head formulation can facilitate adaptation of existing constant-density groundwater flow codes to include density-driven flow by avoiding the need to convert between freshwater head and hydraulic head within the code and by incorporating density-dependent terms as a compartmentalized “correction” to constant-density calculations already performed by the code. The hydraulic-head formulation also accommodates complexities such as unconfined groundwater flow and Newton-Raphson solution schemes more readily than the freshwater-head formulation. Simulation results are presented for four example problems solved using an implementation of the hydraulic-head formulation in MODFLOW.  相似文献   

A Review on Forearc Ophiolite Obduction,Adakite-Like Generation,and Slab Window Development at the Chile Triple Junction Area: Uniformitarian Framework for Spreading-Ridge Subduction     

Jacques Bourgois  Yves Lagabrielle  Hervé Martin  Jérôme Dyment  Jose Frutos  Maria Eugenia Cisternas 《Pure and Applied Geophysics》2016,173(10-11):3217-3246

  相似文献   

Hydrogen emissions from Erebus volcano, Antarctica     

Yves Moussallam  Clive Oppenheimer  Alessandro Aiuppa  Gaetano Giudice  Manuel Moussallam  Philip Kyle 《Bulletin of Volcanology》2012,74(9):2109-2120

The continuous measurement of molecular hydrogen (H₂) emissions from passively degassing volcanoes has recently been made possible using a new generation of low-cost electrochemical sensors. We have used such sensors to measure H₂, along with SO₂, H₂O and CO₂, in the gas and aerosol plume emitted from the phonolite lava lake at Erebus volcano, Antarctica. The measurements were made at the crater rim between December 2010 and January 2011. Combined with measurements of the long-term SO₂ emission rate for Erebus, they indicate a characteristic H₂ flux of 0.03?kg s^–1 (2.8?Mg? day^–1). The observed H₂ content in the plume is consistent with previous estimates of redox conditions in the lava lake inferred from mineral compositions and the observed CO₂/CO ratio in the gas plume (～0.9 log units below the quartz–fayalite–magnetite buffer). These measurements suggest that H₂ does not combust at the surface of the lake, and that H₂ is kinetically inert in the gas/aerosol plume, retaining the signature of the high-temperature chemical equilibrium reached in the lava lake. We also observe a cyclical variation in the H₂/SO₂ ratio with a period of ～10?min. These cycles correspond to oscillatory patterns of surface motion of the lava lake that have been interpreted as signs of a pulsatory magma supply at the top of the magmatic conduit.  相似文献   

Solute and Heat Transport Model of the Henry and Hilleke Laboratory Experiment     

Christian D. Langevin  Alyssa M. Dausman  Michael C. Sukop 《Ground water》2010,48(5):757-770

SEAWAT is a coupled version of MODFLOW and MT3DMS designed to simulate variable-density ground water flow and solute transport. The most recent version of SEAWAT, called SEAWAT Version 4, includes new capabilities to represent simultaneous multispecies solute and heat transport. To test the new features in SEAWAT, the laboratory experiment of Henry and Hilleke (1972) was simulated. Henry and Hilleke used warm fresh water to recharge a large sand-filled glass tank. A cold salt water boundary was represented on one side. Adjustable heating pads were used to heat the bottom and left sides of the tank. In the laboratory experiment, Henry and Hilleke observed both salt water and fresh water flow systems separated by a narrow transition zone. After minor tuning of several input parameters with a parameter estimation program, results from the SEAWAT simulation show good agreement with the experiment. SEAWAT results suggest that heat loss to the room was more than expected by Henry and Hilleke, and that multiple thermal convection cells are the likely cause of the widened transition zone near the hot end of the tank. Other computer programs with similar capabilities may benefit from benchmark testing with the Henry and Hilleke laboratory experiment.  相似文献   

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals     

Tristan Allouis  Jean‐Stéphane Bailly  Yves Pastol  Catherine Le Roux 《地球表面变化过程与地形》2010,35(6):640-650

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

Statement to comment by K. Föllmi on ,,Neubenennungen lithostratigraphischer Einheiten in der helvetischen Kreide“     

Schweizerisches Komitee für Stratigraphie  Thierry Adatte  Jean-Pierre Berger  Hansruedi Bläsi  Reto Burkhalter  Stephan DallAgnolo  Danielle Decrouez  Hanspeter Funk  Heinz Furrer  Yves Gouffon  Hansruedi Graf  Werner Heckendorn  Wilfried Winkler 《Swiss Journal of Geoscience》2009,102(3):519-519

  相似文献   

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

Eric D. Morway  Richard G. Niswonger  Christian D. Langevin  Ryan T. Bailey  Richard W. Healy 《Ground water》2013,51(2):237-251

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.  相似文献   

Dynamique et modélisation des crues dans le bassin du Mono à Nangbéto (Togo/Bénin)     

Ernest Amoussou  Yves Tramblay  Henri S.V. Totin  Gil Mahé  Pierre Camberlin 《水文科学杂志》2013,58(11):2060-2071

Résumé

Les inondations des dernières années en Afrique de l’Ouest rendent nécessaire la prévision hydrologique pour réduire la vulnérabilité des populations. Nous avons étudié la dynamique des crues dans le bassin du Mono en amont du barrage de Nangbéto (Togo et Bénin) par des méthodes statistiques et déterministes. L’augmentation des précipitations au cours de la décennie 1990 et la mise en eau du barrage de Nangbéto en 1988 ont modifié le régime du Mono, dont les débits d’étiages ont augmenté. Les débits entrants et sortants du barrage sont équivalents, même lors des crues. Les maximums annuels des débits entrants peuvent être modélisés par une distribution des valeurs extrêmes (GEV) de type Fréchet, tandis que ceux des débits sortants, suivent une loi GEV de type Weibull (bornée). Le mécanisme dominant de génération des crues dans le bassin est la montée progressive des eaux associée à l’augmentation de la saturation des sols. L’analyse des débits de crue entre 1988 et 2010 n’a pas montré d’augmentation de la magnitude ou de la fréquence des crues. Le modèle GR4J reproduit bien les débits de crues en période calage, mais ces critères de performance se dégradent lorsque l’on s’en éloigne. Dans ces conditions, l’utilisation de ce modèle pour évaluer les risques futurs de crues serait hasardeuse et il est nécessaire d’envisager d’autres types d’approches.  相似文献