A Bayesian spatio-temporal model of COVID-19 spread in England.

Exploring the spatio-temporal variations of COVID-19 transmission and its potential determinants could provide a deeper understanding of the dynamics of disease spread. This study aimed to investigate the spatio-temporal spread of COVID-19 infections in England, and examine its associations with socioeconomic, demographic and environmental risk factors. We obtained weekly reported COVID-19 cases from 7 March 2020 to 26 March 2022 at Middle Layer Super Output Area (MSOA) level in mainland England from publicly available datasets. With these data, we conducted an ecological study to predict the COVID-19 infection risk and identify its associations with socioeconomic, demographic and environmental risk factors using a Bayesian hierarchical spatio-temporal model. The Bayesian model outperformed the ordinary least squares model and geographically weighted regression model in terms of prediction accuracy. The spread of COVID-19 infections over space and time was heterogeneous. Hotspots of infection risk exhibited inconsistent clustering patterns over time. Risk factors found to be positively associated with COVID-19 infection risk were: annual household income [relative risk (RR) = 1.0008, 95% Credible Interval (CI) 1.0005-1.0012], unemployment rate [RR = 1.0027, 95% CI 1.0024-1.0030], population density on the log scale [RR = 1.0146, 95% CI 1.0129-1.0164], percentage of Caribbean population [RR = 1.0022, 95% CI 1.0009-1.0036], percentage of adults aged 45-64 years old [RR = 1.0031, 95% CI 1.0024-1.0039], and particulate matter ( PM 2.5 ) concentrations [RR = 1.0126, 95% CI 1.0083-1.0167]. The study highlights the importance of considering socioeconomic, demographic, and environmental factors in analysing the spatio-temporal variations of COVID-19 infections in England. The findings could assist policymakers in developing tailored public health interventions at a localised level.

Greenland Mass Trends From Airborne and Satellite Altimetry During 2011-2020.

We use satellite and airborne altimetry to estimate annual mass changes of the Greenland Ice Sheet. We estimate ice loss corresponding to a sea-level rise of 6.9 ± 0.4 mm from April 2011 to April 2020, with a highest annual ice loss rate of 1.4 mm/yr sea-level equivalent from April 2019 to April 2020. On a regional scale, our annual mass loss timeseries reveals 10-15 m/yr dynamic thickening at the terminus of Jakobshavn Isbræ from April 2016 to April 2018, followed by a return to dynamic thinning. We observe contrasting patterns of mass loss acceleration in different basins across the ice sheet and suggest that these spatiotemporal trends could be useful for calibrating and validating prognostic ice sheet models. In addition to resolving the spatial and temporal fingerprint of Greenland's recent ice loss, these mass loss grids are key for partitioning contemporary elastic vertical land motion from longer-term glacial isostatic adjustment (GIA) trends at GPS stations around the ice sheet. Our ice-loss product results in a significantly different GIA interpretation from a previous ice-loss product.

Accelerating Ice Loss From Peripheral Glaciers in North Greenland.

In recent decades, Greenland's peripheral glaciers have experienced large-scale mass loss, resulting in a substantial contribution to sea level rise. While their total area of Greenland ice cover is relatively small (4%), their mass loss is disproportionally large compared to the Greenland ice sheet. Satellite altimetry from Ice, Cloud, and land Elevation Satellite (ICESat) and ICESat-2 shows that mass loss from Greenland's peripheral glaciers increased from 27.2 ± 6.2 Gt/yr (February 2003-October 2009) to 42.3 ± 6.2 Gt/yr (October 2018-December 2021). These relatively small glaciers now constitute 11 ± 2% of Greenland's ice loss and contribute to global sea level rise. In the period October 2018-December 2021, mass loss increased by a factor of four for peripheral glaciers in North Greenland. While peripheral glacier mass loss is widespread, we also observe a complex regional pattern where increases in precipitation at high altitudes have partially counteracted increases in melt at low altitude.

Centennial response of Greenland's three largest outlet glaciers.

