Climate change and future flows of Rocky Mountain rivers: converging forecasts from empirical trend projection and down‐scaled global circulation modelling

To predict future river flows, empirical trend projection (ETP) analyses and extends historic trends, while hydroclimatic modelling (HCM) incorporates regional downscaling from global circulation model (GCM) outputs. We applied both approaches to the extensively allocated Oldman River Basin that drains the North American Rocky Mountains and provides an international focus for water sharing. For ETP, we analysed monthly discharges from 1912 to 2008 with non‐parametric regression, and extrapolated changes to 2055. For modelling, we refined the physical models MTCLIM and SNOPAC to provide water inputs into RIVRQ (river discharge), a model that assesses the streamflow regime as involving dynamic peaks superimposed on stable baseflow. After parameterization with 1960–1989 data, we assessed climate forecasts from six GCMs: CGCM1‐A, HadCM3, NCAR‐CCM3, ECHAM4 and 5 and GCM2. Modelling reasonably reconstructed monthly hydrographs (R² about 0·7), and averaging over three decades closely reconstructed the monthly pattern (R² = 0·94). When applied to the GCM forecasts, the model predicted that summer flows would decline considerably, while winter and early spring flows would increase, producing a slight decline in the annual discharge (?3%, 2005–2055). The ETP predicted similarly decreased summer flows but slight change in winter flows and greater annual flow reduction (?9%). The partial convergence of the seasonal flow projections increases confidence in a composite analysis and we thus predict further declines in summer (about ? 15%) and annual flows (about ? 5%). This composite projection indicates a more modest change than had been anticipated based on earlier GCM analyses or trend projections that considered only three or four decades. For other river basins, we recommend the utilization of ETP based on the longest available streamflow records, and HCM with multiple GCMs. The degree of correspondence from these two independent approaches would provide a basis for assessing the confidence in projections for future river flows and surface water supplies. Copyright © 2010 John Wiley & Sons, Ltd.     

The data‐driven approach as an operational real‐time flood forecasting model

Accurate water level forecasts are essential for flood warning. This study adopts a data‐driven approach based on the adaptive network–based fuzzy inference system (ANFIS) to forecast the daily water levels of the Lower Mekong River at Pakse, Lao People's Democratic Republic. ANFIS is a hybrid system combining fuzzy inference system and artificial neural networks. Five ANFIS models were developed to provide water level forecasts from 1 to 5 days ahead, respectively. The results show that although ANFIS forecasts of water levels up to three lead days satisfied the benchmark, four‐ and five‐lead‐day forecasts were only slightly better in performance compared with the currently adopted operational model. This limitation is imposed by the auto‐ and cross‐correlations of the water level time series. Output updating procedures based on the autoregressive (AR) and recursive AR (RAR) models were used to enhance ANFIS model outputs. The RAR model performed better than the AR model. In addition, a partial recursive procedure that reduced the number of recursive steps when applying the AR or the RAR model for multi‐step‐ahead error prediction was superior to the fully recursive procedure. The RAR‐based partial recursive updating procedure significantly improved three‐, four‐ and five‐lead‐day forecasts. Our study further shows that for long lead times, ANFIS model errors are dominated by lag time errors. Although the ANFIS model with the RAR‐based partial recursive updating procedure provided the best results, this method was able to reduce the lag time errors significantly for the falling limbs only. Improvements for the rising limbs were modest. Copyright © 2011 John Wiley & Sons, Ltd.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Knowing wildfire risk: Scientific interactions with risk mitigation policy and practice in Victoria,Australia

Over the past decade, major landscape wildfires (or ‘bushfires’ in Australia) in fire-prone countries have illustrated the seriousness of this global environmental problem. This natural hazard presents a complex mesh of dynamic factors for those seeking to reduce or manage its costs, as ignitions, hazard behaviour, and the reactions of different human and ecological communities during and after hazard events are all extremely uncertain. But while those at risk of wildfire have been subject to significant research, the social dimensions of its management, including the role of science, have received little attention. This paper reports on a case study of the Barwon-Otway area of Victoria in Australia, a high wildfire risk area that has recently been a pilot site for a new risk mitigation strategy utilising the wildfire simulation model PHOENIX RapidFire. Against simple equations between ‘more science’ and ‘less uncertainty,’ this paper presents results from interviews and a workshop with practitioners to investigate how scientific research interacts with and informs both wildfire policy and practice. We suggest that attending to cultural and social specificities of the application of any technical innovation—such as next generation modelling—raises questions for future research about the roles of narrative, performance, and other knowledges in the sedimentation of science.     

The topology of geology 2: Topological uncertainty

Uncertainty is ubiquitous in geology, and efforts to characterise and communicate it are becoming increasingly important. Recent studies have quantified differences between perturbed geological models to gain insight into uncertainty. We build on this approach by quantifying differences in topology, a property that describes geological relationships in a model, introducing the concept of topological uncertainty. Data defining implicit geological models were perturbed to simulate data uncertainties, and the amount of topological variation in the resulting model suite measured to provide probabilistic assessments of specific topological hypotheses, sources of topological uncertainty and the classification of possible model realisations based on their topology. Overall, topology was found to be highly sensitive to small variations in model construction parameters in realistic models, with almost all of the several thousand realisations defining distinct topologies. In particular, uncertainty related to faults and unconformities was found to have profound topological implications. Finally, possible uses of topology as a geodiversity metric and validation filter are discussed, and methods of incorporating topological uncertainty into physical models are suggested.     

