National Oceanography Centre

We present the Met Office Hadley Centre's sea ice and sea surface temperature (SST) data set, HadISST1, and the nighttime marine air temperature (NMAT) data set, HadMAT1. HadISST1 replaces the global sea ice and sea surface temperature (GISST) data sets and is a unique combination of monthly globally complete fields of SST and sea ice concentration on a 1° latitude‐longitude grid from 1871. The companion HadMAT1 runs monthly from 1856 on a 5° latitude‐longitude grid and incorporates new corrections for the effect on NMAT of increasing deck (and hence measurement) heights. HadISST1 and HadMAT1 temperatures are reconstructed using a two‐stage reduced‐space optimal interpolation procedure, followed by superposition of quality‐improved gridded observations onto the reconstructions to restore local detail. The sea ice fields are made more homogeneous by compensating satellite microwave‐based sea ice concentrations for the impact of surface melt effects on retrievals in the Arctic and for algorithm deficiencies in the Antarctic and by making the historical in situ concentrations consistent with the satellite data. SSTs near sea ice are estimated using statistical relationships between SST and sea ice concentration. HadISST1 compares well with other published analyses, capturing trends in global, hemispheric, and regional SST well, containing SST fields with more uniform variance through time and better month‐to‐month persistence than those in GISST. HadMAT1 is more consistent with SST and with collocated land surface air temperatures than previous NMAT data sets.

Distributions of Earth's species are changing at accelerating rates, increasingly driven by human-mediated climate change. Such changes are already altering the composition of ecological communities, but beyond conservation of natural systems, how and why does this matter? We review evidence that climate-driven species redistribution at regional to global scales affects ecosystem functioning, human well-being, and the dynamics of climate change itself. Production of natural resources required for food security, patterns of disease transmission, and processes of carbon sequestration are all altered by changes in species distribution. Consideration of these effects of biodiversity redistribution is critical yet lacking in most mitigation and adaptation strategies, including the United Nation's Sustainable Development Goals.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate – the “global carbon budget” – is important tobetter understand the global carbon cycle, support the development ofclimate policies, and project future climate change. Here we describe andsynthesize data sets and methodology to quantify the five major componentsof the global carbon budget and their uncertainties. Fossil CO2emissions (EFOS) are based on energy statistics and cement productiondata, while emissions from land-use change (ELUC), mainlydeforestation, are based on land use and land-use change data andbookkeeping models. Atmospheric CO2 concentration is measured directlyand its growth rate (GATM) is computed from the annual changes inconcentration. The ocean CO2 sink (SOCEAN) and terrestrialCO2 sink (SLAND) are estimated with global process modelsconstrained by observations. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the lastdecade available (2010–2019), EFOS was 9.6 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.4 ± 0.5 GtC yr−1 when the cement carbonation sink is included), andELUC was 1.6 ± 0.7 GtC yr−1. For the same decade, GATM was 5.1 ± 0.02 GtC yr−1 (2.4 ± 0.01 ppm yr−1), SOCEAN 2.5 ± 0.6 GtC yr−1, and SLAND 3.4 ± 0.9 GtC yr−1, with a budgetimbalance BIM of −0.1 GtC yr−1 indicating a near balance betweenestimated sources and sinks over the last decade. For the year 2019 alone, thegrowth in EFOS was only about 0.1 % with fossil emissions increasingto 9.9 ± 0.5 GtC yr−1 excluding the cement carbonation sink (9.7 ± 0.5 GtC yr−1 when cement carbonation sink is included), and ELUC was 1.8 ± 0.7 GtC yr−1, for total anthropogenic CO2 emissions of 11.5 ± 0.9 GtC yr−1 (42.2 ± 3.3 GtCO2). Also for 2019, GATM was5.4 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEANwas 2.6 ± 0.6 GtC yr−1, and SLAND was 3.1 ± 1.2 GtC yr−1, with a BIM of 0.3 GtC. The global atmospheric CO2concentration reached 409.85 ± 0.1 ppm averaged over 2019. Preliminarydata for 2020, accounting for the COVID-19-induced changes in emissions,suggest a decrease in EFOS relative to 2019 of about −7 % (medianestimate) based on individual estimates from four studies of −6 %, −7 %,−7 % (−3 % to −11 %), and −13 %. Overall, the mean and trend in thecomponents of the global carbon budget are consistently estimated over theperiod 1959–2019, but discrepancies of up to 1 GtC yr−1 persist for therepresentation of semi-decadal variability in CO2 fluxes. Comparison ofestimates from diverse approaches and observations shows (1) no consensusin the mean and trend in land-use change emissions over the last decade, (2)a persistent low agreement between the different methods on the magnitude ofthe land CO2 flux in the northern extra-tropics, and (3) an apparentdiscrepancy between the different methods for the ocean sink outside thetropics, particularly in the Southern Ocean. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set (Friedlingstein et al.,2019; Le Quéré et al., 2018b, a, 2016, 2015b, a, 2014,2013). The data presented in this work are available at https://doi.org/10.18160/gcp-2020 (Friedlingstein et al., 2020).

