Predicting gully erosion using landform evolution models: Insights from mining landforms

Incision as a result of fluvial erosion is an important process to model when simulating landform evolution. For gullies, it is apparent that coupled with the processes that cause incision there must be a range of processes that stop incision. Once started, rills and gullies will grow infinitely without a reduction in support area and/or being arrested by deposition and armouring. Some of these processes have been well studied under the heading of inter-rill erosion. Other limiting processes are related to the shape of the landform and how downstream deposition areas are linked geomorphically to the upstream gullies. Armouring is also an important process that reduces gully incision and extension, where the gully erodes to bedrock and the resistant base limits further development. Post-mining landscapes are new surfaces with new materials and provide the opportunity to examine gully initiation, extension and stabilization. The work presented here has largely been driven by the mining industry, where there has been a need to assess erosion over hazardous wastes like mine tailings and low-level nuclear waste. We demonstrate the usefulness of computer-based landscape evolution models and the more recent soilscape models (that include both surface and subsurface processes) to understand both fluvial and diffusive processes as well as armouring in a digital elevation model framework (as well as landscape evolution). Landscape evolution models provide insights into complex non-linear systems such as gullies. A key need is that of field data to parameterize and validate the models. It is argued that current models have more capability than field data available for parameterization and importantly the validation of model outputs.     

Spatial prediction of temporal soil moisture dynamics using HYDRUS‐1D

This article investigates the soil moisture dynamics within two catchments (Stanley and Krui) in the Goulburn River in NSW during a 3‐year period (2005–2007) using the HYDRUS‐1D soil water model. Sensitivity analyses indicated that soil type, and leaf area index were the key parameters affecting model performance. The model was satisfactorily calibrated on the Stanley microcatchment sites with a single point rainfall record from this microcatchment for both surface 30 cm and full‐profile soil moisture measurements. Good correlations were obtained between observed and simulated soil water storage when calibrations for one site were applied to the other sites. We extended the predictions of soil moisture to a larger spatial scale using the calibrated soil and vegetation parameters to the sites in the Krui catchment where soil moisture measurement sites were up to 30 km distant from Stanley. Similarly good results show that it is possible to use a calibrated soil moisture model with measurements at a single site to extrapolate the soil moisture to other sites for a catchment with an area of up to 1000 km² given similar soils and vegetation and local rainfall data. Site predictions were effectively improved by our simple data assimilation method using only a few sample data collected from the site. This article demonstrates the potential usefulness of continuous time, point‐scale soil moisture data (typical of that measured by permanently installed TDR probes) and simulations for predicting the soil wetness status over a catchment of significant size (up to 1000 km²). Copyright © 2012 John Wiley & Sons, Ltd.     

Sub-ice geology inland of the Transantarctic Mountains in light of new aerogeophysical data

The Transantarctic Mountains are a major geologic boundary that bisects the Antarctic continent, separating the low-lying, tectonically active terrains of West Antarctica from the East Antarctic craton. A new comprehensive aerogeophysical data set, extending 1150 km from the Ross Sea into the interior of East Antarctica provides insights into the complex structure inland of the Transantarctic Mountains. Geophysical maps, compiled from 21 000 km of gravity, magnetic and subglacial topography data, outline the boundaries of several geologic and tectonic segments within the survey area. The coherent pattern in magnetic data and mesa topography suggests a subglacial extent of the Transantarctic Mountains 400–500 km inland the last exposed rock outcrops. We estimate the maximum thickness of a potential sediment infill in the Wilkes Subglacial Basin to be less than 1 km, based on gravity modeling and source depth estimates from magnetic data. The coherent nature of the potential field and topography data, together with the northwest–southeast trends, define the Adventure Subglacial Trench and the Resolution Subglacial Highlands as a tectonic unit. The crustal structure and the strong similarity of the observed gravity with fold-and-thrust belts suggest a compressional scenario for the origin of the Adventure Subglacial Trench and the Resolution Subglacial Highlands. The complexity and apparent structural control of the Wilkes Subglacial Basin raise the issue of what influence pre-existing structures may have played in the formation of the Transantarctic Mountains system. The previous hypothesis of a thermal boundary beneath the mountains is difficult to reconcile with our new gravity data. The apparent difficulties to match our new data with certain key aspects of previous models suggests that a reassessment of the existing uplift models is necessary. We have modeled the prominent gravity anomaly over the Transantarctic Mountains with thicker crust.     

Stratified airflow over one or two hills

