Search NASA⌕ Search

SEARCH · Search NASA

Results for “Global ocean state estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

The Development of the New GEOS-MITgcm Atmosphere-Ocean Model for Coupled Data Assimilation System

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8 deg) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12 deg) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

A New Atmosphere-Ocean Model for Studying Air-Sea Interactions and Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8°) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12°) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

The Development of the New GEOS-MITgcm Atmosphere-Ocean Model for Coupled Data Assimilation System

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8°) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12°) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

Development of the GEOS-MIT gcm Atmosphere-Ocean Model for Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MIT-GCM (General Circulation Model)) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (approximately 1/8 degrees) configuration of the GEOS atmosphere to a latitude-longitude-cap-1080 (approximately 1/12 degrees) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗

Synthesis and Assimilation Systems - Essential Adjuncts to the Global Ocean Observing System

Ocean assimilation systems synthesize diverse in situ and satellite data streams into four-dimensional state estimates by combining the various observations with the model. Assimilation is particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Developments in assimilation and in the observing system have advanced our understanding and prediction of ocean variations at mesoscale and climate scales. Use of these systems for assessing the observing system helps identify the strengths of each observation type. Results indicate that the ocean remains under-sampled and that further improvements in the observing system are needed. Prospects for future advances lie in improved models and better estimates of error statistics for both models and observations. Future developments will be increasingly towards consistent analyses across components of the Earth system. However, even today ocean synthesis and assimilation systems are providing products that are useful for many applications and should be considered an integral part of the global ocean observing and information system.

Rienecker, Michele M.↗

ResORR: A Globally Scalable and Satellite Data-Driven Algorithm for River Flow Regulation Due to Reservoir Operations

We propose a globally scalable algorithm, ResORR (Reservoir Operations driven River Regulation), to predict regulated river flow and tested it over the heavily regulated basin of the Cumberland River in the US. ResORR was found able to model regulated river flow due to upstream reservoir operations of the Cumberland River. Over a mountainous basin dominated by high rainfall, ResORR was effective in capturing extreme flooding modified by upstream hydropower dam operations. On average, ResORR improved regulated river flow simulation by more than 50% across all performance metrics when compared to a hydrologic model without a regulation module. ResORR is a timely software algorithm for understanding human regulation of surface water as satellite-estimated reservoir state is expected to improve globally with the recently launched Surface Water and Ocean Topography (SWOT) mission.

River Regulation↗

The number of past and future regenerations of iron in the oceanand its intrinsic fertilization efficiency

Iron fertilization is explored by tracking dissolved iron (DFe) through its life cycle from injection by aeolian, sedimentary, and hydrothermal sources (birth) to burial in the sediments (death). We develop new diagnostic equations that count iron and phosphate regenerations with each passage through the biological pump and partition the ocean's DFe concentration according to the number of its past or future regenerations. We apply these diagnostics to a family of data-constrained estimates of the iron cycle with sources σ tot in the range 1.9 – 41 Gmol yr -1 . We find that for states with σ tot > 7 Gmol yr -1 , 50 % or more of the DFe inventory has not been regenerated in the past and 85 % or more will not be regenerated in the future. The globally averaged mean number of past or future regenerations scales with the bulk iron lifetime τ ~ σ tot - 1 and has a range of 0.05 – 2.2 for past and 0.01 – 1.4 for future regenerations. Memory of birth location fades rapidly with each regeneration, and DFe regenerated more than approximately five times is found in a pattern shaped by Southern Ocean nutrient trapping. We quantify the intrinsic fertilization efficiency of the unperturbed system at any point r in the ocean as the global export production resulting from the DFe at r per iron molecule. We show that this efficiency is closely related to the mean number of future regenerations that the iron will experience. At the surface, the intrinsic fertilization efficiency has a global mean in the range 0.7 – 7 mol P (mmol Fe) -1 across our family of state estimates and is largest in the central tropical Pacific, with the Southern Ocean having comparable importance only for high-iron-source scenarios.

