An observational and theoretical study of the structure and propagation of the Madden-Julian Oscillation

Thesis (Ph.D.)--University of Washington, 2016-06

Longitudinal localization of tropical intraseasonal variability

Intraseasonal time-scales play an important role in tropical variability. Two modes that contribute significantly to tropical intraseasonal variability (ISV) are the eastward-propagating MaddenJulian Oscillation (MJO), and westward-moving moist equatorial Rossby waves. This note reports on a correspondence between the longitudinal gradient of mean tropical precipitable water (PW), and the geographical regions of genesis, and convective activity, of both these large-scale tropical systems. Our finding is based on an analysis of PW from the MERRA reanalysis product. The data indicate that the mean tropical PW has a dominant wavenumber two (three) structure in longitude in the Northern (Southern) Hemisphere. Departures from a longitudinally homogeneous state are attributed to the influence of subtropical anticyclones, and are accentuated by monsoonal regions of both hemispheres. This mean structure results in a sharply localized longitudinal gradient of PW. Remarkably, regions with positive gradients (such as the Northern and Southern Hemisphere western Indian Ocean), i.e. they have larger PW to the east, are the very zones that are implicated in the formation, and show high levels of convective activity, of the eastward-moving MJO. On the other hand, regions with negative gradients (such as the Southern Hemisphere central Pacific) are the very regions where genesis, and maxima in variance, of westward-moving moist equatorial Rossby waves are known to occur. Apart from providing a first-order longitudinal footprint of the convective phase of these systems, this correspondence reinforces the role of the mean climatic state in tropical ISV. Copyright (c) 2012 Royal Meteorological Society

Dynamics of intraseasonal sea level and thermocline variability in the equatorial Atlantic during 2002-03

Satellite and in situ observations in the equatorial Atlantic Ocean during 2002-03 show dominant spectral peaks at 40-60 days and secondary peaks at 10-40 days in sea level and thermocline within the intraseasonal period band (10-80 days). A detailed investigation of the dynamics of the intraseasonal variations is carried out using an ocean general circulation model, namely, the Hybrid Coordinate Ocean Model (HYCOM). Two parallel experiments are performed in the tropical Atlantic Ocean basin for the period 2000-03: one is forced by daily scatterometer winds from the Quick Scatterometer (QuikSCAT) satellite together with other forcing fields, and the other is forced by the low-passed 80-day version of the above fields. To help in understanding the role played by the wind-driven equatorial waves, a linear continuously stratified ocean model is also used. Within 3 degrees S-3 degrees N of the equatorial region, the strong 40-60-day sea surface height anomaly (SSHA) and thermocline variability result mainly from the first and second baroclinic modes equatorial Kelvin waves that are forced by intraseasonal zonal winds, with the second baroclinic mode playing a more important role. Sharp 40-50-day peaks of zonal and meridional winds appear in both the QuikSCAT and Pilot Research Moored Array in the Tropical Atlantic (PIRATA) data for the period 2002-03, and they are especially strong in 2002. Zonal wind anomaly in the central-western equatorial basin for the period 2000-06 is significantly correlated with SSHA across the equatorial basin, with simultaneous/ lag correlation ranging from-0.62 to 0.74 above 95% significance. Away from the equator (3 degrees-5 degrees N), however, sea level and thermocline variations in the 40-60-day band are caused largely by tropical instability waves (TIWs). On 10-40-day time scales and west of 10 degrees W, the spectral power of sea level and thermocline appears to be dominated by TIWs within 5 degrees S-5 degrees N of the equatorial region. The wind-driven circulation, however, also provides a significant contribution. Interestingly, east of 10 W, SSHA and thermocline variations at 10 40- day periods result almost entirely from wind-driven equatorial waves. During the boreal spring of 2002 when TIWs are weak, Kelvin waves dominate the SSHA across the equatorial basin (2 degrees S-2 degrees N). The observed quasi-biweekly Yanai waves are excited mainly by the quasi-biweekly meridional winds, and they contribute significantly to the SSHA and thermocline variations in 1 degrees-5 degrees N and 1 degrees-5 degrees S regions.

An extended variable-grid global ocean circulation model and its preliminary results of the equatorial Pacific circulation

