Downstream Fining and Sorting of Gravel Clasts in the Braided Rivers of mid-Canterbury, New Zealand

Gravel clast size dimensions have been determined in the Rakaia, Ashburton, and Rangitata rivers by measuring 100 clasts at representative sample locations along each river. In all rivers, gravel size decreases and sorting improves downstream for mean, D₅₀, and D₉₀ fractions of the bed material. Clast size entering the sea is similar in all rivers (30–40 mm b‐axis dimension), despite large variations in transport distance, input size of clasts at their gorges, and discharge. The greatest size reduction occurs in the Rangitata River which has the shortest transport distance and steepest gradient. Rates of downstream clast size reduction are greater than would be assumed from Sternberg's Law, suggesting that additional factors, other than physical abrasion, such as sorting and selective entrainment operate.     

Effect of iron and carbonation on the diffusion of iodine and rhenium in waste encasement concrete and soil fill material under hydraulically unsaturated conditions

Assessing long-term performance of Category 3 cement wasteforms and accurate prediction for radionuclide encasement requires knowledge of the radionuclide–cement interactions and mechanisms of retention (i.e. sorption or precipitation). A set of sediment-concrete half-cell diffusion experiments was conducted under unsaturated conditions (4% and 7% by weight moisture content) using carbonated and non-carbonated concrete–soil half cells. Results indicate the behavior of Re and I release was comparable within a given half-cell test. Diffusivity in soil is a function of moisture content; a 3% increase in moisture content affords a one to two order of magnitude increase in diffusivity. Release of I and Re was 1–3 orders of magnitude less from non-carbonated, relative to carbonated, concrete monoliths. Inclusion of Fe in non-carbonate monoliths resulted in the lowest concrete diffusivity values for both I and Re. This suggests that in the presence of Fe, I and Re are converted to reduced species, which are less soluble and better retained within the concrete monolith. The release of I and Re was greatest from Fe-bearing, carbonated concrete monoliths, suggesting carbonation negates the effect of Fe on the retention of I and Re within concrete monoliths. This is likely due to enhanced formation of microcracks in the presence of Fe, which provide preferential paths for contaminant migration. Although the release of I and Re were greatest from carbonated concrete monoliths containing Fe, the migration of I and Re within a given half cell is dependent on the moisture content, soil diffusivity, and diffusing species.     

Throughfall and its spatial variability in a temperate rainforest of simple structure

There are few throughfall data from southern hemisphere closed-forest, and none from Tasmanian callidendrous cool-temperate rainforest, which has a simpler structure than most primary rainforests. We determined throughfall, measured its local spatial variation, and tested its relationships with rainfall, rainfall intensity, wind speed, canopy dryness, canopy cover, and other structural variables in a cool-temperate callidendrous rainforest in Tasmania. Eighty-two percent of the precipitation was measured as throughfall, which occurred after 2.3 mm of rain fell on dry canopies. The cumulative rainfall in 25 randomly located funnel rain gauges on the forest floor varied from 160 to 567 mm. Canopy cover and other structural variables did not predict the spatial pattern of throughfall. While throughfall in rainfall events was related to rainfall amount and intensity, wind speed did not affect throughfall as a percentage of rainfall. Percentages of throughfall to rainfall over 100 for many low rainfall events may indicate a contribution of fog drip to precipitation on the forest floor. The high local spatial variability in throughfall indicates the mean moisture conditions on the forest floor may not be a good indicator of the potential for localised fire damage.     

Contributions and unrealized potential contributions of cosmogenic-nuclide exposure dating to glacier chronology, 1990–2010

