Transient and Transition Factors in Modeling Permafrost Thaw and Groundwater Flow

Permafrost covers approximately 24% of the Northern Hemisphere, and much of it is degrading, which causes infrastructure failures and ecosystem transitions. Understanding groundwater and heat flow processes in permafrost environments is challenging due to spatially and temporarily varying hydraulic connections between water above and below the near-surface discontinuous frozen zone. To characterize the transitional period of permafrost degradation, a three-dimensional model of a permafrost plateau that includes the supra-permafrost zone and surrounding wetlands was developed. The model is based on the Scotty Creek basin in the Northwest Territories, Canada. FEFLOW groundwater flow and heat transport modeling software is used in conjunction with the piFreeze plug-in, to account for phase changes between ice and water. The Simultaneous Heat and Water (SHAW) flow model is used to calculate ground temperatures and surface water balance, which are then used as FEFLOW boundary conditions. As simulating actual permafrost evolution would require hundreds of years of climate variations over an evolving landscape, whose geomorphic features are unknown, methodologies for developing permafrost initial conditions for transient simulations were investigated. It was found that a model initialized with a transient spin-up methodology, that includes an unfrozen layer between the permafrost table and ground surface, yields better results than with steady-state permafrost initial conditions. This study also demonstrates the critical role that variations in land surface and permafrost table microtopography, along with talik development, play in permafrost degradation. Modeling permafrost dynamics will allow for the testing of remedial measures to stabilize permafrost in high value infrastructure environments.     

Morphology and dynamics of the intertidal floodplain along the Amazon tidal river

Depositional environments along the tidal river downstream of Óbidos have been proposed as important sinks for up to one third of the sediment discharge from the Amazon River. However, the morphology of the intertidal floodplain and the dynamics of sediment exchange along this reach have yet to be described. River-bank surveys in five regions along the Amazon tidal river reveal a distinct transition in bank morphology between the upper, central and lower reaches of the tidal river. The upper tidal-river floodplain is defined by prominent natural levees that control the transfer of water and sediment between the mainstem Amazon River and its floodplain. Greater tidal influence in the central tidal river suppresses levee development, and tidal currents increase sediment transport into the distal parts of the floodplain. In the lower tidal river, the floodplain morphology closely resembles marine intertidal environments (e.g. mud flats, salt marshes), with dendritic tidal channels incising elevated vegetated flats. Theory, morphology and geochronology suggest that the dynamics of sediment delivery to the intertidal floodplain of the Amazon tidal river vary along its length due to the relative influence and coupling of fluvial and tidal dynamics. © 2018 John Wiley & Sons, Ltd.     

Integrated Deep Soil N and Groundwater Isotope Investigation of N Sources Captured by Municipal Wells

Identification of major nitrate sources that adversely impact groundwater quality in municipal well capture zones in areas of emerging nitrate contamination is essential to minimize leaching and prevent exceedance of the nitrate drinking water standard. Vertical profiles of nitrate leachate in deep soils provide an estimate of the amount of nitrate in transit beneath irrigated, row-cropped fields; depths of peak leachate; and the approximate rate of downward movement. Profiles of pore-water soil-nitrate concentrations in thick 60-feet (~18 m), fine-textured soils near Hastings, Nebraska clearly indicate that considerably more nitrate leached beneath furrow-irrigated than center-pivot irrigated fields. Peak leaching appeared to correlate with recorded periods of poor weather conditions during some growing seasons and may best be controlled by “spoon feeding” fertilizer to the crop through the sprinkler irrigation system at times of nutrient need. The presence of trace levels of atrazine and deethylatrazine to 60 feet (18 m) in core samples indicates that larger, more complex anthropogenic molecules also leach through the fine-textured soils. The light δ¹⁵N_NO3 values in the surficial groundwater beneath fertilized and irrigated cropland indicate that ammonium fertilizer is a major N source and suggest that the natural soil-N contribution is negligible. δ¹⁵N_NO3 values were most enriched in irrigation wells located within municipal well capture zones downgradient of a large feedlot. Dual isotope method (DIM) δ¹⁵N_NO3 and δ¹⁸O_NO3 values suggest that the Hastings’ municipal wells farther downgradient are contaminated with a mixture of nitrate from manure and commercial ammonium-based fertilizer. DIM values indicate an absence of denitrification, which has implications for long-term management of the water resources.     

