Search NASA⌕ Search

SEARCH · Search NASA

Results for “subsurface flow”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Structure and Shape of Surface-Mediated Assembly of Surfactants

Achieving controls on the self-assembly of cationic and block copolymer micelles has significant implications for advancing novel material systems for energy, environmental, and biological applications and for tuning fluid flow behavior in subsurface geologic environments. While it is known that the aggregation of micelles is influenced by the composition of the surrounding fluid, the influence of a solid interface has been less explored. In this study, we probe the organization of cetyltrimethylammonium bromide (CTAB) micelles in the absence and presence of the Pluronic P123 block copolymer and quartz substrate using transmission and grazing-incidence small-angle X-ray scattering measurements and classical molecular dynamics (MD) simulations. In the absence of the quartz interface, CTAB with and without P123 molecules assemble as ellipsoid core–shell micelles. Densification of the core on adding P123 is noted from the decrease in the core radius from 15.9 to 14.6 Å. Here, the presence of a quartz substrate causes the micelles to elongate, which is noted by the emergence of a power-law slope in the low q region (<0.02 Å –1 ). Moreover, in the presence of both P123 and the quartz substrate, the micelle shape changes from ellipsoid core–shell to cylindrical core–shell, and a significantly higher number of aggregates are formed. The higher number of aggregates, and faster aggregation kinetics are linked to the organization of the solvent structures as noted from the MD simulations. These findings demonstrate the effects of adding a block copolymer and including a quartz substrate on the self-assembly of cationic (CTAB) micelles, which can be used to inform the design of functional material systems that harness surfactant chemistries.

42 ENGINEERING↗
DOE OSTI2023

DPFEHM: a differentiable subsurface physics simulator

The Earth’s subsurface is a key resource that provides energy via fossil fuels and geothermal sources, stores drinking water, and is used in the fight against climate change via carbon sequestration. Simulating the physical processes that occur in the Earth’s subsurface with computers enables better use of this resource. DPFEHM is a Julia package that includes computer models with a focus on the Earth’s subsurface, especially fluid flow, which is critical for the aforementioned applications. DPFEHM is able to solve the groundwater flow equations (single phase flow), Richards equation (air/water), the advection-dispersion equation, and the 2d wave equation. One of the key features of DPFEHM is that it supports automatic differentiation, so it can be integrated into machine learning workflows using frameworks such as Flux or PyTorch. The automatic differentiation capabilities give it the same performance as adjoint methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

CO2 Saline Storage Talk at the University of Wyoming

These are slides presented virtually to a class on Carbon Capture and Storage (CCS) at the University of Wyoming in Laramie WY. The talk was on the basics of carbon storage in saline formations and reviewed basic geology, important physical properties of CO2 and brine, how fluids flow in the subsurface, environmental risks associated with CO2 storage, the regulations governing CO2 saline storage, and high-level design considerations. The FECM/NETL CO2 Saline Storage Cost Model was used to illustrate how geologic properties vary in different regions of the US (PA, IL and WY) and how the performance of CO2 storage and the cost of CO2 CO2 saline storage depends on geology (October 23, 2024).

Morgan, David↗
NASA NTRS1990

A rough-surface thermophysical model for airless planets

A model for determining diurnal temperatures in spherical-section depressions and which encompasses both subsurface heat-flow and direct and scattered sunlight effects is presently applied to the disk-integrated thermal emission of a rough planetary surface with nonzero thermal inertia. Attention is given to the variation with roughness and thermal inertia of the beaming parameter eta, which characterizes zero-phase thermal emission by comparison with a smooth, nonrotating body and is almost independent of albedo for a given surface roughness. The thermal phase curve of Ceres is noted to be well matched by the model features of (1) prograde rotation, (2) 44-deg rms surface slope, and (3) a thermal inertia that is 30 percent of the lunar value.

Spencer, John R.↗
NASA NTRS2003

Global Characterization of Polygonally Fractured Terrain on Venus and Implications for a Climate Change Origin

On Earth, polygonal fractures with spacings of 1-100 cm form in lava flows that cool slowly under isotropic stress conditions. On Venus, polygonal fractures are observed with typical fracture spacings of 1-2 km, covering regions 10s to 100s of km across. These features have been proposed to form via cooling of lava flows or above subsurface intrusions, or due to surface cooling caused by climate change. A watershed algorithm was modified to search the entire Magellan image database autonomously for terrains with polygonal fractures. Here we characterize the 204 polygon terrains found using this method (including over 170 newly identified regions) with respect to size, morphology, stratigraphy, and geologic setting. These results are used to examine different models of origin.

