Search

We investigate the application of spectral induced polarization (SIP) monitoring to understand seasonal and annual variations in the freeze–thaw processes in permafrost by examining the frequency dependence of subsurface electrical properties. We installed a permanent SIP monitoring profile at a high-mountain permafrost site in the Italian Alps in 2019 and collected SIP data in the frequency range between 0.1–75 Hz over 3 years. The SIP imaging results were interpreted in conjunction with complementary seismic and borehole data sets. In particular, we investigated the phase frequency effect (ϕFE), i.e., the change in the resistivity phase with frequency. We observe that this parameter (ϕFE) is strongly sensitive to temperature changes and might be used as a proxy to delineate spatial and temporal changes in the ice content in the subsurface, providing information not accessible through electrical resistivity tomography (ERT) or single-frequency IP measurements. Temporal changes in ϕFE are validated through laboratory SIP measurements on samples from the site in controlled freeze–thaw experiments. We demonstrate that SIP is capable of resolving temporal changes in the thermal state and the ice / water ratio associated with seasonal freeze–thaw processes. We investigate the consistency between the ϕFE observed in field data and groundwater and ice content estimates derived from petrophysical modeling of ERT and seismic data.

The subglacial landscape of Antarctica records and influences the behaviour of its overlying ice sheet. However, in many places, the evolution of the landscape and its control on ice sheet behaviour have not been investigated in detail. Using recently released radio-echo sounding data, we investigate the subglacial landscape of the Evans–Rutford region of West Antarctica. Following quantitative analysis of the landscape morphology under ice-loaded and ice-unloaded conditions, we identify 10 flat surfaces distributed across the region. Across these 10 surfaces, we identify two distinct populations based on clustering of elevations, which potentially represent remnants of regionally coherent pre-glacial surfaces underlying the West Antarctic Ice Sheet (WAIS). The surfaces are bounded by deeply incised glacial troughs, some of which have potential tectonic controls. We assess two hypotheses for the evolution of the regional landscape: (1) passive-margin evolution associated with the break-up of the Gondwana supercontinent or (2) an extensive planation surface that may have been uplifted in association with either the West Antarctic Rift System or cessation of subduction at the base of the Antarctic Peninsula. We suggest that passive-margin evolution is the most likely of these two mechanisms, with the erosion of glacial troughs adjacent to, and incising, the flat surfaces likely having coincided with the growth of the WAIS. These flat surfaces also demonstrate similarities to other identified surfaces, indicating that a similar formational process may have been acting more widely around the Weddell Sea embayment. The subsequent fluctuations of ice flow, basal thermal regime, and erosion patterns of the WAIS are therefore controlled by the regional tectonic structures.

This work assesses a recently produced 21-member climate model large ensemble (LE) based on the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 2 (E3SM2). The ensemble spans the historical era (1850 to 2014) and 21st century (2015 to 2100), using the SSP370 pathway, allowing for an evaluation of the model's forced response. A companion 500-year preindustrial control simulation is used to initialize the ensemble and estimate drift. Characteristics of the LE are documented and compared against other recently produced ensembles using the E3SM version 1 (E3SM1) and Community Earth System Model (CESM) versions 1 and 2. Simulation drift is found to be smaller, and model agreement with observations is higher in versions 2 of E3SM and CESM versus their version 1 counterparts. Shortcomings in E3SM2 include a lack of warming from the mid to late 20th century, likely due to excessive cooling influence of anthropogenic sulfate aerosols, an issue also evident in E3SM1. Associated impacts on the water cycle and energy budgets are also identified. Considerable model dependence in the response to both aerosols and greenhouse gases is documented and E3SM2's sensitivity to variable prescribed biomass burning emissions is demonstrated. Various E3SM2 and CESM2 model benchmarks are found to be on par with the highest-performing recent generation of climate models, establishing the E3SM2 LE as an important resource for estimating climate variability and responses, though with various caveats as discussed herein. As an illustration of the usefulness of LEs in estimating the potential influence of internal variability, the observed CERES-era trend in net top-of-atmosphere flux is compared to simulated trends and found to be much larger than the forced response in all LEs, with only a few members exhibiting trends as large as observed, thus motivating further study.

The analysis of englacial layers using radio-echo sounding data enables the characterisation and reconstruction of current and past ice-sheet flow. Despite the Lambert Glacier catchment being one of the largest in Antarctica, discharging ∼16 % of East Antarctica's ice, its englacial architecture has been little analysed. Here, we present a comprehensive analysis of Lambert Glacier's englacial architecture using radio-echo sounding data collected by Antarctica's Gamburtsev Province Project (AGAP) North survey. We used an internal layering continuity index (ILCI) to characterise the internal architecture of the ice and identify four macro-scale ILCI zones with distinct glaciological contexts. Whilst the catchment is dominated by continuous englacial layering, disrupted or discontinuous layering is highlighted by the ILCI at both the onset of enhanced ice flow (defined here as >15 m a−1) and along the shear margin, suggesting a transition in englacial deformation conditions and converging ice flow. These zones are characterised by buckled and folded englacial layers which have fold axes aligned with the current ice-flow regime. These folds suggest that the flow direction of the Lambert Glacier trunk has changed little, if at all, during the Holocene. Disturbed englacial layers that do not correspond to modern ice-flow routing found within a deep subglacial channel, however, suggest that ice-flow change has occurred in a former tributary that fed Lambert Glacier from grid north. As large outlet systems such as Lambert Glacier are likely to play a vital role in the future drainage of the East Antarctic Ice Sheet, constraining their englacial architecture to reconstruct their past ice flow and determine basal conditions is important for refining projections of future sea-level change.