Copper demand,supply, and associated energy use to 2050

To a set of well-regarded international scenarios (UNEP’s GEO-4), we have added consideration of the demand, supply, and energy implications related to copper production and use over the period 2010–2050. To our knowledge, these are the first comprehensive metal supply and demand scenarios to be developed. We find that copper demand increases by between 275 and 350% by 2050, depending on the scenario. The scenario with the highest prospective demand is not Market First (a “business as usual” vision), but Equitability First, a scenario of transition to a world of more equitable values and institutions. These copper demands exceed projected copper mineral resources by mid-century and thereafter. Energy demand for copper production also demonstrates strong increases, rising to as much as 2.4% of projected 2050 overall global energy demand. We investigate possible policy responses to these results, concluding that improving the efficiency of the copper cycle and encouraging the development of copper-free energy distribution on the demand side, and improving copper recycling rates on the supply side are the most promising of the possible options. Improving energy efficiency in primary copper production would lead to a reduction in the energy demand by 0.5% of projected 2050 overall global energy demand. In addition, encouraging the shift towards renewable technologies is important to minimize the impacts associated with copper production.     

面向矿床三维动态建模的地质勘探数据库增量更新方法

针对深部找矿过程中地质勘探数据增加和动态三维地质建模需要,提出了一种基于版本管理的矿床地质勘探数据库增量更新方法。在分析矿床地质勘探工程数据与三维地质模型映射关系的基础上,建立时序版本和建模版本作为矿床勘探数据增量更新版本标识。针对矿床地质勘探数据库增量更新的版本管理,提出了一种扩展的有向无环图版本管理模型,设计了基于关系数据库的版本管理方法并探讨了其实现过程。以安徽铜陵凤凰山矿床地质勘探数据为例,基于SQL Server数据库系统建立了实例矿床地质勘探版本数据库,以C#编程实现了实例矿床勘探数据库的版本管理,包括版本的建立、查询与显示等功能,通过勘探工程增量更新版本数据建立了三维矿体的动态修正模型,证实了方法的可行性与有效性。     

Analysis of soil nail walls performance - Case study

This study surveys the performance of temporary soil nail walls with a maximum height of 29.3 m in the Yas project, located in Tehran. Some numerical models, with various modelling approaches, were developed using finite element software and the proper modelling approaches were specified and verified.The modelling results were analysed and compared with the in situ monitored data. The results demonstrate the model’s horizontal deformations are generally greater than the in situ monitored values. To decrease the differences between models and monitored results, the effect of variations of the soil parameters had been surveyed and the limits of needed variations were specified.     

Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

Sensitivity assessment of strontium isotope as indicator of polluted groundwater for hydraulic fracturing flowback fluids produced in the Dameigou Shale of Qaidam Basin

Hydrogeochemical processes that would occur in polluted groundwater and aquifer system, may reduce the sensitivity of Sr isotope being the indicator of hydraulic fracturing flowback fluids(HFFF) in groundwater. In this paper, the Dameigou shale gas field in the northern Qaidam Basin was taken as the study area, where the hydrogeochemical processes affecting Sr isotope was analysed. Then, the model for Sr isotope in HFFF-polluted groundwater was constructed to assess the sensitivity of Sr isotope as HFFF indicator. The results show that the dissolution can release little Sr to polluted groundwater and cannot affect the εSr(the deviation of the 87 Sr/86 Sr ratio) of polluted groundwater. In the meantime, cation exchange can considerably affect Sr composition in the polluted groundwater. The Sr with low εSr is constantly released to groundwater from the solid phase of aquifer media by cation exchange with pollution of Quaternary groundwater by the HFFF and it accounts for 4.6% and 11.0% of Sr in polluted groundwater when the HFFF flux reaches 10% and 30% of the polluted groundwater, respectively. However, the Sr from cation exchange has limited impact on Sr isotope in polluted groundwater. Addition of Sr from cation exchange would only cause a 0.2% and 1.2% decrease in εSr of the polluted groundwater when the HFFF flux reaches 10% and 30% of the polluted groundwater, respectively. These results demonstrate that hydrogeochemical processes have little effect on the sensitivity of Sr isotope being the HFFF indicator in groundwater of the study area. For the scenario of groundwater pollution by HFFF, when the HFFF accounts for 5%(in volume percentage) of the polluted groundwater, the HFFF can result in detectable shifts of εSr(ΔεSr=0.86) in natural groundwater. Therefore, after consideration of hydrogeochemical processes occurred in aquifer with input of the HFFF, Sr isotope is still a sensitive indicator of the Quaternary groundwater pollution by the HFFF produced in the Dameigou shale of Qaidam Basin.