Using seasonal hindcasts to understand the origin of the equatorial cold tongue bias in CGCMs and its impact on ENSO

The cold equatorial SST bias in the tropical Pacific that is persistent in many coupled OAGCMs severely impacts the fidelity of the simulated climate and variability in this key region, such as the ENSO phenomenon. The classical bias analysis in these models usually concentrates on multi-decadal to centennial time series needed to obtain statistically robust features. Yet, this strategy cannot fully explain how the models errors were generated in the first place. Here, we use seasonal re-forecasts (hindcasts) to track back the origin of this cold bias. As such hindcasts are initialized close to observations, the transient drift leading to the cold bias can be analyzed to distinguish pre-existing errors from errors responding to initial ones. A time sequence of processes involved in the advent of the final mean state errors can then be proposed. We apply this strategy to the ENSEMBLES-FP6 project multi-model hindcasts of the last decades. Four of the five AOGCMs develop a persistent equatorial cold tongue bias within a few months. The associated systematic errors are first assessed separately for the warm and cold ENSO phases. We find that the models are able to reproduce either El Niño or La Niña close to observations, but not both. ENSO composites then show that the spurious equatorial cooling is maximum for El Niño years for the February and August start dates. For these events and at this time of the year, zonal wind errors in the equatorial Pacific are present from the beginning of the simulation and are hypothesized to be at the origin of the equatorial cold bias, generating too strong upwelling conditions. The systematic underestimation of the mixed layer depth in several models can also amplify the growth of the SST bias. The seminal role of these zonal wind errors is further demonstrated by carrying out ocean-only experiments forced by the AOCGCMs daily 10-meter wind. In a case study, we show that for several models, this forcing is sufficient to reproduce the main SST error patterns seen after 1 month in the AOCGCM hindcasts.

Interdecadal enhancement of the Walker circulation over the Tropical Pacific in the late 1990s

The Walker circulation is one of the major components of the large-scale tropical atmospheric circulation and variations in its strength are critical to equatorial Pacific Ocean circulation. It has been argued in the literature that during the 20th century the Walker circulation weakened, and that this weakening was attributable to anthropogenic climate change. By using updated observations, we show that there has been a rapid interdecadal enhancement of the Walker circulation since the late 1990s. Associated with this enhancement is enhanced precipitation in the tropical western Pacific, anomalous westerlies in the upper troposphere, descent in the central and eastern tropical Pacific, and anomalous surface easterlies in the western and central tropical Pacific. The characteristics of associated oceanic changes are a strengthened thermocline slope and an enhanced zonal SST gradient across the tropical Pacific. Many characteristics of these changes are similar to those associated with the mid-1970s climate shift with an opposite sign. We also show that the interdecadal variability of the Walker circulation in the tropical Pacific is inversely correlated to the interdecadal variability of the zonal circulation in the tropical Atlantic. An enhancement of the Walker circulation in the tropical Pacific is associated with a weakening zonal circulation in the tropical Atlantic and vise versa, implying an inter-Atlantic-Pacific connection of the zonal overturning circulation variation. Whether these recent changes will be sustained is not yet clear, but our research highlights the importance of understanding the interdecadal variability, as well as the long-term trends, that influence tropical circulation.

Dynamical and observational constraints on tropical Pacific sea surface temperatures at the last glacial maximum

Asynchronously coupled atmosphere and ocean general circulation model simulations are used to examine the consequences of changes in the west/east sea-surface temperature (SST) gradient across the equatorial Pacific at the last glacial maximum (LGM). Simulations forced by the CLIMAP SST for the LGM, where the west/east SST gradient across the Pacific is reduced compared to present, produce a reduction in the strength of the trade winds and a decrease in the west/east slope of the equatorial thermocline that is incompatible with thermocline depths newly inferred from foraminiferal assemblages. Stronger-than-present trade winds, and a more realistic simulation of the thermocline slope, are produced when eastern Pacific SSTs are 2°C cooler than western Pacific SSTs. Our study highlights the importance of spatial heterogeneity in tropical SSTs in determining key features of the glacial climate.

Future changes in the Western North Pacific tropical cyclone activity projected by a multidecadal simulation with a 16-km global atmospheric GCM

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Niño–Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.

Interdecadal changes on the seasonal prediction of the western North Pacific summer climate around the late 1970s and early 1990s

