A global ocean temperature and altimeter data assimilation system for studies of climate variability

 An ocean data assimilation (ODA) system which can assimilate both temperature and altimeter observations has been applied to the global ocean and tested between January 1993–October 1996. A statistical method has been used to convert sea surface height (SSH) anomalies observations from TOPEX/POSEIDON into synthetic temperature profiles. The innovative aspect of this method is the introduction of time dependency in the correlations used to transform the altimeter observations into temperature corrections. The assimilation system is based on a univariate variational optimal interpolation scheme applied to assimilate both in situ and synthetic temperature profiles. In addition, a longer global analysis for the upper-ocean temperature starting from January 1979 and ending November 1997, has been produced to examine the skill of sea temperature assimilation with a rather simple and practical method. The temperature analysis shows encouraging improvement over a corresponding ocean simulation when compared to independent (not assimilated) temperature data both at seasonal and interannual time scales. However, the univariate data assimilation of hydrographic data does not result in an improvement of the velocity field. In fact the assimilation of sparse in situ data can introduce unrealistic spatial variability in the temperature field which affects the velocity field in a negative way. This deficiency is partially overcome when we also assimilate altimeter observations since the coverage is complete and uniform for this data. In particular, our study shows that temperature corrections due to the altimeter signal have a positive impact on the current system in the tropical Pacific. Received: 28 May 2000 / Accepted: 6 November 2000     

Interannual variability of North Pacific eastern subtropical mode water formation in the 1990s derived from a 4-dimensional variational ocean data assimilation experiment

The interannual variability of eastern subtropical mode water (ESTMW) formation in the North Pacific is examined using a new ocean dataset constructed by a 4-dimensional variational data assimilation experiment covering the decade of the 1990s. The volume of newly formed ESTMW varies due to interannual variability in the following three physical processes taking place in the surface layer: (1) convergence in the transport of surface saline water induced by Ekman flow in the vicinity of the formation region, (2) thermal stratification in the preconditioning phase in association with the insolation anomaly induced largely by low-level cloud coverage, and (3) wintertime surface cooling in the eastern subtropics. We find that, in addition to the surface forcing, the properties of both the ESTMW and the upper mixed-layer water are broadly controlled by the volume of the new ESTMW component, and that the variations in the upper mixed-layer water affect the properties of ESTMW formed in the following winter. Due to the combined effect of these processes, the ESTMW subducts down to subsurface layers with a wide range of σ_θ values lying between 24.8 and 25.4 and with significant interannual variation in water mass formation.     

An ocean data assimilation system in the Indian Ocean and west Pacific Ocean

The development and application of a regional ocean data assimilation system are among the aims of the Global Ocean Data Assimilation Experiment. The ocean data assimilation system in the regions including the Indian and West Pacific oceans is an endeavor motivated by this goal. In this study, we describe the system in detail. Moreover, the reanalysis in the joint area of Asia, the Indian Ocean, and the western Pacific Ocean(hereafter AIPOcean) constructed using multi-year model integration with data assimilation is used to test the performance of this system. The ocean model is an eddy-resolving,hybrid coordinate ocean model. Various types of observations including in-situ temperature and salinity profiles(mechanical bathythermograph, expendable bathythermograph, Array for Real-time Geostrophic Oceanography, Tropical Atmosphere Ocean Array, conductivity–temperature–depth, station data), remotely-sensed sea surface temperature, and altimetry sea level anomalies, are assimilated into the reanalysis via the ensemble optimal interpolation method. An ensemble of model states sampled from a long-term integration is allowed to change with season, rather than remaining stationary. The estimated background error covariance matrix may reasonably reflect the seasonality and anisotropy. We evaluate the performance of AIPOcean during the period 1993–2006 by comparisons with independent observations, and some reanalysis products. We show that AIPOcean reduces the errors of subsurface temperature and salinity, and reproduces mesoscale eddies. In contrast to ECCO and SODA products, AIPOcean captures the interannual variability and linear trend of sea level anomalies very well. AIPOcean also shows a good consistency with tide gauges.     

基于弱耦合资料同化的冬季北大西洋涛动年际变率预测

