Abstract

Global estimates of oceanic evaporation (Evp) from 1958 to 2005 have been recently developed by the Objectively Analyzed Air–Sea Fluxes (OAFlux) project at the Woods Hole Oceanographic Institution (WHOI). The nearly 50-yr time series shows that the decadal change of the global oceanic evaporation (Evp) is marked by a distinct transition from a downward trend to an upward trend around 1977–78. Since the transition, the global oceanic Evp has been up about 11 cm yr⁻¹ (∼10%), from a low at 103 cm yr⁻¹ in 1977 to a peak at 114 cm yr⁻¹ in 2003. The increase in Evp was most dramatic during the 1990s. The uncertainty of the estimates is about ±2.74 cm yr⁻¹. By utilizing the newly developed datasets of Evp and related air–sea variables, the study investigated the cause of the decadal change in oceanic Evp. The decadal differences between the 1990s and the 1970s indicates that the increase of Evp in the 1990s occurred over a global scale and had spatially coherent structures. Larger Evp is most pronounced in two key regions—one is the paths of the global western boundary currents and their extensions, and the other is the tropical Indo-Pacific warm water pools. It is also found that Evp was enlarged primarily during the hemispheric wintertime (defined as the mean of December–February for the northern oceans and June–August for the southern oceans). Despite the dominant upward tendency over the global basins, a slight reduction in Evp appeared in such regions as the subtropical centers of the Evp maxima as well as the eastern equatorial Pacific and Atlantic cold tongues.

An empirical orthogonal function (EOF) analysis was performed for the yearly winter-mean time series of Evp and the related air–sea variables [i.e., wind speed (U) and air–sea humidity differences (dq)]. The analysis suggested a dominant role of the wind forcing in the decadal change of both Evp and dq. It is hypothesized that wind impacts Evp in two ways. The first way is direct: the greater wind speed induces more evaporation by carrying water vapor away from the evaporating surface to allow the air–sea humidity gradients to be reestablished at a faster pace. The second way is indirect: the enhanced surface wind strengthens the wind-driven subtropical gyre, which in turn drives a greater heat transport by the western boundary currents, warms up SST along the paths of the currents and extensions, and causes more evaporation by enlarging the air–sea humidity gradients. The EOF analysis performed for the time series of the global annual-mean Evp fields showed that the first three EOF modes account for nearly 50% of the total variance. The mode 1 variability represents the upward trend in Evp after 1978 and is attributable to the increased U, and the mode 2 variability explains much of the downward trend in Evp before 1978 and is correlated to the global dq variability. The EOF mode 3 of Evp captures the interannual variability of Evp on time scales of the El Niño–Southern Oscillation, with the center of action over the eastern equatorial Pacific.

Abstract

The correlation between parameters characterizing observed westerly wind bursts (WWBs) in the equatorial Pacific and the large-scale SST is analyzed using singular value decomposition. The WWB parameters include the amplitude, location, scale, and probability of occurrence for a given SST distribution rather than the wind stress itself. This approach therefore allows for a nonlinear relationship between the SST and the wind signal of the WWBs. It is found that about half of the variance of the WWB parameters is explained by only two large-scale SST modes. The first mode represents a developed El Niño event, while the second mode represents the seasonal cycle. More specifically, the central longitude of WWBs, their longitudinal extent, and their probability seem to be determined to a significant degree by the ENSO-driven signal. The amplitude of the WWBs is found to be strongly influenced by the phase of the seasonal cycle. It is concluded that the WWBs, while partially stochastic, seem an inherent part of the large-scale deterministic ENSO dynamics. Implications for ENSO predictability and prediction are discussed.

Abstract