54 ENVIRONMENTAL SCIENCES↗

Climate and atmospheric deposition drive the inter-annual variability and long-term trend of dissolved organic carbon flux in the conterminous United States

The lateral flux of dissolved organic carbon (DOC) from soils to inland waters and ultimately to the ocean represents a fundamental component of the global carbon cycle. To estimate the DOC flux, we developed an empirical terrestrial-aquatic DOC fluxes model (TAF-DOC). TAF-DOC incorporates various environmental factors (e.g., meteorology, sulfur, and nitrogen deposition) that to-date have not been comprehensively considered or well-represented in existing modeling frameworks. TAF-DOC was applied to estimate spatial-temporal patterns of DOC flux and potential fates across the conterminous United States during the 1985 to 2018 time period. Our results suggest that TAF-DOC successfully characterized spatial-temporal of DOC flux. As expected, the interannual pattern of DOC flux was strongly regulated by precipitation, but the long-term trend was significantly influenced by the rate of atmospheric wet sulfur deposition. From 1985 to 2018, TAF-DOC estimated DOC loading from terrestrial to aquatic ecosystems in the conterminous United States to be 33.5 ± 2.2 TgC per year, which was roughly 0.39–0.49% of total soil organic carbon stock estimates. The dominant fate of terrestrially-derived DOC was delivery to the coastal ocean in riverine export (41%), with another 21% buried in sediment and the remaining 12.8 ± 0.4 TgC per year (38%) returned to the atmosphere through outgassing from inland waters. Finally, assuming the quantities of DOC sediment burial and export to the ocean as an annual sink of terrestrially-derived carbon, budget inventories and models that do not account for DOC flux in the conterminous United States will underestimate the net annual carbon sink by as much as 5.5–6.4%.

54 ENVIRONMENTAL SCIENCES↗

Applications of the GPCP Satellite-Gauge Merged Precipitation Data Set to the Goals of PACS and GCIP

Recently, under the direction of GEWEX/GCIP there has been an attempt to link North American monsoon rainfall to precipitation in the Great Plains for the purposes of predictability on the seasonal time scale. Rain gauge analyses have been used to show continental-scale relationships. However, the interannual variability of warm season precipitation over the United States is linked to the east Pacific climate system. The Global Precipitation Climatology Project (GPCP) satellite estimates of precipitation over oceans add a vital piece of information towards the goals of PACS and GCIP. This preliminary study will begin by focusing on the East Pacific sector. El-Nino-Southern Ossilation (ENSO) modulates the position and strength of the Inter-Tropical Convergence Zone (ITCZ), which affects the precipitation source for the Mexican monsoon. For example, the five driest summers in coastal Mexico occurred during the five strongest El Ninos (as determined by positive values of Nino 3.4) since 1979. Previous work has investigated lag-lead relationships between the North American monsoon and continental rainfall. GPCP has the potential to extend these relationships to a global scale. Finally data sets with a finer time and space scale, including Tropical Rainfall Measurement Mission (TRMM) products, will be used to examine wave dynamics. Kelvin waves, in particular, may be trigger mechanisms for the onset of the North American monsoon.

Curtis, Scott↗

Evaluating the carbon capture potential of industrial waste as a feedstock for enhanced weathering

