Hydrologic Modeling of a Groundwater Dominated Watershed Using a Loosely Coupled Modeling Approach

Total Page:16

File Type:pdf, Size:1020Kb

Hydrologic Modeling of a Groundwater Dominated Watershed Using a Loosely Coupled Modeling Approach Hydrologic Modeling of a Groundwater Dominated Watershed Using a Loosely Coupled Modeling Approach 1Barth, C. , 1P. Krause, 2D.P. Boyle and 3S.L. Markstrom 1Friedrich Schiller University Jena, Department of Geoinformatics, Hydrology,and Modeling, Loebdergraben 32, 07743 Jena, Germany, 2 Desert Research Institute, 2215 Raggio Parkway, Reno, NV 89512, USA, . 3U.S. Geological Survey, Box 25046, Mail Stop 412, Denver Federal Center, Denver, CO 80225, USA E-Mail: [email protected] Keywords: surface-groundwater interaction; loosely coupled model; MODFLOW; PRMS. basin features a strong surface-groundwater EXTENDED ABSTRACT interaction of unknown degree. Since the late 1960s, hydrologic modeling focused This research paper describes the concept of an more and more on interactions between surface inexpensive and simple method of combining the and groundwater flow systems, leading to the advantages of two established hydrologic development of integrated modeling approaches modeling approaches by using output information simulating flow in these two systems and the of the surface water model as input for the simulation of interactions between them. These groundwater model. The basic assumption is made interactions are of variable character as well as that groundwater recharge rates provided by spatial and temporal heterogeneity. PRMS are reasonable because the estimation procedure benefits from the fairly high resolution The primary objective of this study is the of meteorological, land use and soil parameters of application and investigation of a loosely coupled the surface water model. At present, the approach modeling approach, combining two well-known presented in this paper is limited to a models of the United States Geological Survey, the unidirectional coupling of the surface and Precipitation Runoff Modeling System (PRMS) groundwater model, which is only accounting for and the Modular Finite-Difference Groundwater groundwater estimates provided by PRMS as input Flow Model (MODFLOW), in order to simulate for MODFLOW; whereas possible upward complex hydrological processes in a mesoscale movement of water, i.e. groundwater seepage to watershed. Within the scope of the hydrologic the surface is neglected. Performance and modeling groundwater recharge, the parameter efficiency of this approach are assessed in the unidirectionally linking the surface water and scope of applying the model in steady-state as well groundwater system is used to couple these two as transient mode to the Esperstedter Ried basin. models loosely. Considering factors that have a great influence on groundwater recharge (e.g. The results of this case study, which produced a climate, topography, soil, and land use), the number of reasonable as well as some unsatisfying potential of modeling the hydrologic system with a simulations, demonstrate the potential of coupling distributed surface water model is explored. In a surface water model and a groundwater model to addition, the option of constraining the PRMS obtain more complex and accurate analyses and parameter determining groundwater recharge, and simulations of hydrologic systems. However, the therefore the loosely coupled surface-groundwater study also shows the need for further research in model, with observed head levels only, is regard to constraining the parameter linking both investigated. models, i.e. groundwater recharge. Furthermore, the necessity for a more dynamic and complex This modeling approach is applied to the method of integrating surface and subsurface Esperstedter Ried basin, an ungauged, mesoscale, interactions is revealed, which accounts for groundwater-dominated catchment in central downward as well as upward movement of water Germany. The study represents a good example of in the subsurface system. problems frequently occurring in hydrologic modeling. Due to the lack of observed streamflow and surface runoff data, the calibration and validation of any hydrologic model for this watershed proves to be difficult. Furthermore, the 601 Model (MODFLOW), which are both available 1. INTRODUCTION free of charge on the internet. The factor unidirectionally linking both systems (i.e. the Traditional hydrologic modeling approaches parameter coupling the surface and the generally concentrate on either the surface or the groundwater model), namely groundwater subsurface system. Surface water models are recharge, is investigated regarding the performance designed and calibrated with emphasis on surface of this loosely coupled approach. Furthermore, the water movement, evapotranspiration and possibility of constraining the PRMS parameter infiltration; whereas the subsurface is commonly determining the amount of recharge into the viewed as a simple, lumped system and subsurface reservoir by merely using observed groundwater processes are not physically based. groundwater head levels is assessed. Groundwater models, on the other hand, focus on obtaining exact groundwater head levels in respect 2. MODEL DESCRIPTION to occurring stresses while surface water processes are oversimplified. 2.1. The groundwater model MODFLOW For the last three decades, a major objective of The Modular Three-Dimensional Finite-Difference hydrologists has been the coupling of both Groundwater Flow Model (MODFLOW) first approaches in order to model the entire hydrologic released by the U.S. Geological Survey in 1984, system and, furthermore, to investigate the has become the most widely used groundwater interaction between surface and subsurface model (McDonald & Harbaugh 1984) processes. Particularly basic methods of modeling MODFLOW is a mathematical model simulating surface-groundwater interactions based on the groundwater flow through heterogeneous, numerical solutions of differential equations and a porous media. In MODFOW three-dimensional variety of integrated models have been developed groundwater flow is described by the partial- since the late 1960s. Some of the first studies of differential equation surface-subsurface interactions were conducted by ∂ ∂ ∂ ∂ ∂ ∂ ∂ Smith & Woolhiser (1971), who coupled a one- ⎛ h ⎞ + ⎛ h ⎞ + ⎛ h ⎞ − = h (1) ⎜ K xx ⎟ ⎜ K yy ⎟ ⎜ K zz ⎟ W SS dimensional flow vertical infiltration model to an ∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠ ∂z ⎝ ∂z ⎠ ∂t overland flow model, Pinder & Sauer (1971) linked a one-dimensional streamflow model to a where Kxx, Kyy, and Kzz = the hydraulic conductivity two-dimensional aquifer, and Freeze (1972) along the x, y, and z axes that are assumed to be designed a model simulating three-dimensional parallel to the major axes of hydraulic conductivity transient subsurface flow as well as one- (L/T), h = the potentiometric head (L), dimensional streamflow. In recent years, a number W = volumetric flux per unit value representing of integrated models has been developed, sources and/or sinks of water ( W < 0.0 for outflow including SHE, MIKE-SHE, and FHM. SHE of the groundwater system, W > 0.0 for inflow (T-1), SS = specific storage of porous material (L- (SYSTÈME HYDROLOGIQUE EUROPÈEN) is a physically based modeling system, simulating 1), and T = time (McDonald & Harbaugh 1984). various components of water movement (such as evapotranspiration, overland and channel flow, Spatial discretization in MODFLOW is achieved flow in the unsaturated and saturated zones, and by considering real-world aquifers a system of grid many more) based on finite-difference or empirical cells which are characterized in terms of rows (i), equations. In addition to representing the entire columns (j), and layers (k). Every cell represents a land phase of the hydrologic cycle (interception, unit of homogenous properties. In the current evapotranspiration, surface runoff, river routing version MODFLOW-2000, the finite-difference etc.), MIKE-SHE also describes water movement grid is assumed to be of rectangular shape in the unsaturated soil zone and groundwater flow. horizontally, while the vertical dimension is The coupled surface and groundwater model FHM, distorted. Temporal discretization is based on time which incorporates components of the surface steps, which are grouped into stress periods. The water model HSPF and the groundwater model length of particular time steps is user-defined MODFLOW, is another example for an integrated during the model setup. model. In this study MODFLOW-2000 is used, which is This paper describes an alternative approach to operated in the comprehensive graphical user modeling a surface-groundwater system by loosely environment GMS (Department of Defense coupling two well-established models of the Groundwater Modeling System). United States Geological Survey, the Precipitation Runoff Modeling System (PRMS) and the Modular Finite-Difference Groundwater Flow 602 2.2. The surface water model PRMS recharge rates for each grid cell, the proportion of area of each HRU per cell is computed; hence, the The U.S. Geological Survey’s Precipitation Runoff amount of recharge can be determined for every Modeling System (PRMS), a deterministic, active cell in the MODFLOW grid model based on physical-process based watershed model was this proportion. The conversion of temporal scales developed by Leavesley et al. in 1983 (Leavesley needed for the transient MODFLOW simulation is et al. 1983). PRMS which operates either as performed by calculating the mean monthly lumped- or distributed parameter model, is amount of recharge for each HRU based on daily designed to analyze the dynamics and responses of output of ssr2gwrate values. The monthly recharge a hydrologic system regarding the
Recommended publications
  • A GIS-Based Software to Simulate Groundwater Nitrate Load from Septic Systems to Surface Water Bodies
    Computers & Geosciences 52 (2013) 108–116 Contents lists available at SciVerse ScienceDirect Computers & Geosciences journal homepage: www.elsevier.com/locate/cageo ArcNLET: A GIS-based software to simulate groundwater nitrate load from septic systems to surface water bodies J. Fernando Rios a, Ming Ye b,n, Liying Wang b, Paul Z. Lee c, Hal Davis d, Rick Hicks c a Department of Geography, State University of New York at Buffalo, Buffalo, NY, USA b Department of Scientific Computing, Florida State University, Tallahassee, FL, USA c Florida Department of Environmental Protection, Tallahassee, FL, USA d 2625 Vergie Court, Tallahassee, FL 32303, USA article info abstract Article history: Onsite wastewater treatment systems (OWTS), or septic systems, can be a significant source of nitrates Received 17 June 2012 in groundwater and surface water. The adverse effects that nitrates have on human and environmental Received in revised form health have given rise to the need to estimate the actual or potential level of nitrate contamination. 4 October 2012 With the goal of reducing data collection and preparation costs, and decreasing the time required to Accepted 5 October 2012 produce an estimate compared to complex nitrate modeling tools, we developed the ArcGIS-based Available online 16 October 2012 Nitrate Load Estimation Toolkit (ArcNLET) software. Leveraging the power of geographic information Keywords: systems (GIS), ArcNLET is an easy-to-use software capable of simulating nitrate transport in ground- Nitrate transport water and estimating long-term nitrate loads from groundwater to surface water bodies. Data Screening model requirements are reduced by using simplified models of groundwater flow and nitrate transport which Denitrification consider nitrate attenuation mechanisms (subsurface dispersion and denitrification) as well as spatial Septic system Nitrogen variability in the hydraulic parameters and septic tank distribution.
    [Show full text]
  • Climate Models and Their Evaluation
    8 Climate Models and Their Evaluation Coordinating Lead Authors: David A. Randall (USA), Richard A. Wood (UK) Lead Authors: Sandrine Bony (France), Robert Colman (Australia), Thierry Fichefet (Belgium), John Fyfe (Canada), Vladimir Kattsov (Russian Federation), Andrew Pitman (Australia), Jagadish Shukla (USA), Jayaraman Srinivasan (India), Ronald J. Stouffer (USA), Akimasa Sumi (Japan), Karl E. Taylor (USA) Contributing Authors: K. AchutaRao (USA), R. Allan (UK), A. Berger (Belgium), H. Blatter (Switzerland), C. Bonfi ls (USA, France), A. Boone (France, USA), C. Bretherton (USA), A. Broccoli (USA), V. Brovkin (Germany, Russian Federation), W. Cai (Australia), M. Claussen (Germany), P. Dirmeyer (USA), C. Doutriaux (USA, France), H. Drange (Norway), J.-L. Dufresne (France), S. Emori (Japan), P. Forster (UK), A. Frei (USA), A. Ganopolski (Germany), P. Gent (USA), P. Gleckler (USA), H. Goosse (Belgium), R. Graham (UK), J.M. Gregory (UK), R. Gudgel (USA), A. Hall (USA), S. Hallegatte (USA, France), H. Hasumi (Japan), A. Henderson-Sellers (Switzerland), H. Hendon (Australia), K. Hodges (UK), M. Holland (USA), A.A.M. Holtslag (Netherlands), E. Hunke (USA), P. Huybrechts (Belgium), W. Ingram (UK), F. Joos (Switzerland), B. Kirtman (USA), S. Klein (USA), R. Koster (USA), P. Kushner (Canada), J. Lanzante (USA), M. Latif (Germany), N.-C. Lau (USA), M. Meinshausen (Germany), A. Monahan (Canada), J.M. Murphy (UK), T. Osborn (UK), T. Pavlova (Russian Federationi), V. Petoukhov (Germany), T. Phillips (USA), S. Power (Australia), S. Rahmstorf (Germany), S.C.B. Raper (UK), H. Renssen (Netherlands), D. Rind (USA), M. Roberts (UK), A. Rosati (USA), C. Schär (Switzerland), A. Schmittner (USA, Germany), J. Scinocca (Canada), D. Seidov (USA), A.G.
    [Show full text]
  • Vibrational Spectra of Light and Heavy Water with Application to Neutron Cross Section Calculations J
    CORE Metadata, citation and similar papers at core.ac.uk Provided by CONICET Digital Vibrational spectra of light and heavy water with application to neutron cross section calculations J. I. Marquez Damian, D. C. Malaspina, and J. R. Granada Citation: J. Chem. Phys. 139, 024504 (2013); doi: 10.1063/1.4812828 View online: http://dx.doi.org/10.1063/1.4812828 View Table of Contents: http://jcp.aip.org/resource/1/JCPSA6/v139/i2 Published by the AIP Publishing LLC. Additional information on J. Chem. Phys. Journal Homepage: http://jcp.aip.org/ Journal Information: http://jcp.aip.org/about/about_the_journal Top downloads: http://jcp.aip.org/features/most_downloaded Information for Authors: http://jcp.aip.org/authors Downloaded 11 Jul 2013 to 200.0.233.52. This article is copyrighted as indicated in the abstract. Reuse of AIP content is subject to the terms at: http://jcp.aip.org/about/rights_and_permissions THE JOURNAL OF CHEMICAL PHYSICS 139, 024504 (2013) Vibrational spectra of light and heavy water with application to neutron cross section calculations J. I. Marquez Damian,1,a) D. C. Malaspina,2 and J. R. Granada1,b) 1Neutron Physics Department and Instituto Balseiro, Centro Atómico Bariloche, CNEA, Argentina 2Department of Biomedical Engineering and Chemistry of Life Processes Institute, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, USA (Received 4 June 2013; accepted 19 June 2013; published online 11 July 2013) The design of nuclear reactors and neutron moderators require a good representation of the interac- tion of low energy (E < 1 eV) neutrons with hydrogen and deuterium containing materials. These models are based on the dynamics of the material, represented by its vibrational spectrum.
    [Show full text]
  • The Incredible Lightness of Water Vapor
    1 The Incredible Lightness of Water Vapor ∗ 2 Da Yang and Seth Seidel 3 University of California, Davis 4 Lawrence Berkeley National Laboratory, Berkeley ∗ 5 Corresponding author address: Da Yang, 253 Hoagland Hall, Davis, CA 95616. 6 E-mail: [email protected] Generated using v4.3.2 of the AMS LATEX template 1 ABSTRACT 7 The molar mass of water vapor is significantly less than that of dry air. This 8 makes a moist parcel lighter than a dry parcel of the same temperature and 9 pressure. This effect is referred to as the vapor buoyancy effect and has of- 10 ten been overlooked in climate studies. We propose that this effect increases 11 Earth’s outgoing longwave radiation (OLR) and stabilizes Earth’s climate. 12 We illustrate this mechanism in an idealized tropical atmosphere, where there 13 is no horizontal buoyancy gradient in the free troposphere. To maintain the 14 uniform buoyancy distribution, temperature increases toward dry atmosphere 15 columns to compensate reduction of vapor buoyancy. The temperature differ- 16 ence between moist and dry columns would increase with climate warming 17 due to increasing atmospheric water vapor, leading to enhanced OLR and 18 thereby stabilizing Earth’s climate. We estimate that this feedback strength 2 19 is about O(0.2 W/m /K), which compares with cloud feedbacks and surface 20 albedo feedbacks in current climate. 2 21 1. Introduction 22 How fast would Earth’s climate respond to increasing CO2 (Manabe and Wetherald 1975; Flato 23 et al. 2013; Collins et al. 2013)? Why is tropical climate more stable than extratropical climate 24 (Holland and Bitz 2003; Polyakov et al.
    [Show full text]
  • Computational Investigation of Surface Freezing in a Molecular Model Of
    Computational investigation of surface freezing in a molecular model of water Amir Haji-Akbari ( )a,1 and Pablo G. Debenedettia,2 aDepartment of Chemical and Biological Engineering, Princeton University, Princeton, NJ 08544 Contributed by Pablo G. Debenedetti, February 17, 2017 (sent for review December 21, 2016; reviewed by Christoph Dellago and Angelos Michaelides) Water freezes in a wide variety of low-temperature environ- is therefore one of the major open challenges of atmospheric ments, from meteors and atmospheric clouds to soil and bio- sciences. logical cells. In nature, ice usually nucleates at or near inter- One major difficulty in constructing such predictive frame- faces, because homogenous nucleation in the bulk can only be works is the inability to determine the proper scaling of observed at deep supercoolings. Although the effect of proximal volumetric nucleation rates with the size of microdroplets and surfaces on freezing has been extensively studied, major gaps in nanodroplets that constitute clouds. This is practically impor- understanding remain regarding freezing near vapor–liquid inter- tant because atmospheric clouds are composed of droplets of faces, with earlier experimental studies being mostly inconclu- different sizes. This scaling is determined by whether freezing is sive. The question of how a vapor–liquid interface affects freez- enhanced or suppressed near vapor–liquid interfaces. If freezing ing in its vicinity is therefore still a major open question in ice is suppressed at a free interface, the effective volumetric nucle- physics. Here, we address this question computationally by using ation rate will be independent of r, the droplet radius, and bulk the forward-flux sampling algorithm to compute the nucleation freezing will be dominant at all sizes.
    [Show full text]
  • Case Study: Water and Ice
    Case Study: Water and Ice Timothy A. Isgro, Marcos Sotomayor, and Eduardo Cruz-Chu 1 The Universal Solvent Water is essential for sustaining life on Earth. Almost 75% of the Earth’s surface is covered by it. It composes roughly 70% of the human body by mass [1]. It is the medium associated with nearly all microscopic life pro- cesses. Much of the reason that water can sustain life is due to its unique properties. Among the most essential and extreme properties of water is its capabil- ity to absorb large amounts of heat. The heat capacity of water, which is the highest for compounds of its type in the liquid state, measures the amount of heat which needs to be added to water to change its temperature by a given amount. Water, thus, is able to effectively maintain its temperature even 1 when disturbed by great amounts of heat. This property serves to main- tain ocean temperatures, as well as the atmospheric temperature around the oceans. For example, when the sun rises in the morning and a large amount of heat strikes the surface of the Earth, the vast majority of it is absorbed by the ocean. The ocean water, however, does not exhibit a drastic increase in temperature that could make it inhospitable to life. In constrast, when the sun sets in the evening and that heat is taken away, the oceans do not become too cold to harbor life. Water also has high latent heats of vapor- ization and melting, which measure the amount of heat needed to change a certain amount of liquid water to vapor and ice to liquid, respectively.
    [Show full text]
  • Hydrated Sulfate Clusters SO4^2
    Article Cite This: J. Phys. Chem. B 2019, 123, 4065−4069 pubs.acs.org/JPCB 2− n − Hydrated Sulfate Clusters SO4 (H2O)n ( =1 40): Charge Distribution Through Solvation Shells and Stabilization Maksim Kulichenko,† Nikita Fedik,† Konstantin V. Bozhenko,‡,§ and Alexander I. Boldyrev*,† † Department of Chemistry and Biochemistry, Utah State University, 0300 Old Main Hill, Logan, Utah 84322-0300, United States ‡ Department of Physical and Colloid Chemistry, Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St, Moscow 117198, Russian Federation § Institute of Problems of Chemical Physics, Russian Academy of Sciences, Chernogolovka 142432, Moscow Region, Russian Federation *S Supporting Information 2− ABSTRACT: Investigations of inorganic anion SO4 interactions with water are crucial 2− for understanding the chemistry of its aqueous solutions. It is known that the isolated SO4 dianion is unstable, and three H2O molecules are required for its stabilization. In the 2− current work, we report our computational study of hydrated sulfate clusters SO4 (H2O)n (n =1−40) in order to understand the nature of stabilization of this important anion by fi 2− water molecules. We showed that the most signi cant charge transfer from dianion SO4 ≤ 2− to H2O takes place at a number of H2O molecules n 7. The SO4 directly donates its charge only to the first solvation shell and surprisingly, a small amount of electron density of 0.15|e| is enough to be transferred in order to stabilize the dianion. Upon further addition of ff ≤ H2O molecules, we found that the cage e ect played an essential role at n 12, where the fi 2− | | rst solvation shell closes.
    [Show full text]
  • Use of Models for Water Resources Management, Planning, and Policy
    Use of Models for Water Resources Management, Planning, and Policy August 1982 NTIS order #PB83-103655 Library of Congress Catalog Card Number 82-600556 For sale by the Superintendent of Documents, U.S. Government Printing Office, Washington, D.C. 20402 Foreword. The Nation'. water resource policies, affect many problem. in the ~$tates today-food production, energy, regional economic development, enviro:,"qual­ ity, and even our international balance of trade. As the country grows, and ~. ,. _"'I water supplies diminish, it becomes increasingly important to manage existing sup~th the greatest possible efficiency. In re~ent years, successfull11anagement andplanniag of water resources has increasingly been based on the results ,of mathematical~~.: ~~; Leaving aside the mystique of computers and ~pl~_ mathe~atics, mltematical models are simply tools used for understaIJ.ding w~~f~8odIc. and ",at~ re_e man­ agement activities. This part of water resource manaFm.~~ii thou~ not as apparent as dams and reservoirs or pipes and sewers, is a vital cOlripcmetif inIDetltiJ;lg tlriW'l ation' s water resource needs. Sophisticated analysis, through the ute of models, can improve our understanding of water resources and wat~r resource activities, and helpprev~ wasting both water and money. - This assessment of water resource models is therefore 'DOt an assessment Of_'$tical equations or computers, but of the Nation's ability t~ Use models to more e.·' tly and effectively analyze and solve our water resource prob~,m.. :, I •. ~rn-.' .. :", ',.,e ~ a:'~e, ssment :',.,' •.,:. rs not only the usefulness of the technology-the models-but _ability of Feder -,.,.' State water resource agencies to apply these analytic too1a ~~ly ....
    [Show full text]
  • Ice–Water Interfacial Free Energy for the TIP4P, TIP4P/2005, TIP4P/Ice
    Article pubs.acs.org/JPCC Ice−Water Interfacial Free Energy for the TIP4P, TIP4P/2005, TIP4P/ Ice, and mW Models As Obtained from the Mold Integration Technique Jorge R. Espinosa, Carlos Vega, and Eduardo Sanz* Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain ABSTRACT: The freezing of water is greatly influenced by the ice−water interfacial free energy. Yet, no consistent experimental measures of this thermodynamic parameter can be found. In this work we provide estimates for the ice Ih− water interfacial free energy at the normal melting temperature for different crystal planes (basal, primary prismatic, and secondary prismatic) using some widely used water models: TIP4P, TIP4P/2005, TIP4P/Ice, and mW. To compute the interfacial free energy, we use the mold integration method. It consists in calculating the work needed to induce the formation of a crystal slab in the fluid at coexistence conditions with the aid of a mold of potential energy wells whose structure is that of the crystal plane under study. The basal plane has the lowest interfacial free energy in all models of the TIP4P family. For the mW model we could not resolve differences in interfacial free energy between different orientations. The interfacial free energies averaged over all crystal orientations we obtain are 27.2(8), 28.9(8), 29.8(8), and 34.9(8) mJ/m2 for the TIP4P, TIP4P/2005, TIP4P/Ice, and mW models, respectively. The averaged interfacial free energy increases with both the melting temperature and melting enthalpy of the model. Moreover, we compute the interfacial free energy for several crystal orientation of ice Ic using the TIP4P/Ice model and obtain within the accuracy of our calculations the same orientationally averaged interfacial free energy as that of ice Ih.
    [Show full text]
  • Integrated Surface-Subsurface Water Flow Modelling of the Laxemar Area Application of the Hydrological Model ECOFLOW
    R-07-07 Integrated surface-subsurface water flow modelling of the Laxemar area Application of the hydrological model ECOFLOW Nikolay Sokrut, Kent Werner, Johan Holmén Golder Associates AB January 2007 Svensk Kärnbränslehantering AB Swedish Nuclear Fuel and Waste Management Co Box 5864 SE-102 40 Stockholm Sweden Tel 08-459 84 00 +46 8 459 84 00 Fax 08-661 57 19 +46 8 661 57 19 CM Gruppen AB, Bromma, 2007 ISSN 1402-3091 Tänd ett lager: SKB Rapport R-07-07 P, R eller TR. Integrated surface-subsurface water flow modelling of the Laxemar area Application of the hydrological model ECOFLOW Nikolay Sokrut, Kent Werner, Johan Holmén Golder Associates AB January 2007 This report concerns a study which was conducted for SKB. The conclusions and viewpoints presented in the report are those of the authors and do not necessarily coincide with those of the client. A pdf version of this document can be downloaded from www.skb.se Abstract Since 2002, the Swedish Nuclear Fuel and Waste Management Co (SKB) performs site investigations in the Simpevarp area, for the siting of a deep geological repository for spent nuclear fuel. The site descriptive modelling includes conceptual and quantitative modelling of surface-subsurface water interactions, which are key inputs to safety assessment and environmental impact assessment. Such modelling is important also for planning of continued site investigations. In this report, the distributed hydrological model ECOFLOW is applied to the Laxemar subarea to test the ability of the model to simulate surface water and near-surface groundwater flow, and to illustrate ECOFLOW’s advantages and drawbacks.
    [Show full text]
  • Signature Properties of Water: Their Molecular Electronic Origins
    Signature properties of water: Their molecular electronic origins Vlad P. Sokhana,1, Andrew P. Jonesb, Flaviu S. Cipciganb, Jason Craina,b, and Glenn J. Martynab,c aNational Physical Laboratory, Teddington, Middlesex TW11 0LW, United Kingdom; bSchool of Physics and Astronomy, The University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom; and cIBM Thomas J. Watson Research Center, Yorktown Heights, NY 10598 Edited by Paul Madden, University of Oxford, Oxford, United Kingdom, and accepted by the Editorial Board March 31, 2015 (received for review October 1, 2014) Water challenges our fundamental understanding of emergent combination of multiscale descriptions and statistical sampling have materials properties from a molecular perspective. It exhibits a led to insights into physical phenomena across biology, chemistry, uniquely rich phenomenology including dramatic variations in physics, materials science, and engineering (14). Significant progress behavior over the wide temperature range of the liquid into water’s now requires novel predictive models with reduced empirical input crystalline phases and amorphous states. We show that many-body that are rich enough to embody the essential physics of emergent responses arising from water’s electronic structure are essential mech- systems and yet simple enough to retain intuitive features (14). anisms harnessed by the molecule to encode for the distinguishing Typically, atomistic models of materials are derived from a features of its condensed states. We treat the complete set of these common strategy. Interactions are described via a fixed functional many-body responses nonperturbatively within a coarse-grained elec- form with long-range terms taken from low orders of perturbation tronic structure derived exclusively from single-molecule properties.
    [Show full text]
  • Groundwater Modeling System Linkage with the Framework for Risk Analysis in Multimedia Environmental Systems
    PNNL- 15654 Groundwater Modeling System Linkage with the Framework for Risk Analysis in Multimedia Environmental Systems G. Whelan K.J. Castleton February 2006 Prepared for the U.S. Nuclear Regulatory Commission Office of Nuclear Regulatory Research Division of Systems Analysis & Regulatory Effectiveness Rockville, Maryland 20852 under Contract DE-AC05-76RL01830 DISCLAIMER This report was prepared as an account of work sponsored by an agency of the United States Government. Neither the United States Government nor any agency thereof, nor Battelle Memorial Institute, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. PACIFIC NORTHWEST NATIONAL LABORATORY operated by BATTELLE for the UNITED STATES DEPARTMENT OF ENERGY under Contract DE-AC05-76RL01830 Printed in the United States of America Available to DOE and DOE contractors from the Office of Scientific and Technical Information, P.O. Box 62, Oak Ridge, TN 37831-0062; ph: (865) 576-8401 fax: (865) 576-5728 email: [email protected] Available to the public from the National Technical Information Service, U.S.
    [Show full text]