The ability of detainment bunds to decrease sediments transported from pastoral catchments in surface runoff

Erosion leading to sedimentation in surface water may disrupt aquatic habitats and deliver sediment-bound nutrients that contribute to eutrophication. Land use changes causing loss of native vegetation have accelerated already naturally high erosion rates in New Zealand and increased sedimentation in streams and lakes. Sediment-bound phosphorus (P) makes up 71–79% of the 17–19 t P y⁻¹ delivered from anthropogenic sources to Lake Rotorua in New Zealand. Detainment bunds (DBs) were first implemented in the Lake Rotorua catchment in 2010 as a strategy to address P losses from pastoral agriculture. The bunds are 1.5–2 m high earthen stormwater retention structures constructed across the flow path of targeted low-order ephemeral streams with the purpose of temporarily ponding runoff on productive pastures. The current DB design protocol recommends a minimum pond volume of 120 m³ ha⁻¹ of contributing catchment with a maximum pond storage capacity of 10 000 m³. No previous study has investigated the ability of DBs to decrease annual suspended sediment (SS) loads leaving pastoral catchments. Annual SS yields delivered to two DBs with 20 ha and 55 ha catchments were 109 and 28 kg SS ha⁻¹, respectively, during this 12-month study. The DBs retained 1280 kg (59%) and 789 kg (51%) of annual SS loads delivered from the catchments as a result of the bunds' ability to impede stormflow and facilitate soil infiltration and sediment deposition. The results of this study highlight the ability of DBs to decrease SS loads transported from pastures in surface runoff, even during large storm events, and suggests DBs are able to reduce P loading in Lake Rotorua.     

gLucifer: next generation visualization framework for high-performance computational geodynamics

High-performance computing provides unprecedented capabilities to produce higher resolution 4-D models in a fraction of time. Thus, the need exists for a new generation of visualization systems able to maintain parity with the enormous volume of data generated. In attempting to write this much data to disk, each computational step introduces a significant performance bottleneck, yet most existing visualization software packages inherently rely on reading data in from a dump file. Available packages make this assumption of postprocessing at quite a fundamental level and are not very well suited for plotting very large numbers of specialized particles. This necessitates the creation of a new visualization system that meets the needs of large-scale geodynamic modeling. We have developed such a system, gLucifer, using a software framework approach that allows efficient reuse of our efforts in other areas of research. gLucifer is capable of producing movies of a 4-D data set “on the fly” (simultaneously with running the parallel scientific application) without creating a performance bottleneck. By eliminating most of the human efforts involved in visualizing results through postprocessing, gLucifer reconnects the scientist to the numerical experiment as it unfolds. Data sets that were previously very difficult to even manage may be efficiently explored and interrogated without writing to disk, and because this approach is based entirely on memory distributed across as many processors as are being utilized by the scientific application, the visualization solution is scalable into terabytes of data being rendered in real time.     

Interactive visualization of 3D mantle convection

The current availability of thousands of processors at many high performance computing centers has made it feasible to carry out, in near real time, interactive visualization of 3D mantle convection temperature fields, using grid configurations having 10–100 million unknowns. We will describe the technical details involved in carrying out this endeavor, using the facilities available at the Laboratory of Computational Science and Engineering (LCSE) at the University of Minnesota. These technical details involve the modification of a parallel mantle convection program, ACuTEMan; the usage of client–server socket based programs to transfer upwards of a terabyte of time series scientific model data using a local network; a rendering system containing multiple nodes; a high resolution PowerWall display, and the interactive visualization software, DSCVR. We have found that working in an interactive visualizastion mode allows for fast and efficient analysis of mantle convection results. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Increasing the usability and accessibility of geodynamic modelling tools to the geoscience community: UnderworldGUI

Geoscientists are faced with a number of complexities that represent obstacles to the development of realistic simulation of deep earth processes. Realistic 4D thermo-mechanical simulation using software packages like Underworld and Gale, when combined appropriately with geoscientific expertise, can lead to novel insights into the deformation of geological structures at a wide range of time and spatial scales. The challenge for end-user geoscientists lies in applying their knowledge within the framework of the software’s input specification, including initial, internal, and boundary conditions and output visualization parameters. We have built a Graphical User Interface (GUI) to remove many of the difficulties related to editing the Extensible Markup Language (XML) encoded input files of Underworld/Gale geomodels and therefore, to greatly broaden the user base of these software packages. By helping Underworld/Gale to meet a large audience, we provide a tool to the geoscience community that helps to move from untested conceptual models to physically valid, properly scaled modelling. Furthermore, the UnderworldGUI offers a mechanism for storing and retrieving experimental models in a centralised database, thus providing the geoscience community with a means to share the outcomes of its experimental research. Further details of the UnderworldGUI are available at the web site .