北大西洋涛动作为冬季北大西洋地区大气环流的主模态之一,其年际变率对全球许多地区气候变率具有重要影响,但目前其预测技巧并不高。采用降维投影四维变分同化方法,在耦合模式中建立了基于全球大气资料的弱耦合资料同化系统,直接同化月平均再分析资料,并进行了年代际后报试验。结果表明,通过耦合资料同化的手段,可以显著提升耦合模式对冬季北大西洋涛动年际变率及其相关的欧洲北部、美国东部、欧亚大陆北部的冬季近地面温度年际变率的后报效果,相关系数均超过0.05显著水平t检验。该后报效果的改进主要与在耦合同化过程中通过耦合模式中自由发展的海-气相互作用将大气的观测信息储存在耦合模式的海洋分量中,改进了冬季北大西洋地区海表温度“三极”型分布的时空变率及其时间序列的后报效果有关。该研究强调了耦合模式初始状态的准确度对提升冬季北大西洋涛动年际变率的后报技巧具有重要作用。     

Interdecadal climate variability in the subpolar North Atlantic

The statistical relationships between various components of the subpolar North Atlantic air-sea-ice climate system are reexamined in order to investigate potential processes involved in interdecadal climate variability. It is found that sea surface temperature anomalies concentrated in the Labrador Sea region have a strong impact upon atmospheric sea level pressure anomalies over Greenland, which in turn influence the transport of freshwater and ice anomalies out of the Arctic Ocean, via Fram Strait. These freshwater and ice anomalies are advected around the subpolar gyre into the Labrador Sea affecting convection and the formation of Labrador Sea Water. This has an impact upon the transport of North Atlantic Current water into the subpolar gyre and thus, also upon sea surface temperatures in the region. An interdecadal negative feedback loop is therefore proposed as an internal source of climate variability within the subpolar North Atlantic. Through the lags associated with the correlations between different climatic components, observed horizontal advection time scales, and the use of Boolean delay equation models, the time scale for one cycle of this feedback loop is determined to have a period of about 21 years.     

Three inflow pathways of the Indonesian throughflow as seen from the simple ocean data assimilation

The SODA product is used to investigate three Indonesian throughflow (ITF) branches: the flow through the Makassar Strait; through the South China Sea; and through the eastern Indonesian basins. The results reveal strong interannual variation in the Makassar Strait and the eastern Indonesian basins throughflow. Inspection of vertically integrated dynamic height (0–1000 db), a proxy of transport function, suggests that this interannual variation can be traced to the New Guinea Coastal Current, indicative of a strong influence of the South Pacific. The vertically integrated dynamic height along the south Java coast is related to variation in the North Pacific and in particular near the east coast of Mindanao Island, whereas the vertically integrated dynamic height along the coast of West Australia is related to variation in the South Pacific, and in particular near the coast of New Guinea. The integrated dynamic height difference between the Java and New Guinea coast appears to be a good proxy of ITF transport on the interannual time scale. Regression analysis shows a phase dependence of the three ITF pathways on the Nino3.4 index. Decoupling of current anomalies between the surface and subsurface layers is identified in the developing and mature phase of El Nino, reflecting different effects of local and remote forcing through oceanic pathways at the Makassar Strait and eastern Indonesian basins.     

The relative importance of tropical variability forced from the North Pacific through ocean pathways

To what extent is tropical variability forced from the North Pacific through ocean pathways relative to locally generated variability and variability forced through the atmosphere? To address this question, in this study we use an anomaly-coupled model, consisting of a global, atmospheric general circulation model and a 4½-layer, reduced-gravity, Pacific-Ocean model. Three solutions are obtained; with coupling over the entire basin (CNT), with coupling confined to the tropics and wind stress and heat fluxes in the North and South Pacific specified by climatology (TP), and with coupling confined to the Tropics and wind stress and heat fluxes in the North Pacific specified by output from CNT (NPF). It is found that there are two distinct signals forced in the North Pacific that can impact the tropics through ocean pathways. These two signals are forced by wind stress and surface heat flux anomalies in the subtropical North Pacific. The first signal is relatively fast, impacts tropical variability less than a year after forcing, is triggered from November to March, and propagates as a first-mode baroclinic Rossby wave. The second signal is only triggered during springtime when buoyancy forcing can effectively generate higher-order baroclinic modes through subduction anomalies into the permanent thermocline, and it reaches the equator 4–5 years after forcing. The slow signal is found to initiate tropical variability more efficiently than the fast signal with one standard deviation in subtropical zonal wind stress forcing tropical SST anomalies centered on the equator at 135°W of approximately 0.5°C. Allowing extratropically forced tropical variability is found to shift primarily 2-year ENSO variability in a tropics-alone simulation to a more realistic range of 2–6 years.     

