Representing and communicating deep uncertainty in climate-change assessments

IPCC reports provide a synthesis of the state of the science in order to inform the international policy process. This task is made difficult by the presence of deep uncertainty in the climate problem that results from long time scales and complexity. This paper focuses on how deep uncertainty can be effectively communicated. We argue that existing schemes do an inadequate job of communicating deep uncertainty and propose a simple approach that distinguishes between various levels of subjective understanding in a systematic manner. We illustrate our approach with two examples. To cite this article: M. Kandlikar et al., C. R. Geoscience 337 (2005).     

An approach to model the mechanical behavior of transversely isotropic materials

The present work proposes an approach to adapt existing isotropic models to transversely isotropic materials. The main idea is to introduce equivalence relations between the real material and a fictitious isotropic one on which one can take all the advantages of the well‐established isotropic theory. Two applications of this approach are presented here: a failure criterion and a damage model that takes into account the load‐induced anisotropy. In both cases, theoretical predictions are in agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials; nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2015 John Wiley & Sons, Ltd.     

Computational model coupling mode II discrete fracture propagation with continuum damage zone evolution

We propose a numerical method that couples a cohesive zone model (CZM) and a finite element‐based continuum damage mechanics (CDM) model. The CZM represents a mode II macro‐fracture, and CDM finite elements (FE) represent the damage zone of the CZM. The coupled CZM/CDM model can capture the flow of energy that takes place between the bulk material that forms the matrix and the macroscopic fracture surfaces. The CDM model, which does not account for micro‐crack interaction, is calibrated against triaxial compression tests performed on Bakken shale, so as to reproduce the stress/strain curve before the failure peak. Based on a comparison with Kachanov's micro‐mechanical model, we confirm that the critical micro‐crack density value equal to 0.3 reflects the point at which crack interaction cannot be neglected. The CZM is assigned a pure mode II cohesive law that accounts for the dependence of the shear strength and energy release rate on confining pressure. The cohesive shear strength of the CZM is calibrated by calculating the shear stress necessary to reach a CDM damage of 0.3 during a direct shear test. We find that the shear cohesive strength of the CZM depends linearly on the confining pressure. Triaxial compression tests are simulated, in which the shale sample is modeled as an FE CDM continuum that contains a predefined thin cohesive zone representing the idealized shear fracture plane. The shear energy release rate of the CZM is fitted in order to match to the post‐peak stress/strain curves obtained during experimental tests performed on Bakken shale. We find that the energy release rate depends linearly on the shear cohesive strength. We then use the calibrated shale rheology to simulate the propagation of a meter‐scale mode II fracture. Under low confining pressure, the macroscopic crack (CZM) and its damaged zone (CDM) propagate simultaneously (i.e., during the same loading increments). Under high confining pressure, the fracture propagates in slip‐friction, that is, the debonding of the cohesive zone alternates with the propagation of continuum damage. The computational method is applicable to a range of geological injection problems including hydraulic fracturing and fluid storage and should be further enhanced by the addition of mode I and mixed mode (I+II+III) propagation. Copyright © 2016 John Wiley & Sons, Ltd.     

Chemo‐mechanics of cemented granular solids subjected to precipitation and dissolution of mineral species

This paper studies the chemo‐mechanics of cemented granular solids in the context of continuum thermodynamics for fluid‐saturated porous media. For this purpose, an existing constitutive model formulated in the frame of the Breakage Mechanics theory is augmented to cope with reactive processes. Chemical state variables accounting for the reactions between the solid constituents and the solutes in the pore fluid are introduced to enrich the interactions among the microstructural units simulated by the model (i.e., grains and cement bonds). Two different reactive processes are studied (i.e., grain dissolution and cement precipitation), using the chemical variables to describe the progression of the reactions and track changes in the size of grains and bonds. Finally, a homogenization strategy is used to derive the energy potentials of the solid mixture, adopting probability density functions that depend on both mechanical and chemical indices. It is shown that the connection between the statistics of the micro‐scale attributes and the continuum properties of the solid enables the mathematical capture of numerous mechanical effects of lithification and chemical deterioration, such as changes in stiffness, expansion/contraction of the elastic domain, and development of inelastic strains during reaction. In particular, the model offers an interpretation of the plastic strains generated by aggressive environments, which are here interpreted as an outcome of chemically driven debonding and comminution. As a result, the model explains widely observed macroscopic signatures of geomaterial degradation by reconciling the energetics of the deformation/reaction processes with the evolving geometry of the microstructural attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Anisotropic damage of rock joints during cyclic loading: constitutive framework and numerical integration