To investigate the interaction between the tropical Pacific and China seas a variable-grid global ocean circulation model with fine grid[(1/6)degrees] covering the area from 20degreesS to 50degreesN and from 99degrees to 150degreesE is developed. Numerical computation of the annually cyclic circulation fields is performed. The results of the annual mean zonal currents and deep to abyssal western boundary currents in the equatorial Pacific Ocean are reported. The North Equatorial Current,the North Equatorial Countercurrent, the South Equatorial Current and the Equatorial Undercurrent are fairly well simulated. The model well reproduces the northward flowing abyssal western boundary current. From the model results a lower deep western boundary current east of the Bismarck-Solomon-New Hebrides Island chain at depths around 2 000 in has been found. The model results also show that the currents in the equatorial Pacific Ocean have multi-layer structures both in zonal currents and western boundary currents, indicating that the global ocean overturning thermohaline circulation appears of multi-layer pattern.

POST-MONSOON SEA SURFACE TEMPERATURE AND CONVECTION ANOMALIES OVER INDIAN AND PACIFIC OCEANS

We have studied sea surface temperature (SST) anomalies over the Indian and Pacific Oceans (domain 25 °S to 25°N and 40 °E to 160 °W) during the three seasons following the Indian summer monsoon for wet monsoons and also for dry monsoons accompanied or not by El Ni˜no. A dry monsoon is followed by positive SST anomalies in the longitude belt 40 to 120 °E, negative anomalies in 120 to 160 °E and again positive anomalies east of 160 °E. In dry monsoons accompanied by El Ni˜no the anomalies have the same sign, but are much stronger. Wet monsoons have weak anomalies of opposite sign in all three of the longitude belts. Thus El Ni˜no and a dry monsoon have the same types of association with the Indian and Pacific Ocean SSTs. In the sector 40 to 120 °E SST anomalies first appear over the western part of the Indian Ocean (June to September) followed by the same sign of anomalies over its eastern part and China Sea (October to March). By March after a dry monsoon or El Ni˜no the Indian Ocean between 10 °N and 10 °S has a spatially large warm SST anomaly. Anomalies in deep convection tend to follow the SST anomalies, with warm SST anomalies producing positive convection anomalies around the seasonal location of the intertropical convergence zone

The ENSO-Asian monsoon interaction in a coupled ocean-atmosphere GCM

In this study, the authors evaluate the (El Niño–Southern Oscillation) ENSO–Asian monsoon interaction in a version of the Hadley Centre coupled ocean–atmosphere general circulation model (CGCM) known as HadCM3. The main focus is on two evolving anomalous anticyclones: one located over the south Indian Ocean (SIO) and the other over the western North Pacific (WNP). These two anomalous anticyclones are closely related to the developing and decaying phases of the ENSO and play a crucial role in linking the Asian monsoon to ENSO. It is found that the HadCM3 can well simulate the main features of the evolution of both anomalous anticyclones and the related SST dipoles, in association with the different phases of the ENSO cycle. By using the simulated results, the authors examine the relationship between the WNP/SIO anomalous anticyclones and the ENSO cycle, in particular the biennial component of the relationship. It is found that a strong El Niño event tends to be followed by a more rapid decay and is much more likely to become a La Niña event in the subsequent winter. The twin anomalous anticyclones in the western Pacific in the summer of a decaying El Niño are crucial for the transition from an El Niño into a La Niña. The El Niño (La Niña) events, especially the strong ones, strengthen significantly the correspondence between the SIO anticyclonic (cyclonic) anomaly in the preceding autumn and WNP anticyclonic (cyclonic) anomaly in the subsequent spring, and favor the persistence of the WNP anomaly from spring to summer. The present results suggest that both El Niño (La Niña) and the SIO/WNP anticyclonic (cyclonic) anomalies are closely tied with the tropospheric biennial oscillation (TBO). In addition, variability in the East Asian summer monsoon, which is dominated by the internal atmospheric variability, seems to be responsible for the appearance of the WNP anticyclonic anomaly through an upper-tropospheric meridional teleconnection pattern over the western and central Pacific.

The role of air–sea coupling in the simulation of the Madden–Julian oscillation in the Hadley Centre model