Regional magma plumbing and emplacement mechanisms of the Faroe‐Shetland Sill Complex: implications for magma transport and petroleum systems within sedimentary basins

The movement of magma through the shallow crust and the impact of subsurface sill complexes on the hydrocarbon systems of prospective sedimentary basins has long been an area of interest and debate. Based on 3D seismic reflection and well data, we present a regional analysis of the emplacement and magmatic plumbing system of the Palaeogene Faroe‐Shetland Sill Complex (FSSC), which is intruded into the Mesozoic and Cenozoic sequences of the Faroe‐Shetland Basin (FSB). Identification of magma flow directions through detailed seismic interpretation of approximately 100 sills indicates that the main magma input zones into the FSB were controlled primarily by the NE–SW basin structure that compartmentalise the FSB into its constituent sub‐basins. An analysis of well data shows that potentially up to 88% of sills in the FSSC are <40 m in thickness, and thus below the vertical resolution limit of seismic data at depths at which most sills occur. This resolution limitation suggests that caution needs to be exercised when interpreting magmatic systems from seismic data alone, as a large amount of intrusive material could potentially be missed. The interaction of the FSSC with the petroleum systems of the FSB is not well understood. Given the close association between the FSSC and potential petroleum migration routes into some of the oil/gas fields (e.g. Tormore), the role the intrusions may have played in compartmentalisation of basin fill needs to be taken fully into account to further unlock the future petroleum potential of the FSB.     

An investigation of the multi-scale temporal variability of beach profiles at Duck using wavelet packet transforms

In an earlier paper a particular discrete wavelet transform (DWT) was used to study the complex variation of beach profile changes. However, use of the DWT requires that the sequence of spatial and temporal resolution is fixed as a dyadic sequence, which means that the variability over longer intervals is not characterised well. Here we introduce the discrete wavelet packet transform (DWPT) that uses an adaptive scaling to partition the data variance, according to an entropy cost function. The advantages of this approach are demonstrated by its application to the study of temporal variability of a 22 year record of beach profile data from the Field Research Facility (FRF) at Duck, North Carolina, USA. Time series of beach elevations at three locations across a particular profile are investigated in detail. We conclude that the DWPT provides a superior analysis of non-stationary time series to that of the DWT, with improved resolution of the scale intervals of the variability. The beach elevation around the shoreline is shown to respond at both sub-annual and interannual scales, but variability at an annual scale is weak. Moving seaward into deeper water, the variance is partitioned into fewer and longer scales. It is confirmed that elevation changes around the inner bar at Duck exhibit a strong interannual variation consistent with Plant et al. (Plant, N.G., Holman, R.A. and Freilich, M.H., 1999. A simple model for interannual sandbar behaviour. Journal of Geophysical Research 104(C7), 15755–15776). Around 23% of the variance around the inner bar is explained at the temporal scale of 64–128 months, which is consistent with the bar behaviour of 6 years found by Ruessink et al. (Ruessink, B. G., Wijnberg, K. M., Holman, R. A., Kuriyama, Y. and Van Enckevort, I. M. J., 2003. Intersite comparison of interannual nearshore bar behaviour. Journal of Geophysical Research, 108 (C8): 1–12). A significant new finding is, however, that about 26% of the variance is attributable to temporal scales of 16–21.3 months. Reconstruction of the wavelet packet components for individual temporal scales is shown to provide a means for identifying the impact and scale of non-stationary events, such as storms, on the beach response. This provides further information that can be used to interpret the morphological changes in terms of the forcing processes and also serves to inform morphodynamic modelling.     

Nutrient contributions to the Santa Barbara Channel, California, from the ephemeral Santa Clara River

