MODFLOW as a Configurable Multi-Model Hydrologic Simulator

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.     

Evaluation of building separation distance required to prevent pounding during strong earthquakes

Presented herein is an analytical procedure for evaluating the separation distance required between two buildings to prevent pounding during strong earthquakes. The procedure is based on equivalent linearization of non-linear hysteretic behaviour and application of the well-known CQC method of weighting normal mode responses. The numerical results obtained are compared with corresponding results obtained using the well-known SRSS and ABS methods. © 1997 John Wiley & Sons, Ltd.     

Marine inundated archaeological sites and paleofluvial systems: examples from a karst-controlled continental shelf setting in Apalachee Bay,Northeastern Gulf of Mexico

Underwater geoarchaeological research in Apalachee Bay, in the northeastern Gulf of Mexico off northwest Florida, has enabled the reconstruction of portions of the karst-controlled paleodrainage system, the discovery of several inundated prehistoric archaeological sites, and the exposure of sediments accumulated during the drowning of the continental shelf. Diagnostic artifacts discovered at the sites included chipped stone tools and debitage indicating Paleoindian, Early Archaic, and Middle Archaic occupation. A geoarchaeological model using terrestrial analogs was used to locate and investigate inundated sites. Methods employed include seismic profiling, vibracoring, diver tow surveys, diver collection transects, and induction dredge excavations. We document evidence for sea-level rise, related environmental succession and site formation processes for indundated prehistoric sites in the Apalachee Bay region from approximately 8000 to 6000 yr B.P. © 1997 John Wiley & Sons, Inc.     

Artificial neural network simulation for prediction of suspended sediment concentration in the River Ramganga,Ganges Basin,India

The relation between the water discharge (Q) and suspended sediment concentration (SSC) of the River Ramganga at Bareilly, Uttar Pradesh, in the Himalayas, has been modeled using Artificial Neural Networks (ANNs). The current study validates the practical capability and usefulness of this tool for simulating complex nonlinear, real world, river system processes in the Himalayan scenario. The modeling approach is based on the time series data collected from January to December (2008-2010) for Q and SSC. Three ANNs (T1-T3) with different network configurations have been developed and trained using the Levenberg Marquardt Back Propagation Algorithm in the Matlab routines. Networks were optimized using the enumeration technique, and, finally, the best network is used to predict the SSC values for the year 2011. The values thus obtained through the ANN model are compared with the observed values of SSC. The coefficient of determination (R2), for the optimal network was found to be 0.99. The study not only provides insight into ANN modeling in the Himalayan river scenario, but it also focuses on the importance of understanding a river basin and the factors that affect the SSC, before attempting to model it. Despite the temporal variations in the study area, it is possible to model and successfully predict the SSC values with very simplistic ANN models.     

Patterns and drivers of peat topographic changes determined from Structure-from-Motion photogrammetry at field plot and laboratory scales

Little is known about the spatial and temporal variability of peat erosion nor some of its topographic and weather-related drivers. We present field and laboratory observations of peat erosion using Structure-from-Motion (SfM) photogrammetry. Over a 12 month period, 11 repeated SfM surveys were conducted on four geomorphological sites of 18–28 m² (peat hagg, gully wall, riparian area and gully head) in a blanket peatland in northern England. A net topographic change of –14 to +30 mm yr^–1 for the four sites was observed during the whole monitoring period. Cold conditions in the winter of 2016 resulted in highly variable volume change (net surface topographic rise first and lowering afterwards) via freeze–thaw processes. Long periods of dry conditions in the summer of 2017 led to desiccation and drying and cracking of the peat surface and a corresponding surface lowering. Topographic changes were mainly observed over short-term intervals when intense rainfall, flow wash, needle-ice production or surface desiccation was observed. In the laboratory, we applied rainfall simulations on peat blocks and compared the peat losses quantified by traditional sediment flux measurements with SfM derived topographic data. The magnitude of topographic change determined by SfM (mean value: 0.7 mm, SD: 4.3 mm) was very different to the areal average determined by the sediment yield from the blocks (mean value: –0.1 mm, SD: 0.1 mm). Topographic controls on spatial patterns of topographic change were illustrated from both field and laboratory surveys. Roughness was positively correlated to positive topographic change and was negatively correlated to negative topographic change at field plot scale and laboratory macroscale. Overall, the importance of event-scale change and the direct relationship between surface roughness and the rate of topographic change are important characteristics which we suggest are generalizable to other environments. © 2018 John Wiley & Sons, Ltd.     

Sediment and fluvial particulate carbon flux from an eroding peatland catchment