Methods for calculating, interpolating and idealising air flow in complex terrain are reviewed. Then the general structure of stratified airflow over a single hill of height H and length L₁ is studied in detail and shown to be determined by the upwind velocity profile, the magnitude of a characteristic Froude number and the dimensions of the hill. Let N(L₁) be the buoyancy frequency upwind at a height L₁, and u_* and U₀ be the upwind friction and mean velocity respectively; then the flow is effectively neutral if u^*/NL₁>1. But if u_*/NL₁>1 and u₀/NL₁>1, the stratification is weak enough to affect the upwind turbulence and velocity profile but not the dynamics of the flow over the hill. If U₀/NL₁ <1 but U₀/NH>1 the buoyancy forces are strong enough to affect the mean flow over the hill but not strong enough to prevent it passing over the top. In this regime the flow is very sensitive to the form of the upwind temperature profile. If U₀/NH>1, much of the flow passes round the hill. A similar classification, with different flow patterns, is appropriate for unstably stratified flows. When the wind is weak enough, local slope winds can dominate. Results from the analysis of these different regimes are described and compared with laboratory of field measurements where possible. It is shown how some of these results can be extended to groups of hills.     

To what extent have laterites contributed to the geochemical,surface reflectance and magnetic properties of adjacent tropical soils? Evidence from Niger and Burkina Faso

The present study is based on a suite of surface samples from exposures of eroded laterite, considered to be Tertiary in age, and nearby soils in the Sahelian region of SW Niger and Burkina Faso. X‐ray fluorescence, X‐ray diffraction, computer‐controlled scanning electron microscopy, energy‐dispersive X‐ray spectroscopy, diffuse UV–visible reflectance spectroscopy and a suite of magnetic measurements have been used to shed light on the origin of the soils and their possible derivation from the adjacent, eroded laterite outcrops. On the basis of the wide range of data obtained, we conclude that the mineralogy and magnetic properties of the soils preclude direct derivation from the laterites without further weathering and modification. Nor does the evidence support the view that the soils have evolved entirely independently, uninfluenced by input from the laterites. The only conclusion that is consistent with all the lines of evidence is that the erosion of the laterites provided at least a significant part of the material upon which soil formation took place. This must have occurred at a time early enough to permit a long period of subsequent soil development during which the iron oxides, specifically haematite and ferrimagnetic minerals, were significantly modified. From this, we infer that eroded material from discontinuous laterite exposures has contributed significantly to the remotely sensed, distinctive reflectance characteristics of the Sahel surfaces. The magnetic properties of the soils provide evidence for the in situ neo‐formation of fine, secondary, pedogenic magnetite/maghemite grains typical of those found in many soils across the Sahel region and elsewhere in both temperate and tropical environments. Copyright © 2017 John Wiley & Sons, Ltd.     

An examination of the steady‐state assumption in soil development models with application to landscape evolution

Soil depth and soil production are highly complicated phenomena, generated from a complex interaction of physical, biological and chemical processes. It has, nevertheless, become increasingly clear that soil formation rates are closely related to chemical weathering rates. Somewhat paradoxically, it is likewise becoming apparent that such biogeochemical reactions as slowly transform rock to soil are limited by physical processes, such as flowing water and the formation of fractures. We have formulated a theoretical approach that relates soil formation rates to chemical weathering rates, and those, likewise, to solute transport rates. For such a theoretical framework to be relevant, the solute transport rates cannot equal those of the flowing water, as is the case in Gaussian solute transport. Rather, solute transport must be slowed in accordance with heavy‐tailed solute arrival time distributions. The inference is that the traditional advection–dispersion equation formulation for solute transport is inadequate in the typically heterogeneous geological media that weather to form soils. Here we examine the implications of this soil production model on the assumption of the approach to steady state. Particularly at slow erosion rates we find that many soil columns are not in equilibrium. This tendency may be accentuated in dry climates. Copyright © 2017 John Wiley & Sons, Ltd.     

Are Models Too Simple? Arguments for Increased Parameterization

The idea that models should be as simple as possible is often accepted without question. However, too much simplification and parsimony may degrade a model's utility. Models are often constructed to make predictions; yet, they are commonly parameterized with a focus on calibration, regardless of whether (1) the calibration data can constrain simulated predictions or (2) the number and type of calibration parameters are commensurate with the hydraulic property details on which key predictions may depend. Parameterization estimated through the calibration process is commonly limited by the necessity that the number of calibration parameters be smaller than the number of observations. This limitation largely stems from historical restrictions in calibration and computing capability; we argue here that better methods and computing capabilities are now available and should become more widely used. To make this case, two approaches to model calibration are contrasted: (1) a traditional approach based on a small number of homogeneous parameter zones defined by the modeler a priori and (2) regularized inversion, which includes many more parameters than the traditional approach. We discuss some advantages of regularized inversion, focusing on the increased insight that can be gained from calibration data. We present these issues using reasoning that we believe has a common sense appeal to modelers; knowledge of mathematics is not required to follow our arguments. We present equations in an Appendix, however, to illustrate the fundamental differences between traditional model calibration and a regularized inversion approach.     

