Probabilistic reasoning about measurements of equilibrium climate sensitivity: combining disparate lines of evidence

Where policy and science intersect, there are always issues of ambiguous and conflicting lines of evidence. Combining disparate information sources is mathematically complex; common heuristics based on simple statistical models easily lead us astray. Here, we use Bayesian Nets (BNs) to illustrate the complexity in reasoning under uncertainty. Data from joint research at Resources for the Future and NASA Langley are used to populate a BN for predicting equilibrium climate sensitivity (ECS). The information sources consist of measuring the rate of decadal temperature rise (DTR) and measuring the rate of percentage change in cloud radiative forcing (CRF), with both the existing configuration of satellites and with a proposed enhanced measuring system. The goal of all measurements is to reduce uncertainty in equilibrium climate sensitivity. Subtle aspects of probabilistic reasoning with concordant and discordant measurements are illustrated. Relative to the current prior distribution on ECS, we show that after 30 years of observing with the current systems, the 2σ uncertainty band for ECS would be shrunk on average to 73% of its current value. With the enhanced systems over the same time, it would be shrunk to 32% of its current value. The actual shrinkage depends on the values actually observed. These results are based on models recommended by the Social Cost of Carbon methodology and assume a Business as Usual emissions path.     

Land use and climate change on the Yamal Peninsula of north-west Siberia: some ecological and socio-economic implications

Compared to climate, land use change is expected to comprise a more important component of global change in the coming decades. However, climate is anticipated to supass land use as a factor later in the next century, particularly in the Arctic. Discussed here are the implications of land use and climate change on the Yamal peninsula of north-west Siberia, homeland of the Yamal Nenets. Since the discovery of super-giant natural gas fields in the 1960s, extensive exploration has resulted in direct withdrawal of large areas for infrastructure development and associated disturbance regimes have led to cumulative impacts on thousands of additional hectares of land. The land withdrawals have pushed a relatively consistent or increasing number of reindeer onto progressively smaller parcels of pasture. This has led to excessive grazing and trampling of lichens, bryophytes and shrubs and, in many areas, erosion of sandy soils via deflation. The low Arctic tundra lies entirely within the continuous permafrost zone and ice-rich substrates are widespread. One implication of this is that both anthropogenic and zoogenic distubance regimes may easily initiate thermokarst and aeolian rosion, leading to significant further losses of pastures. Even without industrial distubance, a slight change of the climate would result in massive thermokarst erosion. This would have negative consequences equal to or greater than the mechanical distubances described above. The synergistic effects of land use coupled with climate, change therefore have profound implication for the ecosystems of Yamal, as well as the future of the Nenets culture, society and economy.     

Mapping Submerged Aquatic Vegetation with GIS in the Caloosahatchee Estuary: Evaluation of Different Interpolation Methods

This article evaluates different spatial interpolation methods for mapping submerged aquatic vegetation (SAV) in the Caloosahatchee Estuary, Florida. Data used for interpolation were collected by the Submersed Aquatic Vegetation Early Warning System (SAVEWS). The system consists of hydro-acoustic equipment, which operates from a slow-moving boat and records bottom depth, seagrass height, and seagrass density. This information is coupled with geographic location coordinates from a Global Positioning System (GPS) and stored together in digital files, representing SAV status at points along transect lines. Adequate spatial interpolation is needed to present the SAV information, including density, height, and water depth, as spatially continuous data for mapping and for comparison between seasons and years. Interpolation methods examined in this study include ordinary kriging with five different semivariance models combined with a variable number of neighboring points, the inverse distance weighted (IDW) method with different parameters, and the triangulated irregular network (TIN) method with linear and quintic options. Interpolation results were compared with survey data at selected calibration transects to examine the suitability of different interpolation methods. Suitability was quantified by the determination coefficient (R2) and the root-mean-square error (RMSE) between interpolated and observed values. The most suitable interpolation method was identified as the one yielding the highest R2 value and/or the lowest RMSE value. For different geographic conditions, seasons, and SAV parameters, different interpolation methods were recommended. This study identified that kriging was more suitable than the IDW or TIN method for spatial interpolation of all SAV parameters measured. It also suggested that transect data with irregular spatial distribution patterns such as SAV parameters are sensitive to interpolation methods. An inappropriate interpolation method such as TIN can lead to erroneous spatial representation of the SAV status. With a functional geographic system and adequate computing power, the evaluation and selection of interpolation methods can be automated and quantitative, leading to a more efficient and accurate decision.     

Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

Soil water dynamics under different forest vegetation cover: Implications for hillslope stability

Though it is well known that vegetation affects the water balance of soils through canopy interception and evapotranspiration, its hydrological contribution to soil hydrology and stability is yet to be fully quantified. To improve understanding of this hydrological process, soil water dynamics have been monitored at three adjacent hillslopes with different vegetation covers (deciduous tree cover, coniferous tree cover, and grass cover), for nine months from December 2014 to September 2015. The monitored soil moisture values were translated into soil matric suction (SMS) values to facilitate the analysis of hillslope stability. Our observations showed significant seasonal variations in SMS for each vegetation cover condition. However, a significant difference between different vegetation covers was only evident during the winter season where the mean SMS under coniferous tree cover (83.6 kPa) was significantly greater than that under grass cover (41 kPa). The hydrological reinforcing contribution due to matric suction was highest for the hillslope with coniferous tree cover, while the hillslope with deciduous tree cover was second and the hillslope with grass cover was third. The greatest contributions for all cover types were during the summer season. During the winter season, the wettest period of the monitoring study, the additional hydrological reinforcing contributions provided by the deciduous tree cover (1.5 to 6.5 kPa) or the grass cover (0.9 to 5.4 kPa) were insufficient to avoid potential slope failure conditions. However, the additional hydrological reinforcing contribution from the coniferous tree cover (5.8 to 10.4 kPa) was sufficient to provide potentially stable hillslope conditions during the winter season. Our study clearly suggests that during the winter season the hydrological effects from both deciduous tree and grass covers are insufficient to promote slope stability, while the hydrological reinforcing effects from the coniferous tree cover are sufficient even during the winter season. Copyright © 2018 John Wiley & Sons, Ltd.     