The Greenland Ice Sheet is the largest land ice contributor to sea level rise. This will continue in the future but at an uncertain rate and observational estimates are limited to the last few decades. Understanding the long-term glacier response to external forcing is key to improving projections. Here we use historical photographs to calculate ice loss from 1880-2012 for Jakobshavn, Helheim, and Kangerlussuaq glacier. We estimate ice loss corresponding to a sea level rise of 8.1 ± 1.1 millimetres from these three glaciers. Projections of mass loss for these glaciers, using the worst-case scenario, Representative Concentration Pathways 8.5, suggest a sea level contribution of 9.1-14.9 mm by 2100. RCP8.5 implies an additional global temperature increase of 3.7 °C by 2100, approximately four times larger than that which has taken place since 1880. We infer that projections forced by RCP8.5 underestimate glacier mass loss which could exceed this worst-case scenario.

Ice sheet contributions to future sea-level rise from structured expert judgment.

Despite considerable advances in process understanding, numerical modeling, and the observational record of ice sheet contributions to global mean sea-level rise (SLR) since the Fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change, severe limitations remain in the predictive capability of ice sheet models. As a consequence, the potential contributions of ice sheets remain the largest source of uncertainty in projecting future SLR. Here, we report the findings of a structured expert judgement study, using unique techniques for modeling correlations between inter- and intra-ice sheet processes and their tail dependences. We find that since the AR5, expert uncertainty has grown, in particular because of uncertain ice dynamic effects. For a +2 °C temperature scenario consistent with the Paris Agreement, we obtain a median estimate of a 26 cm SLR contribution by 2100, with a 95th percentile value of 81 cm. For a +5 °C temperature scenario more consistent with unchecked emissions growth, the corresponding values are 51 and 178 cm, respectively. Inclusion of thermal expansion and glacier contributions results in a global total SLR estimate that exceeds 2 m at the 95th percentile. Our findings support the use of scenarios of 21st century global total SLR exceeding 2 m for planning purposes. Beyond 2100, uncertainty and projected SLR increase rapidly. The 95th percentile ice sheet contribution by 2200, for the +5 °C scenario, is 7.5 m as a result of instabilities coming into play in both West and East Antarctica. Introducing process correlations and tail dependences increases estimates by roughly 15%.

Decreasing cloud cover drives the recent mass loss on the Greenland Ice Sheet.

The Greenland Ice Sheet (GrIS) has been losing mass at an accelerating rate since the mid-1990s. This has been due to both increased ice discharge into the ocean and melting at the surface, with the latter being the dominant contribution. This change in state has been attributed to rising temperatures and a decrease in surface albedo. We show, using satellite data and climate model output, that the abrupt reduction in surface mass balance since about 1995 can be attributed largely to a coincident trend of decreasing summer cloud cover enhancing the melt-albedo feedback. Satellite observations show that, from 1995 to 2009, summer cloud cover decreased by 0.9 ± 0.3% per year. Model output indicates that the GrIS summer melt increases by 27 ± 13 gigatons (Gt) per percent reduction in summer cloud cover, principally because of the impact of increased shortwave radiation over the low albedo ablation zone. The observed reduction in cloud cover is strongly correlated with a state shift in the North Atlantic Oscillation promoting anticyclonic conditions in summer and suggests that the enhanced surface mass loss from the GrIS is driven by synoptic-scale changes in Arctic-wide atmospheric circulation.

An assessment of forward and inverse GIA solutions for Antarctica.

In this work we assess the most recent estimates of glacial isostatic adjustment (GIA) for Antarctica, including those from both forward and inverse methods. The assessment is based on a comparison of the estimated uplift rates with a set of elastic-corrected GPS vertical velocities. These have been observed from an extensive GPS network and computed using data over the period 2009-2014. We find systematic underestimations of the observed uplift rates in both inverse and forward methods over specific regions of Antarctica characterized by low mantle viscosities and thin lithosphere, such as the northern Antarctic Peninsula and the Amundsen Sea Embayment, where its recent ice discharge history is likely to be playing a role in current GIA. Uplift estimates for regions where many GIA models have traditionally placed their uplift maxima, such as the margins of Filchner-Ronne and Ross ice shelves, are found to be overestimated. GIA estimates show large variability over the interior of East Antarctica which results in increased uncertainties on the ice-sheet mass balance derived from gravimetry methods.

Geodetic measurements reveal similarities between post-Last Glacial Maximum and present-day mass loss from the Greenland ice sheet.