Erosion and the associated loss of carbon is a major environmental concern in many peatlands and remains difficult to accurately quantify beyond the plot scale. Erosion was measured in an upland blanket peatland catchment (0.017 km²) in northern England using structure-from-motion (SfM) photogrammetry, sediment traps and stream sediment sampling at different spatial scales. A net median topographic change of –27 mm yr^–1 was recorded by SfM over the 12-month monitoring period for the entire surveyed area (598 m²). Within the entire surveyed area there were six nested catchments where both SfM and sediment traps were used to measure erosion. Substantial amounts of peat were captured in sediment traps during summer storm events after two months of dry weather where desiccation of the peat surface occurred. The magnitude of topographic change for the six nested catchments determined by SfM (mean value: 5.3 mm, standard deviation: 5.2 mm) was very different to the areal average derived from sediment traps (mean value: –0.3 mm, standard deviation: 0.1 mm). Thus, direct interpolation of peat erosion from local net topographic change into sediment yield at the catchment outlet appears problematic. Peat loss measured at the hillslope scale was not representative of that at the catchment scale. Stream sediment sampling at the outlet of the research catchment (0.017 km²) suggested that the yields of suspended sediment and particulate organic carbon were 926.3 t km^–2 yr^–1 and 340.9 t km^–2 yr^–1, respectively, with highest losses occurring during the autumn. Both freeze–thaw during winter and desiccation during long periods of dry weather in spring and summer were identified as important peat weathering processes during the study. Such weathering was a key enabler of subsequent fluvial peat loss from the catchment. © 2019 John Wiley & Sons, Ltd.     

Radiocarbon Ages and Environments of Deposition of the Wono and Trego Hot Springs Tephra Layers in the Pyramid Lake Subbasin, Nevada

Uncalibrated radiocarbon data from core PLC92B taken from Wizards Cove in the Pyramid Lake subbasin indicate that the Trego Hot Springs and Wono tephra layers were deposited 23,200 ± 300 and 27,300 ± 300¹⁴C yr B.P. (uncorrected for reservoir effect). Sedimentological data from sites in the Pyramid Lake and Smoke Creek–Black Rock Desert subbasins indicate that the Trego Hot Springs tephra layer was deposited during a relatively dry period when Pyramid Lake was at or below its spill point (1177 m) to the Winnemucca Lake subbasin. The Wono tephra layer was deposited when lake depth was controlled by spill across Emerson Pass sill (1207 m) to the Smoke Creek–Black Rock Desert subbasin.¹⁸O data from core PLC92B also support the concept that the Trego Hot Springs tephra fell into a relatively shallow Pyramid Lake and that the Wono tephra fell into a deeper spilling lake.     

Predicting sediment transport by interrill overland flow on rough surfaces

Modelling soil erosion requires an equation for predicting the sediment transport capacity by interrill overland flow on rough surfaces. The conventional practice of partitioning total shear stress into grain and form shear stress and predicting transport capacity using grain shear stress lacks rigour and is prone to underestimation. This study therefore explores the possibility that inasmuch as surface roughness affects flow hydraulic variables which, in turn, determine transport capacity, there may be one or more hydraulic variables which capture the effect of surface roughness on transport capacity suffciently well for good predictions of transport capacity to be achieved from data on these variables alone. To investigate this possibility, regression analyses were performed on data from 1506 flume experiments in which discharge, slope, water temperature, rainfall intensity, and roughness size, shape and concentration were varied. The analyses reveal that 89·8 per cent of the variance in transport capacity can be accounted for by excess flow power and flow depth. Including roughness size and concentration in the regression improves that explained variance by only 3·5 per cent. Evidently, flow depth, when used in combination with excess flow power, largely captures the effect of surface roughness on transport capacity. This finding promises to simplify greatly the task of developing a general sediment equation for interrill overland flow on rough surfaces. Copyright © 1998 John Wiley & Sons, Ltd.     

Wetland Loss Patterns and Inundation-Productivity Relationships Prognosticate Widespread Salt Marsh Loss for Southern New England

Tidal salt marsh is a key defense against, yet is especially vulnerable to, the effects of accelerated sea level rise. To determine whether salt marshes in southern New England will be stable given increasing inundation over the coming decades, we examined current loss patterns, inundation-productivity feedbacks, and sustaining processes. A multi-decadal analysis of salt marsh aerial extent using historic imagery and maps revealed that salt marsh vegetation loss is both widespread and accelerating, with vegetation loss rates over the past four decades summing to 17.3 %. Landward retreat of the marsh edge, widening and headward expansion of tidal channel networks, loss of marsh islands, and the development and enlargement of interior depressions found on the marsh platform contributed to vegetation loss. Inundation due to sea level rise is strongly suggested as a primary driver: vegetation loss rates were significantly negatively correlated with marsh elevation (r ²?=?0.96; p?=?0.0038), with marshes situated below mean high water (MHW) experiencing greater declines than marshes sitting well above MHW. Growth experiments with Spartina alterniflora, the Atlantic salt marsh ecosystem dominant, across a range of elevations and inundation regimes further established that greater inundation decreases belowground biomass production of S. alterniflora and, thus, negatively impacts organic matter accumulation. These results suggest that southern New England salt marshes are already experiencing deterioration and fragmentation in response to sea level rise and may not be stable as tidal flooding increases in the future.