Search for: All records

The 1991–2020 climate normals for sea surface temperature (SST) are computed based on the NOAA Daily Optimum Interpolation SST dataset. This is the first time that high‐resolution SST normals with global coverage can be achieved in the satellite SST era. Normals are one of the fundamental parameters in describing and understanding weather and climate and provide decision‐making information to industry, public, and scientific communities. This product suite includes SST mean, standard deviation, count and extreme parameters at daily, monthly, seasonal and annual time scales on 0.25° spatial grids. The main feature of the SST mean state revealed by the normals is that in the Tropics, the Indo‐Pacific Ocean is dominated by the warm pool (SST ≥ 28°C) while the eastern Pacific is characterized by the cold tongue (SST ≤ 24°C); in the midlatitudes, SSTs are in zonal patterns with high meridional gradients. Daily SST standard deviations are generally small (<1.0°C) except in frontal zones (>1.5°C) mostly associated with ocean currents such as the Gulf Stream, Kuroshio and Equatorial Currents. Compared to the 1982–2011 climatology, the 1991–2020 mean SSTs increased over most global areas but obvious cooling is seen in the Southern Ocean, eastern tropical South Pacific Ocean and North Atlantic warming hole. The Indo‐Pacific warm pool (IPWP) is found to have strengthened in both intensity and coverage since 1982–2011. By a count parameter criterion of ≥300 days annually with SST ≥ 28°C, the IPWP coverage increased 33% from 1982–2011 to 1991–2020. The global mean SST of 1991–2020 is warmer than that of 1982–2011, and the warming rate over 1991–2020 doubles that over 1901–2020.

 

This study extends initial work by Sun and Penny and Sun et al. to explore the inclusion of path information from surface drifters using an augmented-state Lagrangian data assimilation based on the local ensemble transform Kalman filter (LETKF-LaDA) with vertical localization to improve analysis of the ocean. The region of interest is the Gulf of Mexico during the passage of Hurricane Isaac in the summer of 2012. Results from experiments with a regional ocean model at eddy-permitting and eddy-resolving model resolutions are used to quantify improvements to the analysis of sea surface velocity, sea surface temperature, and sea surface height in a data assimilation system. The data assimilation system assimilates surface drifter positions, as well as vertical profiles of temperature and salinity. Data were used from drifters deployed as a part of the Grand Lagrangian Deployment beginning 20 July 2012. Comparison of experiment results shows that at both eddy-permitting and eddy-resolving horizontal resolutions Lagrangian assimilation of drifter positions significantly improves analysis of the ocean state responding to hurricane conditions. These results, which should be applicable to other tropical oceans such as the Bay of Bengal, open new avenues for estimating ocean initial conditions to improve tropical cyclone forecasting.

 

Our study shows that the intercomparison among sea surface temperature (SST) products is influenced by the choice of SST reference, and the interpolation of SST products. The influence of reference SST depends on whether the reference SSTs are averaged to a grid or in pointwise in situ locations, including buoy or Argo observations, and filtered by first-guess or climatology quality control (QC) algorithms. The influence of the interpolation depends on whether SST products are in their original grids or preprocessed into common coarse grids. The impacts of these factors are demonstrated in our assessments of eight widely used SST products (DOISST, MUR25, MGDSST, GAMSSA, OSTIA, GPB, CCI, CMC) relative to buoy observations: (i) when the reference SSTs are averaged onto 0.25° × 0.25° grid boxes, the magnitude of biases is lower in DOISST and MGDSST (<0.03°C), and magnitude of root-mean-square differences (RMSDs) is lower in DOISST (0.38°C) and OSTIA (0.43°C); (ii) when the same reference SSTs are evaluated at pointwise in situ locations, the standard deviations (SDs) are smaller in DOISST (0.38°C) and OSTIA (0.39°C) on 0.25° × 0.25° grids; but the SDs become smaller in OSTIA (0.34°C) and CMC (0.37°C) on products’ original grids, showing the advantage of those high-resolution analyses for resolving finer-scale SSTs; (iii) when a loose QC algorithm is applied to the reference buoy observations, SDs increase; and vice versa; however, the relative performance of products remains the same; and (iv) when the drifting-buoy or Argo observations are used as the reference, the magnitude of RMSDs and SDs become smaller, potentially due to changes in observing intervals. These results suggest that high-resolution SST analyses may take advantage in intercomparisons.