This paper reviews the application of cosmogenic-nuclide exposure dating to glacier chronology. Exposure dating of glacial landforms has made an outsize impact on this field because the technique filled an obvious need that had already been recognized by glacial geologists. By now, hundreds of studies have used cosmogenic-nuclide exposure dating to date glacial deposits, and in fact it is rare to find a study of glacial geology or glacier chronology, or any paleoclimate synthesis that makes use of such studies, that does not involve exposure dating. These developments have resulted in major contributions to glacier chronology and paleoclimate, in particular i) reconstructing Antarctic ice sheet change, ii) establishing the chronology of late Pleistocene and Holocene glacier change in mountain regions where it was previously unknown; iii) establishing the broad chronological outlines of mountain glaciations prior to the Last Glacial Maximum; and iv) gaining insight into subglacial erosional processes through the observation that many glaciated surfaces preserve cosmogenic-nuclide inventories from long past ice-free periods as well as the present one. An important potential future contribution will be the application of the large data set of exposure-dated glacier chronologies to better understand global and regional climate dynamics during Lateglacial and Holocene millennial-scale climate changes. However, this contribution cannot be realized without significant progress in two areas: i) understanding and accounting for geologic processes that cause apparent exposure ages on glacial landforms to differ from the true age of the landform, and ii) minimizing systematic uncertainties in exposure ages that stem from cosmogenic-nuclide production-rate estimates and scaling schemes. At present there exists an enormous data set of exposure ages on glacial deposits, but these data cannot be used to their full potential in paleoclimate syntheses due to an inadequate understanding of geologic scatter and production-rate uncertainties. The intent of this paper is to highlight this situation and suggest some strategies for realizing this potential.     

Exploring the limits of identifying sub-pixel thermal features using ASTER TIR data

Understanding the characteristics of volcanic thermal emissions and how they change with time is important for forecasting and monitoring volcanic activity and potential hazards. Satellite instruments view volcanic thermal features across the globe at various temporal and spatial resolutions. Thermal features that may be a precursor to a major eruption, or indicative of important changes in an on-going eruption can be subtle, making them challenging to reliably identify with satellite instruments. The goal of this study was to explore the limits of the types and magnitudes of thermal anomalies that could be detected using satellite thermal infrared (TIR) data. Specifically, the characterization of sub-pixel thermal features with a wide range of temperatures is considered using ASTER multispectral TIR data. First, theoretical calculations were made to define a “thermal mixing detection threshold” for ASTER, which quantifies the limits of ASTER's ability to resolve sub-pixel thermal mixing over a range of hot target temperatures and % pixel areas. Then, ASTER TIR data were used to model sub-pixel thermal features at the Yellowstone National Park geothermal area (hot spring pools with temperatures from 40 to 90 °C) and at Mount Erebus Volcano, Antarctica (an active lava lake with temperatures from 200 to 800 °C). Finally, various sources of uncertainty in sub-pixel thermal calculations were quantified for these empirical measurements, including pixel resampling, atmospheric correction, and background temperature and emissivity assumptions.     

A Model-Averaging Method for Assessing Groundwater Conceptual Model Uncertainty

This study evaluates alternative groundwater models with different recharge and geologic components at the northern Yucca Flat area of the Death Valley Regional Flow System (DVRFS), USA. Recharge over the DVRFS has been estimated using five methods, and five geological interpretations are available at the northern Yucca Flat area. Combining the recharge and geological components together with additional modeling components that represent other hydrogeological conditions yields a total of 25 groundwater flow models. As all the models are plausible given available data and information, evaluating model uncertainty becomes inevitable. On the other hand, hydraulic parameters (e.g., hydraulic conductivity) are uncertain in each model, giving rise to parametric uncertainty. Propagation of the uncertainty in the models and model parameters through groundwater modeling causes predictive uncertainty in model predictions (e.g., hydraulic head and flow). Parametric uncertainty within each model is assessed using Monte Carlo simulation, and model uncertainty is evaluated using the model averaging method. Two model-averaging techniques (on the basis of information criteria and GLUE) are discussed. This study shows that contribution of model uncertainty to predictive uncertainty is significantly larger than that of parametric uncertainty. For the recharge and geological components, uncertainty in the geological interpretations has more significant effect on model predictions than uncertainty in the recharge estimates. In addition, weighted residuals vary more for the different geological models than for different recharge models. Most of the calibrated observations are not important for discriminating between the alternative models, because their weighted residuals vary only slightly from one model to another.     

Investigation of an early Pleistocene marine osmium isotope record from the eastern equatorial Pacific