The role of air–sea coupling in the simulation of the Madden–Julian oscillation (MJO) is explored using two configurations of the Hadley Centre atmospheric model (AGCM), GA3.0, which differ only in F, a parameter controlling convective entrainment and detrainment. Increasing F considerably improves deficient MJO-like variability in the Indian and Pacific Oceans, but variability in and propagation through the Maritime Continent remains weak. By coupling GA3.0 in the tropical Indo-Pacific to a boundary-layer ocean model, KPP, and employing climatological temperature corrections, well resolved air–sea interactions are simulated with limited alterations to the mean state. At default F, when GA3.0 has a poor MJO, coupling produces a stronger MJO with some eastward propagation, although both aspects remain deficient. These results agree with previous sensitivity studies using AGCMs with poor variability. At higher F, coupling does not affect MJO amplitude but enhances propagation through the Maritime Continent, resulting in an MJO that resembles observations. A sensitivity experiment with coupling in only the Indian Ocean reverses these improvements, suggesting coupling in the Maritime Continent and West Pacific is critical for propagation. We hypothesise that for AGCMs with a poor MJO, coupling provides a “crutch” to artificially augment MJO-like activity through high-frequency SST anomalies. In related experiments, we employ the KPP framework to analyse the impact of air–sea interactions in the fully coupled GA3.0, which at default F shows a similar MJO to uncoupled GA3.0. This is due to compensating effects: an improvement from coupling and a degradation from mean-state errors. Future studies on the role of coupling should carefully separate these effects.

A multi-model analysis of the role of the ocean on the African and Indian monsoon during mid-Holocene

We investigate the role of the ocean feedback on the climate in response to insolation forcing during the mid-Holocene (6,000 year BP) using results from seven coupled ocean–atmosphere general circulation models. We examine how the dipole in late summer sea-surface temperature (SST) anomalies in the tropical Atlantic increases the length of the African monsoon, how this dipole structure is created and maintained, and how the late summer SST warming in the northwest Indian Ocean affects the monsoon retreat in this sector. Similar mechanisms are found in all of the models, including a strong wind evaporation feedback and changes in the mixed layer depth that enhance the insolation forcing, as well as increased Ekman transport in the Atlantic that sharpens the Atlantic dipole pattern. We also consider changes in interannual variability over West Africa and the Indian Ocean. The teleconnection between variations in SST and Sahelian precipitation favor a larger impact of the Atlantic dipole mode in this region. In the Indian Ocean, the strengthening of the Indian dipole structure in autumn has a damping effect on the Indian dipole mode at the interannual time scale

Propagation of the Madden-Julian Oscillation through the Maritime Continent and scale interaction with the diurnal cycle of precipitation

The convectively active part of the Madden-Julian Oscillation (MJO) propagates eastward through the warm pool, from the Indian Ocean through the Maritime Continent (the Indonesian archipelago) to the western Pacific. The Maritime Continent's complex topography means the exact nature of the MJO propagation through this region is unclear. Model simulations of the MJO are often poor over the region, leading to local errors in latent heat release and global errors in medium-range weather prediction and climate simulation. Using 14 northern winters of TRMM satellite data it is shown that, where the mean diurnal cycle of precipitation is strong, 80% of the MJO precipitation signal in the Maritime Continent is accounted for by changes in the amplitude of the diurnal cycle. Additionally, the relationship between outgoing long-wave radiation (OLR) and precipitation is weakened here, such that OLR is no longer a reliable proxy for precipitation. The canonical view of the MJO as the smooth eastward propagation of a large-scale precipitation envelope also breaks down over the islands of the Maritime Continent. Instead, a vanguard of precipitation (anomalies of 2.5 mm day^-1 over 10^6 km^2) jumps ahead of the main body by approximately 6 days or 2000 km. Hence, there can be enhanced precipitation over Sumatra, Borneo or New Guinea when the large-scale MJO envelope over the surrounding ocean is one of suppressed precipitation. This behaviour can be accommodated into existing MJO theories. Frictional and topographic moisture convergence and relatively clear skies ahead of the main convective envelope combine with the low thermal inertia of the islands, to allow a rapid response in the diurnal cycle which rectifies onto the lower-frequency MJO. Hence, accurate representations of the diurnal cycle and its scale interaction appear to be necessary for models to simulate the MJO successfully.

An assessment of Indian monsoon seasonal forecasts and mechanisms underlying monsoon interannual variability in the Met Office GloSea5-GC2 system