One of the key components of this research has been the mapping of Antarctic bed topography and ice thickness parameters that are crucial for modelling ice flow and hence for predicting future ice loss and the ensuing sea level rise. Supported by the Scientific Committee on Antarctic Research (SCAR), the Bedmap3 Action Group aims not only to produce new gridded maps of ice thickness and bed topography for the international scientific community, but also to standardize and make available all the geophysical survey data points used in producing the Bedmap gridded products. Here, we document the survey data used in the latest iteration, Bedmap3, incorporating and adding to all of the datasets previously used for Bedmap1 and Bedmap2, including ice bed, surface and thickness point data from all Antarctic geophysical campaigns since the 1950s. More specifically, we describe the processes used to standardize and make these and future surveys and gridded datasets accessible under the Findable, Accessible, Interoperable, and Reusable (FAIR) data principles. With the goals of making the gridding process reproducible and allowing scientists to re-use the data freely for their own analysis, we introduce the new SCAR Bedmap Data Portal (https://bedmap.scar.org, last access: 1 March 2023) created to provide unprecedented open access to these important datasets through a web-map interface. We believe that this data release will be a valuable asset to Antarctic research and will greatly extend the life cycle of the data held within it. Data are available from the UK Polar Data Centre: https://data.bas.ac.uk (last access: 5 May 2023​​​​​​​). See the Data availability section for the complete list of datasets.

The last decade has seen the development of a growing tendency to use the most modern technologies, in the field of Cultural Heritage, with the aim of digitizing and facilitating protection and conservation activities. Much research has focused on the development of innovative methods such as photogrammetry or Terrestrial Laser Scanners, in terms of reliability, precision, time and costs. In this research, however, the use of the smartphone was investigated by comparing the point clouds obtained via videogrammetry from smartphones, with those generated by different digital survey techniques, such as Terrestrial Laser Scanners and photogrammetry via SLR camera. Specifically, a smartphone was used and the comparison between point clouds was conducted based on four criteria: point clouds fitting, density evaluation, profiling and texture quality with the aim, therefore, of verifying the geometric reliability of the data and the quality of polygonal lines and the mesh/texture derived. The case study selected was the sepulchral monument of the Pascopepe Lambertini family (14th century), located in the Crypt of Santa Maria della Scala, the Cathedral of Trani (South of Italy). Finally, the research demonstrates how this methodology, in the documentation of heritage, allows for even greater portability and accessibility compared to other methodologies, maintaining widely acceptable standards of accuracy and therefore going to constitute a valid alternative in the documentation of historical heritage.

The presence of subglacial lakes and subglacial hydrological networks under the East and West Antarctic ice sheets is now relatively well understood, whilst their influence on ice dynamics is the subject of ongoing research. In contrast, little is known about subglacial lakes and hydrological networks under the Antarctic Peninsula Ice Sheet and how these are influencing glacier behaviour. Here we describe the rapid drainage and slow refill of a subglacial lake under Mars Glacier using remote sensing and aerogeophysics. Results suggest drainage of the subglacial lake occurred prior to 2011, resulting in the collapse of the overlying ice into the newly formed subglacial cavity. The cavity has been refilling since this time, with peak rates of infilling associated with seasonal surface meltwater activity. We review evidence for similar features elsewhere in the Antarctic Peninsula and discuss whether their appearance marks a threshold shift in glacier thermal regimes and the activation or enhancement of their subglacial hydrological networks by surface meltwater. Collectively, these features show coupling of climate processes and the bed of the region's glaciers highlighting their ongoing vulnerability to climate change.

Solar climate intervention using stratospheric aerosol injection is a proposed method of reducing global mean temperatures to reduce the worst consequences of climate change. A detailed assessment of responses and impacts of such an intervention is needed with multiple global models to support societal decisions regarding the use of these approaches to help address climate change. We present a new modeling protocol aimed at simulating a plausible deployment of stratospheric aerosol injection and reproducibility of simulations using other Earth system models: Assessing Responses and Impacts of Solar climate intervention on the Earth system with stratospheric aerosol injection (ARISE-SAI). The protocol and simulations are aimed at enabling community assessment of responses of the Earth system to solar climate intervention. ARISE-SAI simulations are designed to be more policy-relevant than existing large ensembles or multi-model simulation sets. We describe in detail the first set of ARISE-SAI simulations, ARISE-SAI-1.5, which utilize a moderate emissions scenario, introduce stratospheric aerosol injection at ∼21.5 km in the year 2035, and keep global mean surface air temperature near 1.5 ∘C above the pre-industrial value utilizing a feedback or control algorithm. We present the detailed setup, aerosol injection strategy, and preliminary climate analysis from a 10-member ensemble of these simulations carried out with the Community Earth System Model version 2 with the Whole Atmosphere Community Climate Model version 6 as its atmospheric component.

The potential for multiyear prediction of impactful Earth system change remains relatively underexplored compared to shorter (subseasonal to seasonal) and longer (decadal) timescales. In this study, we introduce a new initialized prediction system using the Community Earth System Model version 2 (CESM2) that is specifically designed to probe potential and actual prediction skill at lead times ranging from 1 month out to 2 years. The Seasonal-to-Multiyear Large Ensemble (SMYLE) consists of a collection of 2-year-long hindcast simulations, with four initializations per year from 1970 to 2019 and an ensemble size of 20. A full suite of output is available for exploring near-term predictability of all Earth system components represented in CESM2. We show that SMYLE skill for El Niño–Southern Oscillation is competitive with other prominent seasonal prediction systems, with correlations exceeding 0.5 beyond a lead time of 12 months. A broad overview of prediction skill reveals varying degrees of potential for useful multiyear predictions of seasonal anomalies in the atmosphere, ocean, land, and sea ice. The SMYLE dataset, experimental design, model, initial conditions, and associated analysis tools are all publicly available, providing a foundation for research on multiyear prediction of environmental change by the wider community.

The mid-Pliocene (∼3 Ma) is one of the most recent warm periods with high CO2 concentrations in the atmosphere and resulting high temperatures, and it is often cited as an analog for near-term future climate change. Here, we apply a moisture budget analysis to investigate the response of the large-scale hydrological cycle at low latitudes within a 13-model ensemble from the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2). The results show that increased atmospheric moisture content within the mid-Pliocene ensemble (due to the thermodynamic effect) results in wetter conditions over the deep tropics, i.e., the Pacific intertropical convergence zone (ITCZ) and the Maritime Continent, and drier conditions over the subtropics. Note that the dynamic effect plays a more important role than the thermodynamic effect in regional precipitation minus evaporation (PmE) changes (i.e., northward ITCZ shift and wetter northern Indian Ocean). The thermodynamic effect is offset to some extent by a dynamic effect involving a northward shift of the Hadley circulation that dries the deep tropics and moistens the subtropics in the Northern Hemisphere (i.e., the subtropical Pacific). From the perspective of Earth's energy budget, the enhanced southward cross-equatorial atmospheric transport (0.22 PW), induced by the hemispheric asymmetries of the atmospheric energy, favors an approximately 1∘ northward shift of the ITCZ. The shift of the ITCZ reorganizes atmospheric circulation, favoring a northward shift of the Hadley circulation. In addition, the Walker circulation consistently shifts westward within PlioMIP2 models, leading to wetter conditions over the northern Indian Ocean. The PlioMIP2 ensemble highlights that an imbalance of interhemispheric atmospheric energy during the mid-Pliocene could have led to changes in the dynamic effect, offsetting the thermodynamic effect and, hence, altering mid-Pliocene hydroclimate.