Accurate quantification of the millennial-scale mass balance of the Greenland ice sheet (GrIS) and its contribution to global sea-level rise remain challenging because of sparse in situ observations in key regions. Glacial isostatic adjustment (GIA) is the ongoing response of the solid Earth to ice and ocean load changes occurring since the Last Glacial Maximum (LGM; ~21 thousand years ago) and may be used to constrain the GrIS deglaciation history. We use data from the Greenland Global Positioning System network to directly measure GIA and estimate basin-wide mass changes since the LGM. Unpredicted, large GIA uplift rates of +12 mm/year are found in southeast Greenland. These rates are due to low upper mantle viscosity in the region, from when Greenland passed over the Iceland hot spot about 40 million years ago. This region of concentrated soft rheology has a profound influence on reconstructing the deglaciation history of Greenland. We reevaluate the evolution of the GrIS since LGM and obtain a loss of 1.5-m sea-level equivalent from the northwest and southeast. These same sectors are dominating modern mass loss. We suggest that the present destabilization of these marine-based sectors may increase sea level for centuries to come. Our new deglaciation history and GIA uplift estimates suggest that studies that use the Gravity Recovery and Climate Experiment satellite mission to infer present-day changes in the GrIS may have erroneously corrected for GIA and underestimated the mass loss by about 20 gigatons/year.

Spatial and temporal Antarctic Ice Sheet mass trends, glacio-isostatic adjustment, and surface processes from a joint inversion of satellite altimeter, gravity, and GPS data.

We present spatiotemporal mass balance trends for the Antarctic Ice Sheet from a statistical inversion of satellite altimetry, gravimetry, and elastic-corrected GPS data for the period 2003-2013. Our method simultaneously determines annual trends in ice dynamics, surface mass balance anomalies, and a time-invariant solution for glacio-isostatic adjustment while remaining largely independent of forward models. We establish that over the period 2003-2013, Antarctica has been losing mass at a rate of -84 ± 22 Gt yr-1, with a sustained negative mean trend of dynamic imbalance of -111 ± 13 Gt yr-1. West Antarctica is the largest contributor with -112 ± 10 Gt yr-1, mainly triggered by high thinning rates of glaciers draining into the Amundsen Sea Embayment. The Antarctic Peninsula has experienced a dramatic increase in mass loss in the last decade, with a mean rate of -28 ± 7 Gt yr-1 and significantly higher values for the most recent years following the destabilization of the Southern Antarctic Peninsula around 2010. The total mass loss is partly compensated by a significant mass gain of 56 ± 18 Gt yr-1 in East Antarctica due to a positive trend of surface mass balance anomalies.

Spatial and temporal distribution of mass loss from the Greenland Ice Sheet since AD 1900.

The response of the Greenland Ice Sheet (GIS) to changes in temperature during the twentieth century remains contentious, largely owing to difficulties in estimating the spatial and temporal distribution of ice mass changes before 1992, when Greenland-wide observations first became available. The only previous estimates of change during the twentieth century are based on empirical modelling and energy balance modelling. Consequently, no observation-based estimates of the contribution from the GIS to the global-mean sea level budget before 1990 are included in the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Here we calculate spatial ice mass loss around the entire GIS from 1900 to the present using aerial imagery from the 1980s. This allows accurate high-resolution mapping of geomorphic features related to the maximum extent of the GIS during the Little Ice Age at the end of the nineteenth century. We estimate the total ice mass loss and its spatial distribution for three periods: 1900-1983 (75.1 ± 29.4 gigatonnes per year), 1983-2003 (73.8 ± 40.5 gigatonnes per year), and 2003-2010 (186.4 ± 18.9 gigatonnes per year). Furthermore, using two surface mass balance models we partition the mass balance into a term for surface mass balance (that is, total precipitation minus total sublimation minus runoff) and a dynamic term. We find that many areas currently undergoing change are identical to those that experienced considerable thinning throughout the twentieth century. We also reveal that the surface mass balance term shows a considerable decrease since 2003, whereas the dynamic term is constant over the past 110 years. Overall, our observation-based findings show that during the twentieth century the GIS contributed at least 25.0 ± 9.4 millimetres of global-mean sea level rise. Our result will help to close the twentieth-century sea level budget, which remains crucial for evaluating the reliability of models used to predict global sea level rise.