We assess Indian summer monsoon seasonal forecasts in GloSea5-GC2, the Met Office fully coupled subseasonal to seasonal ensemble forecasting system. Using several metrics, GloSea5-GC2 shows similar skill to other state-of-the-art forecast systems. The prediction skill of the large-scale South Asian monsoon circulation is higher than that of Indian monsoon rainfall. Using multiple linear regression analysis we evaluate relationships between Indian monsoon rainfall and five possible drivers of monsoon interannual variability. Over the time period studied (1992-2011), the El Nino-Southern Oscillation (ENSO) and the Indian Ocean dipole (IOD) are the most important of these drivers in both observations and GloSea5-GC2. Our analysis indicates that ENSO and its teleconnection with the Indian rainfall are well represented in GloSea5-GC2. However, the relationship between the IOD and Indian rainfall anomalies is too weak in GloSea5-GC2, which may be limiting the prediction skill of the local monsoon circulation and Indian rainfall. We show that this weak relationship likely results from a coupled mean state bias that limits the impact of anomalous wind forcing on SST variability, resulting in erroneous IOD SST anomalies. Known difficulties in representing convective precipitation over India may also play a role. Since Indian rainfall responds weakly to the IOD, it responds more consistently to ENSO than in observations. Our assessment identifies specific coupled biases that are likely limiting GloSea5-GC2 prediction skill, providing targets for model improvement.

Impacts of interruption of the Agulhas leakage on the tropical Atlantic in coupled ocean-atmosphere simulations

In this paper we use a coupled ocean-atmosphere model to investigate the impact of the interruption of Agulhas leakage of Indian ocean water on the tropical Atlantic, a region where strong coupled ocean-atmosphere interactions occur. The effect of a shut down of leakage of Indian ocean water is isolated from the effect of a collapse of the MOC. In our experiments, the ocean model is forced with boundary conditions in the southeastern corner of the domain that correspond to no interocean exchange of Indian ocean water into the Atlantic. The southern boundary condition is taken from the Levitus data and ensures an MOC in the Atlantic. Within this configuration, instead of warm and salty Indian ocean water temperature (cold) and salinity (fresh) anomalies of southern ocean origin propagate into the South Atlantic and eventually reach the equatorial region, mainly in the thermocline. This set up mimics the closure of the ""warm water path"" in favor of the ""cold water path"". As part of the atmospheric response, there is a northward shift of the intertropical convergence zone (ITCZ). The changes in trade winds lead to reduced Ekman pumping in the equatorial region. This leads to a freshening and warming of the surface waters along the equator. Especially in the Cold Tongue region, the cold and fresh subsurface anomalies do not reach the surface due to the reduced upwelling. The anomaly signals are transported by the equatorial undercurrent and spread away from the equator within the thermocline. Part of the anomaly eventually reaches the Tropical North Atlantic, where it affects the Guinea Dome. Surprisingly, the main effect at the surface is small on the equator and relatively large at the Guinea Dome. In the atmosphere, the northward shift of the ITCZ is associated with a band of negative precipitation anomalies and higher salinities over the Tropical South Atlantic. An important implication of these results is that the modified water characteristics due to a shut down of the Agulhas leakage remain largely unaffected when crossing the equatorial Atlantic and therefore can affect the deepwater formation in the North Atlantic. This supports the hypothesis that the Agulhas leakage is an important source region for climate change and decadal variability of the Atlantic.

Ocean dynamics, not dust, have controlled equatorial Pacific productivity over the past 500,000 years

Biological productivity in the equatorial Pacific is relatively high compared with other low-latitude regimes, especially east of the dateline, where divergence driven by the trade winds brings nutrient rich waters of the Equatorial Undercurrent to the surface. The equatorial Pacific is one of the three principal high-nutrient low chlorophyll ocean regimes where biological utilization of nitrate and phosphate is limited, in part, by the availability of iron. Throughout most of the equatorial Pacific, upwelling of water from the Equatorial Undercurrent supplies far more dissolved iron than is delivered by dust, by as much as two orders of magnitude. Nevertheless, recent studies have inferred that the greater supply of dust during ice ages stimulated greater utilization of nutrients within the region of upwelling on the equator, thereby contributing to the sequestration of carbon in the ocean interior. Here we present proxy records for dust and for biological productivity over the past 500 ky at three sites spanning the breadth of the equatorial Pacific Ocean to test the dust fertilization hypothesis. Dust supply peaked under glacial conditions, consistent with previous studies, whereas proxies of export production exhibit maxima during ice age terminations. Temporal decoupling between dust supply and biological productivity indicates that other factors, likely involving ocean dynamics, played a greater role than dust in regulating equatorial Pacific productivity.

Mg/Ca ratios of Cenozoic benthic foraminifera from the Indian and South Atlantic Ocean