Much of our understanding of Earth's past climate comes from the measurement of oxygen and carbon isotope variations in deep-sea benthic foraminifera. Yet, long intervals in existing records lack the temporal resolution and age control needed to thoroughly categorize climate states of the Cenozoic era and to study their dynamics. Here, we present a new, highly resolved, astronomically dated, continuous composite of benthic foraminifer isotope records developed in our laboratories. Four climate states-Hothouse, Warmhouse, Coolhouse, Icehouse-are identified on the basis of their distinctive response to astronomical forcing depending on greenhouse gas concentrations and polar ice sheet volume. Statistical analysis of the nonlinear behavior encoded in our record reveals the key role that polar ice volume plays in the predictability of Cenozoic climate dynamics.

Abstract. Accurate assessment of anthropogenic carbon dioxide (CO2) emissions andtheir redistribution among the atmosphere, ocean, and terrestrial biospherein a changing climate is critical to better understand the global carboncycle, support the development of climate policies, and project futureclimate change. Here we describe and synthesize data sets and methodologies toquantify the five major components of the global carbon budget and theiruncertainties. Fossil CO2 emissions (EFOS) are based on energystatistics and cement production data, while emissions from land-use change(ELUC), mainly deforestation, are based on land use and land-use changedata and bookkeeping models. Atmospheric CO2 concentration is measureddirectly, and its growth rate (GATM) is computed from the annualchanges in concentration. The ocean CO2 sink (SOCEAN) is estimatedwith global ocean biogeochemistry models and observation-baseddata products. The terrestrial CO2 sink (SLAND) is estimated withdynamic global vegetation models. The resulting carbon budget imbalance(BIM), the difference between the estimated total emissions and theestimated changes in the atmosphere, ocean, and terrestrial biosphere, is ameasure of imperfect data and understanding of the contemporary carboncycle. All uncertainties are reported as ±1σ. For the year 2021, EFOS increased by 5.1 % relative to 2020, withfossil emissions at 10.1 ± 0.5 GtC yr−1 (9.9 ± 0.5 GtC yr−1 when the cement carbonation sink is included), and ELUC was 1.1 ± 0.7 GtC yr−1, for a total anthropogenic CO2 emission(including the cement carbonation sink) of 10.9 ± 0.8 GtC yr−1(40.0 ± 2.9 GtCO2). Also, for 2021, GATM was 5.2 ± 0.2 GtC yr−1 (2.5 ± 0.1 ppm yr−1), SOCEAN was 2.9 ± 0.4 GtC yr−1, and SLAND was 3.5 ± 0.9 GtC yr−1, with aBIM of −0.6 GtC yr−1 (i.e. the total estimated sources were too low orsinks were too high). The global atmospheric CO2 concentration averaged over2021 reached 414.71 ± 0.1 ppm. Preliminary data for 2022 suggest anincrease in EFOS relative to 2021 of +1.0 % (0.1 % to 1.9 %)globally and atmospheric CO2 concentration reaching 417.2 ppm, morethan 50 % above pre-industrial levels (around 278 ppm). Overall, the meanand trend in the components of the global carbon budget are consistentlyestimated over the period 1959–2021, but discrepancies of up to 1 GtC yr−1 persist for the representation of annual to semi-decadalvariability in CO2 fluxes. Comparison of estimates from multipleapproaches and observations shows (1) a persistent large uncertainty in theestimate of land-use change emissions, (2) a low agreement between thedifferent methods on the magnitude of the land CO2 flux in the northernextratropics, and (3) a discrepancy between the different methods on thestrength of the ocean sink over the last decade. This living data updatedocuments changes in the methods and data sets used in this new globalcarbon budget and the progress in understanding of the global carbon cyclecompared with previous publications of this data set. The data presented inthis work are available at https://doi.org/10.18160/GCP-2022 (Friedlingstein et al., 2022b).

