Browse
Abstract
Vertical profiles of shortwave and longwave irradiances computed with satellite-derived cloud properties and temperature and humidity profiles from reanalysis are used to estimate entropy production. Entropy production by shortwave radiation is computed by the absorbed irradiance within layers in the atmosphere and by the surface divided by their temperatures. Similarly, entropy production by longwave radiation is computed by emitted irradiance to space from layers in the atmosphere and surface divided by their temperatures. Global annual mean entropy production by shortwave absorption and longwave emission to space are, respectively, 0.852 and 0.928 W m−2 K−1. With a steady-state assumption, entropy production by irreversible processes within the Earth system is estimated to be 0.076 W m−2 K−1 and by nonradiative irreversible processes to be 0.049 W m−2 K−1. Both global annual mean entropy productions by shortwave absorption and longwave emission to space increase with increasing shortwave absorption (i.e., with decreasing the planetary albedo). The increase of entropy production by shortwave absorption is, however, larger than the increase of entropy production by longwave emission to space. The result implies that global annual mean entropy production by irreversible processes decreases with increasing shortwave absorption. Input and output temperatures derived by dividing the absorbed shortwave irradiance and emitted longwave irradiance to space by respective entropy production are, respectively, 282 and 259 K, which give the Carnot efficiency of the Earth system of 8.5%.
Abstract
Vertical profiles of shortwave and longwave irradiances computed with satellite-derived cloud properties and temperature and humidity profiles from reanalysis are used to estimate entropy production. Entropy production by shortwave radiation is computed by the absorbed irradiance within layers in the atmosphere and by the surface divided by their temperatures. Similarly, entropy production by longwave radiation is computed by emitted irradiance to space from layers in the atmosphere and surface divided by their temperatures. Global annual mean entropy production by shortwave absorption and longwave emission to space are, respectively, 0.852 and 0.928 W m−2 K−1. With a steady-state assumption, entropy production by irreversible processes within the Earth system is estimated to be 0.076 W m−2 K−1 and by nonradiative irreversible processes to be 0.049 W m−2 K−1. Both global annual mean entropy productions by shortwave absorption and longwave emission to space increase with increasing shortwave absorption (i.e., with decreasing the planetary albedo). The increase of entropy production by shortwave absorption is, however, larger than the increase of entropy production by longwave emission to space. The result implies that global annual mean entropy production by irreversible processes decreases with increasing shortwave absorption. Input and output temperatures derived by dividing the absorbed shortwave irradiance and emitted longwave irradiance to space by respective entropy production are, respectively, 282 and 259 K, which give the Carnot efficiency of the Earth system of 8.5%.
Abstract
The NASA Energy and Water Cycle Study (NEWS) climatology is a self-consistent coupled annual and seasonal cycle solution for radiative, turbulent, and water fluxes over Earth’s surface using Earth observation data covering 2000–09. Here we seek to improve the NEWS solution, particularly over the ocean basins, by considering spatial covariances in the observation errors (some evidence for which is found by comparing five turbulent flux products over the oceans) and by introducing additional horizontal transports from ocean reanalyses as weak constraints. By explicitly representing large error covariances between surface heat flux components over the major ocean basins we retain the flux contrasts present in the original data and infer additional heat losses over the North Atlantic Ocean, more consistent with a strong Atlantic overturning. This change does not alter the global flux balance but if only the errors in evaporation and precipitation are correlated then those fluxes experience larger adjustments (e.g., the surface latent heat flux increases to 85 ± 2 W m−2). Replacing SeaFlux v1 with J-OFURO v3 (Japanese Ocean Flux Data Sets with Use of Remote Sensing Observations) ocean fluxes also leads to a considerable increase in the global latent heat loss as well as a larger North Atlantic heat loss. Furthermore, including a weak constraint on the horizontal transports of heat and freshwater from high-resolution ocean reanalyses improves the net fluxes over the North Atlantic, Caribbean Sea, and Arctic Ocean, without any impact on the global flux balances. These results suggest that better characterized flux uncertainties can greatly improve the quality of the optimized flux solution.
Abstract
The NASA Energy and Water Cycle Study (NEWS) climatology is a self-consistent coupled annual and seasonal cycle solution for radiative, turbulent, and water fluxes over Earth’s surface using Earth observation data covering 2000–09. Here we seek to improve the NEWS solution, particularly over the ocean basins, by considering spatial covariances in the observation errors (some evidence for which is found by comparing five turbulent flux products over the oceans) and by introducing additional horizontal transports from ocean reanalyses as weak constraints. By explicitly representing large error covariances between surface heat flux components over the major ocean basins we retain the flux contrasts present in the original data and infer additional heat losses over the North Atlantic Ocean, more consistent with a strong Atlantic overturning. This change does not alter the global flux balance but if only the errors in evaporation and precipitation are correlated then those fluxes experience larger adjustments (e.g., the surface latent heat flux increases to 85 ± 2 W m−2). Replacing SeaFlux v1 with J-OFURO v3 (Japanese Ocean Flux Data Sets with Use of Remote Sensing Observations) ocean fluxes also leads to a considerable increase in the global latent heat loss as well as a larger North Atlantic heat loss. Furthermore, including a weak constraint on the horizontal transports of heat and freshwater from high-resolution ocean reanalyses improves the net fluxes over the North Atlantic, Caribbean Sea, and Arctic Ocean, without any impact on the global flux balances. These results suggest that better characterized flux uncertainties can greatly improve the quality of the optimized flux solution.
Abstract
The ocean surface albedo is responsible for the distribution of solar (shortwave) radiant energy between the atmosphere and ocean and therefore is a key parameter in Earth’s surface energy budget. In situ ocean observations typically do not measure upward reflected solar radiation, which is necessary to compute net solar radiation into the ocean. Instead, the upward component is computed from the measured downward component using an albedo estimate. At two NOAA Ocean Climate Station buoy sites in the North Pacific, the International Satellite Cloud Climatology Project (ISCCP) monthly climatological albedo has been used, while for the NOAA Global Tropical Buoy Array a constant albedo is used. This constant albedo is also used in the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk flux algorithm. This study considers the impacts of using the more recently available NASA Cloud and the Earth’s Radiant Energy System (CERES) albedo product for these ocean surface heat flux products. Differences between albedo estimates in global satellite products like these imply uncertainty in the net surface solar radiation heat flux estimates that locally exceed the target uncertainty of 1.0 W m−2 for the global mean, set by the Global Climate Observing System (GCOS) of the World Meteorological Organization (WMO). Albedo has large spatiotemporal variability on hourly, monthly, and interannual time scales. Biases in high-resolution SWnet (the difference between surface downwelling and upwelling shortwave radiation) can arise if the albedo diurnal cycle is unresolved. As a result, for periods when satellite albedo data are not available it is recommended that an hourly climatology be used when computing high-resolution net surface shortwave radiation.
Abstract
The ocean surface albedo is responsible for the distribution of solar (shortwave) radiant energy between the atmosphere and ocean and therefore is a key parameter in Earth’s surface energy budget. In situ ocean observations typically do not measure upward reflected solar radiation, which is necessary to compute net solar radiation into the ocean. Instead, the upward component is computed from the measured downward component using an albedo estimate. At two NOAA Ocean Climate Station buoy sites in the North Pacific, the International Satellite Cloud Climatology Project (ISCCP) monthly climatological albedo has been used, while for the NOAA Global Tropical Buoy Array a constant albedo is used. This constant albedo is also used in the Coupled Ocean–Atmosphere Response Experiment (COARE) bulk flux algorithm. This study considers the impacts of using the more recently available NASA Cloud and the Earth’s Radiant Energy System (CERES) albedo product for these ocean surface heat flux products. Differences between albedo estimates in global satellite products like these imply uncertainty in the net surface solar radiation heat flux estimates that locally exceed the target uncertainty of 1.0 W m−2 for the global mean, set by the Global Climate Observing System (GCOS) of the World Meteorological Organization (WMO). Albedo has large spatiotemporal variability on hourly, monthly, and interannual time scales. Biases in high-resolution SWnet (the difference between surface downwelling and upwelling shortwave radiation) can arise if the albedo diurnal cycle is unresolved. As a result, for periods when satellite albedo data are not available it is recommended that an hourly climatology be used when computing high-resolution net surface shortwave radiation.
Abstract
The net surface energy flux is computed as a residual of the energy budget using top-of-atmosphere radiation combined with the divergence of the column-integrated atmospheric energy transports, and then used with the vertically integrated ocean heat content tendencies to compute the ocean meridional heat transports (MHTs). The mean annual cycles and 12-month running mean MHTs as a function of latitude are presented for 2000–16. Effects of the Indonesian Throughflow (ITF), associated with a net volume flow around Australia accompanied by a heat transport, are fully included. Because the ITF-related flow necessitates a return current northward in the Tasman Sea that relaxes during El Niño, the reduced ITF during El Niño may contribute to warming in the south Tasman Sea by allowing the East Australian Current to push farther south even as it gains volume from the tropical waters not flowing through the ITF. Although evident in 2015/16, when a major marine heat wave occurred, these effects can be overwhelmed by changes in the atmospheric circulation. Large interannual MHT variability in the Pacific is 4 times that of the Atlantic. Strong relationships reveal influences from the southern subtropics on ENSO for this period. At the equator, northward ocean MHT arises mainly in the Atlantic (0.75 PW), offset by the Pacific (−0.33 PW) and Indian Oceans (−0.20 PW) while the atmosphere transports energy southward (−0.35 PW). The net equatorial MHT southward (−0.18 PW) is enhanced by −0.1 PW that contributes to the greater warming of the southern (vs northern) oceans.
Abstract
The net surface energy flux is computed as a residual of the energy budget using top-of-atmosphere radiation combined with the divergence of the column-integrated atmospheric energy transports, and then used with the vertically integrated ocean heat content tendencies to compute the ocean meridional heat transports (MHTs). The mean annual cycles and 12-month running mean MHTs as a function of latitude are presented for 2000–16. Effects of the Indonesian Throughflow (ITF), associated with a net volume flow around Australia accompanied by a heat transport, are fully included. Because the ITF-related flow necessitates a return current northward in the Tasman Sea that relaxes during El Niño, the reduced ITF during El Niño may contribute to warming in the south Tasman Sea by allowing the East Australian Current to push farther south even as it gains volume from the tropical waters not flowing through the ITF. Although evident in 2015/16, when a major marine heat wave occurred, these effects can be overwhelmed by changes in the atmospheric circulation. Large interannual MHT variability in the Pacific is 4 times that of the Atlantic. Strong relationships reveal influences from the southern subtropics on ENSO for this period. At the equator, northward ocean MHT arises mainly in the Atlantic (0.75 PW), offset by the Pacific (−0.33 PW) and Indian Oceans (−0.20 PW) while the atmosphere transports energy southward (−0.35 PW). The net equatorial MHT southward (−0.18 PW) is enhanced by −0.1 PW that contributes to the greater warming of the southern (vs northern) oceans.
Abstract
This study combines state-of-the-art reanalyses such as the fifth-generation European Re-Analysis (ERA5) and the Ocean Reanalysis System 5 (ORAS5) with novel observational products to present an updated estimate of the coupled atmosphere–ocean–sea ice Arctic energy budget, including flux and storage terms covering 2001–17. Observational products provide independent estimates of crucial budget terms, including oceanic heat transport from unique mooring-derived data, radiative fluxes from satellites, and sea ice volume from merged satellite data. Results show that the time averages of independent estimates of radiative, atmospheric, and oceanic energy fluxes into the Arctic Ocean domain are remarkably consistent in the sense that their sum closely matches the observed rate of regional long-term oceanic heat accumulation of ~1 W m−2. Atmospheric and oceanic heat transports are found to be stronger compared to earlier assessments (~100 and ~16 W m−2, respectively). Data inconsistencies are larger when considering the mean annual cycle of the coupled energy budget, with RMS values of the monthly budget residual between 7 and 15 W m−2, depending on the employed datasets. This nevertheless represents an average reduction of ~72% of the residual compared to earlier work and demonstrates the progress made in data quality and diagnostic techniques. Finally, the budget residual is eliminated using a variational approach to provide a best estimate of the mean annual cycle. The largest remaining sources of uncertainty are ocean heat content and latent heat associated with sea ice melt and freeze, which both suffer from the lack of observational constraints. More ocean in situ observations and reliable sea ice thickness observations and their routinely assimilation into reanalyses are needed to further reduce uncertainty.
Abstract
This study combines state-of-the-art reanalyses such as the fifth-generation European Re-Analysis (ERA5) and the Ocean Reanalysis System 5 (ORAS5) with novel observational products to present an updated estimate of the coupled atmosphere–ocean–sea ice Arctic energy budget, including flux and storage terms covering 2001–17. Observational products provide independent estimates of crucial budget terms, including oceanic heat transport from unique mooring-derived data, radiative fluxes from satellites, and sea ice volume from merged satellite data. Results show that the time averages of independent estimates of radiative, atmospheric, and oceanic energy fluxes into the Arctic Ocean domain are remarkably consistent in the sense that their sum closely matches the observed rate of regional long-term oceanic heat accumulation of ~1 W m−2. Atmospheric and oceanic heat transports are found to be stronger compared to earlier assessments (~100 and ~16 W m−2, respectively). Data inconsistencies are larger when considering the mean annual cycle of the coupled energy budget, with RMS values of the monthly budget residual between 7 and 15 W m−2, depending on the employed datasets. This nevertheless represents an average reduction of ~72% of the residual compared to earlier work and demonstrates the progress made in data quality and diagnostic techniques. Finally, the budget residual is eliminated using a variational approach to provide a best estimate of the mean annual cycle. The largest remaining sources of uncertainty are ocean heat content and latent heat associated with sea ice melt and freeze, which both suffer from the lack of observational constraints. More ocean in situ observations and reliable sea ice thickness observations and their routinely assimilation into reanalyses are needed to further reduce uncertainty.
Abstract
This study revisits the classical problem of quantifying the radiative effects of unique cloud types in the era of spaceborne active observations. The radiative effects of nine cloud types, distinguished based on their vertical structure defined by CloudSat and CALIPSO observations, are assessed at both the top of the atmosphere and the surface. The contributions from single- and multilayered clouds are explicitly diagnosed. The global, annual mean net cloud radiative effect at the top of the atmosphere is found to be −17.1 ± 4.2 W m−2 owing to −44.2 ± 2 W m−2 of shortwave cooling and 27.1 ± 3.7 W m−2 of longwave heating. Leveraging explicit cloud base and vertical structure information, we further estimate the annual mean net cloud radiative effect at the surface to be −24.8 ± 8.7 W m−2 (−51.1 ± 7.8 W m−2 in the shortwave and 26.3 ± 3.8 W m−2 in the longwave). Multilayered clouds are found to exert the strongest influence on the top-of-atmosphere energy balance. However, a strong asymmetry in net cloud radiative cooling between the hemispheres (8.6 W m−2) is dominated by enhanced cooling from stratocumulus over the southern oceans. It is found that there is no corresponding asymmetry at the surface owing to enhanced longwave emission by southern ocean clouds in winter, which offsets a substantial fraction of their impact on solar absorption in summer. Thus the asymmetry in cloud radiative effects is entirely realized as an atmosphere heating imbalance between the hemispheres.
Abstract
This study revisits the classical problem of quantifying the radiative effects of unique cloud types in the era of spaceborne active observations. The radiative effects of nine cloud types, distinguished based on their vertical structure defined by CloudSat and CALIPSO observations, are assessed at both the top of the atmosphere and the surface. The contributions from single- and multilayered clouds are explicitly diagnosed. The global, annual mean net cloud radiative effect at the top of the atmosphere is found to be −17.1 ± 4.2 W m−2 owing to −44.2 ± 2 W m−2 of shortwave cooling and 27.1 ± 3.7 W m−2 of longwave heating. Leveraging explicit cloud base and vertical structure information, we further estimate the annual mean net cloud radiative effect at the surface to be −24.8 ± 8.7 W m−2 (−51.1 ± 7.8 W m−2 in the shortwave and 26.3 ± 3.8 W m−2 in the longwave). Multilayered clouds are found to exert the strongest influence on the top-of-atmosphere energy balance. However, a strong asymmetry in net cloud radiative cooling between the hemispheres (8.6 W m−2) is dominated by enhanced cooling from stratocumulus over the southern oceans. It is found that there is no corresponding asymmetry at the surface owing to enhanced longwave emission by southern ocean clouds in winter, which offsets a substantial fraction of their impact on solar absorption in summer. Thus the asymmetry in cloud radiative effects is entirely realized as an atmosphere heating imbalance between the hemispheres.
Abstract
A diagnostic tool for determining surface and atmospheric contributions to interannual variations in top-of-atmosphere (TOA) reflected shortwave (SW) and net downward SW surface radiative fluxes is introduced. The method requires only upward and downward radiative fluxes at the TOA and surface as input and therefore can readily be applied to both satellite-derived and model-generated radiative fluxes. Observations from the Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Edition 4.0 product show that 81% of the monthly variability in global mean reflected SW TOA flux anomalies is associated with atmospheric variations (mainly clouds), 6% is from surface variations, and 13% is from atmosphere–surface covariability. Over the Arctic Ocean, most of the variability in both reflected SW TOA flux and net downward SW surface flux anomalies is explained by variations in sea ice and cloud fraction alone (r 2 = 0.94). Compared to CERES, variability in two reanalyses—the ECMWF interim reanalysis (ERA-Interim) and NASA’s Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2)—show large differences in the regional distribution of variance for both the atmospheric and surface contributions to anomalies in net downward SW surface flux. For MERRA-2 the atmospheric contribution is 17% too large compared to CERES while ERA-Interim underestimates the variance by 15%. The difference is mainly due to how cloud variations are represented in the reanalyses. The overall surface contribution in both ERA-Interim and MERRA-2 is smaller than CERES EBAF by 15% for ERA-Interim and 58% for MERRA-2, highlighting limitations of the reanalyses in representing surface albedo variations and their influence on SW radiative fluxes.
Abstract
A diagnostic tool for determining surface and atmospheric contributions to interannual variations in top-of-atmosphere (TOA) reflected shortwave (SW) and net downward SW surface radiative fluxes is introduced. The method requires only upward and downward radiative fluxes at the TOA and surface as input and therefore can readily be applied to both satellite-derived and model-generated radiative fluxes. Observations from the Clouds and the Earth’s Radiant Energy System (CERES) Energy Balanced and Filled (EBAF) Edition 4.0 product show that 81% of the monthly variability in global mean reflected SW TOA flux anomalies is associated with atmospheric variations (mainly clouds), 6% is from surface variations, and 13% is from atmosphere–surface covariability. Over the Arctic Ocean, most of the variability in both reflected SW TOA flux and net downward SW surface flux anomalies is explained by variations in sea ice and cloud fraction alone (r 2 = 0.94). Compared to CERES, variability in two reanalyses—the ECMWF interim reanalysis (ERA-Interim) and NASA’s Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2)—show large differences in the regional distribution of variance for both the atmospheric and surface contributions to anomalies in net downward SW surface flux. For MERRA-2 the atmospheric contribution is 17% too large compared to CERES while ERA-Interim underestimates the variance by 15%. The difference is mainly due to how cloud variations are represented in the reanalyses. The overall surface contribution in both ERA-Interim and MERRA-2 is smaller than CERES EBAF by 15% for ERA-Interim and 58% for MERRA-2, highlighting limitations of the reanalyses in representing surface albedo variations and their influence on SW radiative fluxes.
Abstract
The spectrum of reflected solar radiation emerging at the top of the atmosphere is rich with Earth system information. To identify spectral signatures in the reflected solar radiation and directly relate them to the underlying physical properties controlling their structure, over 90 000 solar reflectance spectra are computed over West Africa in 2010 using a fast radiation code employing the spectral characteristics of the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY). Cluster analysis applied to the computed spectra reveals spectral signatures related to distinct surface properties, and cloud regimes distinguished by their spectral shortwave cloud radiative effect (SWCRE). The cloud regimes exhibit a diverse variety of mean broadband SWCREs, and offer an alternative approach to define cloud type for SWCRE applications that does not require any prior assumptions. The direct link between spectral signatures and distinct physical properties extracted from clustering remains robust between spatial scales of 1, 20, and 240 km, and presents an excellent opportunity to understand the underlying properties controlling real spectral reflectance observations. Observed SCIAMACHY spectra are assigned to the calculated spectral clusters, showing that cloud regimes are most frequent during the active West African monsoon season of June–October in 2010, and all cloud regimes have a higher frequency of occurrence during the active monsoon season of 2003 compared with the inactive monsoon season of 2004. Overall, the distinct underlying physical properties controlling spectral signatures show great promise for monitoring evolution of the Earth system directly from solar spectral reflectance observations.
Abstract
The spectrum of reflected solar radiation emerging at the top of the atmosphere is rich with Earth system information. To identify spectral signatures in the reflected solar radiation and directly relate them to the underlying physical properties controlling their structure, over 90 000 solar reflectance spectra are computed over West Africa in 2010 using a fast radiation code employing the spectral characteristics of the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY). Cluster analysis applied to the computed spectra reveals spectral signatures related to distinct surface properties, and cloud regimes distinguished by their spectral shortwave cloud radiative effect (SWCRE). The cloud regimes exhibit a diverse variety of mean broadband SWCREs, and offer an alternative approach to define cloud type for SWCRE applications that does not require any prior assumptions. The direct link between spectral signatures and distinct physical properties extracted from clustering remains robust between spatial scales of 1, 20, and 240 km, and presents an excellent opportunity to understand the underlying properties controlling real spectral reflectance observations. Observed SCIAMACHY spectra are assigned to the calculated spectral clusters, showing that cloud regimes are most frequent during the active West African monsoon season of June–October in 2010, and all cloud regimes have a higher frequency of occurrence during the active monsoon season of 2003 compared with the inactive monsoon season of 2004. Overall, the distinct underlying physical properties controlling spectral signatures show great promise for monitoring evolution of the Earth system directly from solar spectral reflectance observations.
Abstract
Ocean meridional heat transports (MHTs) are deduced as a residual using energy budgets to produce latitude versus time series for the globe, Indo-Pacific, and Atlantic. The top-of-atmosphere (TOA) radiation is combined with the vertically integrated atmospheric energy divergence from atmospheric reanalyses to produce the net surface energy fluxes everywhere. The latter is then combined with estimates of the vertically integrated ocean heat content (OHC) tendency to produce estimates of the ocean heat divergence. Because seasonal sea ice and land runoff effects are not fully considered, the mean annual cycle is incomplete, but those effects are small for interannual variability. However, there is a mismatch between 12-month inferred surface flux and the corresponding OHC changes globally, requiring adjustments to account for the Earth’s global energy imbalance. Estimates are greatly improved by building in the constraint that MHT must go to zero at the northern and southern extents of the ocean basin at all times, enabling biases between the TOA and OHC data to be reconciled. Zonal mean global, Indo-Pacific, and Atlantic basin ocean MHTs are computed and presented as 12-month running means and for the mean annual cycle for 2000–16. For the Indo-Pacific, the tropical and subtropical MHTs feature a strong relationship with El Niño–Southern Oscillation (ENSO), and in the Atlantic, MHT interannual variability is significantly affected by and likely influences the North Atlantic Oscillation (NAO). However, Atlantic and Pacific changes are linked, suggesting that the northern annular mode (as opposed to NAO) is predominant. There is also evidence of decadal variability or trends.
Abstract
Ocean meridional heat transports (MHTs) are deduced as a residual using energy budgets to produce latitude versus time series for the globe, Indo-Pacific, and Atlantic. The top-of-atmosphere (TOA) radiation is combined with the vertically integrated atmospheric energy divergence from atmospheric reanalyses to produce the net surface energy fluxes everywhere. The latter is then combined with estimates of the vertically integrated ocean heat content (OHC) tendency to produce estimates of the ocean heat divergence. Because seasonal sea ice and land runoff effects are not fully considered, the mean annual cycle is incomplete, but those effects are small for interannual variability. However, there is a mismatch between 12-month inferred surface flux and the corresponding OHC changes globally, requiring adjustments to account for the Earth’s global energy imbalance. Estimates are greatly improved by building in the constraint that MHT must go to zero at the northern and southern extents of the ocean basin at all times, enabling biases between the TOA and OHC data to be reconciled. Zonal mean global, Indo-Pacific, and Atlantic basin ocean MHTs are computed and presented as 12-month running means and for the mean annual cycle for 2000–16. For the Indo-Pacific, the tropical and subtropical MHTs feature a strong relationship with El Niño–Southern Oscillation (ENSO), and in the Atlantic, MHT interannual variability is significantly affected by and likely influences the North Atlantic Oscillation (NAO). However, Atlantic and Pacific changes are linked, suggesting that the northern annular mode (as opposed to NAO) is predominant. There is also evidence of decadal variability or trends.