Controls on the lateral channel-migration rate of braided channel systems in coarse non-cohesive sediment

Lateral movements of alluvial river channels control the extent and reworking rates of alluvial fans, floodplains, deltas, and alluvial sections of bedrock rivers. These lateral movements can occur by gradual channel migration or by sudden changes in channel position (avulsions). Whereas models exist for rates of river avulsion, we lack a detailed understanding of the rates of lateral channel migration on the scale of a channel belt. In a two-step process, we develop here an expression for the lateral migration rate of braided channel systems in coarse, non-cohesive sediment. On the basis of photographic and topographic data from laboratory experiments of braided channels performed under constant external boundary conditions, we first explore the impact of autogenic variations of the channel-system geometry (i.e. channel-bank heights, water depths, channel-system width, and channel slope) on channel-migration rates. In agreement with theoretical expectations, we find that, under such constant boundary conditions, the laterally reworked volume of sediment is constant and lateral channel-migration rates scale inversely with the channel-bank height. Furthermore, when channel-bank heights are accounted for, lateral migration rates are independent of the remaining channel geometry parameters. These constraints allow us, in a second step, to derive two alternative expressions for lateral channel-migration rates under different boundary conditions using dimensional analysis. Fits of a compilation of laboratory experiments to these expressions suggest that, for a given channel bank-height, migration rates are strongly sensitive to water discharges and more weakly sensitive to sediment discharges. In addition, external perturbations, such as changes in sediment and water discharges or base level fall, can indirectly affect lateral channel-migration rates by modulating channel-bank heights. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd. © 2019 The Author. Earth Surface Processes and Landforms published by John Wiley & Sons, Ltd.     

Local gravity field modeling using spherical radial basis functions and a genetic algorithm

Spherical Radial Basis Functions (SRBFs) can express the local gravity field model of the Earth if they are parameterized optimally on or below the Bjerhammar sphere. This parameterization is generally defined as the shape of the base functions, their number, center locations, bandwidths, and scale coefficients. The number/location and bandwidths of the base functions are the most important parameters for accurately representing the gravity field; once they are determined, the scale coefficients can then be computed accordingly. In this study, the point-mass kernel, as the simplest shape of SRBFs, is chosen to evaluate the synthesized free-air gravity anomalies over the rough area in Auvergne and GNSS/Leveling points (synthetic height anomalies) are used to validate the results. A two-step automatic approach is proposed to determine the optimum distribution of the base functions. First, the location of the base functions and their bandwidths are found using the genetic algorithm; second, the conjugate gradient least squares method is employed to estimate the scale coefficients. The proposed methodology shows promising results. On the one hand, when using the genetic algorithm, the base functions do not need to be set to a regular grid and they can move according to the roughness of topography. In this way, the models meet the desired accuracy with a low number of base functions. On the other hand, the conjugate gradient method removes the bias between derived quasigeoid heights from the model and from the GNSS/leveling points; this means there is no need for a corrector surface. The numerical test on the area of interest revealed an RMS of 0.48 mGal for the differences between predicted and observed gravity anomalies, and a corresponding 9 cm for the differences in GNSS/leveling points.     

Sensitivity of GRACE-derived estimates of groundwater-level changes in southern Ontario,Canada