Abstract Accurate inventory of tropical peatland is important in order to (a) determine the magnitude of the carbon pool; (b) estimate the scale of transfers of peat‐derived greenhouse gases to the atmosphere resulting from land use change; and (c) support carbon emissions reduction policies. We review available information on tropical peatland area and thickness and calculate peat volume and carbon content in order to determine their best estimates and ranges of variation. Our best estimate of tropical peatland area is 441 025 km 2 (∼11% of global peatland area) of which 247 778 km 2 (56%) is in Southeast Asia. We estimate the volume of tropical peat to be 1758 Gm 3 (∼18–25% of global peat volume) with 1359 Gm 3 in Southeast Asia (77% of all tropical peat). This new assessment reveals a larger tropical peatland carbon pool than previous estimates, with a best estimate of 88.6 Gt (range 81.7–91.9 Gt) equal to 15–19% of the global peat carbon pool. Of this, 68.5 Gt (77%) is in Southeast Asia, equal to 11–14% of global peat carbon. A single country, Indonesia, has the largest share of tropical peat carbon (57.4 Gt, 65%), followed by Malaysia (9.1 Gt, 10%). These data are used to provide revised estimates for Indonesian and Malaysian forest soil carbon pools of 77 and 15 Gt, respectively, and total forest carbon pools (biomass plus soil) of 97 and 19 Gt. Peat carbon contributes 60% to the total forest soil carbon pool in Malaysia and 74% in Indonesia. These results emphasize the prominent global and regional roles played by the tropical peat carbon pool and the importance of including this pool in national and regional assessments of terrestrial carbon stocks and the prediction of peat‐derived greenhouse gas emissions.

Abstract We document the development of the first version of the U.K. Earth System Model UKESM1. The model represents a major advance on its predecessor HadGEM2‐ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well‐resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric‐stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two‐moment, five‐species, modal aerosol; and ocean biogeochemistry with two‐way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre‐industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger‐than‐observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC −1 .

Large, abrupt, and widespread climate changes with major impacts have occurred repeatedly in the past, when the Earth system was forced across thresholds. Although abrupt climate changes can occur for many reasons, it is conceivable that human forcing of climate change is increasing the probability of large, abrupt events. Were such an event to recur, the economic and ecological impacts could be large and potentially serious. Unpredictability exhibited near climate thresholds in simple models shows that some uncertainty will always be associated with projections. In light of these uncertainties, policy-makers should consider expanding research into abrupt climate change, improving monitoring systems, and taking actions designed to enhance the adaptability and resilience of ecosystems and economies.

