Search Results

You are looking at 1 - 10 of 21 items for

  • Author or Editor: Lijing Cheng x
  • Refine by Access: All Content x
Clear All Modify Search
Viktor Gouretski and Lijing Cheng

Abstract

A homogeneous, consistent, high-quality in situ temperature dataset covering some decades in time is crucial for the detection of climate changes in the ocean. For the period from 1940 to the present, this study investigates the data quality of temperature profiles from mechanical bathythermographs (MBT) by comparing these data with reference data obtained from Nansen bottle casts and conductivity–temperature–depth (CTD) profilers. This comparison reveals significant systematic errors in MBT measurements. The MBT bias is as large as 0.2°C before 1980 on the global average and reduces to less than 0.1°C after 1980. A new empirical correction scheme for MBT data is derived, where the MBT correction is country, depth, and time dependent. Comparison of the new MBT correction scheme with three schemes proposed earlier in the literature suggests a better performance of the new schemes. The reduction of the biases increases the homogeneity of the global ocean database being mostly important for climate change–related studies, such as the improved estimation of the ocean heat content changes.

Free access
Lijing Cheng and Jiang Zhu

Abstract

A complete map of the ocean subsurface temperature is essential for monitoring aspects of climate change such as the ocean heat content (OHC) and sea level changes and for understanding the dynamics of the ocean/climate variation. However, global observations have not been available in the past, so a mapping strategy is required to fill the data gaps. In this study, an advanced mapping method is proposed to reconstruct the historical ocean subsurface (0–700 m) temperature field from 1940 to 2014 by using ensemble optimal interpolation with a dynamic ensemble (EnOI-DE) approach and a multimodel ensemble of phase 5 of the Coupled Model Intercomparison Project (CMIP5) historical and representative concentration pathway 4.5 simulations. The reconstructed field is a combination of two parts: a first guess provided by the ensemble mean of CMIP5 models and an adjustment by minimizing the analysis error with the assistance of error covariance determined by the CMIP5 models. The uncertainty of the field can also be assessed. This new approach was evaluated using a series of tests, including subsample tests by using data from the Argo period, idealized tests by specifying a truth field from the models, and withdrawn-data tests by removing 20% of the observations for validation. In addition, the authors showed that the ocean mean state, long-term trends, and interannual and decadal variability are all well represented. Furthermore, the most significant benefit of this method is to provide an improved estimate of the long-term historical OHC changes since 1940, which have important implications for Earth’s energy budget.

Full access
Lijing Cheng and Jiang Zhu

Abstract

The choice of climatology is an essential step in calculating the key climate indicators, such as historical ocean heat content (OHC) change. The anomaly field is required during the calculation and is obtained by subtracting the climatology from the absolute field. The climatology represents the ocean spatial variability and seasonal circle. This study found a considerable weaker long-term trend when historical climatologies (constructed by using historical observations within a long time period, i.e., 45 yr) were used rather than Argo-period climatologies (i.e., constructed by using observations during the Argo period, i.e., since 2004). The change of the locations of the observations (horizontal sampling) during the past 50 yr is responsible for this divergence, because the ship-based system pre-2000 has insufficient sampling of the global ocean, for instance, in the Southern Hemisphere, whereas this area began to achieve full sampling in this century by the Argo system. The horizontal sampling change leads to the change of the reference time (and reference OHC) when the historical-period climatology is used, which weakens the long-term OHC trend. Therefore, Argo-period climatologies should be used to accurately assess the long-term trend of the climate indicators, such as OHC.

Full access
Lijing Cheng and Jiang Zhu

Abstract

Assessment of the upper-ocean (0–700 m) heat content (OHC) is a key task for monitoring climate change. However, irregular spatial and temporal distribution of historical subsurface observations has induced uncertainties in OHC estimation. In this study, a new source of uncertainties in calculating OHC due to the insufficiency of vertical resolution in historical ocean subsurface temperature profile observations was diagnosed. This error was examined by sampling a high-vertical-resolution climatological ocean according to the depth intervals of in situ subsurface observations, and then the error was defined as the difference between the OHC calculated by subsampled profiles and the OHC of the climatological ocean. The obtained resolution-induced error appeared to be cold in the upper 100 m (with a peak of approximately −0.1°C), warm within 100–700 m (with a peak of ~0.1°C near 180 m), and warm when averaged over 0–700-m depths (with a global average of ~0.01°–0.025°C, ~1–2.5 × 1022 J). Geographically, it showed a warm bias within 30°S–30°N and a cold bias at higher latitudes in both hemispheres, the sign of which depended on the concave or convex shape of the vertical temperature profiles. Finally, the authors recommend maintaining an unbiased observation system in the future: a minimal vertical depth bin of 5% of the depth was needed to reduce the vertical-resolution-induced bias to less than 0.005°C on global average (equal to Argo accuracy).