Identifying predictability and the corresponding sources for the western North Pacific (WNP) summer climate in the case of non-stationary teleconnections during recent decades benefits for further improvements of long-range prediction on the WNP and East Asian summers. In the past few decades, pronounced increases on the summer sea surface temperature (SST) and associated interannual variability are observed over the tropical Indian Ocean and eastern Pacific around the late 1970s and over the Maritime Continent and western–central Pacific around the early 1990s. These increases are associated with significant enhancements of the interannual variability for the lower-tropospheric wind over the WNP. In this study, we further assess interdecadal changes on the seasonal prediction of the WNP summer anomalies, using May-start retrospective forecasts from the ENSEMBLES multi-model project in the period 1960–2005. It is found that prediction of the WNP summer anomalies exhibits an interdecadal shift with higher prediction skills since the late 1970s, particularly after the early 1990s. Improvements of the prediction skills for SSTs after the late 1970s are mainly found around tropical Indian Ocean and the WNP. The better prediction of the WNP after the late 1970s may arise mainly from the improvement of the SST prediction around the tropical eastern Indian Ocean. The close teleconnections between the tropical eastern Indian Ocean and WNP summer variability work both in the model predictions and observations. After the early 1990s, on the other hand, the improvements are detected mainly around the South China Sea and Philippines for the lower-tropospheric zonal wind and precipitation anomalies, associating with a better description of the SST anomalies around the Maritime Continent. A dipole SST pattern over the Maritime Continent and the central equatorial Pacific Ocean is closely related to the WNP summer anomalies after the early 1990s. This teleconnection mode is quite predictable, which is realistically reproduced by the models, presenting more predictable signals to the WNP summer climate after the early 1990s.

Influences of deep convective cloud systems on tropospheric trace gases and photochemistry over the tropical West Pacific

A numerical model for studying the influences of deep convective cloud systems on photochemistry was developed based on a non-hydrostatic meteorological model and chemistry from a global chemistry transport model. The transport of trace gases, the scavenging of soluble trace gases, and the influences of lightning produced nitrogen oxides (NOx=NO+NO2) on the local ozone-related photochemistry were investigated in a multi-day case study for an oceanic region located in the tropical western Pacific. Model runs considering influences of large scale flows, previously neglected in multi-day cloud resolving and single column model studies of tracer transport, yielded that the influence of the mesoscale subsidence (between clouds) on trace gas transport was considerably overestimated in these studies. The simulated vertical transport and scavenging of highly soluble tracers were found to depend on the initial profiles, reconciling contrasting results from two previous studies. Influences of the modeled uptake of trace gases by hydrometeors in the liquid and the ice phase were studied in some detail for a small number of atmospheric trace gases and novel aspects concerning the role of the retention coefficient (i.e. the fraction of a dissolved trace gas that is retained in the ice phase upon freezing) on the vertical transport of highly soluble gases were illuminated. Including lightning NOx production inside a 500 km 2-D model domain was found to be important for the NOx budget and caused small to moderate changes in the domain averaged ozone concentrations. A number of sensitivity studies yielded that the fraction of lightning associated NOx which was lost through photochemical reactions in the vicinity of the lightning source was considerable, but strongly depended on assumptions about the magnitude and the altitude of the lightning NOx source. In contrast to a suggestion from an earlier study, it was argued that the near zero upper tropospheric ozone mixing ratios which were observed close to the study region were most probably not caused by the formation of NO associated with lightning. Instead, it was argued in agreement with suggestions from other studies that the deep convective transport of ozone-poor air masses from the relatively unpolluted marine boundary layer, which have most likely been advected horizontally over relatively large distances (both before and after encountering deep convection) probably played a role. In particular, it was suggested that the ozone profiles observed during CEPEX (Central Equatorial Pacific Experiment) were strongly influenced by the deep convection and the larger scale flow which are associated with the intra-seasonal oscillation.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

Intensification of tropical Pacific biological productivity due to volcanic eruptions

Major volcanic eruptions generate widespread ocean cooling, which reduces upper ocean stratification. This effect has the potential to increase nutrient delivery into the euphotic zone and boost biological productivity. Using externally forced last millennium simulations of three climate/Earth System models (Model for Interdisciplinary Research On Climate (MIROC), Community Earth System Model (CESM), and LOch-Vecode-Ecbilt-CLio-agIsm Model (LOVECLIM)), we test the hypothesis that large volcanic eruptions intensify nutrient-driven export production. It is found that strong volcanic radiative forcing enhances the likelihood of eastern Pacific El Niño-like warming in CESM and LOVECLIM. This leads to an initial reduction of nutrients and export production in the eastern equatorial Pacific. However, this initial response reverses after about 3 years in association with La Niña cooling. The resulting delayed enhancement of biological production resembles the multiyear response in MIROC. The model simulations show that volcanic impacts on tropical Pacific dynamics and biogeochemistry persist for several years, thus providing a new source for potential multiyear ecosystem predictability.