The effect of riparian tree roots on the mass‐stability of riverbanks

Plants interact with and modify the processes of riverbank erosion by altering bank hydrology, flow hydraulics and bank geotechnical properties. The physically based slope stability model GWEDGEM was used to assess how changes in bank geotechnical properties due to the roots of native Australian riparian trees affected the stability of bank sections surveyed along the Latrobe River. Modelling bank stability against mass failure with and without the reinforcing effects of River Red Gum (Eucalyptus camaldulensis) or Swamp Paperbark (Melaleuca ericifolia) indicates that root reinforcement of the bank substrate provides high levels of bank protection. The model indicates that the addition of root reinforcement to an otherwise unstable bank section can raise the factor of safety (F _s) from F _s = 1·0 up to about F _s = 1·6. The addition of roots to riverbanks improves stability even under worst‐case hydrological conditions and is apparent over a range of bank geometries, varying with tree position. Trees growing close to potential failure plane locations, either low on the bank or on the floodplain, realize the greatest bank reinforcement. Copyright © 2000 John Wiley & Sons, Ltd.     

Deposition and preservation of fluvio-tidal shallow-marine sandstones: A re-evaluation of the Neoproterozoic Jura Quartzite (western Scotland)

The 2 to 5 km thick, sandstone-dominated (>90%) Jura Quartzite is an extreme example of a mature Neoproterozoic sandstone, previously interpreted as a tide-influenced shelf deposit and herein re-interpreted within a fluvio-tidal deltaic depositional model. Three issues are addressed: (i) evidence for the re-interpretation from tidal shelf to tidal delta; (ii) reasons for vertical facies uniformity; and (iii) sand supply mechanisms to form thick tidal-shelf sandstones. The predominant facies (compound cross-bedded, coarse-grained sandstones) represents the lower parts of metres to tens of metres high, transverse fluvio-tidal bedforms with superimposed smaller bedforms. Ubiquitous erosional surfaces, some with granule–pebble lags, record erosion of the upper parts of those bedforms. There was selective preservation of the higher energy, topographically-lower, parts of channel-bar systems. Strongly asymmetrical, bimodal, palaeocurrents are interpreted as due to associated selective preservation of fluvially-enhanced ebb tidal currents. Finer-grained facies are scarce, due largely to suspended sediment bypass. They record deposition in lower-energy environments, including channel mouth bars, between and down depositional-dip of higher energy fluvio-ebb tidal bars. The lack of wave-formed sedimentary structures and low continuity of mudstone and sandstone interbeds, support deposition in a non-shelf setting. Hence, a sand-rich, fluvial–tidal, current-dominated, largely sub-tidal, delta setting is proposed. This new interpretation avoids the problem of transporting large amounts of coarse sand to a shelf. Facies uniformity and vertical stacking are likely due to sediment oversupply and bypass rather than balanced sediment supply and subsidence rates. However, facies evidence of relative sea level changes is difficult to recognise, which is attributed to: (i) the areally extensive and polygenetic nature of the preserved facies, and (ii) a large stored sediment buffer that dampened response to relative sea-level and/or sediment supply changes. Consideration of preservation bias towards high-energy deposits may be more generally relevant, especially to thick Neoproterozoic and Lower Palaeozoic marine sandstones.     

The impact of storm-triggered landslides on sediment dynamics and catchment-wide denudation rates in the southern Central Range of Taiwan following the extreme rainfall event of Typhoon Morakot

Extreme erosion events can produce large short-term sediment fluxes. Such events complicate erosion rates estimated from cosmogenic nuclide concentrations in river sediment by providing sediment with a concentration different from the long-term basin average. We present a detrital ¹⁰Be study in southern Taiwan, with multiple samples obtained in a time sequence bracketing the 2009 Typhoon Morakot, to assess the impact of landslide sediment on ¹⁰Be concentrations (N¹⁰Be) in river sediment. Sediment samples were collected from 13 major basins, two or three times over the last decade, to observe the temporal variation of N¹⁰Be. Landslide inventories with time intervals of 5–6 years were used to quantify sediment flux changes. A negative correlation between N¹⁰Be and landslide areal density indicates dilution of N¹⁰Be by landslide sediment. Denudation rates estimated from the diluted N¹⁰Be can be up to three times higher than the lowest rate derived from the same basins. Observed increases imply that, 3 years after the passage of Typhoon Morakot, fluvial channels still contain a considerable amount of sediment provided by hillslope landslides during the event. However, higher N¹⁰Be in 2016 samples indicate that the contribution from landslide sediment at the sampled grain size has decreased with time. The correlation between changes in N¹⁰Be and landslide area and volume is not strong, likely resulting from the stochastic and complex nature of sediment transport. To simultaneously evaluate the volume of landslide-derived sediment and estimate the background denudation rate, associated with less impulsive sediment supply, we constructed a sediment-mixing model with the time series of N¹⁰Be and landslide inventories. The spatial pattern of background erosion rate in southern Taiwan is consistent with the regional tectonic framework, indicating that the landscape is evolving mainly in response to the tectonic forcing, and this signal is modified, but not obscured by impulsive sediment supply. © 2019 John Wiley & Sons, Ltd.