While climate change mitigation targets necessarily concern maximum mean state changes, understanding impacts and developing adaptation strategies will be largely contingent on how climate variability responds to increasing anthropogenic perturbations. Thus far Earth system modeling efforts have primarily focused on projected mean state changes and the sensitivity of specific modes of climate variability, such as the El Niño–Southern Oscillation. However, our knowledge of forced changes in the overall spectrum of climate variability and higher-order statistics is relatively limited. Here we present a new 100-member large ensemble of climate change projections conducted with the Community Earth System Model version 2 over 1850–2100 to examine the sensitivity of internal climate fluctuations to greenhouse warming. Our unprecedented simulations reveal that changes in variability, considered broadly in terms of probability distribution, amplitude, frequency, phasing, and patterns, are ubiquitous and span a wide range of physical and ecosystem variables across many spatial and temporal scales. Greenhouse warming in the model alters variance spectra of Earth system variables that are characterized by non-Gaussian probability distributions, such as rainfall, primary production, or fire occurrence. Our modeling results have important implications for climate adaptation efforts, resource management, seasonal predictions, and assessing potential stressors for terrestrial and marine ecosystems.

The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.

Subglacial water plays an important role in ice sheet dynamics and stability. Subglacial lakes are often located at the onset of ice streams and have been hypothesised to enhance ice flow downstream by lubricating the ice–bed interface. The most recent subglacial-lake inventory of Antarctica mapped nearly 400 lakes, of which ∼ 14 % are found in West Antarctica. Despite the potential importance of subglacial water for ice dynamics, there is a lack of detailed subglacial-water characterisation in West Antarctica. Using radio-echo sounding data, we analyse the ice–bed interface to detect subglacial lakes. We report 33 previously uncharted subglacial lakes and present a systematic analysis of their physical properties. This represents a ∼ 40 % increase in subglacial lakes in West Antarctica. Additionally, a new digital elevation model of basal topography of the Ellsworth Subglacial Highlands was built and used to create a hydropotential model to simulate the subglacial hydrological network. This allows us to characterise basal hydrology, determine subglacial water catchments and assess their connectivity. We show that the simulated subglacial hydrological catchments of the Rutford Ice Stream, Pine Island Glacier and Thwaites Glacier do not correspond to their ice surface catchments.