We investigate the evolution of Cenozoic climate and ice volume as evidenced by the oxygen isotopic composition of seawater (delta18Osw) derived from benthic foraminiferal Mg/Ca ratios to constrain the temperature effect contained in foraminiferal delta18O values. We have constructed two benthic foraminiferal Mg/Ca records from intermediate water depth sites (Ocean Drilling Program sites 757 and 689 from the subtropical Indian Ocean and the Weddell Sea, respectively). Together with the previously published composite record of Lear et al. (2002, doi:10.1126/science.287.5451.269) and the Neogene record from the Southern Ocean of Billups and Schrag (2002, doi:10.1029/2000PA000567), we obtain three, almost complete representations of the delta18Osw for the past 52 Myr. We discuss the sensitivity of early Cenozoic Mg/Ca-derived paleotemperatures (and hence the delta18Osw) to assumptions about seawater Mg/Ca ratios. We find that during the middle Eocene (~ 49-40 Ma), modern seawater ratios yield Mg/Ca-derived temperatures that are in good agreement with the oxygen isotope paleothermometer assuming ice-free conditions. Intermediate waters cooled during the middle Eocene reaching minimum temperatures by 40 Ma. The corresponding delta18Osw reconstructions support ice growth on Antarctica beginning by at least 40 Ma. At the Eocene/Oligocene boundary, Mg/Ca ratios (and hence temperatures) from Weddell Sea site 689 display a well-defined maximum. We caution against a paleoclimatic significance of this result and put forth that the partitioning coefficient of Mg in benthic foraminifera may be sensitive to factors other than temperature. Throughout the remainder of the Cenozoic, the temporal variability among delta18Osw records is similar and similar to longer-term trends in the benthic foraminiferal delta18O record. An exception occurs during the Pliocene when delta18Osw minima in two of the three records suggest reductions in global ice volume that are not apparent in foraminiferal delta18O records, which provides a new perspective to the ongoing debate about the stability of the Antarctic ice sheet. Maximum delta18Osw values recorded during the Pleistocene at Southern Ocean site 747 agree well with values derived from the geochemistry of pore waters (Schrag et al., 1996, doi:10.1126/science.272.5270.1930) further highlighting the value of the new Mg/Ca calibrations of Martin et al. (2002, doi:10.1016/S0012-821X(02)00472-7) and Lear et al. (2002, doi:10.1016/S0016-7037(02)00941-9) applied in this study. We conclude that the application of foraminiferal Mg/Ca ratios allows a refined view of Cenozoic ice volume history despite uncertainties related to the geochemical cycling of Mg and Ca on long time scales.

Faunal analysis and Mg/Ca records of planktic foraminifera in sediment core SO189/2_039KL

Climate phenomena like the monsoon system, El Niño Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) are interconnected via various feedback mechanisms and control the climate of the Indian Ocean and its surrounding continents on various timescales. The eastern tropical Indian Ocean is a key area for the interplay of these phenomena and for reconstructing their past changes and forcing mechanisms. Here we present records of upper ocean thermal gradient, thermocline temperatures (TT) and relative abundances of planktic foraminifera in core SO 189-39KL taken off western Sumatra (0°47.400' S, 99°54.510' E) for the last 8 ka that we use as proxies for changes in upper ocean structure. The records suggest a deeper thermocline between 8 ka and ca 3 ka compared to the late Holocene. We find a shoaling of the thermocline after 3 ka, most likely indicating an increased occurrence of upwelling during the late Holocene compared to the mid-Holocene which might represent changes in the IOD-like mean state of the Indian Ocean with a more negative IOD-like mean state during the mid-Holocene and a more positive IOD-like mean state during the past 3 ka. This interpretation is supported by a transient Holocene climate model simulation in which an IOD-like mode is identified that involves an insolation-forced long-term trend of increasing anomalous surface easterlies over the equatorial eastern Indian Ocean.

Age determination and sea surface temperature reconstruction for the Southern Ocean

The isotopic and micropaleontological deglacial records of three deep-sea cores from 44°S to 55°S have been dated by accelerator mass spectrometry. The available records did not allow accurate dating of the initiation of the deglaciation. By 13,000 years B.P., sea surface temperatures reached values similar to the present values. A cool oscillation abruptly interrupted this warm phase between 12,000 and 11,000 years B.P. Initiation of this cooling therefore preceded the northern hemisphere Younger Dryas by approximately 1000 years. Complete warming was reached by 10,000 years B.P., more or less synchronous with the northeast Atlantic Ocean.