On the seesaw in interannual variability of upper ocean heat advection between the North Atlantic Subpolar Gyre and the Nordic Seas

In this study we investigate the interannual variability of the heat content of the upper North Atlantic and the main factors, which influence the observed variability: the ocean-atmosphere heat exchange and the ocean heat transport. The data from the combined in situ and satellite dataset ARMOR-3D, and from the ocean reanalyzes ORAS5 and SODA3 (two versions) show a similar decadal interannual variability of the heat content, as well as of oceanic heat transports the study regions, though the mean values may differ. The observed variations are linked to the North Atlantic Oscillation Index (NAOI). The current velocity of the North Atlantic, East Greenland and Labrador currents of the Subpolar Gyre increases with the NAOI, but the velocity of the Irminger and West Greenland currents decreases. This forms a seesaw of heat advection by the North Atlantic Current between the Subpolar Gyre and Nordic Seas. In the Subpolar gyre during the periods of high NAOI, this negative anomaly of the oceanic heat convergence adds to the intensified sea-surface heat release to the atmosphere, together effectively reducing the upper ocean heat content. The upper ocean heat content of the Norwegian Sea shows practically no link to the NAOI, in spite of a somewhat larger oceanic heat flux across its southern boundary linked to high NAOI.     

Spatiotemporal variability of meridional mass transport in the North Atlantic

Analyzed is the interannual variability of the meridional mass transport ψ_S in the North Atlantic based on the Sverdrup relation. The continuous (1980–2005) monthly wind stress dataset with the spatial resolution of 1 × 1° was used as the initial data. Sverdrup transport analysis performed for different latitudinal transects within the North Atlantic subtropical gyre demonstrated that the maximum long-term Sverdrup transport (?25.2 Sv) can be found at 33°N. Studied is a mechanism of the interaction between the meridional Sverdrup transport and the water flow in the Florida Strait. The significant correlation coefficient (0.5) is revealed for the Florida Strait water discharge and the mass transport at 27°N. Analyzed is the relationship between ψ_S and the North Atlantic Oscillation index and the statistically significant correlation coefficient (0.45) is obtained for the Sverdrup transport at 49°N.     

Impact of freshwater discharge from the Greenland ice sheet on North Atlantic climate variability

Using a coupled ocean–atmosphere general circulation model, we investigated the impact of Greenland ice sheet melting on North Atlantic climate variability. The positive-degree day (PDD) method was incorporated into the model to control continental ice melting (PDD run). Models with and without the PDD method produce a realistic pattern of North Atlantic sea surface temperature (SST) variability that fluctuates from decadal to multidecadal periods. However, the interdecadal variability in PDD run is significantly dominated in the longer time scale compared to that in the run without PDD method. The main oscillatory feature in these experiments likely resembles the density-driven oscillatory mode. A reduction in the ocean density over the subpolar Atlantic results in suppression of the Atlantic Meridional Overturning Circulation (AMOC), leading to a cold SST due to a weakening of northward heat transport. The decreased surface evaporation associated with the cold SST further reduces the ocean density and thus, simultaneously acts as a positive feedback mechanism. The southward meridional current associated with the suppressed AMOC causes a positive tendency in the ocean density through density advection, which accounts for the phase transition of this oscillatory mode. The Greenland ice melting process reduces the mean meridional current and meridional density gradient because of additional fresh water flux, which suppress the delayed negative feedback due to meridional density advection. As a result, the oscillation period becomes longer and the transition is more delayed.     