Smrekar, S. E.↗
NASA NTRS2006

Europa Surface Radiation Environment for Lander Assessment

The Jovian magnetospheric particle environment at Europa's surface is critical to assessment of landed astrobiological experiments in three respects: (1) the landing site must be chosen for the best prospects for detectable organic or inorganic signs of Life, e.g. regions of freshly emergent flows from the subsurface; (2) lander systems must reach the surface through the Jovian magnetospheric environment and operate long enough on the surface to return useful data; (3) lander instrumentation must be capable of detecting signs of life in the context of the local environmental radiation and associated chemistry. The Galileo, Voyager, and Pioneer missions have provided a wealth of data on energetic particle intensities throughout the Jovian magnetosphere including from many flybys of Europa. cumulative radiation dosages for spacecraft enroute to Europa can be well characterized, but knowledge of the surface radiation environment is very limited. Energetic electrons should primarily impact the trailing hemisphere with decreasing intensity towards the center of the leading hemisphere and are the most significant radiation component down to meter depths in the surface regolith due to secondary interactions. Observed surface distribution for sulfates is suggestive of electron irradiation but may have alternative interpretations. Having much-larger magnetic gyroradii than electrons, energetic protons and heavier ions irradiate more of the global surface. The particular orientations of electron, proton, and ion gyromotion would project into corresponding directional (e.g., east-west) anisotropies of particle flu into the surface. Particular topographic features at the landing site may therefore offer shielding from part of the incident radiation.

Cooper, John F.↗
NASA NTRSDate not supplied

A Notional Artemis Lunar Surface Exploration Package (ArLSEP) based on the Gandalf Staff Platform

Introduction: The Artemis program is planning to deliver crew and cargo to the lunar surface, but there is no current package for supporting lunar in-struments and experiments similar to the Apollo Lunar Surface Exploration Package (ALSEP). This abstract provides a possible concept for such a package using the Gandalf Staff Platform as a common core. Gandalf Staff: The Gandalf Staff is an early prototype system developed over FY’21/FY’22 using NASA Science Technology Mission Directorate (STMD) Center Information Fund (CIF) grants to de-sign, build and test “proof-of-concept” components. These components include a 24v battery powered monopole that powers a suite of subsystems, including a Graphical User Interface (GUI) for crew, surface voice and data communications, Lunar Search and Rescue (LunaSAR) navigation and communications, LiDAR, field site external lighting, 360-degree camera, and a geothermal instrument for measuring sub-surface temperature gradient. The staff can be carried independently by an Extra-Vehicular Activity (EVA) astronaut, or can be mounted into a tripod for “hands free” support at a surface site being investigated. The staff can be attached to an external solar array and power storage system for long-duration operations. [1,2] ALSEP: An ASLEP flew on each mission Apollo 12 to Apollo 17. For Apollo 11, a simplified packaged called the Early Apollo Scientific Experiments Pack-age (EASEP) was flown. Each package included a “Central Station” that provided the power and communications connected to a variety of instruments and sensors. The power was provided by a Radioisotope Thermoelectric Generator (RTG) fueled by Plutoni-um-238 generating 70 watts of power (initially, decayed over time) [3]. The communications system provide for direct to Earth data transfer from the lunar surface. Each pack-age was stowed externally in the Lunar Module (LM) Scientific Equipment (SEQ) bay with a mass up to 163 kg (Apollo 17). The crew unloaded the ALSEP from the LM and deployed the instruments on the lunar surface. Although designed to operate for only 1 year, many sites operated for up to 8 years successfully [4]. The Active Seismic Experiment (ASE) included 3 geophones for detecting seismic waves created by mortars and thumpers deployed by the crew. Other active experiments measured the lunar atmosphere, the heat flow in the subsurface, the lunar gravity and potential gravity waves, the lunar magnetic field, the solar wind and plasma interactions in cislunar space. Passive experiments included collectors for dust and cosmic rays, and retroreflectors for precise measurements of distance using a laser from Earth. The ALSEP program continues to generate insights into lunar formation and evolution. ArLSEP Concepts: The lunar surface science package for the Artemis program will hopefully exceed the capability of the ALSEP. There are multiple issues for discussion leading to the design of a new ArLSEP, needing requirements definition from the science community, NASA mission architecture, and NASA budget planners. 1. Delivery Mechanism Two possible projects currently provide capability to deliver scientific cargo to the lunar surface: 1) the Commercial Lunar Payload Services (CLPS) [5] and the Human Landing System (HLS) [6, 7]. Each project is controlled by a different organization within NASA and budgeted with different criteria although both support lunar exploration. The HLS system delivers crew (and potentially cargo) to human landing sites. If an ArLSEP is “predeployed” to such a site, the design must include power (either from the vehicle or independently) to keep the electronics functioning until deployed by the crew. If an ArLSEP is delivered on a vehicle after the crew is present on the lunar surface, safety protocols require adequate distance from the humans for impact from descent propelled sur-face regolith ejecta. This distance can not exceed the capability of the crew to walk (if no rover) to the vehicle for ArLSEP deployment. 2. Overall Guidelines The general design of ArLSEP will likely follow the ALSEP with a common system for communications and power; however, significant architecture differences between Apollo and Artemis exist. Power: The RTG will not be available for early Artemis missions nor likely follow-on Lunar Exploration Transportation Services (LETS) missions [8]. Thus, ArLSEP power must be supplied by solar arrays with sufficient battery capability to “keep alive” necessary electronics during any lunar surface eclipse period. Communication: The Artemis program is developing a series of communications satellites for lunar orbit to provide surface transmission of data and voice to Earth. Called “LunaNET”, this network is component useful for ArLSEP since south polar locations may not always have direct “line-of-sight” to Earth [9]. 3. Concept of Operations (ConOps) The general ConOps for ArLSEP is to deliver the package to lunar surface before the crew arrives, and then have the crew deploy the package after some period of time. This requires coordinated design (for power systems) and launch window (for schedule) on both the cargo and crew missions. Once the ArLSEP is deployed, it will operate autonomously for a number of years. It should be designed to be EVA compatible for crew maintenance and upgrade. 4. Notional Design (for discussion purpose only) The landing site near the South Pole is expected to have no eclipse cycle exceeding 5 days, so the “keep alive” power is 144 hours (6 days to include margin). A 12v ArLSEP will use rechargeable LiFePO4 cells, which are common in the Electric Vehicle (EV) industry. With a current of 5 amps and a 125 watt system, the mass is about 90kg. The comm. system and structure adds another 10kg, thus the “Central Station” is approximately 100kg. The solar power is collected on four arrays (each 2m above the surface), and the entire ArLSEP is designed to stow in a 2m x 1m x 1m volume. The experiment and instrument design will vary for each installation and add mass to the total (although they are expected to fit within the 2m3 volume). Seismic wave generation will likely not be provided with mortars, thus an electric “thumper” will be required. Active instruments such as imaging systems and sensing instruments will benefit from the additional power and communication capability provided by ArLSEP. Passive systems such as retroreflectors, witness plates, and cosmic dust collectors can be added to either the landing vehicle and/or the ArLSEP. With repeated HLS missions to the same human site, the ArLSEP can be expanded and easily maintained for long duration science collection on the lunar surface.