Abstract Limiting anthropogenic global climate warming since the start of the industrial period to less than 2 °C will very likely require both deep and rapid reductions in anthropogenic greenhouse gas emissions and a range of approaches toward carbon dioxide removal (CDR). One prominent CDR approach is enhanced weathering (EW), in which crushed silicate rock is applied on land or in the open ocean to accelerate natural weathering processes that absorb carbon dioxide from Earth’s ocean–atmosphere system. However, in addition to a range of potential environmental, socioeconomic, and ethical issues associated with this pathway, bottlenecks in feedstock sourcing represent a key barrier for deployment of EW at scale. Here, we evaluate the potential of silicate wastes produced from industrial processes—such as steel slag and cement waste—as feedstocks for the EW process. An empirical model that links industrial alkaline waste production to gross domestic product at purchase power parity is developed to forecast waste production in the alternative futures described by the shared socioeconomic pathway (SSP) framework. By incorporating these results into an intermediate-complexity Earth system model, we also explore the impacts of EW using industrial waste on changes to global temperature, ocean pH, and ocean aragonite saturation state, while also quantifying overall CDR efficiency through the end of the century. We estimate a maximum cumulative end-of-century capture potential of ∼400 GtCO 2 for industrial waste, which could represent a significant fraction of the projected CDR requirement of many mitigation scenarios in the SSP framework. However, feedstock-dependent environmental impacts and the technoeconomics of feedstock redistribution may ultimately limit deployment scope.

Xu, Pengxiao (ORCID:0009000633724293)↗

Adaptive Error Estimation in Linearized Ocean General Circulation Models

Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Chechelnitsky, Michael Y.↗

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗

Analyzing the Effect of Ocean Internal Variability on Depth-Integrated Steric Sea-Level Rise Trends Using a Low-Resolution CESM Ensemble

Ocean heat uptake is a key indicator of climate change, in part because it contributes to sea-level rise. Quantifying the uncertainties surrounding ocean heat uptake and sea-level rise are important in assessing climate-related risks. Here, comprehensive global climate model ensembles are used to evaluate uncertainties surrounding decadal trends in depth-integrated global steric sea-level rise due to thermal expansion of the ocean. Results are presented against observational estimates, which are used as a guide to the state of recent literature. The first ensemble uses the Community Earth System Model (CESM), which samples the effects of internal variability within the coupled Earth system including contributions from the sub-surface ocean. We compare and contrast these results with an ensemble based on the Coupled Model Intercomparison Project Phase 5 (CMIP5), which samples the combined effects of structural model differences and internal variability. The effects of both internal variability and structural model differences contribute substantially to uncertainties in modeled steric sea-level trends for recent decades, and the magnitude of these effects varies with depth. The 95% range in total sea-level rise trends across the CESM ensemble is 0.151 mm·year –1 for 1957–2013, while this range is 0.895 mm·year –1 for CMIP5. These ranges increase during the more recent decade of 2005–2015 to 0.509 mm·year –1 and 1.096 mm·year –1 for CESM and CMIP5, respectively. The uncertainties are amplified for regional assessments, highlighting the importance of both internal variability and structural model differences when considering uncertainties surrounding modeled sea-level trends. Results can potentially provide useful constraints on estimations of global and regional sea-level variability, in particular for areas with few observations such as the deep ocean and Southern Hemisphere.

54 ENVIRONMENTAL SCIENCES↗

Long-Term Measurements Show Little Evidence for Large Increases in Total U.S. Methane Emissions Over the Past Decade

Recent studies show conflicting estimates of trends in methane (CH 4 ) emissions from oil and natural gas (ONG) operations in the United States. We analyze atmospheric CH 4 measurements from 20 North American sites in the National Oceanic and Atmospheric Administration Global Greenhouse Gas Reference Network and determined trends for 2006–2015. Using CH 4 vertical gradients as an indicator of regional surface emissions, we find no significant increase in emissions at most sites and modest increases at three sites heavily influenced by ONG activities. Our estimated increases in North American ONG CH 4 emissions (on average approximately 3.4 ± 1.4 %/year for 2006–2015, ±σ) are much smaller than estimates from some previous studies and below our detection threshold for total emissions increases at the east coast sites that are sensitive to U.S. outflows. We also find an increasing trend in ethane/methane emission ratios, which has resulted in major overestimation of oil and gas emissions trends in some previous studies.

54 ENVIRONMENTAL SCIENCES↗

Direct Radiative Effect of Aerosols Based on PARASOL and OMI Satellite Observations