Abstract We assess evidence relevant to Earth's equilibrium climate sensitivity per doubling of atmospheric CO 2 , characterized by an effective sensitivity S . This evidence includes feedback process understanding, the historical climate record, and the paleoclimate record. An S value lower than 2 K is difficult to reconcile with any of the three lines of evidence. The amount of cooling during the Last Glacial Maximum provides strong evidence against values of S greater than 4.5 K. Other lines of evidence in combination also show that this is relatively unlikely. We use a Bayesian approach to produce a probability density function (PDF) for S given all the evidence, including tests of robustness to difficult‐to‐quantify uncertainties and different priors. The 66% range is 2.6–3.9 K for our Baseline calculation and remains within 2.3–4.5 K under the robustness tests; corresponding 5–95% ranges are 2.3–4.7 K, bounded by 2.0–5.7 K (although such high‐confidence ranges should be regarded more cautiously). This indicates a stronger constraint on S than reported in past assessments, by lifting the low end of the range. This narrowing occurs because the three lines of evidence agree and are judged to be largely independent and because of greater confidence in understanding feedback processes and in combining evidence. We identify promising avenues for further narrowing the range in S , in particular using comprehensive models and process understanding to address limitations in the traditional forcing‐feedback paradigm for interpreting past changes.

From measurements of Hf–Yb mixtures, we have found that the correction of isobaric interferences involving accepted Yb isotope ratios and reasonable estimates of mass bias result in a significantly under-corrected 176Hf, which is proportional to the amount of Yb added. This can be explained by (1) a significant difference in the instrumental mass bias between Hf and Yb, and (2) that the accepted values for isotopic ratios within the Yb and/or Hf systems are incorrect. We have evaluated these possibilities by measuring mixed solutions of Yb and Hf on two MC-ICP-MS instruments and undertaking a series of REE fractionation experiments using a thermal ionisation mass spectrometer (TIMS). Our results indicate that the presently accepted abundances of the Yb isotopes are not appropriate. We present new values for Yb isotopic abundances based on the TIMS and MC-ICP-MS results. Using the newly defined Yb values, we demonstrate that Yb and Hf have similar levels of mass bias in plasma ionisation instruments, and that Hf isotope ratios can be used to correct Yb mass bias before subsequent correction of isobaric interference. A laser ablation comparison of Yb and Hf indicates that similar relationships exist, and can be applied to micro-analytical techniques where chemical separation is not possible.

Physiological thermal-tolerance limits of terrestrial ectotherms often exceed local air temperatures, implying a high degree of thermal safety (an excess of warm or cold thermal tolerance). However, air temperatures can be very different from the equilibrium body temperature of an individual ectotherm. Here, we compile thermal-tolerance limits of ectotherms across a wide range of latitudes and elevations and compare these thermal limits both to air and to operative body temperatures (theoretically equilibrated body temperatures) of small ectothermic animals during the warmest and coldest times of the year. We show that extreme operative body temperatures in exposed habitats match or exceed the physiological thermal limits of most ectotherms. Therefore, contrary to previous findings using air temperatures, most ectotherms do not have a physiological thermal-safety margin. They must therefore rely on behavior to avoid overheating during the warmest times, especially in the lowland tropics. Likewise, species living at temperate latitudes and in alpine habitats must retreat to avoid lethal cold exposure. Behavioral plasticity of habitat use and the energetic consequences of thermal retreats are therefore critical aspects of species' vulnerability to climate warming and extreme events.

Increasing quantities of atmospheric anthropogenic fixed nitrogen entering the open ocean could account for up to about a third of the ocean's external (nonrecycled) nitrogen supply and up to approximately 3% of the annual new marine biological production, approximately 0.3 petagram of carbon per year. This input could account for the production of up to approximately 1.6 teragrams of nitrous oxide (N2O) per year. Although approximately 10% of the ocean's drawdown of atmospheric anthropogenic carbon dioxide may result from this atmospheric nitrogen fertilization, leading to a decrease in radiative forcing, up to about two-thirds of this amount may be offset by the increase in N2O emissions. The effects of increasing atmospheric nitrogen deposition are expected to continue to grow in the future.