An assessment is made of the mean and variability of the net air–sea heat flux, Q _net, from four products (ECCO, OAFlux–CERES, ERA-Interim, and NCEP1) over the global ice-free ocean from January 2001 to December 2010. For the 10-yr “hiatus” period, all products agree on an overall net heat gain over the global ice-free ocean, but the magnitude varies from 1.7 to 9.5 W m⁻². The differences among products are particularly large in the Southern Ocean, where they cannot even agree on whether the region gains or loses heat on the annual mean basis. Decadal trends of Q _net differ significantly between products. ECCO and OAFlux–CERES show almost no trend, whereas ERA-Interim suggests a downward trend and NCEP1 shows an upward trend. Therefore, numerical simulations utilizing different surface flux forcing products will likely produce diverged trends of the ocean heat content during this period. The downward trend in ERA-Interim started from 2006, driven by a peculiar pattern change in the tropical regions. ECCO, which used ERA-Interim as initial surface forcings and is constrained by ocean dynamics and ocean observations, corrected the pattern. Among the four products, ECCO and OAFlux–CERES show great similarities in the examined spatial and temporal patterns. Given that the two estimates were obtained using different approaches and based on largely independent observations, these similarities are encouraging and instructive. It is more likely that the global net air–sea heat flux does not change much during the so-called hiatus period.

Abstract

The study examined global variability of air–sea sensible heat flux (SHF) from 1980 to 2009 and the large-scale atmospheric and ocean circulations that gave rise to this variability. The contribution of high-latitude wintertime SHF was identified, and the relative importance of the effect of the sea–air temperature difference versus the effect of wind on decadal SHF variability was analyzed using an empirical orthogonal function (EOF) approach. The study showed that global SHF anomalies are strongly modulated by SHF at high latitudes (poleward of 45°) during winter seasons. Decadal variability of global wintertime SHF can be reasonably represented by the sum of two leading EOF modes, namely, the boreal wintertime SHF in the northern oceans and the austral wintertime SHF in the southern oceans. The study also showed that global wintertime SHF is modulated by the prominent modes of the large-scale atmospheric circulation at high latitudes. The increase of global SHF in the 1990s is attributable to the strengthening of the Southern Hemisphere annular mode index, while the decrease of global SHF after 2000 is due primarily to the downward trend of the Arctic Oscillation index. This study identified the important effects of wind direction and speed on SHF variability. Changes in winds modify the sea–air temperature gradient by advecting cold and dry air from continents and by imposing changes in wind-driven oceanic processes that affect sea surface temperature (SST). The pattern of air temperature anomalies dominates over the pattern of SST anomalies and dictates the pattern of decadal SHF variability.

Abstract

Westerly wind bursts (WWBs) in the equatorial Pacific occur during the development of most El Niño events and are believed to be a major factor in ENSO’s dynamics. Because of their short time scale, WWBs are normally considered part of a stochastic forcing of ENSO, completely external to the interannual ENSO variability. Recent observational studies, however, suggest that the occurrence and characteristics of WWBs may depend to some extent on the state of ENSO components, implying that WWBs, which force ENSO, are modulated by ENSO itself.

Satellite and in situ observations are used here to show that WWBs are significantly more likely to occur when the warm pool is extended eastward. Based on these observations, WWBs are added to an intermediate complexity coupled ocean–atmosphere ENSO model. The representation of WWBs is idealized such that their occurrence is modulated by the warm pool extent. The resulting model run is compared with a run in which the WWBs are stochastically applied. The modulation of WWBs by ENSO results in an enhancement of the slow frequency component of the WWBs. This causes the amplitude of ENSO events forced by modulated WWBs to be twice as large as the amplitude of ENSO events forced by stochastic WWBs with the same amplitude and average frequency. Based on this result, it is suggested that the modulation of WWBs by the equatorial Pacific SST is a critical element of ENSO’s dynamics, and that WWBs should not be regarded as purely stochastic forcing. In the paradigm proposed here, WWBs are still an important aspect of ENSO’s dynamics, but they are treated as being partially stochastic and partially affected by the large-scale ENSO dynamics, rather than being completely external to ENSO.

It is further shown that WWB modulation by the large-scale equatorial SST field is roughly equivalent to an increase in the ocean–atmosphere coupling strength, making the coupled equatorial Pacific effectively self-sustained.

Abstract