(Table 1) Contents of tungsten and other metals in Fe-Mn nodules from the Pacific Ocean

Tungsten contents in iron-manganese nodules and crusts from different parts of the World Ocean, as well as its relationships with a number of chemical elements are under consideration. A trend to correlation of tungsten with Fe, Ti, W, Pb, and Co is noticed. Comparison of tungsten contents in the nodules and host sediments indicates its low geochemical mobility.

Dust flux records from the Subarctic North Pacific

We present a new record of eolian dust flux to the western Subarctic North Pacific (SNP) covering the past 27000 years based on a core from the Detroit Seamount. Comparing the SNP dust record to the NGRIP ice core record shows significant differences in the amplitude of dust changes to the two regions during the last deglaciation, while the timing of abrupt changes is synchronous. If dust deposition in the SNP faithfully records its mobilization in East Asian source regions, then the difference in the relative amplitude must reflect climate-related changes in atmospheric dust transport to Greenland. Based on the synchronicity in the timing of dust changes in the SNP and Greenland, we tie abrupt deglacial transitions in the 230Th-normalized 4He flux record to corresponding transitions in the well-dated NGRIP dust flux record to provide a new chronostratigraphic technique for marine sediments from the SNP. Results from this technique are complemented by radiocarbon dating, which allows us to independently constrain radiocarbon paleoreservoir ages. We find paleoreservoir ages of 745 ± 140 yr at 11653 yr BP, 680 ± 228 yr at 14630 yr BP and 790 ± 498 yr at 23290 yr BP. Our reconstructed paleoreservoir ages are consistent with modern surface water reservoir ages in the western SNP. Good temporal synchronicity between eolian dust records from the Subantarctic Atlantic and equatorial Pacific and the ice core record from Antarctica supports the reliability of the proposed dust tuning method to be used more widely in other global ocean regions.

Nd, Pb, and Be isotopes of ferromanaganese crusts of the South Pacific

The application of radiogenic isotopes to the study of Cenozoic circulation patterns in the South Pacific Ocean has been hampered by the fact that records from only equatorial Pacific deep water have been available. We present new Pb and Nd isotope time series for two ferromanganese crusts that grew from equatorial Pacific bottom water (D137-01, 'Nova', 7219 m water depth) and southwest Pacific deep water (63KD, 'Tasman', 1700 m water depth). The crusts were dated using 10Be/9Be ratios combined with constant Co-flux dating and yield time series for the past 38 and 23 Myr, respectively. The surface Nd and Pb isotope distributions are consistent with the present-day circulation pattern, and therefore the new records are considered suitable to reconstruct Eocene through Miocene paleoceanography for the South Pacific. The isotope time series of crusts Nova and Tasman suggest that equatorial Pacific deep water and waters from the Southern Ocean supplied the dissolved trace metals to both sites over the past 38 Myr. Changes in the isotopic composition of crust Nova are interpreted to reflect development of the Antarctic Circumpolar Current and changes in Pacific deep water circulation caused by the build up of the East Antarctic Ice Sheet. The Nd isotopic composition of the shallower water site in the southwest Pacific appears to have been more sensitive to circulation changes resulting from closure of the Indonesian seaway.

The Fohsella chronocline in middle Miocene sediments from the Ontong Java Plateau, Pacific Ocean