This work describes a constitutive framework for modeling the behavior of rough joints under cyclic loading. Particular attention is paid to the intrinsic links between dilatancy, surface degradation, and mobilized shear strength. The framework also accounts for the important effect of shear‐induced anisotropy. The resulting approach is fully three‐dimensional and is not restricted to plane‐displacement kinematics. Both the governing formulation and an algorithm for implicit numerical integration are presented. While the proposed methods are general, we also postulate a specific model that is compared with experimental data. It employs relatively few free parameters but shows good agreement with laboratory tests. Copyright © 2013 John Wiley & Sons, Ltd.     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Modelling atmospheric and hydrologic processes for assessment of meadow restoration impact on flow and sediment in a sparsely gauged California watershed

The restoration of meadowland using the pond and plug technique of gully elimination was performed in a 9‐mile segment along Last Chance Creek, Feather River Basin, California, in order to rehabilitate floodplain functions such as mitigating floods, retaining groundwater, and reducing sediment yield associated with bank erosion and to significantly alter the hydrologic regime. However, because the atmospheric and hydrological conditions have evolved over the restoration period, it was difficult to obtain a comprehensible evaluation of the impact of restoration activities by means of field measurements. In this paper, a new use of physically based models for environmental assessment is described. The atmospheric conditions over the sparsely gauged Last Chance Creek watershed (which does not have any precipitation or weather stations) during the combined historical critical dry and wet period (1982–1993) were reconstructed over the whole watershed using the atmospheric fifth‐generation mesoscale model driven with the US National Center for Atmospheric Research and US National Center for Environmental Prediction reanalysis data. Using the downscaled atmospheric data as its input, the watershed environmental hydrology (WEHY) model was applied to this watershed. All physical parameters of the WEHY model were derived from the existing geographic information system and satellite‐driven data sets. By comparing the prerestoration and postrestoration simulation results under the identical atmospheric conditions, a more complete environmental assessment of the restoration project was made. Model results indicate that the flood peak may be reduced by 10–20% during the wet year and the baseflow may be enhanced by 10–20% during the following dry seasons (summer to fall) in the postrestoration condition. The model results also showed that the hydrologic impact of the land management associated with the restoration mitigates bank erosion and sediment discharge during winter storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

An original semi‐discrete approach to assess gas conductivity of concrete structures

For civil engineering structures with a tightness role, structural permeability is a key issue. In this context, this paper presents a new proposition of a numerical modelling of leakage rate through a cracked concrete structure undergoing mode I cracking. The mechanical state of the material, considered in the framework of continuum mechanics based on finite element modelling, is described by means of the stress‐based nonlocal damage model which takes into account the stress state and provides realistic local mechanical fields. A semi‐discrete method based on the strong discontinuity approach to estimate crack opening is then considered in the post‐treatment phase. Using a Poiseuille's like relation, the coupling between the mechanical state of the material and its dry gas conductivity is performed. For validation purposes, an original experimental campaign is conducted on a dry concrete disc loaded in a splitting setup. During the loading, gas conductivity and digital image correlation analysis are performed. The comparison with the 3D experimental mechanical global response highlights the performance of the mechanical model. The comparison between crack openings measured by digital image correlation and estimated by the strong discontinuity method shows a good agreement. Finally, the results of the semi‐discrete approach coupled with the gas conductivity compared with experimental data show a good estimation of the structural conductivity. Consequently, if the mechanical problem is well modelled at the global scale, then the proposed approach provides good estimation of gas conductivity. Copyright © 2016 John Wiley & Sons, Ltd.     

Review and enhancement of 3D concrete models for large‐scale numerical simulations of concrete structures

The present paper focuses on selected plasticity and damage‐plasticity models for describing the 3D material behavior of concrete. In particular, a plasticity model and a damage‐plasticity model are reviewed and evaluated. Based on the results of the evaluation, enhancements are proposed, aiming at improving the correspondence between predicted and observed material behavior and aiming at implementing a robust and efficient stress update algorithm in a finite element program for performing large‐scale 3D numerical simulations of concrete structures. The capabilities of the concrete models are demonstrated by 3D numerical simulations of benchmark tests with combined bending and torsional loading and combined compression and shear loading and by a large‐scale 3D finite element analysis of a model test of a concrete arch dam. Copyright © 2011 John Wiley & Sons, Ltd.     

基于决策树分类的大屯矿区地物信息提取及矿区污染分析

本文在对国内外遥感图像分类方法充分研究分析的基础上,选择决策树分类法对大屯矿区的Landsat 8遥感图像进行分类研究。选取样本提取并分析研究区典型地类光谱特征曲线,依据光谱曲线特征和归一化植被指数建立了土地利用分类决策树模型,通过反复试验和修正,筛选出适宜大屯矿区地物分类的决策树最优阈值,对研究区进行分类和精度评价,最后通过分类结果对研究区的水体污染状况进行简要分析。