Landslides in glacial environments are high-magnitude, long-runout events, believed to be increasing in frequency as a paraglacial response to ice retreat and thinning and, arguably, due to warming temperatures and degrading permafrost above current glaciers. However, our ability to test these assumptions by quantifying the temporal sequencing of debris inputs over large spatial and temporal extents is limited in areas with glacier ice. Discrete landslide debris inputs, particularly in accumulation areas, are rapidly “lost”, being reworked by motion and icefalls and/or covered by snowfall. Although large landslides can be detected and located using their seismic signature, smaller (M≤5.0) landslides frequently go undetected because their seismic signature is less than the noise floor, particularly supraglacially deposited landslides, which feature a “quiet” runout over snow. Here, we present GERALDINE (Google Earth Engine supRaglAciaL Debris INput dEtector): a new free-to-use tool leveraging Landsat 4–8 satellite imagery and Google Earth Engine. GERALDINE outputs maps of new supraglacial debris additions within user-defined areas and time ranges, providing a user with a reference map, from which large debris inputs such as supraglacial landslides (>0.05 km2) can be rapidly identified. We validate the effectiveness of GERALDINE outputs using published supraglacial rock avalanche inventories, and then demonstrate its potential by identifying two previously unknown, large (>2 km2) landslide-derived supraglacial debris inputs onto glaciers in the Hayes Range, Alaska, one of which was not detected seismically. GERALDINE is a first step towards a complete global magnitude–frequency of landslide inputs onto glaciers over the 38 years of Landsat Thematic Mapper imagery.

We present a topographic digital elevation model (DEM) for Princess Elizabeth Land (PEL), East Antarctica. The DEM covers an area of ∼900 000 km2 and was built from radio-echo sounding data collected during four campaigns since 2015. Previously, to generate the Bedmap2 topographic product, PEL's bed was characterized from low-resolution satellite gravity data across an otherwise large (>200 km wide) data-free zone. We use the mass conservation (MC) method to produce an ice thickness grid across faster flowing (>30 m yr−1) regions of the ice sheet and streamline diffusion in slower flowing areas. The resulting ice thickness model is integrated with an ice surface model to build the bed DEM. Together with BedMachine Antarctica and Bedmap2, this new bed DEM completes the first-order measurement of subglacial continental Antarctica – an international mission that began around 70 years ago. The ice thickness data and bed DEMs of PEL (resolved horizontally at 500 m relative to ice surface elevations obtained from the Reference Elevation Model of Antarctica – REMA) are accessible from https://doi.org/10.5281/zenodo.4023343 (Cui et al., 2020a) and https://doi.org/10.5281/zenodo.4023393 (Cui et al., 2020b).

With semiconductor technology gradually approaching its physical and thermal limits, recent supercomputers have adopted major architectural changes to continue increasing the performance through more power-efficient heterogeneous many-core systems. Examples include Sunway TaihuLight that has four management processing elements (MPEs) and 256 computing processing elements (CPEs) inside one processor and Summit that has two central processing units (CPUs) and six graphics processing units (GPUs) inside one node. Meanwhile, current high-resolution Earth system models that desperately require more computing power generally consist of millions of lines of legacy code developed for traditional homogeneous multicore processors and cannot automatically benefit from the advancement of supercomputer hardware. As a result, refactoring and optimizing the legacy models for new architectures become key challenges along the road of taking advantage of greener and faster supercomputers, providing better support for the global climate research community and contributing to the long-lasting societal task of addressing long-term climate change. This article reports the efforts of a large group in the International Laboratory for High-Resolution Earth System Prediction (iHESP) that was established by the cooperation of Qingdao Pilot National Laboratory for Marine Science and Technology (QNLM), Texas A&M University (TAMU), and the National Center for Atmospheric Research (NCAR), with the goal of enabling highly efficient simulations of the high-resolution (25 km atmosphere and 10 km ocean) Community Earth System Model (CESM-HR) on Sunway TaihuLight. The refactoring and optimizing efforts have improved the simulation speed of CESM-HR from 1 SYPD (simulation years per day) to 3.4 SYPD (with output disabled) and supported several hundred years of pre-industrial control simulations. With further strategies on deeper refactoring and optimizing for remaining computing hotspots, as well as redesigning architecture-oriented algorithms, we expect an equivalent or even better efficiency to be gained on the new platform than traditional homogeneous CPU platforms. The refactoring and optimizing processes detailed in this paper on the Sunway system should have implications for similar efforts on other heterogeneous many-core systems such as GPU-based high-performance computing (HPC) systems.

It has been hypothesized that complex englacial structures identified within the East Antarctic and Greenland ice sheets are generated by (i) water freezing to the ice sheet base and evolving under ice flow, (ii) deformation of ice of varying rheology, or (iii) entrainment of basal material. Using ice-penetrating radar, we identify a widespread complex of deep-ice facies in West Antarctica that exist in the absence of basal water. These deep-ice units are extensive, thick (>500 m), and incorporate multiple highly reflective englacial layers. At the lateral margin of an enhanced flow tributary of the Institute Ice Stream, these units are heavily deformed and folded by the action of lateral flow convergence. Radar reflectivity analysis demonstrates that the uppermost reflector of the deep-ice package is highly anisotropic, due to abrupt alternations in crystal orientation fabric, and consequently will have a different rheology to the ice above and below it. Deformation and folding of the deep-ice package is an englacial response to the combination of laterally-convergent ice flow and the physical properties of the ice column.