The Santa Clara River delivers nutrient rich runoff to the eastern Santa Barbara Channel during brief (1–3 day) episodic events. Using both river and oceanographic measurements, we evaluate river loading and dispersal of dissolved macronutrients (silicate, inorganic N and P) and comment on the biological implications of these nutrient contributions. Both river and ocean observations suggest that river nutrient concentrations are inversely related to river flow rates. Land use is suggested to influence these concentrations, since runoff from a subwatershed with substantial agriculture and urban areas had much higher nitrate than runoff from a wooded subwatershed. During runoff events, river nutrients were observed to conservatively mix into the buoyant, surface plume immediately seaward of the Santa Clara River mouth. Dispersal of these river nutrients extended 10s of km into the channel. Growth of phytoplankton and nutrient uptake was low during our observations (1–3 days following runoff), presumably due to the very low light levels resulting from high turbidity. However, nutrient quality of runoff (Si:N:P = 16:5:1) was found to be significantly different than upwelling inputs (13:10:1), which may influence different algal responses once sediments settle. Evaluation of total river nitrate loads suggests that most of the annual river nutrient fluxes to the ocean occur during the brief winter flooding events. Wet winters (such as El Niño) contribute nutrients at rates approximately an order-of-magnitude greater than “average” winters. Although total river nitrate delivery is considerably less than that supplied by upwelling, the timing and location of these types of events are very different, with river discharge (upwelling) occurring predominantly in the winter (summer) and in the eastern (western) channel.     

Results from microseismic monitoring,conventional instrumentation,and tomography surveys in the creation and thinning of a burst-prone sill pillar

Located in northern Québec, the Lac Shortt Mine was a small gold mine consisting of a thin subvertical orebody which was mined in three main phases. High stress and rockbursting conditions were experienced when ore was extracted in the upper zone between the surface and a depth of 500 metres during the first two phases of mining. Severe rockbursts were experienced in late 1989 near the shaft and in the footwall development following a deepening of the mine shaft to a depth of 830 m and partial development of footwall drift access for the third phase of mining (the mining of the lower zone starting at a depth of 830 m moving upward toward a depth of 500 m). A 16-channel Electrolab MP250 microseismic system, with a Queen's University Full-Waveform piggy-back system, was installed underground at the site due to these problems.It was expected that the thinning sill would be subjected to an ever-increasing load as the thickness of the 500 m sill pillar decreased in the face of the mining excavation from below. A monitoring program consisting of the microseismic monitoring system, a range of conventional geomechanics monitoring tools as well as the undertaking of periodic seismic tomography surveys to assess the ongoing state of stress and rock mass condition within the sill was therefore warranted.The anomalously high-magnitude stress field and the brittle rockmass created a situation in which rockmass failure was common and violent. In the creation and thinning of the sill pillar, the location of banded microseismic activity was crucial in tracing rockmass failure and the associated ground control problems. Reliable source-location determination enabled the identification of areas of stress increase. The movement of the rockmass failure front could be followed, and was responsible for stope dilution, footwall and orebody development deterioration, and caving.Source-mechanism analyses gave accurate double-couple solutions for approximately forty percent of these events having at least ten recognizable polarities. Results suggested movement along vertical north-south striking or vertical east-west striking features. Underground observation of damaged access points showed that vertical north-south striking joints were experiencing failure.The microseismic activity, which was consistently concentrated close to the southwest and northeast corners of current production stopes, could be explained by a stress field oriented obliquely to the strike of the orebody, as measured prior to shrinkage of the sill pillar byin situ stress measurements and observed borehole overbreaks. The orientations of theP andT axes for the microseismic activity further confirmed that the stress field oriented obliquely to strike.While an increase in compressional-wave velocity of 2.3 percent, corresponding to a measured stress increase of approximately 10 MPa could be measured by repeated tomographic surveys, it was relatively small and only a factor of two or so above the velocity measured uncertainty. The relative insensitivity of thein situ rock mass modulus to the applied stress is believed to be largely due to the rockmass discontinuities being relatively closed prior to stress increase, as substantiated by the small deformations seen by the extensometer and borehole camera. This situation existed because of the very high pre-mining stress level.The experimental demonstration that the rock could not absorb substantially increased load through the mechanism of discontinuity closure or tightening (which would be reflected in the modulus) may be evidence in itself of potentially burst-prone ground, such as encountered at Lac Shortt.     

Malham Cove: splendour and enigma

A dry waterfall that recently came to life, a landform that is not karstic in the middle of a great karst terrain, underground drainage that is not as simple as it seems, and on‐going debate about how it was formed: Malham Cove is more than just a spectacular sight.     

Lötschberg–geology for tunnellers

Switzerland's Lötschberg tunnels present a long story of both failure and success. Tunnelling that began 100 years ago ended in a massive disaster–which is a classic of engineering geology–but subsequent works were very much more successful.