The present study used a new net surface heat flux (Q _net) product obtained from the Objective Analyzed Air–Sea Fluxes (OAFlux) project and the International Satellite Cloud Climatology Project (ISCCP) to examine two specific issues—one is to which degree Q _net controls seasonal variations of sea surface temperature (SST) in the tropical Atlantic Ocean (20°S–20°N, east of 60°W), and the other is whether the physical relation can serve as a measure to evaluate the physical representation of a heat flux product. To better address the two issues, the study included the analysis of three additional heat flux products: the Southampton Oceanographic Centre (SOC) heat flux analysis based on ship reports, and the model fluxes from the National Centers for Environmental Prediction–National Center for Atmospheric Research (NCEP–NCAR) reanalysis and the 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40). The study also uses the monthly subsurface temperature fields from the World Ocean Atlas to help analyze the seasonal changes of the mixed layer depth (h _MLD).

The study showed that the tropical Atlantic sector could be divided into two regimes based on the influence level of Q _net. SST variability poleward of 5°S and 10°N is dominated by the annual cycle of Q _net. In these regions the warming (cooling) of the sea surface is highly correlated with the increased (decreased) Q _net confined in a relatively shallow (deep) h _MLD. The seasonal evolution of SST variability is well predicted by simply relating the local Q _net with a variable h _MLD. On the other hand, the influence of Q _net diminishes in the deep Tropics within 5°S and 10°N and ocean dynamic processes play a dominant role. The dynamics-induced changes in SST are most evident along the two belts, one of which is located on the equator and the other off the equator at about 3°N in the west, which tilts to about 10°N near the northwestern African coast.

The study also showed that if the degree of consistency between the correlation relationships of Q _net, h _MLD, and SST variability serves as a measure of the quality of the Q _net product, then the Q _net from OAFlux + ISCCP and ERA-40 are most physically representative, followed by SOC. The NCEP–NCAR Q _net is least representative. It should be noted that the Q _net from OAFlux + ISCCP and ERA-40 have a quite different annual mean pattern. OAFlux + ISCCP agrees with SOC in that the tropical Atlantic sector gains heat from the atmosphere on the annual mean basis, where the ERA-40 and the NCEP–NCAR model reanalyses indicate that positive Q _net occurs only in the narrow equatorial band and in the eastern portion of the tropical basin. Nevertheless, seasonal variances of the Q _net from OAFlux + ISCCP and ERA-40 are very similar once the respective mean is removed, which explains why the two agree with each other in accounting for the seasonal variability of SST.

In summary, the study suggests that an accurate estimation of surface heat flux is crucially important for understanding and predicting SST fluctuations in the tropical Atlantic Ocean. It also suggests that future emphasis on improving the surface heat flux estimation should be placed more on reducing the mean bias.

Abstract

Daily latent and sensible heat fluxes for the Atlantic Ocean from 1988 to 1999 with 1° × 1° resolution have been recently developed at Woods Hole Oceanographic Institution (WHOI) by using a variational object analysis approach. The present study evaluated the degree of improvement made by the WHOI analysis using in situ buoy/ship measurements as verification data. The measurements were taken from the following field experiments: the five-buoy Subduction Experiment in the eastern subtropical North Atlantic, three coastal field programs in the western Atlantic, two winter cruises by R/V Knorr from the Labrador Sea Deep Convection Experiment, and the Pilot Research Moored Array in the Tropical Atlantic (PIRATA). The differences between the observed and the WHOI-analyzed fluxes and surface meteorological variables were quantified. Comparisons with the outputs from two numerical weather prediction (NWP) models were also conducted.