Micropaleontologists have traditionally recognized the mid-Miocene Fohsella lineage as a flagship for phyletic gradualism within the planktic foraminifera. However, study of a deep-sea record from the western equatorial Pacific (ODP Site 806) reveals that coiling ratios within this clade suddenly (<5 kyr) shift after a prolonged, ancestral state of near randomness (~50%) to a transient phase (13.42-13.43 Ma) of dextral dominance (~75%) immediately following the first common occurrence of keeled fohsellids. This brief period of dextral dominance was abruptly (<5 kyr) succeeded by an irreversible change to sinistral dominance (~96%). Fohsellid abundances decline markedly through the interval in which the sinistral preference is established. The shift to sinistrality (13.42 Ma) predated the deepening of fohsellid depth ecology by ~240-488 kyr, indicating that these two events were unrelated. This view is supported by a lack of delta 18O evidence for depth-habitat differences between the two chiral forms, which refutes the notion that sinistral fohsellids were "pre-adapted" for ensuing hydrographic change because they occupied a deeper depth habitat than their dextral counterparts. Planktic foraminiferal assemblages become strongly oligotrophic in character through the interval in which the fohsellid delta 18O increase is recorded, indicating that the migration to deeper depths was fostered by an expansion of the mixed layer in the western equatorial Pacific. Salient aspects of this brief, but conspicuous faunal change are a marked increase in the abundance of symbiont-bearing globigerinoidids, a concomitant collapse of local Jenkinsella mayeri/siakensis populations, and reduced fohsellid abundances. The rapid and permanent nature of the Fohsella sinistral shift provides a distinct, unequivocal datum that may prove useful for correlating mid-Miocene sections throughout the Caribbean Sea and tropical regions in the western sectors of the Pacific and Atlantic. The coiling ratio changes that occurred during the evolution of the Fohsella chronocline probably reflect changing population dynamics between cryptic genotypes with different coiling preferences.

Pliocene-Pleistocene oxygen isotope record DSDP Site 586, Ontong Java Plateau, Pacific Ocean i

Oceanographic changes in the western equatorial Pacific during the past 6 Ma are inferred from oxygen isotopic analyses of planktic and benthic foraminifera from Ontong Java Plateau (DSDP Site 586). The taxa are Globigerinoides sacculifer, Pulleniatina, Cibicidoides wuellerstorfi, and Oridorsalis umbonatus. Cooling and ice buildup are indicated by an 18O enrichment of 0.3 per mil in the planktic species near 3.4 Ma. This shift apparently is compensated in the benthic data by a warming of the deep waters by between 1° and 2° C. We suggest that the dominant source of upper deep water supply to the Pacific changed from Antarctic to North Atlantic at that time, the North Atlantic-derived water being warmer. Near 2.8 Ma (approximately) the planktic foraminifera again record an enrichment in 18O (Delta delta18O=0.25 per mil). We suggest ice buildup in the northern hemisphere as the cause, because of subsequent sharp increase in fluctuations of the delta18O signal, that is, instability. The enrichment is magnified in the benthic foraminifera (Delta delta18O = 0.5 per mil) by a cooling of the deep water by 1.5° at the time, presumably signalling a glacial-type reduction of North Atlantic Deep Water (NADW) production. Episodic divergence between the signals of G. sacculifer and Pulleniatina in the Pleistocene apparently reflects periods of increased upwelling in the western equatorial Pacific. The amplitude of ice volume fluctuations cannot be reconstructed from delta18O data alone, unless there are constraints on temperature variations. The increase in amplitude of fluctuation of the benthic and planktic signals during the Pleistocene may be attributed either to an increase in maximum ice volume, or to an increase in the fractionation of continental ice, or a combination of both causes.

Part of the global DOC versus AOU (dissolved organic carbon/apparent oxygen utilization) data compilation, Pacific Ocean

Total organic carbon (TOC) was analyzed on four transects along 140°W in 1992 using a high temperature combustion/discrete injection (HTC/DI) analyzer. For two of the transects, the analyses were conducted on-board ship. Mixed-layer concentrations of organic carbon varied from about 80 µM C at either end of the transect (12°N and 12°S) to about 60 µM C at the equator. Total organic carbon concentrations decreased rapidly below the mixed-layer to about 38-40 µM C at 1000 m across the transect. Little variation was observed below this depth; deep water concentrations below 2000 m were virtually monotonic at about 36 µM C. Repeat measurements made on subsequent cruises consistently found the same concentrations at 1000 m or deeper, but substantial variations were observed in the mixed-layer and the upper water column above 400 m depth. Linear mixing models of total organic carbon versus sigmaT exhibited zones of organic carbon formation and consumption. TOC was found to be inversely correlated with apparent oxygen utilization (AOU) in the region between the mixed-layer and the oxygen minimum. In the mixed-layer, TOC concentrations varied seasonally. Part of the variations in TOC at the equator was driven by changes in the upwelling rate in response to variations in physical forcing related to an El Niño and to the passage of tropical instability waves. TOC export fluxes, calculated from simple box models, averaged 8±4 mmol C/m**2/day at the equator and also varied seasonally. These export fluxes account for 50-75% of the total carbon deficit and are consistent with other estimates and model predictions.