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.     

Geophysical models for the tectonic framework of the Lake Vostok region, East Antarctica

Aerogeophysical and seismological data from a geophysical survey in the interior of East Antarctica were used to develop a conceptual tectonic model for the Lake Vostok region. The model is constrained using three independent data sets: magnetic, seismic, and gravimetric. A distinct change in the aeromagnetic anomaly character across Lake Vostok defines a crustal boundary. Depth to magnetic basement estimates image a 400-km-wide and more than 10-km-deep sedimentary basin west of the lake. Analysis of teleseismic earthquakes suggests a relatively thin crust beneath Lake Vostok consistent with predictions from kinematic and flexural gravity modelling. Magnetic, gravity, and subglacial topography data reveal a tectonic boundary within East Antarctica. Based on our kinematic and flexural gravity modelling, this tectonic boundary appears to be the result of thrust sheet emplacement onto an earlier passive continental margin. No data presently exist to date directly either the timing of passive margin formation or the subsequent shortening phase. The preserved thrust sheet thickness is related to the thickness of the passive margin crust. Because a significant amount of time is required to erode the thrust sheet topography, we suggest that these tectonic events are Proterozoic in age. Minor normal reactivation of the thrust sheet offers a simple mechanism to explain the formation of the Lake Vostok Basin. A low level of seismicity exists in the vicinity of this tectonic boundary. The existence of a crustal boundary in the Antarctic interior provides new constraints on the Proterozoic architecture of the East Antarctic craton.     

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.     

Performance characteristics of a new electronic snow water equivalent sensor in different climates

The US Army ERDC CRREL and the US Department of Agriculture Natural Resources Conservation Service developed a square electronic snow water equivalent (e‐SWE) sensor as an alternative to using fluid‐filled snow pillows to measure SWE. The sensors consist of a centre panel to measure SWE and eight outer panels to buffer edge stress concentrations. Seven 3 m square e‐SWE sensors were installed in five different climate zones. During the 2011–2012 winter, 1.8 and 1.2 m square e‐SWE sensors were installed and operated in Oregon. With the exception of New York State and Newfoundland, the e‐SWE sensors accurately measured SWE, with R² values between the sensor and manual SWE measurements of between 0.86 and 0.98. The e‐SWE sensor at Hogg Pass, Oregon, accurately measured SWE during the past 8 years of operations. In the thin, icy snow of New York during midwinter 2008–2009, the e‐SWE sensors overmeasured SWE because of edge stress concentrations associated with strong icy layers and a shallow snow cover. The New York e‐SWE sensors' measurement accuracy improved in spring 2009 and further improved during the 2011–2012 winter with operating experience. At Santiam Junction, measured SWE from the 1.8 and 1.2 m square e‐SWE sensors agreed well with the snow pillow, 3 m square e‐SWE sensor, and manual SWE measurements until February 2013, when dust and gravel blew onto the testing area resulting in anomalous measurements. © 2014 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Carbon and oxygen isotope constraints on fluid sources and fluid–wallrock interaction in regional alteration and iron-oxide–copper–gold mineralisation, eastern Mt Isa Block, Australia

The source of metasomatic fluids in iron-oxide–copper–gold districts is contentious with models for magmatic and other fluid sources having been proposed. For this study, δ ¹⁸O and δ ¹³C ratios were measured from carbonate mineral separates in the Proterozoic eastern Mt Isa Block of Northwest Queensland, Australia. Isotopic analyses are supported by petrography, mineral chemistry and cathodoluminescence imagery. Marine meta-carbonate rocks (ca. 20.5‰ δ ¹⁸O and 0.5‰ δ ¹³C calcite) and graphitic meta-sedimentary rocks (ca. 14‰ δ ¹⁸O and −18‰ δ ¹³C calcite) are the main supracrustal reservoirs of carbon and oxygen in the district. The isotopic ratios for calcite from the cores of Na–(Ca) alteration systems strongly cluster around 11‰ δ ¹⁸O and −7‰ δ ¹³C, with shifts towards higher δ ¹⁸O values and higher and lower δ ¹³C values, reflecting interaction with different hostrocks. Na–(Ca)-rich assemblages are out of isotopic equilibrium with their metamorphic hostrocks, and isotopic values are consistent with fluids derived from or equilibrated with igneous rocks. However, igneous rocks in the eastern Mt Isa Block contain negligible carbon and are incapable of buffering the δ ¹³C signatures of CO₂-rich metasomatic fluids associated with Na–(Ca) alteration. In contrast, plutons in the eastern Mt Isa Block have been documented as having exsolved saline CO₂-rich fluids and represent the most probable fluid source for Na–(Ca) alteration. Intrusion-proximal, skarn-like Cu–Au orebodies that lack significant K and Fe enrichment (e.g. Mt Elliott) display isotopic ratios that cluster around values of 11‰ δ ¹⁸O and −7‰ δ ¹³C (calcite), indicating an isotopically similar fluid source as for Na–(Ca) alteration and that significant fluid–wallrock interaction was not required in the genesis of these deposits. In contrast, K- and Fe-rich, intrusion-distal deposits (e.g. Ernest Henry) record significant shifts in δ ¹⁸O and δ ¹³C towards values characteristic of the broader hostrocks to the deposits, reflecting fluid–wallrock equilibration before mineralisation. Low temperature, low salinity, low δ ¹⁸O (<10‰ calcite) and CO₂-poor fluids are documented in retrograde metasomatic assemblages, but these fluids are paragenetically late and have not contributed significantly to the mass budgets of Cu–Au mineralisation.     

The July 2007 rock and ice avalanches at Mount Steele, St. Elias Mountains, Yukon, Canada

A large rock and ice avalanche occurred on the north face of Mount Steele, southwest Yukon Territory, Canada, on July 24, 2007. In the days and weeks preceding the landslide, several smaller avalanches initiated from the same slope. The ice and rock debris traveled a maximum horizontal distance 5.76 km with a maximum vertical descent of 2,160 m, leaving a deposit 3.66 km² in area on Steele Glacier. The seismic magnitude estimated from long-period surface waves (M _s) is 5.2. Modeling of the waveforms suggests an estimated duration of approximately 100 s and an average velocity of between 35 and 65 m/s. This landslide is one of 18 large rock avalanches known to have occurred since 1899 on slopes adjacent to glaciers in western Canada. We describe the setting, reconstruct the event chronology and present a preliminary characterization of the Mount Steele ice and rock avalanches based on field reconnaissance, analysis of seismic records and an airborne LiDAR survey. We also present the results of a successful dynamic simulation for the July 24 event.     

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.