Plant nutrients can be recycled through microbial decomposition of organic matter but replacement of base cations and phosphorus, lost through harvesting of biomass/biofuels or leaching, requires de novo supply of fresh nutrients released through weathering of soil parent material (minerals and rocks). Weathering involves physical and chemical processes that are modified by biological activity of plants, microorganisms and animals. This article reviews recent progress made in understanding biological processes contributing to weathering. A perspective of increasing spatial scale is adopted, examining the consequences of biological activity for weathering from nanoscale interactions, through in vitro and in planta microcosm and mesocosm studies, to field experiments, and finally ecosystem and global level effects. The topics discussed include the physical alteration of minerals and mineral surfaces; the composition, amounts, chemical properties, and effects of plant and microbial secretions; and the role of carbon flow (including stabilisation and sequestration of C in organic and inorganic forms). Although the predominant focus is on the effects of fungi in forest ecosystems, the properties of biofilms, including bacterial interactions, are also discussed. The implications of these biological processes for modelling are discussed, and we attempt to identify some key questions and knowledge gaps, as well as experimental approaches and areas of research in which future studies are likely to yield useful results. A particular focus of this article is to improve the representation of the ways in which biological processes complement physical and chemical processes that mobilise mineral elements, making them available for plant uptake. This is necessary to produce better estimates of weathering that are required for sustainable management of forests in a post-fossil-fuel economy. While there are abundant examples of nanometre- and micrometre-scale physical interactions between microorganisms and different minerals, opinion appears to be divided with respect to the quantitative significance of these observations for overall weathering. Numerous in vitro experiments and microcosm studies involving plants and their associated microorganisms suggest that the allocation of plant-derived carbon, mineral dissolution and plant nutrient status are tightly coupled, but there is still disagreement about the extent to which these processes contribute to field-scale observations. Apart from providing dynamically responsive pathways for the allocation of plant-derived carbon to power dissolution of minerals, mycorrhizal mycelia provide conduits for the long-distance transportation of weathering products back to plants that are also quantitatively significant sinks for released nutrients. These mycelial pathways bridge heterogeneous substrates, reducing the influence of local variation in C:N ratios. The production of polysaccharide matrices by biofilms of interacting bacteria and/or fungi at interfaces with mineral surfaces and roots influences patterns of production of antibiotics and quorum sensing molecules, with concomitant effects on microbial community structure, and the qualitative and quantitative composition of mineral-solubilising compounds and weathering products. Patterns of carbon allocation and nutrient mobilisation from both organic and inorganic substrates have been studied at larger spatial and temporal scales, including both ecosystem and global levels, and there is a generally wider degree of acceptance of the “systemic” effects of microorganisms on patterns of nutrient mobilisation. Theories about the evolutionary development of weathering processes have been advanced but there is still a lack of information connecting processes at different spatial scales. Detailed studies of the liquid chemistry of local weathering sites at the micrometre scale, together with upscaling to soil-scale dissolution rates, are advocated, as well as new approaches involving stable isotopes.

We report three active subglacial lakes within 2 km of the lateral margin of Isunguata Sermia, West Greenland, identified by differencing time-stamped ArcticDEM strips. Each lake underwent one drainage–refill event between 2009 and 2017, with two lakes draining in < 1 month in August 2014 and August 2015. The 2015 drainage caused a ∼ 1-month down-glacier slowdown in ice flow and flooded the foreland, aggrading the proglacial channel by 8 m. The proglacial flooding confirms the ice-surface elevation anomalies as subglacial water bodies and demonstrates how their drainage can significantly modify proglacial environments. These subglacial lakes offer accessible targets for geophysical investigations and exploration.

The Community Atmosphere Model version 6 (CAM6), released in 2018 as part of the Community Earth System Model version 2 (CESM2), is a major upgrade over the previous CAM5 that has been used in numerous global and regional climate studies. Since CESM2–CAM6 will participate in the upcoming Coupled Model Intercomparison Project phase 6 (CMIP6) and is likely to be adopted in many future studies, its simulation fidelity needs to be thoroughly examined. Here we evaluate the performance of a developmental version of the Community Atmosphere Model with parameterizations that will be used in version 6 (CAM6α), with a default 1∘ horizontal resolution (0.9∘×1.25∘, CAM6α-1∘) and a high-resolution configuration (approximately 0.25∘, CAM6α-0.25∘), against various observational and reanalysis datasets of precipitation over Asia. CAM6α performance is compared with CAM5 at default 1∘ horizontal resolution (CAM5-1∘) and a high-resolution configuration at 0.25∘ (CAM5-0.25∘). With the prognostic treatment of precipitation processes and the new microphysics module, CAM6α is able to better simulate climatological mean and extreme precipitation over Asia, better capture the heaviest precipitation events, better reproduce the diurnal cycle of precipitation rates over most of Asia, and better simulate the probability density distributions of daily precipitation over Tibet, Korea, Japan and northern China. Higher horizontal resolution in CAM6α improves the simulation of mean and extreme precipitation over northern China, but the performance degrades over the Maritime Continent. Moisture budget diagnosis suggests that the physical processes leading to model improvement are different over different regions. Both upgraded physical parameterizations and higher horizontal resolution affect the simulated precipitation response to the internal variability of the climate system (e.g., Asian monsoon variability, El Niño–Southern Oscillation – ENSO, Pacific Decadal Oscillation – PDO), but the effects vary across different regions. For example, higher horizontal resolution degrades the model performance in simulating precipitation variability over southern China associated with the East Asian summer monsoon. In contrast, precipitation variability associated with ENSO improves with upgraded physical parameterizations and higher horizontal resolution. CAM6α-0.25∘ and CAM6α-1∘ show an opposite response to the PDO over southern China. Basically, the response to increases in horizontal resolution is dependent on the CAM version.