Amidst changing climates, understanding the world’s water resources is of increasing importance. In Ontario, Canada, low water conditions are currently assessed using only precipitation and watershed-based stream gauges by the Conservation Authorities in Ontario and the Ministry of Natural Resources and Forestry. Regional groundwater-storage changes in Ontario are not currently measured using satellite data by research institutes. In this study, contributions from the Gravity Recovery and Climate Experiment (GRACE) data are compared to a hydrogeological database covering southern Ontario from 2003 to 2013, to determine the suitability of GRACE total water storage estimates for monitoring groundwater storage in this location. Terrestrial water storage data from GRACE were used to determine monthly groundwater storage (GWS) anomaly values. GWS values were also determined by multiplying groundwater-level elevations (from the Provincial Groundwater Monitoring Network wells) by specific yield. Comparisons of GRACE-derived GWS to well-based GWS data determined that GRACE is sufficiently sensitive to obtain a meaningful signal in southern Ontario. Results show that GWS values produced by GRACE are useful for identifying regional changes in groundwater storage in areas with limited available hydrogeological characterization data. Results also indicate that GRACE may have an ability to forecast changes in groundwater storage, which will become useful when monitoring climate shifts in the near future.     

A review of definitions of the Himalayan Main Central Thrust

International Journal of Earth Sciences - Most workers regard the Main Central Thrust (MCT) as one of the key high strain zones in the Himalaya because it accommodated at least 90 km of...     

A Paddock to reef monitoring and modelling framework for the Great Barrier Reef: Paddock and catchment component

Targets for improvements in water quality entering the Great Barrier Reef (GBR) have been set through the Reef Water Quality Protection Plan (Reef Plan). To measure and report on progress towards the targets set a program has been established that combines monitoring and modelling at paddock through to catchment and reef scales; the Paddock to Reef Integrated Monitoring, Modelling and Reporting Program (Paddock to Reef Program). This program aims to provide evidence of links between land management activities, water quality and reef health. Five lines of evidence are used: the effectiveness of management practices to improve water quality; the prevalence of management practice adoption and change in catchment indicators; long-term monitoring of catchment water quality; paddock & catchment modelling to provide a relative assessment of progress towards meeting targets; and finally marine monitoring of GBR water quality and reef ecosystem health. This paper outlines the first four lines of evidence.     

Spatial dynamics of overbank sedimentation in floodplain systems

Floodplains provide valuable social and ecological functions, and understanding the rates and patterns of overbank sedimentation is critical for river basin management and rehabilitation. Channelization of alluvial systems throughout the world has altered hydrological and sedimentation processes within floodplain ecosystems. In the loess belt region of the Lower Mississippi Alluvial Valley of the United States, channelization, the geology of the region, and past land-use practices have resulted in the formation of dozens of valley plugs in stream channels and the formation of shoals at the confluence of stream systems. Valley plugs completely block stream channels with sediment and debris and can result in greater deposition rates on floodplain surfaces. Presently, however, information is lacking on the rates and variability of overbank sedimentation associated with valley plugs and shoals.We quantified deposition rates and textures in floodplains along channelized streams that contained valley plugs and shoals, in addition to floodplains occurring along an unchannelized stream, to improve our understanding of overbank sedimentation associated with channelized streams. Feldspar clay marker horizons and marker poles were used to measure floodplain deposition from 2002 to 2005 and data were analyzed with geospatial statistics to determine the spatial dynamics of sedimentation within the floodplains.Mean sediment deposition rates ranged from 0.09 to 0.67 cm/y at unchannelized sites, 0.16 to 2.27 cm/y at shoal sites, and 3.44 to 6.20 cm/y at valley plug sites. Valley plug sites had greater rates of deposition, and the deposited sediments contained more coarse sand material than either shoal or unchannelized sites. A total of 59 of 183 valley plug study plots had mean deposition rates > 5 cm/y. The geospatial analyses showed that the spatial dynamics of sedimentation can be influenced by the formation of valley plugs and shoals on channelized streams; however, responses can vary. Restoration efforts in the region need to have basinwide collaboration with landowners and address catchment-scale processes, including the geomorphic instability of the region, to be successful.