Arctic circulation regimes.

Between 1948 and 1996, mean annual environmental parameters in the Arctic experienced a well-pronounced decadal variability with two basic circulation patterns: cyclonic and anticyclonic alternating at 5 to 7 year intervals. During cyclonic regimes, low sea-level atmospheric pressure (SLP) dominated over the Arctic Ocean driving sea ice and the upper ocean counterclockwise; the Arctic atmosphere was relatively warm and humid, and freshwater flux from the Arctic Ocean towards the subarctic seas was intensified. By contrast, during anticylonic circulation regimes, high SLP dominated driving sea ice and the upper ocean clockwise. Meanwhile, the atmosphere was cold and dry and the freshwater flux from the Arctic to the subarctic seas was reduced. Since 1997, however, the Arctic system has been under the influence of an anticyclonic circulation regime (17 years) with a set of environmental parameters that are atypical for this regime. We discuss a hypothesis explaining the causes and mechanisms regulating the intensity and duration of Arctic circulation regimes, and speculate how changes in freshwater fluxes from the Arctic Ocean and Greenland impact environmental conditions and interrupt their decadal variability.

Paleofluvial mega-canyon beneath the central Greenland ice sheet.

Subglacial topography plays an important role in modulating the distribution and flow of basal water. Where topography predates ice sheet inception, it can also reveal insights into former tectonic and geomorphological processes. Although such associations are known in Antarctica, little consideration has been given to them in Greenland, partly because much of the ice sheet bed is thought to be relatively flat and smooth. Here, we present evidence from ice-penetrating radar data for a 750-km-long subglacial canyon in northern Greenland that is likely to have influenced basal water flow from the ice sheet interior to the margin. We suggest that the mega-canyon predates ice sheet inception and will have influenced basal hydrology in Greenland over past glacial cycles.

Aerial photographs reveal late-20th-century dynamic ice loss in northwestern Greenland.

Global warming is predicted to have a profound impact on the Greenland Ice Sheet and its contribution to global sea-level rise. Recent mass loss in the northwest of Greenland has been substantial. Using aerial photographs, we produced digital elevation models and extended the time record of recent observed marginal dynamic thinning back to the mid-1980s. We reveal two independent dynamic ice loss events on the northwestern Greenland Ice Sheet margin: from 1985 to 1993 and 2005 to 2010, which were separated by limited mass changes. Our results suggest that the ice mass changes in this sector were primarily caused by short-lived dynamic ice loss events rather than changes in the surface mass balance. This finding challenges predictions about the future response of the Greenland Ice Sheet to increasing global temperatures.

Reassessment of the potential sea-level rise from a collapse of the West Antarctic Ice Sheet.

Theory has suggested that the West Antarctic Ice Sheet may be inherently unstable. Recent observations lend weight to this hypothesis. We reassess the potential contribution to eustatic and regional sea level from a rapid collapse of the ice sheet and find that previous assessments have substantially overestimated its likely primary contribution. We obtain a value for the global, eustatic sea-level rise contribution of about 3.3 meters, with important regional variations. The maximum increase is concentrated along the Pacific and Atlantic seaboard of the United States, where the value is about 25% greater than the global mean, even for the case of a partial collapse.

Integrating satellite observations with modelling: basal shear stress of the Filcher-Ronne ice streams, Antarctica.

Using inverse methods constrained by recent satellite observations, we have produced a comprehensive estimate of the basal shear stress beneath the Filchner-Ronne ice streams. The inversions indicate that a weak bed (approx. 4-20kPa) underlies much of these ice streams. Compared to the Ross ice streams, the distribution of weak subglacial till is more heterogeneous, with 'sticky spots' providing much of the resistance to flow. A weak bed beneath Recovery ice stream extends several hundred kilometres inland with flow. Along this ice stream, discrepancies between thickness measurements and flux estimates suggest the existence of a deep (-1400m) trough not resolved by existing maps of subglacial topography. We hypothesize that the presence of this and other deep troughs is a major influence on this sector of the ice sheet that is not fully incorporated in current models of ice-sheet evolution.