The North Atlantic Oscillation and oceanic precipitation variability

Global North Atlantic Oscillation (NAO) oceanic precipitation features in the latter half of the twentieth century are documented based on the intercomparison of multiple state-of-the-art precipitation datasets and the analysis of the NAO atmospheric circulation and SST anomalies. Most prominent precipitation anomalies occur over the ocean in the North Atlantic, where in winter a “quadrupole-like” pattern is found with centers in the western tropical Atlantic, sub-tropical Atlantic, high-latitude eastern Atlantic and over the Labrador Sea. The extent of the sub-tropical and high-latitude center and the amount of explained variance (over 50%) are quite remarkable. However, the tropical Atlantic center is probably the most intriguing feature of this pattern apparently linking the NAO with ITCZ variability. In summer, the pattern is “tripole-like” with centers in the eastern Mediterranean Sea, the North Sea/Baltic Sea and in the sub-polar Atlantic. In the eastern Indian Ocean, the correlation is positive in winter and negative in summer, with some link to ENSO variability. The sensitivity of these patterns to the choice of the NAO index is minor in winter while quite important in summer. Interannual NAO precipitation anomalies have driven similar fresh water variations in these “key” regions. In the sub-tropical and high-latitude Atlantic in winter precipitation anomalies have been roughly 15 and 10% of climatology per unit change of the NAO, respectively. Decadal changes of the NAO during the last 50 years have also influenced precipitation and fresh water flux at these time-scales, with values lower (higher) than usual in the high-latitude eastern North Atlantic (Labrador Sea) in the 1960s and the late 1970s, and an opposite situation since the early 1980s; in summer the North Sea/Baltic region has been drier than usual during the period 1965–1975 when the NAO was generally positive.     

Evaluation of ocean data assimilation in CAS-ESM-C: Constraining the SST field

A weakly coupled assimilation system, in which SST observations are assimilated into a coupled climate model(CASESM-C) through an ensemble optimal interpolation scheme, was established. This system is a useful tool for historical climate simulation, showing substantial advantages, including maintaining the atmospheric feedback, and keeping the oceanic fields from drifting far away from the observation, among others. During the coupled model integration, the bias of both surface and subsurface oceanic fields in the analysis can be reduced compared to unassimilated fields. Based on 30 model years of output from the system, the climatology and interannual variability of the climate system were evaluated. The results showed that the system can reasonably reproduce the climatological global precipitation and SLP, but it still suffers from the double ITCZ problem. Besides, the ENSO footprint, which is revealed by ENSO-related surface air temperature, geopotential height and precipitation during El Ni ?no evolution, is basically reproduced by the system. The system can also simulate the observed SST–rainfall relationships well on both interannual and intraseasonal timescales in the western North Pacific region, in which atmospheric feedback is crucial for climate simulation.     

A 20-year coupled ocean-sea ice-atmosphere variability mode in the North Atlantic in an AOGCM

In order to understand potential predictability of the ocean and climate at the decadal time scales, it is crucial to improve our understanding of internal variability at this time scale. Here, we describe a 20-year mode of variability found in the North Atlantic in a 1,000-year pre-industrial simulation of the IPSL-CM5A-LR climate model. This mode involves the propagation of near-surface temperature and salinity anomalies along the southern branch of the subpolar gyre, leading to anomalous sea-ice melting in the Nordic Seas, which then forces sea-level pressure anomalies through anomalous surface atmospheric temperatures. The wind stress associated to this atmospheric structure influences the strength of the East Greenland Current across the Denmark Strait, which, in turn, induces near-surface temperature and salinity anomalies of opposite sign at the entrance of the Labrador Sea. This starts the second half of the cycle after approximatively 10 years. The time scale of the cycle is thus essentially set by advection of tracers along the southern branch of the subpolar gyre, and by the time needed for anomalous East Greenland Current to accumulate heat and freshwater anomalies at the entrance of the Labrador Sea. The Atlantic meridional overturning circulation (AMOC) does not play a dominant role in the mode that is confined in the subpolar North Atlantic, but it also has a 20-year preferred timescale. This is due to the influence of the propagating salinity anomalies on the oceanic deep convection. The existence of this preferred timescale has important implications in terms of potential predictability of the North Atlantic climate in the model, although its realism remains questionable and is discussed.     