ALSEP↗
NASA NTRS1985

A cool anomaly off northern California - An investigation using IR imagery and in situ data

The OPTOMA (Ocean Prediction Through Observation, Modeling and Analysis) program is developing an ocean descriptive-predictive system for four-dimensional data assimilation. It is presently concerned with the mesoscale variability in the California Current System (CCS). The present paper has the objective to assess the relationship of surface temperature structure to subsurface temperature structure and flow fields. Surface temperature anomalies are related to the mesoscale horizontal advective and subsurface mass fields. The cool anomaly off northern California in summer 1982 is discussed, taking into account the temperature structure along transects, T-S variations and inferred water masses, variation of horizontal temperature patterns with depth, and synoptic information from IR data.

Rienecker, M. M.↗
DOE OSTI2023

A Perturbative Solution for Nonlinear Stratified Upwelling over a Frictional Slope

Abstract A perturbative solution of simplified primitive equations for nonlinear weakly stratified upwelling over a frictional slope is found that resolves the vertical structure of velocity fields and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The solution uses assumptions consistent with the model proposed by Lentz and Chapman, including a steady-state, constant cross-shore density gradient, no alongshore gradients, laterally inviscid, and consideration of cross-shore advection of alongshore momentum. The solution resolves the vertical structure of velocity fields (including subsurface maxima of compensational flow, not resolved by Lentz and Chapman) and can satisfy Ertel’s potential vorticity conservation in the stratified inviscid interior. The dynamics are similar to Lentz and Chapman; bottom stress balances alongshore wind stress in a homogeneous density ocean and is replaced by nonlinear cross-shore transport of alongshore momentum as the Burger number (S=αN/f, whereα,N, andfare the bottom slope, buoyancy frequency, Coriolis frequency, respectively) increases. When the solution uses the empirical relation between cross-shore and vertical density gradients proposed by Lentz and Chapman, vorticity conservation is not satisfied and the nonlinear momentum transport estimated by the solution linearly increases withS, asymptotically matching Lentz and Chapman forS< 1. When the solution conserves interior potential vorticity, the momentum transport is proportional toS 2 forS< 1 and is in better agreement with numerical simulations.