Ocean acidification may have severe consequences for marine ecosystems; however, assessing its future impact is difficult because laboratory experiments and field observations are limited by their reduced ecologic complexity and sample period, respectively. In contrast, the geological record contains long-term evidence for a variety of global environmental perturbations, including ocean acidification plus their associated biotic responses. We review events exhibiting evidence for elevated atmospheric CO(2), global warming, and ocean acidification over the past ~300 million years of Earth's history, some with contemporaneous extinction or evolutionary turnover among marine calcifiers. Although similarities exist, no past event perfectly parallels future projections in terms of disrupting the balance of ocean carbonate chemistry-a consequence of the unprecedented rapidity of CO(2) release currently taking place.

The response of coastal wetlands to sea-level rise during the twentyfirst century remains uncertain. Global-scale projections suggest that between 20 and 90 per cent (for low and high sea-level rise scenarios, respectively) of the present-day coastal wetland area will be lost, which will in turn result in the loss of biodiversity and highly valued ecosystem services 1-3 . These projections do not necessarily take into account all essential geomorphological 4-7 and socio-economic system feedbacks 8 . Here we present an integrated global modelling approach that considers both the ability of coastal wetlands to build up vertically by sediment accretion, and the accommodation space, namely, the vertical and lateral space available for fine sediments to accumulate and be colonized by wetland vegetation. We use this approach to assess global-scale changes in coastal wetland area in response to global sea-level rise and anthropogenic coastal occupation during the twenty-first century. On the basis of our simulations, we find that, globally, rather than losses, wetland gains of up to 60 per cent of the current area are possible, if more than 37 per cent (our upper estimate for current accommodation space) of coastal wetlands have sufficient accommodation space, and sediment supply remains at present levels. In contrast to previous studies 1-3 , we project that until 2100, the loss of global coastal wetland area will range between 0 and 30 per cent, assuming no further accommodation space in addition to current levels. Our simulations suggest that the resilience of global wetlands is primarily driven by the availability of accommodation space, which is strongly influenced by the building of anthropogenic infrastructure in the coastal zone and such infrastructure is expected to change over the twenty-first century. Rather than being an inevitable consequence of global sea-level rise, our findings indicate that large-scale loss of coastal wetlands might be avoidable, if sufficient additional accommodation space can be created through careful nature-based adaptation solutions to coastal management.

Historical reanalyses that span more than a century are needed for a wide range of studies, from understanding large‐scale climate trends to diagnosing the impacts of individual historical extreme weather events. The Twentieth Century Reanalysis (20CR) Project is an effort to fill this need. It is supported by the National Oceanic and Atmospheric Administration (NOAA), the Cooperative Institute for Research in Environmental Sciences (CIRES), and the U.S. Department of Energy (DOE), and is facilitated by collaboration with the international Atmospheric Circulation Reconstructions over the Earth initiative. 20CR is the first ensemble of sub‐daily global atmospheric conditions spanning over 100 years. This provides a best estimate of the weather at any given place and time as well as an estimate of its confidence and uncertainty. While extremely useful, version 2c of this dataset (20CRv2c) has several significant issues, including inaccurate estimates of confidence and a global sea level pressure bias in the mid‐19th century. These and other issues can reduce its effectiveness for studies at many spatial and temporal scales. Therefore, the 20CR system underwent a series of developments to generate a significant new version of the reanalysis. The version 3 system (NOAA‐CIRES‐DOE 20CRv3) uses upgraded data assimilation methods including an adaptive inflation algorithm; has a newer, higher‐resolution forecast model that specifies dry air mass; and assimilates a larger set of pressure observations. These changes have improved the ensemble‐based estimates of confidence, removed spin‐up effects in the precipitation fields, and diminished the sea‐level pressure bias. Other improvements include more accurate representations of storm intensity, smaller errors, and large‐scale reductions in model bias. The 20CRv3 system is comprehensively reviewed, focusing on the aspects that have ameliorated issues in 20CRv2c. Despite the many improvements, some challenges remain, including a systematic bias in tropical precipitation and time‐varying biases in southern high‐latitude pressure fields.