Osmium isotope composition (¹⁸⁷Os/¹⁸⁸Os) and concentrations of Os, Ir and Pt are reported for an early Pleistocene section from the ODP Site 849 in the eastern equatorial Pacific. Using the data obtained in this study, the contributions from detrital and extraterrestrial particulate matter to Os concentration and ¹⁸⁷Os/¹⁸⁸Os of sediment are estimated. Our calculations show that detrital contributions to sedimentary Os are too small (<2%) to significantly shift measured bulk sediment ¹⁸⁷Os/¹⁸⁸Os away from seawater values. A moderate but significant negative correlation between ¹⁸⁷Os/¹⁸⁸Os and ³He/¹⁸⁸Os indicate that the average particulate extraterrestrial Os flux to this site is 1.21 ± 0.47 pg cm⁻² kyr⁻¹, which constitutes ?3% of total Os burial flux. The estimates of detrital and extraterrestrial Os are used to calculate the seawater ¹⁸⁷Os/¹⁸⁸Os in the early Pleistocene. The most notable features of this early Pleistocene ¹⁸⁷Os/¹⁸⁸Os record are: (1) glacial-interglacial ¹⁸⁷Os/¹⁸⁸Os differences are insignificant within errors of estimates, (2) glacial ¹⁸⁷Os/¹⁸⁸Os values are higher compared to those reported for the late Pleistocene glacials. Comparison of ¹⁸⁷Os/¹⁸⁸Os values at Site 849 to the late Pleistocene records suggests that average seawater ¹⁸⁷Os/¹⁸⁸Os change has been modest (∼5%) since the early Pleistocene. Assuming that ¹⁸⁷Os/¹⁸⁸Os difference between the glacial periods of the late and the early Pleistocene results solely from temperature dependence of weathering rates, it has been calculated that average surface temperature during the late Pleistocene glacials was 0.8 ± 0.2 °C lower than glacials in the early Pleistocene. This inference is consistent with temperature estimates based on a recent study of pCO₂ reconstruction in the Pleistocene. This observation based on limited studies of marine ¹⁸⁷Os/¹⁸⁸Os records seems to suggest that temperature played an important role in influencing chemical weathering during the Pleistocene glacials. However, more studies are needed to confirm if this temperature-weathering feedback was operational throughout the Pleistocene. A significant down core Ir-³He co-variation coupled with similar burial fluxes of Ir at Site 849 and at LL44 GPC-3 in the north Pacific point to the utility of Ir concentration as a point paleoflux tracer. However, a twofold difference in Ir burial fluxes between the eastern and the western equatorial Pacific suggests that calibration in space and time is required to use Ir concentration as a robust indicator of paleoflux through time. Significant co-variation of concentrations of Os and total alkenone during the glacials coupled with lighter δ¹³C of benthic foraminifera indicates that productivity and carbon burial played a dominant control on scavenging of Os at Site 849. In a broader context, this data set encourages future investigation of response of PGE behavior to paleoceanographic processes.     

Evaluating a new landform evolution model: A case study using a proposed mine rehabilitation landform

Landform evolution models are powerful tools for determining long-term erosional stability and denudation rates spanning geological timescales. SIBERIA, CAESAR and CHILD are examples of these model. The newly developed State Space Soil Production and Assessment Model (SSSPAM) coupled soilscape-landform evolution model has the ability to assess overall erosion rates of catchment scale landforms either using short-term precipitation events, variable precipitation or time-averaged precipitation (annual average). In addition, SSSPAM has the capability of developing the subsurface soil profile through weathering and armouring. In SSSPAM, physical processes of pedogenesis such as erosion and armouring, diffusion, sediment deposition and weathering are modelled using a state space matrix approach. In this article we simulate the short-term evolution (100 years) of a proposed post-mining landform using both SIBERIA and SSSPAM and compare the erosion and sediment output results. For the short-term simulations SSSPAM's armouring capability was disabled. The models were then used to simulate the evolution of the catchment for 10,000 years. Results demonstrate that the short-term SSSPAM simulation results compare well with the results from the established landform evolution model SIBERIA. The long-term armouring disabled SSSPAM simulations produces simulated erosion rates comparable with SIBERIA simulations both of which are similar to upper limit of field measured denudation rates. The SSSPAM simulation using armouring demonstrated that armouring reduced the erosion rate of the catchment by a factor of 4 which is comparable with the lower limit of field measured denudation rates. This observation emphasizes the importance of armouring in long-term evolution of landforms. Soil profile cross-sections developed from the same results show that SSSPAM can also reproduce subsurface soil evolution and stratification and spatial variability of soil profile characteristics typically observed in the field.     

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.