The current assessments of the carbon turnover in the Arctic tundra are subject to large uncertainties. This problem can (inter alia) be ascribed to both the general shortage of flux data from the vast and sparsely inhabited Arctic region, as well as the typically high spatiotemporal variability of carbon fluxes in tundra ecosystems. Addressing these challenges, carbon dioxide fluxes on an active flood plain situated in the Siberian Lena River Delta were studied during two growing seasons with the eddy covariance method. The footprint exhibited a heterogeneous surface, which generated mixed flux signals that could be partitioned in such a way that both respiratory loss and photosynthetic gain were obtained for each of two vegetation classes. This downscaling of the observed fluxes revealed a differing seasonality in the net uptake of bushes (−0.89 µmol m−2 s−1) and sedges (−0.38 µmol mm−2 s−1) in 2014. That discrepancy, which was concealed in the net signal, resulted from a comparatively warm spring in conjunction with an early snowmelt and a varying canopy structure. Thus, the representativeness of footprints may adversely be affected in response to prolonged unusual weather conditions. In 2015, when air temperatures on average corresponded to climatological means, both vegetation-class-specific flux rates were of similar magnitude (−0.69 µmol m−2 s−1). A comprehensive set of measures (e.g. phenocam) corroborated the reliability of the partitioned fluxes and hence confirmed the utility of flux decomposition for enhanced flux data analysis. This scrutiny encompassed insights into both the phenological dynamic of individual vegetation classes and their respective functional flux to flux driver relationships with the aid of ecophysiologically interpretable parameters. For comparison with other sites, the decomposed fluxes were employed in a vegetation class area-weighted upscaling that was based on a classified high-resolution orthomosaic of the flood plain. In this way, robust budgets that take the heterogeneous surface characteristics into account were estimated. In relation to the average sink strength of various Arctic flux sites, the flood plain constitutes a distinctly stronger carbon dioxide sink. Roughly 42 % of this net uptake, however, was on average offset by methane emissions lowering the sink strength for greenhouse gases. With growing concern about rising greenhouse gas emissions in high-latitude regions, providing robust carbon budgets from tundra ecosystems is critical in view of accelerating permafrost thaw, which can impact the global climate for centuries.

We present a new digital elevation model (DEM) of the bed, with a 1 km gridding, of the Weddell Sea (WS) sector of the West Antarctic Ice Sheet (WAIS). The DEM has a total area of ∼ 125 000 km2 covering the Institute, Möller and Foundation ice streams, as well as the Bungenstock ice rise. In comparison with the Bedmap2 product, our DEM includes new aerogeophysical datasets acquired by the Center for Remote Sensing of Ice Sheets (CReSIS) through the NASA Operation IceBridge (OIB) program in 2012, 2014 and 2016. We also improve bed elevation information from the single largest existing dataset in the region, collected by the British Antarctic Survey (BAS) Polarimetric radar Airborne Science Instrument (PASIN) in 2010–2011, from the relatively crude measurements determined in the field for quality control purposes used in Bedmap2. While the gross form of the new DEM is similar to Bedmap2, there are some notable differences. For example, the position and size of a deep subglacial trough (∼ 2 km below sea level) between the ice-sheet interior and the grounding line of the Foundation Ice Stream have been redefined. From the revised DEM, we are able to better derive the expected routing of basal water and, by comparison with that calculated using Bedmap2, we are able to assess regions where hydraulic flow is sensitive to change. Given the potential vulnerability of this sector to ocean-induced melting at the grounding line, especially in light of the improved definition of the Foundation Ice Stream trough, our revised DEM will be of value to ice-sheet modelling in efforts to quantify future glaciological changes in the region and, from this, the potential impact on global sea level. The new 1 km bed elevation product of the WS sector can be found at https://doi.org/10.5281/zenodo.1035488.

In the paper an example for the application of kriging methods to estimate damage to buildings in crisis scenarios is introduced. Furthermore, the Java implementations for Ordinary and Universal Kriging on mobile GIS are presented. As variogram models an exponential, a Gaussian and a spherical variogram are tested in detail. Different test constellations are introduced with various information densities. As test data set, public data from the analysis of the 2010 Haiti earthquake by satellite images are pre-processed and visualized in a Geographic Information System. As buildings, topography and other external influences cannot be seen as being constant for the whole area under investigation, semi variograms are calculated by consulting neighboured classified buildings using the so called moving window method. The evaluation of the methods shows that the underlying variogram model is the determining factor for the quality of the interpolation rather than the choice of the kriging method or increasing the information density of a random sample. The implementation is completely realized with the programming language Java. Thereafter, the implemented software component is integrated into GeoTech Mobile, a mobile GIS Android application based on the processing of standardized spatial data representations defined by the Open Geospatial Consortium (OGC). As a result the implemented methods can be used on mobile devices, i.e. they may be transferred to other application fields. That is why we finally point out further research with new applications in the Dubai region.