Full access
John P. Abraham, Rebecca Cowley, and Lijing Cheng

Abstract

A very large portion of the historical information on ocean temperatures has been measured using expendable bathythermograph (XBT) devices. For decades, these devices provided the majority of global information. It is, therefore, important to quantify their accuracy and identify biases in this important historical dataset. Here, calculations are made of the influence of water temperature on the rate of descent of the XBT devices into the ocean waters. In colder regions, the larger viscosity of the water is expected to cause a greater drag force on the device, which would slow the descent. It was found through computational fluid dynamic models that the impact of temperature and viscosity on the probe descent is approximately 2.2% for water temperatures that range from 0° to 27°C. Probe-specific temperature-dependent fall rate equations were applied to 269 collocated XBT/conductivity–temperature–depth (CTD) measurements from two different research cruises. It was found that the probe-specific descent equations were an improvement over the uncorrected method. Next, in an effort to automate the process, the fall rate coefficients were related to the topmost measured temperature in the water column. With this relationship, comparisons were made between the probe-specific descent calculations and 2937 high-resolution XBT–CTD pairs. It was found that again, the new methodology outperformed the standard fall rate equation. The new method was also compared with an independent correction method that was previously published. It was found that both new methods were improvements upon the industry-standard fall rate calculation. Subsequent calculations using the top-100-m water temperature were performed and were found to be statistically insignificant compared to the proposed simplified method.

Full access
Andrea Storto, Lijing Cheng, and Chunxue Yang

Abstract

Recent estimates of the global warming rates suggest that approximately 9% of the Earth’s excess heat is cumulated in the deep and abyssal oceans (below 2000 m depth) during the last two decades. Such estimates assume stationary trends deducted as long-term rates. In order to reassess the deep ocean warming and potentially shed light on its inter-annual variability, we formulate the balance between the Earth’s Energy Imbalance (EEI), the steric sea level and the ocean heat content (OHC), at yearly time scales during the 2003-2018 period, as a variational problem. The solution is achieved through variational minimization, merging together observational data from top-of-atmosphere EEI, inferred from CERES, steric sea level estimates from altimetry minus gravimetry, and upper ocean heat content estimates from in-situ platforms (mostly Argo floats). Global ocean reanalyses provide background error covariances for the OHC analysis. The analysis indicates a 2000m-bottom warming of 0.08 ± 0.04 W m−2 for the period 2003-2018, equal to 13% of the total ocean warming (0.62 ± 0.08 W m−2), slightly larger than previous estimates but consistent within the error bars. The analysis provides a fully consistent optimized solution also for the steric sea level and EEI. Moreover, the simultaneous use of the different heat budget observing networks is able to decrease the analysis uncertainty with respect to the observational one, for all observation types and especially for the 0-700m OHC and steric sea level (more than 12% reduction). The sensitivity of the analysis to the choice of the background timeseries is proved insignificant.

Restricted access
Lijing Cheng, Jiang Zhu, Franco Reseghetti, and Qingping Liu

Abstract

A new technique to estimate three major biases of XBT probes (improper fall rate, start-up transient, and pure temperature error) has been developed. Different from the well-known and standard “temperature error free” differential method, the new method analyses temperature profiles instead of vertical gradient temperature profiles. Consequently, it seems to be more noise resistant because it uses the integral property over the entire vertical profile instead of gradients. Its validity and robustness have been checked in two ways. In the first case, the new integral technique and the standard differential method have been applied to a set of simulated XBT profiles having a known fall-rate equation to which various combinations of pure temperature errors, random errors, and spikes have been added for the sake of this simulation. Results indicated that the single pure temperature error has little impact on the fall-rate coefficients for both methods, whereas with the added random error and spikes the simulation leads to better results with the new integral technique than with the standard differential method. In the second case, two sets of profiles from actual XBT versus CTD comparisons, collected near Barbados in 1990 and in the western Mediterranean (2003–04 and 2008–09), have been used. The individual fall-rate coefficients and start-up transient for each XBT profile, along with the overall pure temperature correction, have been calculated for the XBT profiles. To standardize procedures and to improve the terms of comparison, the individual start-up transient estimated by the integral method was also assigned and included in calculations with the differential method. The new integral method significantly reduces both the temperature difference between XBT and CTD profiles and the standard deviation. Finally, the validity of the mean fall-rate coefficients and the mean start-up transient, respectively, for DB and T7 probes as precalculated equations was verified. In this case, the temperature difference is reduced to less than 0.1°C for both datasets, and it randomly distributes around the null value. In addition, the standard deviation on depth values is largely reduced, and the maximum depth error computed with the datasets near Barbados is within 1.1% of its real value. Results also indicate that the integral method has a good performance mainly when applied to profiles in regions with either a very large temperature gradient, at the thermocline or a very small one, toward the bottom.