Ground water discharge and nitrate flux to the Gulf of Mexico

Ground water samples (37 to 186 m depth) from Baldwin County, Alabama, are used to define the hydrogeology of Gulf coastal aquifers and calculate the subsurface discharge of nutrients to the Gulf of Mexico. The ground water flow and nitrate flux have been determined by linking ground water concentrations to 3H/3He and 4He age dates. The middle aquifer (A2) is an active flow system characterized by postnuclear tritium levels, moderate vertical velocities, and high nitrate concentrations. Ground water discharge could be an unaccounted source for nutrients in the coastal oceans. The aquifers annually discharge 1.1 +/- 0.01 x 10(8) moles of nitrate to the Gulf of Mexico, or 50% and 0.8% of the annual contributions from the Mobile-Alabama River System and the Mississippi River System, respectively. In southern Baldwin County, south of Loxley, increasing reliance on ground water in the deeper A3 aquifer requires accurate estimates of safe ground water withdrawal. This aquifer, partially confined by Pliocene clay above and Pensacola Clay below, is tritium dead and contains elevated 4He concentrations with no nitrate and estimated ground water ages from 100 to 7000 years. The isotopic composition and concentration of natural gas diffusing from the Pensacola Clay into the A3 aquifer aids in defining the deep ground water discharge. The highest 4He and CH4 concentrations are found only in the deepest sample (Gulf State Park), indicating that ground water flow into the Gulf of Mexico suppresses the natural gas plume. Using the shape of the CH4-He plume and the accumulation of 4He rate (2.2 +/- 0.8 microcc/kg/1000 years), we estimate the natural submarine discharge and the replenishment rate for the A3 aquifer.     

Subsurface Injection of Treated Sewage into a Saline-Water Aquifer at St. Petersburg, Florida — Water-Quality Changes and Potential for Recovery of Injected Sewage

The city of St. Petersburg is testing subsurface injection of treated sewage into the Floridan aquifer as a means of eliminating discharge of sewage to surface waters and as a means of storing treated sewage for future nonpotable reuse. The injection zone at the test site at the start of injection contained saline water with chloride concentrations ranging from 14,000 to 20,000 milligrams per liter (mg/1). Treated sewage with a mean chloride concentration of 170 mg/1 was injected through a single well for 12 months at a mean rate of 4.7 × 10⁵ cubic feet per day. The volume of water injected during the year was 1.7 × 10⁸ cubic feet. Dissolved oxygen was contained in the sewage prior to injection. Water removed from the injection zone during injection was essentially free of oxygen. Probable growth of denitrifying bacteria and, thus, microbial denitri-fication, was suggested by bacterial counts in water from two observation wells that were close to the injection well. The volume fraction of treated sewage in water from wells located 35 feet and 733 feet from the injection well and open to the upper part of the injection zone stabilized at about 0.9 and 0.75, respectively. Chloride concentrations stabilized at about 1,900 mg/1 in water from the well that was 35 feet from the injection well and stabilized at about 4,000 mg/1 in water from the well that was 733 feet from the injection well. These and other data suggest that very little near injection-quality treated sewage would be recoverable from storage in the injection zone.     

Citywide distribution of lead and other element in soils and indoor dusts in Syracuse,NY

Associations between Pb, Zn, Fe and Mn in soils and indoor dusts in urban Syracuse, NY have been investigated at different levels of spatial aggregation. The concentrations of these elements of interest (EOI) in 3566 soil samples were mapped across the city to investigate variations in concentration potentially associated with specific city locations. Indoor floor dust loadings for the EOI determined at 433 residences were mapped in a similar fashion. Pearson product correlation coefficients at different levels of soil sample aggregation (individual sampling points, block group averages, and census tract (CT) averages), consistently demonstrated a strong correlation between soil Pb and Zn concentrations. Correlations between Pb and both Fe and Mn soil concentrations were also significant. However, the correlation between the Fe and Mn floor dust loadings was much stronger than that for soil, as were the correlations between the Pb and both the Fe and Mn dust loadings. The correlation between the Fe and Zn floor dust loadings was far less significant. Surprisingly, most of the correlations between the paired EOI in the dusts and soils at the individual homes were mostly not statistically significant. The two correlations of any significance, and these were between the Pb in soil and the Pb in the dust, and between Pb in dust and Zn in soil.