We demonstrate a powerful and practical method of characterizing interactions in fine magnetic particle systems utilizing a class of hysteresis curves known as first order reversal curves. This method is tested on samples of highly dispersed magnetic particles, where it leads to a more detailed understanding of interactions than has previously been possible. In a quantitative comparison between this method and the δM method, which is based on the Wohlfarth relation, our method provides a more precise measure of the strength of the interactions. Our method also has the advantage that it can be used to decouple the effects of the mean interaction field from the effects of local interaction field variance.

The rapid anthropogenic climate change that is being experienced in the early twenty-first century is intimately entwined with the health and functioning of the biosphere. Climate change is impacting ecosystems through changes in mean conditions and in climate variability, coupled with other associated changes such as increased ocean acidification and atmospheric carbon dioxide concentrations. It also interacts with other pressures on ecosystems, including degradation, defaunation and fragmentation. There is a need to understand the ecological dynamics of these climate impacts, to identify hotspots of vulnerability and resilience and to identify management interventions that may assist biosphere resilience to climate change. At the same time, ecosystems can also assist in the mitigation of, and adaptation to, climate change. The mechanisms, potential and limits of such nature-based solutions to climate change need to be explored and quantified. This paper introduces a thematic issue dedicated to the interaction between climate change and the biosphere. It explores novel perspectives on how ecosystems respond to climate change, how ecosystem resilience can be enhanced and how ecosystems can assist in addressing the challenge of a changing climate. It draws on a Royal Society-National Academy of Sciences Forum held in Washington DC in November 2018, where these themes and issues were discussed. We conclude by identifying some priorities for academic research and practical implementation, in order to maximize the potential for maintaining a diverse, resilient and well-functioning biosphere under the challenging conditions of the twenty-first century. This article is part of the theme issue 'Climate change and ecosystems: threats, opportunities and solutions'.

Abstract Submarine sediment density flows are one of the most important processes for moving sediment across our planet, yet they are extremely difficult to monitor directly. The speed of long run‐out submarine density flows has been measured directly in just five locations worldwide and their sediment concentration has never been measured directly. The only record of most density flows is their sediment deposit. This article summarizes the processes by which density flows deposit sediment and proposes a new single classification for the resulting types of deposit. Colloidal properties of fine cohesive mud ensure that mud deposition is complex, and large volumes of mud can sometimes pond or drain‐back for long distances into basinal lows. Deposition of ungraded mud (T E‐3 ) most probably finally results from en masse consolidation in relatively thin and dense flows, although initial size sorting of mud indicates earlier stages of dilute and expanded flow. Graded mud (T E‐2 ) and finely laminated mud (T E‐1 ) most probably result from floc settling at lower mud concentrations. Grain‐size breaks beneath mud intervals are commonplace, and record bypass of intermediate grain sizes due to colloidal mud behaviour. Planar‐laminated (T D ) and ripple cross‐laminated (T C ) non‐cohesive silt or fine sand is deposited by dilute flow, and the external deposit shape is consistent with previous models of spatial decelerating (dissipative) dilute flow. A grain‐size break beneath the ripple cross‐laminated (T C ) interval is common, and records a period of sediment reworking (sometimes into dunes) or bypass. Finely planar‐laminated sand can be deposited by low‐amplitude bed waves in dilute flow (T B‐1 ), but it is most likely to be deposited mainly by high‐concentration near‐bed layers beneath high‐density flows (T B‐2 ). More widely spaced planar lamination (T B‐3 ) occurs beneath massive clean sand (T A ), and is also formed by high‐density turbidity currents. High‐density turbidite deposits (T A , T B‐2 and T B‐3 ) have a tabular shape consistent with hindered settling, and are typically overlain by a more extensive drape of low‐density turbidite (T D and T C ,). This core and drape shape suggests that events sometimes comprise two distinct flow components. Massive clean sand is less commonly deposited en masse by liquefied debris flow (D CS ), in which case the clean sand is ungraded or has a patchy grain‐size texture. Clean‐sand debrites can extend for several tens of kilometres before pinching out abruptly. Up‐current transitions suggest that clean‐sand debris flows sometimes form via transformation from high‐density turbidity currents. Cohesive debris flows can deposit three types of ungraded muddy sand that may contain clasts. Thick cohesive debrites tend to occur in more proximal settings and extend from an initial slope failure. Thinner and highly mobile low‐strength cohesive debris flows produce extensive deposits restricted to distal areas. These low‐strength debris flows may contain clasts and travel long distances (D M‐2 ), or result from more local flow transformation due to turbulence damping by cohesive mud (D M‐1 ). Mapping of individual flow deposits (beds) emphasizes how a single event can contain several flow types, with transformations between flow types. Flow transformation may be from dilute to dense flow, as well as from dense to dilute flow. Flow state, deposit type and flow transformation are strongly dependent on the volume fraction of cohesive fine mud within a flow. Recent field observations show significant deviations from previous widely cited models, and many hypotheses linking flow type to deposit type are poorly tested. There is much still to learn about these remarkable flows.