Mineral dust plays an important role in the climate system by interacting with radiation, clouds, and biogeochemical cycles. In addition, natural archives show that the dust cycle experienced variability in the past in response to global and local climate change. The compilation of the DIRTMAP (Dust Indicators and Records from Terrestrial and MArine Palaeoenvironments) paleodust data sets in the last 2 decades provided a benchmark for paleoclimate models that include the dust cycle, following a time slice approach. We propose an innovative framework to organize a paleodust data set that builds on the positive experience of DIRTMAP and takes into account new scientific challenges by providing a concise and accessible data set of temporally resolved records of dust mass accumulation rates and particle grain size distributions. We consider data from ice cores, marine sediments, loess–paleosol sequences, lake sediments, and peat bogs for this compilation, with a temporal focus on the Holocene period. This global compilation allows the investigation of the potential, uncertainties, and confidence level of dust mass accumulation rate reconstructions and highlights the importance of dust particle size information for accurate and quantitative reconstructions of the dust cycle. After applying criteria that help to establish that the data considered represent changes in dust deposition, 45 paleodust records have been identified, with the highest density of dust deposition data occurring in the North Atlantic region. Although the temporal evolution of dust in the North Atlantic appears consistent across several cores and suggests that minimum dust fluxes are likely observed during the early to mid-Holocene period (6000–8000 years ago), the magnitude of dust fluxes in these observations is not fully consistent, suggesting that more work needs to be done to synthesize data sets for the Holocene. Based on the data compilation, we used the Community Earth System Model to estimate the mass balance of and variability in the global dust cycle during the Holocene, with dust loads ranging from 17.2 to 20.8 Tg between 2000 and 10 000 years ago and with a minimum in the early to mid-Holocene (6000–8000 years ago).

We present ice-penetrating radar evidence for ancient (pre-glacial) and extensive erosion surfaces preserved beneath the upstream Institute and Möller ice streams, West Antarctica. Radar data reveal a smooth, laterally continuous, gently sloping topographic block, comprising two surfaces separated by a distinct break in slope. The erosion surfaces are preserved in this location due to the collective action of the Pirrit and Martin–Nash hills on ice sheet flow, resulting in a region of slow flowing, cold-based ice downstream of these major topographic barriers. Our analysis reveals that smooth, flat subglacial topography does not always correspond to regions of either present or former fast ice flow, as has previously been assumed. We discuss the potential origins of the erosion surfaces. Erosion rates across the surfaces are currently low, precluding formation via present-day glacial erosion. We suggest that fluvial or marine processes are most likely to have resulted in the formation of these surfaces, but we acknowledge that distinguishing between these processes with certainty requires further data.

A recent ocean modelling study indicates that possible changes in circulation may bring warm deep-ocean water into direct contact with the grounding lines of the Filchner–Ronne ice streams, suggesting the potential for future ice losses from this sector equivalent to ~0.3 m of sea-level rise. Significant advancements have been made in our knowledge of both the basal topography and ice velocity in the Weddell Sea sector, and the ability to accurately model marine ice sheet dynamics, thus enabling an assessment to be made of the relative sensitivities of the diverse collection of ice streams feeding the Filchner–Ronne Ice Shelf. Here we use the BISICLES ice sheet model, which employs adaptive-mesh refinement to resolve grounding line dynamics, to carry out such an assessment. The impact of realistic perturbations to the surface and sub-shelf mass balance forcing fields from our 2000-year "reference" model run indicate that both the Institute and Möller ice streams are highly sensitive to changes in basal melting either near to their respective grounding lines, or in the region of the ice rises within the Filchner–Ronne Ice Shelf. These same perturbations have little impact, however, on the Rutford, Carlson or Foundation ice streams, while the Evans Ice Stream is found to enter a phase of unstable retreat only after melt at its grounding line has increased by 50% of likely present-day values.

Repeat-pass ICESat altimetry has revealed 124 discrete surface height changes across the Antarctic Ice Sheet, interpreted to be caused by subglacial lake discharges (surface lowering) and inputs (surface uplift). Few of these active lakes have been confirmed by radio-echo sounding (RES) despite several attempts (notable exceptions are Lake Whillans and three in the Adventure Subglacial Trench). Here we present targeted RES and radar altimeter data from an "active lake" location within the upstream Institute Ice Stream, into which at least 0.12 km3 of water was previously calculated to have flowed between October 2003 and February 2008. We use a series of transects to establish an accurate depiction of the influences of bed topography and ice surface elevation on water storage potential. The location of surface height change is downstream of a subglacial hill on the flank of a distinct topographic hollow, where RES reveals no obvious evidence for deep (> 10 m) water. The regional hydropotential reveals a sink coincident with the surface change, however. Governed by the location of the hydrological sink, basal water will likely "drape" over topography in a manner dissimilar to subglacial lakes where flat strong specular RES reflections are measured. The inability of RES to detect the active lake means that more of the Antarctic ice sheet bed may contain stored water than is currently appreciated. Variation in ice surface elevation data sets leads to significant alteration in calculations of the local flow of basal water indicating the value of, and need for, high-resolution altimetry data in both space and time to establish and characterise subglacial hydrological processes.