The mean and daily variability of the latent and sensible heat fluxes from the WHOI analysis are an improvement over the NWP fluxes at all of the measurement sites. The improved flux representation is due to the use of not only a better flux algorithm but also the improved estimates for flux-related variables. The mean differences from the observations in latent heat flux and sensible heat flux, respectively, range from 2.9 (3% of the corresponding mean measurement value) and 1.0 W m⁻² (13%) at the Subduction Experiment site, to 11.9 (13%) and 0.7 W m⁻² (11%) across the PIRATA array, to 15.9 (20%) and 10.5 W m⁻² (34%) at the coastal buoy sites, to 8.7 (7%) and 9.7 W m⁻² (6%) along the Knorr cruise tracks. The study also suggests that further improvement in the accuracy of latent and sensible heat fluxes will depend on the availability of high-quality SST observations and improved representation/observations of air humidity in the tropical Atlantic.

Abstract

Surface meteorological variables and turbulent heat fluxes in the National Centers for Environmental Prediction–National Center for Atmospheric Research reanalyses 1 and 2 (NCEP1 and NCEP2) and the analysis from the operational system of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared with high-quality moored buoy observations in regions of the Atlantic including the eastern North Atlantic, the coastal regions of the western North Atlantic, and the Tropics. The buoy latent and sensible heat fluxes are determined from buoy measurements using the recently improved Tropical Ocean Global Atmosphere Coupled Ocean–Atmosphere Response Experiment (TOGA COARE) flux algorithm.

The time mean oceanic heat loss from the model analyses is systematically overestimated in all the regions. The overestimation in latent heat loss ranges from about 14 W m⁻² (13%) in the eastern subtropical North Atlantic to about 29 W m⁻² (30%) in the Tropics to about 30 W m⁻² (49%) in the midlatitude coastal areas, where the overestimation in sensible heat flux reaches about 20 W m⁻² (60%). Depending upon the region and the NWP model, these systematic overestimations are either reduced, or change to underestimations, or remain unchanged when the TOGA COARE flux algorithm is used to recalculate the fluxes. The bias in surface meteorological variables, one of the major factors related to the biases in the revised NWP heat fluxes, varies with region and NWP analysis. Generally the temperature and humidity biases in the coastal regions are much larger than other regions. In the extratropical regions, NCEP1 and NCEP2 generally show a wet bias, which is mainly responsible for the underestimation in the revised NWP latent heat loss. In the Tropics a dry bias is found in the NWP analyses, particularly in ECMWF and NCEP2, which contributes to the overestimation in the revised NWP latent heat loss. Compared to NCEP1, NCEP2 shows less cold bias in 2-m air temperature and thus less biased sensible heat flux; NCEP2 also shows less humid bias in 2-m humidity in the extratropical regions but more dry bias in 2-m humidity in the Tropics, either of which leads to a more biased latent heat flux in NCEP2.

Despite the significant biases in the NWP surface fields and the poor representation of short-time sea surface temperature variability, the NWP models are able to represent the dominant short-time variability in other basic variables and thus the variability in heat fluxes in the wintertime coastal regions of the western North Atlantic (on timescales of 3–4 days and 1 week) and the northern and southern subtropical regions (on a timescale of about 2 weeks), but ECMWF and particularly the NCEP analyses do not represent well the 2–3-week variability in the tropical Atlantic.

Abstract

A new daily latent and sensible flux product developed at the Woods Hole Oceanographic Institution (WHOI) with 1° × 1° resolution for the Atlantic Ocean (65°S–65°N) for the period from 1988 to 1999 was presented. The flux product was developed by using a variational objective analysis approach to obtain best estimates of the flux-related basic surface meteorological variables (e.g., wind speed, air humidity, air temperature, and sea surface temperature) through synthesizing satellite data and outputs of numerical weather prediction (NWP) models at the National Centers for Environmental Prediction (NCEP) and the European Centre for Medium-Range Weather Forecasts (ECMWF). The state-of-the-art bulk flux algorithm 2.6a, developed from the field experiments of the Coupled Ocean–Atmosphere Response Experiment (COARE), was applied to compute the flux fields.