Oceanography↗
DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗
DOE OSTI2021

Multilevel Graph Partitioning for Three-Dimensional Discrete Fracture Network Flow Simulations

We present a topology-based method for mesh-partitioning in three-dimensional discrete fracture network (DFN) simulations that takes advantage of the intrinsic multi-level nature of a DFN. DFN models are used to simulate flow and transport through low-permeability fractured media in the subsurface by explicitly representing fractures as discrete entities. The governing equations for flow and transport are numerically integrated on computational meshes generated on the interconnected fracture networks. Modern high-fidelity DFN simulations require high-performance computing on multiple processors where performance and scalability depends partially on obtaining a high-quality partition of the mesh to balance work-loads and minimize communication across all processors. The discrete structure of a DFN naturally lends itself to various graph representations, which can be thought of as coarse-scale representations of the computational mesh. Using this concept, we develop two applications of the multilevel graph partitioning algorithm to partition the mesh of a DFN. In the first, we project a partition of the graph based on the DFN topology onto the mesh of the DFN and in the second, this DFN-based projection is used as the initial condition for further partitioning refinement of the mesh. We compare the performance of these methods with standard multi-level graph partitioning using graph-based metrics (cut, imbalance, partitioning time), computational-based metrics (FLOPS, iterations, solver time), and total run time. The DFN-based and the mesh-based partitioning methods are comparable in terms of the graph-based metrics, but the time required to obtain the partition is several orders of magnitude faster using the DFN-based partitions. The computation-based metrics show comparable performance between both methods so, in combination, the DFN-based partitions are several orders of magnitude faster than the mesh-based partition. Furthermore, the method which uses the DFN-partition solution as the initial condition of the mesh partition provided cut and imbalance values that were close to the mesh-based partition but in a fraction of the time. In turn, this hybrid method outperformed both of the other methods in terms of the total run time.

58 GEOSCIENCES↗
NASA NTRS1976

Revised lunar heat-flow values

The 3.5- and 2-year subsurface temperature histories at the Apollo 15 and 17 heat-flow sites have been analyzed, and the results yield significantly lower thermal conductivity determinations than the results of previous short-term experiments. The thermal conductivity determined by probes at a depth of about 150 cm and 250 cm lies in the range 0.9-1.3 times 10 to the -4th W/cm K. On the basis of measurements of variations of surface thorium abundance and inferred crustal thicknesses, the average global heat flux is estimated to be about 1.8 microwatts/sq cm. This requires a uranium concentration of 46 ppb.

Langseth, M. G.↗
NASA NTRS2003

The Broken Belt: Meteorite Concentrations on Stranded Ice

Since the first Antarctic meteorite concentrations were discovered more than 25 years ago, many theories regarding the role of iceflow in the production of meteorite concentrations have been put forward, and most agree on the basic principles. These models suggest that as the East Antarctic icesheet flows toward the margins of the continent, meteorites randomly located within the volume of ice are transported toward the icesheet margin. Where mountains or subsurface obstructions block glacial flow, diversion of ice around or over an obstruction reduces horizontal ice movement rates adjacent to the barriers and creates a vertical (upward) component of movement. If local mechanisms for ice loss (ablation) exist at such sites, an equilibrium surface will develop according to the balance between ice supply and loss, and the cargo of meteorites is exhumed on a blue ice surface. The result is a conceptual conveyor belt bringing meteorite-bearing volumes of ice from the interior of the continent to stagnant or slowmoving surfaces where ice is then lost and a precious cargo is left as a lag deposit. Cassidy et al. provides an excellent overview of how this model has been adapted to several Antarctic stranding surfaces.

Harvey, R. P.↗
DOE OSTI2025

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

From Fluid Flow to Coupled Processes in Fractured Rock: Recent Advances and New Frontiers

