Search Results
You are looking at 1 - 10 of 20 items for
- Author or Editor: Tim Boyer x
- Refine by Access: All Content x
On the Consistency of the Bottle and CTD Profile DataAbstract
Nansen bottle casts served as the main oceanographic instrumentation type for more than a century since the establishing of the technique in the late 1890s. Between the end of the 1960s and the end of the 1990s Nansen cast technique has been gradually replaced by electronic sensor profilers (CTD). Both instrumentation types are considered as the most accurate among other oceanographic instruments and are often used as the unbiased reference. We conducted a comprehensive investigation of the consistency of the temperature data from Nansen casts and CTD profilers analyzing the quasi-collocated bottle and CTD data between the 1960s and the 1990s when both instrumentation types overlap. We found that Nansen casts tend to overestimate the sample depth with reversing mercury-in-glass thermometer temperatures being on average slightly lower compared to CTD data. Respectively, depth and temperature corrections are provided. Further, we estimated the ocean heat content changes between 1955 and 1990 using (along with all other instrumentation types) corrected and uncorrected Nansen cast data. These calculations show that for the upper 2 km layer the global average warming trend for this time period increases from 0.20 ± 0.05 W m−2 for the uncorrected data to 0.28 ± 0.06 W m−2 for the corrected data at the 90% confidence level. Finally, we suggest that the Nansen bottle cast profiles be put into a separate instrumentation group within the World Ocean Database.
Abstract
Nansen bottle casts served as the main oceanographic instrumentation type for more than a century since the establishing of the technique in the late 1890s. Between the end of the 1960s and the end of the 1990s Nansen cast technique has been gradually replaced by electronic sensor profilers (CTD). Both instrumentation types are considered as the most accurate among other oceanographic instruments and are often used as the unbiased reference. We conducted a comprehensive investigation of the consistency of the temperature data from Nansen casts and CTD profilers analyzing the quasi-collocated bottle and CTD data between the 1960s and the 1990s when both instrumentation types overlap. We found that Nansen casts tend to overestimate the sample depth with reversing mercury-in-glass thermometer temperatures being on average slightly lower compared to CTD data. Respectively, depth and temperature corrections are provided. Further, we estimated the ocean heat content changes between 1955 and 1990 using (along with all other instrumentation types) corrected and uncorrected Nansen cast data. These calculations show that for the upper 2 km layer the global average warming trend for this time period increases from 0.20 ± 0.05 W m−2 for the uncorrected data to 0.28 ± 0.06 W m−2 for the corrected data at the 90% confidence level. Finally, we suggest that the Nansen bottle cast profiles be put into a separate instrumentation group within the World Ocean Database.
Abstract
We estimate ocean heat content (OHC) change in the upper 2000 m in the Gulf of Mexico (GOM) from 1950 to 2020 to improve understanding of regional warming. Our estimates are based on 192 890 temperature profiles from the World Ocean Database. Warming occurs at all depths and in most regions except for a small region at northeastern GOM between 200 and 600 m. GOM OHC in the upper 2000 m increases at a rate of 0.38 ± 0.13 ZJ decade−1 between 1970 and 2020, which is equivalent to 1.21 ± 0.41 terawatts (TW). The GOM sea surface temperature (SST) increased ∼1.0° ± 0.25°C between 1970 and 2020, equivalent to a warming rate of 0.19° ± 0.05°C decade−1. Although SST in the GOM increases at a rate approximately twice that for the global ocean, the full-depth ocean heat storage rate in the GOM (0.86 ± 0.26 W m−2) applied to the entire GOM surface is comparable to that for the global ocean (0.82–1.11 W m−2). The upper-1000-m layer accounts for approximately 80%–90% of the total warming and variations in the upper 2000 m in the GOM. The Loop Current advective net heat flux is estimated to be 40.7 ± 6.3 TW through the GOM. A heat budget analysis shows the difference between the advective heat flux and the ocean heat storage rate (1.76 ± 1.36 TW, 1992–2017) can be roughly balanced with the annual net surface heat flux from ECCO (−37.9 TW).
Abstract
We estimate ocean heat content (OHC) change in the upper 2000 m in the Gulf of Mexico (GOM) from 1950 to 2020 to improve understanding of regional warming. Our estimates are based on 192 890 temperature profiles from the World Ocean Database. Warming occurs at all depths and in most regions except for a small region at northeastern GOM between 200 and 600 m. GOM OHC in the upper 2000 m increases at a rate of 0.38 ± 0.13 ZJ decade−1 between 1970 and 2020, which is equivalent to 1.21 ± 0.41 terawatts (TW). The GOM sea surface temperature (SST) increased ∼1.0° ± 0.25°C between 1970 and 2020, equivalent to a warming rate of 0.19° ± 0.05°C decade−1. Although SST in the GOM increases at a rate approximately twice that for the global ocean, the full-depth ocean heat storage rate in the GOM (0.86 ± 0.26 W m−2) applied to the entire GOM surface is comparable to that for the global ocean (0.82–1.11 W m−2). The upper-1000-m layer accounts for approximately 80%–90% of the total warming and variations in the upper 2000 m in the GOM. The Loop Current advective net heat flux is estimated to be 40.7 ± 6.3 TW through the GOM. A heat budget analysis shows the difference between the advective heat flux and the ocean heat storage rate (1.76 ± 1.36 TW, 1992–2017) can be roughly balanced with the annual net surface heat flux from ECCO (−37.9 TW).
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.
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.
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.
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.
Abstract
Historical bathythermograph datasets are known to be biased, and there have been several efforts to model this bias. Three different correction models of temperature bias in the historical bathythermograph dataset are compared here: the steady model of Hanawa et al. and the time-dependent models of Levitus et al. and Wijffels et al. The impact of these different models is examined in the context of global analysis experiments using the Simple Ocean Data Assimilation system. The results show that the two time-dependent bias models significantly reduce warm bias in global heat content, notably in the 10 years starting in the early 1970s and again in the early 1990s. Overall, the Levitus et al. model has its greatest impact near the surface and the Wijffels et al. model has its greatest impact at subtropical thermocline depths. Examination of the vertical structure of temperature error shows that at thermocline depths the Wijffels et al. model overcompensates, leading to a slight cool bias, while at shallow levels the same model causes a slight warm bias in the central and eastern subtropics and at thermocline depths on the equator in the Pacific Ocean as a result of reduced vertical entrainment. The results also show that the bias-correction models may alter the representation of interannual variability. During the 1997/98 El Niño and the subsequent La Niña the Levitus et al. model, which has its main impact at shallow depths, reduces the 50-m temperature anomalies in the eastern equatorial Pacific by 10%–20% and strengthens the zonal currents by up to 50%. The Wijffels et al. correction, which has its main impact at deeper levels, has much less effect on the oceanic expression of ENSO.
Abstract
Historical bathythermograph datasets are known to be biased, and there have been several efforts to model this bias. Three different correction models of temperature bias in the historical bathythermograph dataset are compared here: the steady model of Hanawa et al. and the time-dependent models of Levitus et al. and Wijffels et al. The impact of these different models is examined in the context of global analysis experiments using the Simple Ocean Data Assimilation system. The results show that the two time-dependent bias models significantly reduce warm bias in global heat content, notably in the 10 years starting in the early 1970s and again in the early 1990s. Overall, the Levitus et al. model has its greatest impact near the surface and the Wijffels et al. model has its greatest impact at subtropical thermocline depths. Examination of the vertical structure of temperature error shows that at thermocline depths the Wijffels et al. model overcompensates, leading to a slight cool bias, while at shallow levels the same model causes a slight warm bias in the central and eastern subtropics and at thermocline depths on the equator in the Pacific Ocean as a result of reduced vertical entrainment. The results also show that the bias-correction models may alter the representation of interannual variability. During the 1997/98 El Niño and the subsequent La Niña the Levitus et al. model, which has its main impact at shallow depths, reduces the 50-m temperature anomalies in the eastern equatorial Pacific by 10%–20% and strengthens the zonal currents by up to 50%. The Wijffels et al. correction, which has its main impact at deeper levels, has much less effect on the oceanic expression of ENSO.
An Algorithm for Classifying Unknown Expendable Bathythermograph (XBT) Instruments Based on Existing MetadataAbstract
Time-varying biases in expendable bathythermograph (XBT) instruments have emerged as a key uncertainty in estimates of historical ocean heat content variability and change. One of the challenges in the development of XBT bias corrections is the lack of metadata in ocean profile databases. Approximately 50% of XBT profiles in the World Ocean database (WOD) have no information about manufacturer or probe type. Building on previous research efforts, this paper presents a deterministic algorithm for assigning missing XBT manufacturer and probe type for individual temperature profiles based on 1) the reporting country, 2) the maximum reported depth, and 3) the record date. The criteria used are based on bulk analysis of known XBT profiles in the WOD for the period 1966–2015. A basic skill assessment demonstrates a 77% success rate at correctly assigning manufacturer and probe type for profiles where this information is available. The skill rate is lowest during the early 1990s, which is also a period when metadata information is particularly poor. The results suggest that substantive improvements could be made through further data analysis and that future algorithms may benefit from including a larger number of predictor variables.
Abstract
Time-varying biases in expendable bathythermograph (XBT) instruments have emerged as a key uncertainty in estimates of historical ocean heat content variability and change. One of the challenges in the development of XBT bias corrections is the lack of metadata in ocean profile databases. Approximately 50% of XBT profiles in the World Ocean database (WOD) have no information about manufacturer or probe type. Building on previous research efforts, this paper presents a deterministic algorithm for assigning missing XBT manufacturer and probe type for individual temperature profiles based on 1) the reporting country, 2) the maximum reported depth, and 3) the record date. The criteria used are based on bulk analysis of known XBT profiles in the WOD for the period 1966–2015. A basic skill assessment demonstrates a 77% success rate at correctly assigning manufacturer and probe type for profiles where this information is available. The skill rate is lowest during the early 1990s, which is also a period when metadata information is particularly poor. The results suggest that substantive improvements could be made through further data analysis and that future algorithms may benefit from including a larger number of predictor variables.
Abstract
Long time series of XBT data in the Bay of Bengal and the Arabian Sea are valuable datasets for exploring and understanding climate variability. However, such studies of interannual and longer-scale variability of temperature require an understanding, and, if possible, a correction of errors introduced by biases in the XBT and expendable conductivity–temperature–depth (XCTD) data. Two cruises in each basin were undertaken in 2008/09 on which series of tests of XBTs and XCTDs dropped simultaneously with CTD casts were performed. The XBT and XCTD depths were corrected by comparison with CTD data using a modification of an existing algorithm. Significant probe-to-probe fall-rate equation (FRE) velocity and deceleration coefficient variability was found within a cruise, as well as cruise-to-cruise variability. A small (∼0.01°C) temperature bias was also identified for XBTs on each cruise. No new FRE can be proposed for either the Bay of Bengal or the Arabian Sea for XBTs. For the more consistent XCTD, basin-specific FREs are possible for the Bay of Bengal, but not for the Arabian Sea. The XCTD FRE velocity coefficients are significantly higher than the XCTD manufacturers’ FRE coefficient or those from previous tests, possibly resulting from the influence of temperature on XCTD FRE. Mean temperature anomalies versus a long-term mean climatology for XBT data using the present default FRE have a bias (which is positive for three cruises and negative for one cruise) compared to the mean temperature anomalies for CTD data. Some improvement is found when applying newly calculated cruise-specific FREs. This temperature error must be accounted for in any study of temperature change in the regions.
Abstract
Long time series of XBT data in the Bay of Bengal and the Arabian Sea are valuable datasets for exploring and understanding climate variability. However, such studies of interannual and longer-scale variability of temperature require an understanding, and, if possible, a correction of errors introduced by biases in the XBT and expendable conductivity–temperature–depth (XCTD) data. Two cruises in each basin were undertaken in 2008/09 on which series of tests of XBTs and XCTDs dropped simultaneously with CTD casts were performed. The XBT and XCTD depths were corrected by comparison with CTD data using a modification of an existing algorithm. Significant probe-to-probe fall-rate equation (FRE) velocity and deceleration coefficient variability was found within a cruise, as well as cruise-to-cruise variability. A small (∼0.01°C) temperature bias was also identified for XBTs on each cruise. No new FRE can be proposed for either the Bay of Bengal or the Arabian Sea for XBTs. For the more consistent XCTD, basin-specific FREs are possible for the Bay of Bengal, but not for the Arabian Sea. The XCTD FRE velocity coefficients are significantly higher than the XCTD manufacturers’ FRE coefficient or those from previous tests, possibly resulting from the influence of temperature on XCTD FRE. Mean temperature anomalies versus a long-term mean climatology for XBT data using the present default FRE have a bias (which is positive for three cruises and negative for one cruise) compared to the mean temperature anomalies for CTD data. Some improvement is found when applying newly calculated cruise-specific FREs. This temperature error must be accounted for in any study of temperature change in the regions.
Abstract
Ocean heat content (HC) is one of the key indicators of climate variability and also provides ocean memory critical for seasonal and decadal predictions. The availability of multiple operational ocean analyses (ORAs) now routinely produced around the world is an opportunity for estimation of uncertainties in HC analysis and development of ensemble-based operational HC climate indices. In this context, the spread across the ORAs is used to quantify uncertainties in HC analysis and the ensemble mean of ORAs to identify, and to monitor, climate signals. Toward this goal, this study analyzed 10 ORAs, two objective analyses based on in situ data only, and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability, and long-term trend of HC in the upper 300 m (HC300) from 1980 to 2009 are compared.
The spread across HC300 analyses generally decreased with time and reached a minimum in the early 2000s when the Argo data became available. There was a good correspondence between the increase of data counts and reduction of the spread. The agreement of HC300 anomalies among different ORAs, measured by the signal-to-noise ratio (S/N), is generally high in the tropical Pacific, tropical Indian Ocean, North Pacific, and North Atlantic but low in the tropical Atlantic and extratropical southern oceans where observations are very sparse. A set of climate indices was derived as HC300 anomalies averaged over the areas where the covariability between SST and HC300 represents the major climate modes such as ENSO, Indian Ocean dipole, Atlantic Niño, Pacific decadal oscillation, and Atlantic multidecadal oscillation.
Abstract
Ocean heat content (HC) is one of the key indicators of climate variability and also provides ocean memory critical for seasonal and decadal predictions. The availability of multiple operational ocean analyses (ORAs) now routinely produced around the world is an opportunity for estimation of uncertainties in HC analysis and development of ensemble-based operational HC climate indices. In this context, the spread across the ORAs is used to quantify uncertainties in HC analysis and the ensemble mean of ORAs to identify, and to monitor, climate signals. Toward this goal, this study analyzed 10 ORAs, two objective analyses based on in situ data only, and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability, and long-term trend of HC in the upper 300 m (HC300) from 1980 to 2009 are compared.
The spread across HC300 analyses generally decreased with time and reached a minimum in the early 2000s when the Argo data became available. There was a good correspondence between the increase of data counts and reduction of the spread. The agreement of HC300 anomalies among different ORAs, measured by the signal-to-noise ratio (S/N), is generally high in the tropical Pacific, tropical Indian Ocean, North Pacific, and North Atlantic but low in the tropical Atlantic and extratropical southern oceans where observations are very sparse. A set of climate indices was derived as HC300 anomalies averaged over the areas where the covariability between SST and HC300 represents the major climate modes such as ENSO, Indian Ocean dipole, Atlantic Niño, Pacific decadal oscillation, and Atlantic multidecadal oscillation.
Abstract
The difficulty in effectively evaluating sea surface temperature (SST) analyses is finding independent observations, since most available observations have been used in the SST analyses. In this study, the ocean profile measurements [from reverse thermometer, CTD, mechanical bathythermograph (MBT), and XBT] above 5-m depth over 1950–2016 from the World Ocean Database (WOD) are used (data labeled pSSTW). The biases of MBT and XBT are corrected based on currently available algorithms. The bias-corrected pSSTW over 1950–2016 and satellite-based SST from the European Space Agency (ESA) Climate Change Initiative (CCI) over 1992–2010 are used to evaluate commonly available SST analyses. These SST analyses are the Extended Reconstructed SST (ERSST), versions 5, 4, and 3b, the Met Office Hadley Centre Sea Ice and SST dataset (HadISST), and the Japan Meteorological Administration (JMA) Centennial In Situ Observation-Based Estimates of SST version 2.9.2 (COBE-SST2). Our comparisons show that the SST from COBE-SST2 is the closest to pSSTW and CCI in most of the Pacific, Atlantic, and Southern Oceans, which may result from its unique bias correction to ship observations. The SST from ERSST version 5 is more consistent with pSSTW than its previous versions over 1950–2016, and is more consistent with CCI than its previous versions over 1992–2010. The better performance of ERSST version 5 over its previous versions is attributed to its improved bias correction applied to ship observations with a baseline of buoy observations, and is seen in most of the Pacific and Atlantic.
Abstract
The difficulty in effectively evaluating sea surface temperature (SST) analyses is finding independent observations, since most available observations have been used in the SST analyses. In this study, the ocean profile measurements [from reverse thermometer, CTD, mechanical bathythermograph (MBT), and XBT] above 5-m depth over 1950–2016 from the World Ocean Database (WOD) are used (data labeled pSSTW). The biases of MBT and XBT are corrected based on currently available algorithms. The bias-corrected pSSTW over 1950–2016 and satellite-based SST from the European Space Agency (ESA) Climate Change Initiative (CCI) over 1992–2010 are used to evaluate commonly available SST analyses. These SST analyses are the Extended Reconstructed SST (ERSST), versions 5, 4, and 3b, the Met Office Hadley Centre Sea Ice and SST dataset (HadISST), and the Japan Meteorological Administration (JMA) Centennial In Situ Observation-Based Estimates of SST version 2.9.2 (COBE-SST2). Our comparisons show that the SST from COBE-SST2 is the closest to pSSTW and CCI in most of the Pacific, Atlantic, and Southern Oceans, which may result from its unique bias correction to ship observations. The SST from ERSST version 5 is more consistent with pSSTW than its previous versions over 1950–2016, and is more consistent with CCI than its previous versions over 1992–2010. The better performance of ERSST version 5 over its previous versions is attributed to its improved bias correction applied to ship observations with a baseline of buoy observations, and is seen in most of the Pacific and Atlantic.