A predictability study of simulated North Atlantic multidecadal variability

 The North Atlantic is one of the few places on the globe where the atmosphere is linked to the deep ocean through air–sea interaction. While the internal variability of the atmosphere by itself is only predictable over a period of one to two weeks, climate variations are potentially predictable for much longer periods of months or even years because of coupling with the ocean. This work presents details from the first study to quantify the predictability for simulated multidecadal climate variability over the North Atlantic. The model used for this purpose is the GFDL coupled ocean-atmosphere climate model used extensively for studies of global warming and natural climate variability. This model contains fluctuations of the North Atlantic and high-latitude oceanic circulation with variability concentrated in the 40–60 year range. Oceanic predictability is quantified through analysis of the time-dependent behavior of large-scale empirical orthogonal function (EOF) patterns for the meridional stream function, dynamic topography, 170 m temperature, surface temperature and surface salinity. The results indicate that predictability in the North Atlantic depends on three main physical mechanisms. The first involves the oceanic deep convection in the subpolar region which acts to integrate atmospheric fluctuations, thus providing for a red noise oceanic response as elaborated by Hasselmann. The second involves the large-scale dynamics of the thermohaline circulation, which can cause the oceanic variations to have an oscillatory character on the multidecadal time scale. The third involves nonlocal effects on the North Atlantic arising from periodic anomalous fresh water transport advecting southward from the polar regions in the East Greenland Current. When the multidecadal oscillatory variations of the thermohaline circulation are active, the first and second EOF patterns for the North Atlantic dynamic topography have predictability time scales on the order of 10–20 y, whereas EOF-1 of SST has predictability time scales of 5–7 y. When the thermohaline variability has weak multidecadal power, the Hasselmann mechanism is dominant and the predictability is reduced by at least a factor of two. When the third mechanism is in an extreme phase, the North Atlantic dynamic topography patterns realize a 10–20 year predictability time scale. Additional analysis of SST in the Greenland Sea, in a region associated with the southward propagating fresh water anomalies, indicates the potential for decadal scale predictability for this high latitude region as well. The model calculations also allow insight into regional variations of predictability, which might be useful information for the design of a monitoring system for the North Atlantic. Predictability appears to break down most rapidly in regions of active convection in the high-latitude regions of the North Atlantic. Received: 28 October 1996 / Accepted: 21 March 1997     

Quantifying atmosphere and ocean origins of North American precipitation variability

How atmospheric and oceanic processes control North American precipitation variability has been extensively investigated, and yet debates remain. Here we address this question in a 50 km-resolution flux-adjusted global climate model. The high spatial resolution and flux adjustment greatly improve the model’s ability to realistically simulate North American precipitation, the relevant tropical and midlatitude variability and their teleconnections. Comparing two millennium-long simulations with and without an interactive ocean, we find that the leading modes of North American precipitation variability on seasonal and longer timescales exhibit nearly identical spatial and spectral characteristics, explained fraction of total variance and associated atmospheric circulation. This finding suggests that these leading modes arise from internal atmospheric dynamics and atmosphere-land coupling. However, in the fully coupled simulation, North American precipitation variability still correlates significantly with tropical ocean variability, consistent with observations and prior literature. We find that tropical ocean variability does not create its own type of atmospheric variability but excites internal atmospheric modes of variability in midlatitudes. This oceanic impact on North American precipitation is secondary to atmospheric impacts based on correlation. However, relative to the simulation without an interactive ocean, the fully coupled simulation amplifies precipitation variance over southwest North America (SWNA) during late spring to summer by up to 90%. The amplification is caused by a stronger variability in atmospheric moisture content that is attributed to tropical Pacific sea surface temperature variability. Enhanced atmospheric moisture variations over the tropical Pacific are transported by seasonal mean southwesterly winds into SWNA, resulting in larger precipitation variance.

     

A new North Atlantic Oscillation index and its variability