Accurate portrayal of the aerosol characteristics is crucial to determine aerosol contribution to the Earth's radiation budget. We employ novel satellite retrievals to make a new measurement-based estimate of the shortwave direct radiative effect of aerosols (DREA), both over land and ocean. Global satellite measurements of aerosol optical depth, single-scattering albedo (SSA), and phase function from PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) are used in synergy with OMI (Ozone Monitoring Instrument) SSA. Aerosol information is combined with land-surface bidirectional reflectance distribution function and cloud characteristics from MODIS (Moderate Resolution Imaging Spectroradiometer) satellite products. Eventual gaps in observations are filled with the state-of-the-art global aerosol model ECHAM5-HAM2. It is found that our estimate of DREA is largely insensitive to model choice. Radiative transfer calculations show that DREA at top-of-atmosphere is −4.6 +/- 1.5 W/sq m for cloud-free and −2.1 +/- 0.7 W/sq m for all-sky conditions, during year 2006. These fluxes are consistent with, albeit generally less negative over ocean than, former assessments. Unlike previous studies, our estimate is constrained by retrievals of global coverage SSA, which may justify different DREA values. Remarkable consistency is found in comparison with DREA based on CERES (Clouds and the Earth's Radiant Energy System) and MODIS observations.

Globalsatellite measurements↗

Predictability of Ocean Heat Content From Electrical Conductance

Ocean heat content (OHC) is a key climate variable that needs to be monitored to know how Earth's energy imbalance is changing, yet observing OHC remains a challenge. The present study examines whether a depth integral of the ocean's electrical conductivity (“conductance”), which may be inferred from both in situ methods and satellite magnetometers over the global ocean, could help monitor OHC. The ocean's electrical conductivity locally depends on temperature, salinity, and pressure, but it is not as well known how the conductance depends on OHC and ocean salt content. By examining the output of an ocean state estimate shown to agree well with observations that have not been assimilated, this study evaluates the fundamental limitations of using perfectly known ocean conductance to predict OHC, rather than the challenges associated with accounting for observational error. It is found that the ocean's conductance and OHC fields are nonlinearly related but nevertheless highly correlated. A statistical framework tends to predict OHC more accurately than ocean salt content from ocean conductance in regions where conductivity is more sensitive to salinity than temperature. The annually (bidecadally) averaged OHC can be predicted from a combination of conductance and depth‐averaged conductivity ocean fields to within nearly 0.1% (1%) error globally and even more accurately in many poorly observed (e.g., ice‐covered) regions. Practical application of this statistical framework to monitor OHC requires examination of the effect of uncertainties in the observed bathymetry and ocean conductance, which vary with application.

Trossman, D. S.↗

Improvements in the Uncertainty Model in the Goddard Institute for Space Studies Surface Temperature (GISTEMP) Analysis

We outline a new and improved uncertainty analysis for the Goddard Institute of Space Studies (GISS) Surface Temperature product version 4 (GISTEMP v4). Historical spatial variations in surface temperature anomalies are derived from historical weather station data and ocean data from ships, buoys and other sensors. Uncertainties arise from measurement uncertainty, changes in spatial coverage of the station record, and systematic biases due to technology shifts and land cover changes. Previously published uncertainty estimates for GISTEMP included only the effect of incomplete station coverage. Here, we update this term using currently available spatial distributions of source data, state‐of‐the‐art reanalyses and incorporate independently derived estimates for ocean data processing, station homogenization and other structural biases. The resulting 95% uncertainties are near 0.05° C in the global annual mean for the last 50 years, and increase going back further in time reaching 0.15° C in 1880. In addition, we quantify the benefits and inherent uncertainty due to the GISTEMP interpolation and averaging method. We use the total uncertainties to estimate the probability for each record year in the GISTEMP to actually be the true record year (to that date), and conclude with 86% likelihood that 2016 was indeed the hottest year of the instrumental period (so far).

Nathan J L Lenssen↗