The study focused on analyzing the mean field properties of the WHOI daily latent and sensible heat fluxes and their comparisons with the ship-based climatology from the Southampton Oceanography Centre (SOC) and NWP outputs. It is found that the WHOI yearly mean fluxes are consistent with the SOC climatology in both structure and amplitude, but the WHOI yearly mean basic variables are not always consistent with SOC; the better agreement in the fluxes is due to the effects of error compensation during variable combinations. Both ECMWF and NCEP–Department of Energy (DOE) Atmospheric Model Intercomparison Project (AMIP) Reanalysis-2 (NCEP2) model data have larger turbulent heat loss (∼20 W m⁻²) than the WHOI product. Nevertheless, the WHOI fluxes agree well with the NCEP-2 Reanalysis fluxes in structure and the trend of year-to-year variations, but not with the ECMWF operational outputs; the latter have a few abrupt changes coinciding with the modifications in the model forecast–analysis system. The degree of impact of the model changes on the basic variables is not as dramatic, a factor that justifies the inclusion of the basic variables, not the fluxes, from the ECMWF operational model in the synthesis. The flux algorithms of the two NWP models give a larger latent and sensible heat loss. Recalculating the NWP fluxes using the COARE algorithm considerably reduces the strength but does not replicate the WHOI results. The present analysis could not quantify the degree of improvement in the mean aspect of the WHOI daily flux fields as accurate basinwide verification data are lacking.

This study is the first to demonstrate that the synthesis approach is a useful tool for combining the NWP and satellite data sources and improving the mean representativeness of daily basic variable fields and, hence, the daily flux fields. It is anticipated that such an approach may become increasingly relied upon in the preparation of future high-quality flux products.

Abstract

This study investigated the accuracy and physical representation of air–sea surface heat flux estimates for the Indian Ocean on annual, seasonal, and interannual time scales. Six heat flux products were analyzed, including the newly developed latent and sensible heat fluxes from the Objectively Analyzed Air–Sea Heat Fluxes (OAFlux) project and net shortwave and longwave radiation results from the International Satellite Cloud Climatology Project (ISCCP), the heat flux analysis from the Southampton Oceanography Centre (SOC), the National Centers for Environmental Prediction reanalysis 1 (NCEP1) and reanalysis-2 (NCEP2) datasets, and the European Centre for Medium-Range Weather Forecasts operational (ECMWF-OP) and 40-yr Re-Analysis (ERA-40) products.

This paper presents the analysis of the six products in depicting the mean, the seasonal cycle, and the interannual variability of the net heat flux into the ocean. Two time series of in situ flux measurements, one taken from a 1-yr Arabian Sea Experiment field program and the other from a 1-month Joint Air–Sea Monsoon Interaction Experiment (JASMINE) field program in the Bay of Bengal were used to evaluate the statistical properties of the flux products over the measurement periods. The consistency between the six products on seasonal and interannual time scales was investigated using a standard deviation analysis and a physically based correlation analysis.

The study has three findings. First of all, large differences exist in the mean value of the six heat flux products. Part of the differences may be attributable to the bias in the numerical weather prediction (NWP) models that underestimates the net heat flux into the Indian Ocean. Along the JASMINE ship tracks, the four NWP modeled mean fluxes all have a sign opposite to the observations, with NCEP1 being underestimated by 53 W m⁻² (the least biased) and ECMWF-OP by 108 W m⁻² (the most biased). At the Arabian Sea buoy site, the NWP mean fluxes also have an underestimation bias, with the smallest bias of 26 W m⁻² (ERA-40) and the largest bias of 69 W m⁻² (NCEP1). On the other hand, the OAFlux+ISCCP has the best comparison at both measurement sites. Second, the bias effect changes with the time scale. Despite the fact that the mean is biased significantly, there is no major bias in the seasonal cycle of all the products except for ECMWF-OP. The latter does not have a fixed mean due to the frequent updates of the model platform. Finally, among the four products (OAFlux+ISCCP, ERA-40, NCEP1, and NCEP2) that can be used for studying interannual variability, OAFlux+ISCCP and ERA-40 Q _net have good consistency as judged from both statistical and physical measures. NCEP1 shows broad agreement with the two products, with varying details. By comparison, NCEP2 is the least representative of the Q _net variabilities over the basin scale.