A new North Atlantic Oscillation (NAO) index, the NAOI, is defined as the differences of normalized sea level pressures regionally zonal-averaged over a broad range of longitudes 80°W-30°E. A comprehensive comparison of six NAO indices indicates that the new NAOI provides a more faithful representation of the spatial-temporal variability associated with the NAO on all timescales. A very high signal-to-noise ratio for the NAOI exists for all seasons, and the life cycle represented by the NAOI describes well the seasonal migration for action centers of the NAO. The NAOI captures a larger fraction of the variance of sea level pressure over the North Atlantic sector (20°-90°N, 80°W-30°E), on average 10% more than any other NAO index. There are quite different relationships between the NAOI and surface air temperature during winter and summer. A novel feature, however, is that the NAOI is significantly negative correlated with surface air temperature over the North Atlantic Ocean between 10°-25°N and     

Formation and pathways of North Atlantic Deep Water in a coupled ice–ocean model of the Arctic–North Atlantic Oceans

We investigate the formation process and pathways of deep water masses in a coupled ice–ocean model of the Arctic and North Atlantic Oceans. The intent is to determine the relative roles of these water masses from the different source regions (Arctic Ocean, Nordic Seas, and Subpolar Atlantic) in the meridional overturning circulation. The model exhibits significant decadal variability in the deep western boundary current and the overturning circulation. We use detailed diagnostics to understand the process of water mass formation in the model and the resulting effects on the North Atlantic overturning circulation. Particular emphasis is given to the multiple sources of North Atlantic Deep Water, the dominant deep water masses of the world ocean. The correct balance of Labrador Sea, Greenland Sea and Norwegian Sea sources is difficult to achieve in climate models, owing to small-scale sinking and convection processes. The global overturning circulation is described as a function of potential temperature and salinity, which more clearly signifies dynamical processes and clarifies resolution problems inherent to the high latitude oceans. We find that fluxes of deep water masses through various passages in the model are higher than observed estimates. Despite the excessive volume flux, the Nordic Seas overflow waters are diluted by strong mixing and enter the Labrador Sea at a lighter density. Through strong subpolar convection, these waters along with other North Atlantic water masses are converted into the densest waters [similar density to Antarctic Bottom Water (AABW)] in the North Atlantic. We describe the diminished role of salinity in the Labrador Sea, where a shortage of buoyant surface water (or excess of high salinity water) leads to overly strong convection. The result is that the Atlantic overturning circulation in the model is very sensitive to the surface heat flux in the Labrador Sea and hence is correlated with the North Atlantic Oscillation. As strong subpolar convection is found in other models, we discuss broader implications.     

Stochastically-forced multidecadal variability in the North Atlantic: a model study

Observations show a multidecadal signal in the North Atlantic ocean, but the underlying mechanism and cause of its timescale remain unknown. Previous studies have suggested that it may be driven by the North Atlantic Oscillation (NAO), which is the dominant pattern of winter atmospheric variability. To further address this issue, the global ocean general circulation model, Nucleus for European Modelling of the Ocean (NEMO), is driven using a 2,000 years long white noise forcing associated with the NAO. Focusing on key ocean circulation patterns, we show that the Atlantic Meridional Overturning Circulation (AMOC) and Sub-polar gyre (SPG) strength both have enhanced power at low frequencies but no dominant timescale, and thus provide no evidence for a oscillatory ocean-only mode of variability. Instead, both indices respond linearly to the NAO forcing, but with different response times. The variability of the AMOC at 30°N is strongly enhanced on timescales longer than 90 years, while that of the SPG strength starts increasing at 15 years. The different response characteristics are confirmed by constructing simple statistical models that show AMOC and SPG variability can be related to the NAO variability of the previous 53 and 10 winters, respectively. Alternatively, the AMOC and the SPG strength can be reconstructed with Auto-regressive (AR) models of order seven and five, respectively. Both statistical models reconstruct interannual and multidecadal AMOC variability well, while on the other hand, the AR(5) reconstruction of the SPG strength only captures multidecadal variability. Using these methods to reconstruct ocean variables can be useful for prediction and model intercomparision.