Full access
Lijing Cheng, Jiang Zhu, Rebecca Cowley, Tim Boyer, and Susan Wijffels

Abstract

Systematic biases in historical expendable bathythermograph (XBT) data are examined using two datasets: 4151 XBT–CTD side-by-side pairs from 1967 to 2011 and 218 653 global-scale XBT–CTD pairs (within one month and 1°) extracted from the World Ocean Database 2009 (WOD09) from 1966 to 2010. Using the side-by-side dataset, it was found that both the pure thermal bias and the XBT fall rate (from which the depth of observation is calculated) increase with water temperature. Correlations between the terminal velocity A and deceleration B terms of the fall-rate equation (FRE) and between A and the offset from the surface terms are obtained, with A as the dominant term in XBT fall-rate behavior. To quantify the time variation of the XBT fall-rate and pure temperature biases, global-scale XBT–CTD pairs are used. Based on the results from the two datasets, a new correction scheme for historical XBT data is proposed for nine independent probe-type groups. The scheme includes corrections for both temperature and depth records, which are all variable with calendar year, water temperature, and probe type. The results confirm those found in previous studies: a slowing in fall rate during the 1970s and 2000s and the large pure thermal biases during 1970–85. The performance of nine different correction schemes is compared. After the proposed corrections are applied to the XBT data in the WOD09 dataset, global ocean heat content from 1967 to 2010 is reestimated.

Full access
Rebecca Cowley, Susan Wijffels, Lijing Cheng, Tim Boyer, and Shoichi Kizu

Abstract

Because they make up 56% of ocean temperature profile data between 1967 and 2001, quantifying the biases in expendable bathythermograph (XBT) data is fundamental to understanding the evolution of the planetary energy and sea level budgets over recent decades. The nature and time history of these biases remain in dispute and dominate differences in analyses of the history of ocean warming. A database of over 4100 side-by-side deployments of XBTs and conductivity–temperature–depth (CTD) data has been assembled, and this unique resource is used to characterize and separate out the pure temperature bias from depth error in a way that was not previously possible. Two independent methods of bias extraction confirm that the results are robust to bias model and fitting method. It was found that there is a pure temperature bias in Sippican probes of ~0.05°C, independent of depth. The temperature bias has a time dependency, being larger (~0.1°C) in the earlier analog acquisition era and being likely due to changes in recorder type. Large depth errors are found in the 1970s–80s in shallower-measuring Sippican T4/T6 probe types, but the deeper-measuring Sippican T7/Deep Blue (DB) types have no error during this time. The Sippican T7/DB fall rate slows from ~1990 onward. It is found that year-to-year variations in fall rate have a bigger effect on corrections to the global XBT database than do any small effects of ocean temperature on fall rate. This study has large implications for the future development of better schemes to correct the global historical XBT archive.

Full access
Lijing Cheng, Kevin E. Trenberth, John T. Fasullo, Michael Mayer, Magdalena Balmaseda, and Jiang Zhu

Abstract

As the strongest interannual perturbation to the climate system, El Niño–Southern Oscillation (ENSO) dominates the year-to-year variability of the ocean energy budget. Here we combine ocean observations, reanalyses, and surface flux data with Earth system model simulations to obtain estimates of the different terms affecting the redistribution of energy in the Earth system during ENSO events, including exchanges between ocean and atmosphere and among different ocean basins, and lateral and vertical rearrangements. This comprehensive inventory allows better understanding of the regional and global evolution of ocean heat related to ENSO and provides observational metrics to benchmark performance of climate models. Results confirm that there is a strong negative ocean heat content tendency (OHCT) in the tropical Pacific Ocean during El Niño, mainly through enhanced air–sea heat fluxes Q into the atmosphere driven by high sea surface temperatures. In addition to this diabatic component, there is an adiabatic redistribution of heat both laterally and vertically (0–100 and 100–300 m) in the tropical Pacific and Indian oceans that dominates the local OHCT. Heat is also transported and discharged from 20°S–5°N into off-equatorial regions within 5°–20°N during and after El Niño. OHCT and Q changes outside the tropical Pacific Ocean indicate the ENSO-driven atmospheric teleconnections and changes of ocean heat transport (i.e., Indonesian Throughflow). The tropical Atlantic and Indian Oceans warm during El Niño, partly offsetting the tropical Pacific cooling for the tropical oceans as a whole. While there are distinct regional OHCT changes, many compensate each other, resulting in a weak but robust net global ocean cooling during and after El Niño.

Open access