Autonomous Underwater Vehicles (AUVs) have a wide range of applications in marine geoscience, and are increasingly being used in the scientific, military, commercial, and policy sectors. Their ability to operate autonomously of a host vessel makes them well suited to exploration of extreme environments, from the world's deepest hydrothermal vents to beneath polar ice sheets. They have revolutionised our ability to image the seafloor, providing higher resolution seafloor mapping data than can be achieved from surface vessels, particularly in deep water. This contribution focuses on the major advances in marine geoscience that have resulted from AUV data. The primary applications are i) submarine volcanism and hydrothermal vent studies, ii) mapping and monitoring of low-temperature fluid escape features and chemosynthetic ecosystems, iii) benthic habitat mapping in shallow- and deep-water environments, and iv) mapping of seafloor morphological features (e.g. bedforms generated beneath ice or sediment-gravity flows). A series of new datasets is presented that highlight the growing versatility of AUVs for marine geoscience studies, including i) multi-frequency acoustic imaging of trawling impacts on deep-water coral mounds, iii) collection of high-resolution seafloor photomosaics at abyssal depths, and iii) velocity measurements of active submarine density flows. Future developments in AUV technology of potential relevance to marine geoscience include new vehicles with enhanced hovering, long endurance, extreme depth, or rapid response capabilities, while development of new sensors will further expand the range of geochemical parameters that can be measured.

We present a new surface-atmospheric dataset for driving ocean–sea-ice models based on Japanese 55-year atmospheric reanalysis (JRA-55), referred to here as JRA55-do. The JRA55-do dataset aims to replace the CORE interannual forcing version 2 (hereafter called the CORE dataset), which is currently used in the framework of the Coordinated Ocean-ice Reference Experiments (COREs) and the Ocean Model Intercomparison Project (OMIP). A major improvement in JRA55-do is the refined horizontal grid spacing ( ∼ 55 km) and temporal interval (3 hr). The data production method for JRA55-do essentially follows that of the CORE dataset, whereby the surface fields from an atmospheric reanalysis are adjusted relative to reference datasets. To improve the adjustment method, we use high-quality products derived from satellites and from several other atmospheric reanalysis projects, as well as feedback on the CORE dataset from the ocean modelling community. Notably, the surface air temperature and specific humidity are adjusted using multi-reanalysis ensemble means. In JRA55-do, the downwelling radiative fluxes and precipitation, which are affected by an ambiguous cloud parameterisation employed in the atmospheric model used for the reanalysis, are based on the reanalysis products. This approach represents a notable change from the CORE dataset, which imported independent observational products. Consequently, the JRA55-do dataset is more self-contained than the CORE dataset, and thus can be continually updated in near real-time. The JRA55-do dataset extends from 1958 to the present, with updates expected at least annually. This paper details the adjustments to the original JRA-55 fields, the scientific rationale for these adjustments, and the evaluation of JRA55-do. The adjustments successfully corrected the biases in the original JRA-55 fields. The globally averaged features are similar between the JRA55-do and CORE datasets, implying that JRA55-do can suitably replace the CORE dataset for use in driving global ocean–sea-ice models.