Abstract Quantitative predictions of natural and induced phenomena in fractured rock is one of the great challenges in the Earth and Energy Sciences with far‐reaching economic and environmental impacts. Fractures occupy a very small volume of a subsurface formation but often dominate fluid flow, solute transport and mechanical deformation behavior. They play a central role in CO 2 sequestration, nuclear waste disposal, hydrogen storage, geothermal energy production, nuclear nonproliferation, and hydrocarbon extraction. These applications require predictions of fracture‐dependent quantities of interest such as CO 2 leakage rate, hydrocarbon production, radionuclide plume migration, and seismicity; to be useful, these predictions must account for uncertainty inherent in subsurface systems. Here, we review recent advances in fractured rock research covering field‐ and laboratory‐scale experimentation, numerical simulations, and uncertainty quantification. We discuss how these have greatly improved the fundamental understanding of fractures and one's ability to predict flow and transport in fractured systems. Dedicated field sites provide quantitative measurements of fracture flow that can be used to identify dominant coupled processes and to validate models. Laboratory‐scale experiments fill critical knowledge gaps by providing direct observations and measurements of fracture geometry and flow under controlled conditions that cannot be obtained in the field. Physics‐based simulation of flow and transport provide a bridge in understanding between controlled simple laboratory experiments and the massively complex field‐scale fracture systems. Finally, we review the use of machine learning‐based emulators to rapidly investigate different fracture property scenarios and accelerate physics‐based models by orders of magnitude to enable uncertainty quantification and near real‐time analysis.

58 GEOSCIENCES↗
DOE OSTI2020

A higher-order finite element reactive transport model for unstructured and fractured grids

Abstract This work presents a new reactive transport framework that combines a powerful geochemistry engine with advanced numerical methods for flow and transport in subsurface fractured porous media. Specifically, the PhreeqcRM interface (developed by the USGS) is used to take advantage of a large library of equilibrium and kinetic aqueous and fluid-rock reactions, which has been validated by numerous experiments and benchmark studies. Fluid flow is modeled by the Mixed Hybrid Finite Element (FE) method, which provides smooth velocity fields even in highly heterogenous formations with discrete fractures. A multilinear Discontinuous Galerkin FE method is used to solve the multicomponent transport problem. This method is locally mass conserving and its second order convergence significantly reduces numerical dispersion. In terms of thermodynamics, the aqueous phase is considered as a compressible fluid and its properties are derived from a Cubic Plus Association (CPA) equation of state. The new simulator is validated against several benchmark problems (involving, e.g., Fickian and Nernst-Planck diffusion, isotope fractionation, advection-dispersion transport, and rock-fluid reactions) before demonstrating the expanded capabilities offered by the underlying FE foundation, such as high computational efficiency, parallelizability, low numerical dispersion, unstructured 3D gridding, and discrete fraction modeling.

58 GEOSCIENCES↗
DOE OSTI2024

Calibration of V-Notch and Compound Weirs for Subsurface Drainage Water Level Control Structures

Highlights Accurate discharge estimation is important when evaluating edge-of-field conservation practices. V-notch weir equations were developed for three sizes of subsurface drainage water level control structures. Compound weir equations were developed for subsurface drainage water level control structures. The compound weir equation accurately estimates discharge for flows within and overtopping the V-notch. Abstract.Numerous edge-of-field conservation practices use subsurface drainage water level control structures to monitor water levels and estimate discharge. In a control structure, procedures for calculating discharge when flow depth (head) exceeds the V-notch depth and overflows in the rectangular portion of the compound weir (CW) are ambiguous. In this study, we developed calibration equations for V-notch weirs in Agri Drain inline water level control structures of different sizes for flows within the V-notch and overtopping flow events. The discharge equation for overtopping events (Q CW , L s -1 ) was determined as: Q CW = a 1 (h b1 -h 1 b1 )+a 2 (W e -W v )h 1 b2 , where h and h 1 are heads above vertex/bottom and top of V-notch (cm), respectively, W is the effective crest width of rectangular weir (cm), W v is the top width of V-notch (cm), a 1 and b 1 are parameters for V-notch weir obtained by calibration, and a 2 and b 2 are calibration parameters for rectangular weir obtained from literature. Results were compared with a weir equation available in the literature (Q V+R ), which combines a V-notch equation with a head equal to V-depth and a rectangular weir equation for flow above V-depth. Discharge at overflow was estimated with high accuracy with Q CW, whereas Q V+R underestimated discharge (e.g., PBIAS of 0.67% vs. 17.82% for a 15.2 cm structure). An example using Q V+R resulted in a 14% lower annual estimation of nitrate-N load diverted to a saturated buffer than Q CW due to underestimation of drainage discharge during overflow events. Results suggest that the developed equation (Q CW ) accurately estimates discharge and will thus improve the estimated N load compared to Q V+R . Keywords: Compound weir, Flow monitoring, Subsurface drainage, V-notch weir, Water level control structure, Weir calibration.

Agriculture↗