The last interglaciation (~130 to 116 ka) is a time period with a strong astronomically induced seasonal forcing of insolation compared to the present. Proxy records indicate a significantly different climate to that of the modern, in particular Arctic summer warming and higher eustatic sea level. Because the forcings are relatively well constrained, it provides an opportunity to test numerical models which are used for future climate prediction. In this paper we compile a set of climate model simulations of the early last interglaciation (130 to 125 ka), encompassing a range of model complexities. We compare the simulations to each other and to a recently published compilation of last interglacial temperature estimates. We show that the annual mean response of the models is rather small, with no clear signal in many regions. However, the seasonal response is more robust, and there is significant agreement amongst models as to the regions of warming vs cooling. However, the quantitative agreement of the model simulations with data is poor, with the models in general underestimating the magnitude of response seen in the proxies. Taking possible seasonal biases in the proxies into account improves the agreement, but only marginally. However, a lack of uncertainty estimates in the data does not allow us to draw firm conclusions. Instead, this paper points to several ways in which both modelling and data could be improved, to allow a more robust model–data comparison.

We present Bedmap2, a new suite of gridded products describing surface elevation, ice-thickness and the seafloor and subglacial bed elevation of the Antarctic south of 60° S. We derived these products using data from a variety of sources, including many substantial surveys completed since the original Bedmap compilation (Bedmap1) in 2001. In particular, the Bedmap2 ice thickness grid is made from 25 million measurements, over two orders of magnitude more than were used in Bedmap1. In most parts of Antarctica the subglacial landscape is visible in much greater detail than was previously available and the improved data-coverage has in many areas revealed the full scale of mountain ranges, valleys, basins and troughs, only fragments of which were previously indicated in local surveys. The derived statistics for Bedmap2 show that the volume of ice contained in the Antarctic ice sheet (27 million km3) and its potential contribution to sea-level rise (58 m) are similar to those of Bedmap1, but the mean thickness of the ice sheet is 4.6% greater, the mean depth of the bed beneath the grounded ice sheet is 72 m lower and the area of ice sheet grounded on bed below sea level is increased by 10%. The Bedmap2 compilation highlights several areas beneath the ice sheet where the bed elevation is substantially lower than the deepest bed indicated by Bedmap1. These products, along with grids of data coverage and uncertainty, provide new opportunities for detailed modelling of the past and future evolution of the Antarctic ice sheets.

Vegetation-climate interactions are thought to have amplified polar warmth during past warm periods. Here, we explore the vegetation-climate interactions in the mid-Cretaceous using a fully coupled ocean-atmosphere general circulation model with a dynamic vegetation component. We run simulations with 1x, 10x and 16x pre-industrial atmospheric CO2. Results show that forests expand from mid-latitudes to high latitudes as CO2 increases from 1x to 10x and 16x, mainly due to the CO2-induced warming. This expansion of mid-to-high latitude forests are largely supported by the distribution of mid-Cretaceous fossil woods and coal deposits. Globally, the presence of vegetation increases mean annual temperature and precipitation by 0.9 °C and 0.11 mm day−1 relative to bare ground. High-latitude warming induced by the presence of vegetation (∼1.9 °C) is less than half of that reported in previous studies. The weaker warming here is mainly due to less pronounced albedo feedbacks, and to a less extent, reduced poleward heat transport via weakening of the meridional overturning circulation. Our results suggest that other mechanisms in addition to high atmospheric CO2 and high-latitude vegetation are required to maintain the polar warmth.

We present results from new geophysical data allowing modelling of the water flow within Subglacial Lake Ellsworth (SLE), West Antarctica. Our simulations indicate that this lake has a novel temperature distribution due to significantly thinner ice than other surveyed subglacial lakes. The critical pressure boundary (tipping depth), established from the semi-empirical Equation of State, defines whether the lake's flow regime is convective or stratified. It passes through SLE and separates different temperature (and flow) regimes on either side of the lake. Our results have implications for the location of proposed access holes into SLE, the choice of which will depend on scientific or operational priorities. If an understanding of subglacial lake water properties and dynamics is the priority, holes are required in a basal freezing area at the North end of the lake. This would be the preferred priority suggested by this paper, requiring temperature and salinity profiles in the water column. A location near the Southern end, where bottom currents are lowest, is optimum for detecting the record of life in the bed sediments; to minimise operational risk and maximise the time span of a bed sediment core, a location close to the middle of the lake, where the basal interface is melting and the lake bed is at its deepest, remains the best choice. Considering potential lake-water salinity and ice-density variations, we estimate the critical tipping depth, separating different temperature regimes within subglacial lakes, to be in about 2900 to 3045 m depth.

In 2008 the temporal focus of the Palaeoclimate Modelling Intercomparison Project was expanded to include a model intercomparison for the mid-Pliocene warm period (3.29–2.97 million years ago). This project is referred to as PlioMIP (Pliocene Model Intercomparison Project). Two experiments have been agreed upon and comprise phase 1 of PlioMIP. The first (Experiment 1) will be performed with atmosphere-only climate models. The second (Experiment 2) will utilise fully coupled ocean-atmosphere climate models. The aim of this paper is to provide a detailed model intercomparison project description which documents the experimental design in a more detailed way than has previously been done in the literature. Specifically, this paper describes the experimental design and boundary conditions that will be utilised for Experiment 1 of PlioMIP.

No abstract available.