Intercomparisons of gridded SST products be affected by how the products are compared with in situ observations: whether the products are in coarse (0.25°) or original (0.05°–0.10°) grids, whether the in situ SSTs are in their reported locations or gridded and how they are quality controlled, and whether the biases of satellite SSTs are corrected by localized matchups or large-scale patterns. By taking all these factors into account, our analyses indicate that the NOAA DOISST is among the best SST products for the long period (1981–present) and relatively coarse (0.25°) resolution that it was designed for.

 

This paper describes the new Regional Arctic Ocean/sea ice Reanalysis (RARE) with a domain that spans a subpolar/polar cap poleward of 45°N. Sequential data assimilation constrains temperature and salinity using World Ocean Database profiles as well as in situ and satellite SST, and PIOMAS sea ice thickness estimates. The 41-yr (1980–2020) RARE1.15.2 reanalysis with resolution varying between 2 and 5 km horizontally and 1–10 m vertically in the upper 100 m is examined. To explore the impact of resolution RARE1.15.2 is compared to a coarser-resolution SODA3.15.2, which uses the same modeling and data assimilation system. Improving resolution in the reanalysis system improves agreement with observations. It produces stronger more compact currents, enhances eddy kinetic energy, and strengthens along-isopycnal heat and salt transports, but reduces vertical exchanges and thus strengthens upper ocean haline stratification. RARE1.15.2 and SODA3.15.2 are also compared to the Hadley Center EN4.2.2 statistical objective analysis. In regions of reasonable data coverage such as the Nordic seas the three products produce similar time-mean distributions of temperature and salinity. But in regions of poor coverage and in regions where the coverage changes in time EN4.2.2 suffers more from those inhomogeneities. Finally, the impact on the Arctic of interannual temperature fluctuations in the subpolar gyres on the Arctic Ocean is compared. The influence of the subpolar North Pacific is limited to a region surrounding Bering Strait. The influence of the subpolar North Atlantic, in contrast, spreads throughout the Nordic seas and Barents Sea in all three products within two years.

The Arctic Ocean/sea ice system plays crucial roles in climate variability and change by controlling the northern end of the oceanic overturning circulation, the equator to pole air pressure gradient, and Earth’s energy balance. Yet the historical ocean observation set is sparse and inhomogeneous, while ocean dynamics has challengingly fine horizontal and vertical scales. This paper introduces a new Regional Arctic Ocean/sea ice Reanalysis (RARE) whose goal is to use the combined constraints of mesoscale ocean dynamics, historical observations, surface meteorology, and continental runoff in a data assimilation framework to reconstruct historical variability. RARE is used to produce a 41-yr ocean/sea ice reanalysis 1980–2020 whose results are described here.

 

We estimate the effect of model deficiencies in the Global Forecast System that lead to systematic forecast errors, as a first step toward correcting them online (i.e., within the model) as in Danforth & Kalnay (2008a, 2008b). Since the analysis increments represent the corrections that new observations make on the 6 h forecast in the analysis cycle, we estimate the model bias corrections from the time average of the analysis increments divided by 6 h, assuming that initial model errors grow linearly and first ignoring the impact of observation bias. During 2012–2016, seasonal means of the 6 h model bias are generally robust despite changes in model resolution and data assimilation systems, and their broad continental scales explain their insensitivity to model resolution. The daily bias dominates the submonthly analysis increments and consists primarily of diurnal and semidiurnal components, also requiring a low dimensional correction. Analysis increments in 2015 and 2016 are reduced over oceans, which we attribute to improvements in the specification of the sea surface temperatures. These results provide support for future efforts to make online correction of the mean, seasonal, and diurnal and semidiurnal model biases of Global Forecast System to reduce both systematic and random errors, as suggested by Danforth & Kalnay (2008a, 2008b). It also raises the possibility that analysis increments could be used to provide guidance in testing new physical parameterizations.