Joint modeling of lithosphere and mantle dynamics: Evaluation of constraints from global tomography models

With the advances in technology, seismological theory, and data acquisition, a number of high-resolution seismic tomography models have been published. However, discrepancies between tomography models often arise from different theoretical treatments of seismic wave propagation, different inversion strategies, and different data sets. Using a fixed velocity-to-density scaling and a fixed radial viscosity profile, we compute global mantle flow models associated with the different tomography models and test the impact of these for explaining surface geophysical observations (geoid, dynamic topography, stress, and strain rates). We use the joint modeling of lithosphere and mantle dynamics approach of Ghosh and Holt (2012) to compute the full lithosphere stresses, except that we use HC for the mantle circulation model, which accounts for the primary flow-coupling features associated with density-driven mantle flow. Our results show that the seismic tomography models of S40RTS and SAW642AN provide a better match with surface observables on a global scale than other models tested. Both of these tomography models have important similarities, including upwellings located in Pacific, Eastern Africa, Iceland, and mid-ocean ridges in the Atlantic and Indian Ocean and downwelling flows mainly located beneath the Andes, the Middle East, and central and Southeast Asia.

On the roles of the northeast cold surge, the Borneo vortex, the Madden-Julian Oscillation, and the Indian Ocean Dipole during the extreme 2006/2007 flood in southern Peninsular Malaysia

The mid-December 2006 to late January 2007 flood in southern Peninsular Malaysia was the worst flood in a century and was caused by three extreme precipitation episodes. These extreme precipitation events were mainly associated with strong northeasterly winds over the South China Sea. In all cases, the northeasterlies penetrated anomalously far south and followed almost a straight trajectory. The elevated terrain over Sumatra and southern Peninsular Malaysia caused low-level convergence. The strong easterly winds near Java associated with the Rossby wave-type response to Madden-Julian Oscillation (MJO) inhibited the counter-clockwise turning of the northeasterlies and the formation of the Borneo vortex, which, in turn, enhanced the low-level convergence over the region. The abrupt termination of the Indian Ocean Dipole (IOD) in December 2006 played a secondary role as warmer equatorial Indian Ocean helped in the MJO formation.

Atmospheric boundary-layer across Hadley and Ferrel cells over the Indian Ocean

A Southern Ocean Pilot cruise covering the latitudes from 10 degrees N to 56 degrees S in the open Indian Ocean was carried out during January February 2004. Surface and upper air data collected during this cruise are reported here. It is shown that the broad features of the atmosphere, in particular that of temperature, follow the tropical and mid-latitude weather expected during January February in this region. However, the atmospheric boundary-layer shows large variations, both in its height and structure between tropics and high latitudes. Strong influence of the surface heat flux on boundary layer structure is clearly seen. Humidity field reveals several local maxima and minima, suggesting a laminated atmosphere with air from different sources moving almost unmixed in adjacent layers.

Dynamics of the Indian-Ocean oxygen minimum zones

In the Indian Ocean, mid-depth oxygen minimum zones (OMZs) occur in the Arabian Sea and the Bay of Bengal. The lower part of the Arabian-Sea OMZ (ASOMZ; below 400 m) intensifies northward across the basin; in contrast, its upper part (above 400 m) is located in the central/eastern basin, well east of the most productive regions along the western boundary. The Bay-of-Bengal OMZ (BBOMZ), although strong, is weaker than the ASOMZ. To investigate the processes that maintain the Indian-Ocean OMZs, we obtain a suite of solutions to a coupled biological/physical model. Its physical component is a variable-density, 6 1/2-layer model, in which each layer corresponds to a distinct dynamical regime or water-mass type. Its biological component has six compartments: nutrients, phytoplankton, zooplankton, two size classes of detritus, and oxygen. Because the model grid is non-eddy resolving (0.5 degrees), the biological model also includes a parameterization of enhanced mixing based on the eddy kinetic energy derived from satellite observations. To explore further the impact of local processes on OMZs, we also obtain analytic solutions to a one-dimensional, simplified version of the biological model. Our control run is able to simulate basic features of the oxygen, nutrient, and phytoplankton fields throughout the Indian Ocean. The model OMZs result from a balance, or lack thereof, between a sink of oxygen by remineralization and subsurface oxygen sources due primarily to northward spreading of oxygenated water from the Southern Hemisphere, with a contribution from Persian-Gulf water in the northern Arabian Sea. The northward intensification of the lower ASOMZ results mostly from horizontal mixing since advection is weak in its depth range. The eastward shift of the upper ASOMZ is due primarily to enhanced advection and vertical eddy mixing in the western Arabian Sea, which spread oxygenated waters both horizontally and vertically. Advection carries small detritus from the western boundary into the central/eastern Arabian Sea, where it provides an additional source of remineralization that drives the ASOMZ to suboxic levels. The model BBOMZ is weaker than the ASOMZ because the Bay lacks a remote source of detritus from the western boundary. Although detritus has a prominent annual cycle, the model OMZs do not because there is not enough time for significant remineralization to occur.

The 2014 M-w 6.1 Bay of Bengal, Indian Ocean, Earthquake: A Possible Association with the 85 degrees E Ridge

Oceanic intraplate earthquakes are known to occur either on active ridge-transform structures or by reactivation of their inactive counterparts, generally referred to as fossil ridges or transforms. The Indian Ocean, one of the most active oceanic intraplate regions, has generated large earthquakes associated with both these types of structures. The moderate earthquake that occurred on 21 May 2014 (M-w 6.1) in the northern Bay of Bengal followed an alternate mechanism, as it showed no clear association either with active or extinct ridge-transform structures. Its focal depth of >50 km is uncommon but not improbable, given the similar to 90 Ma age of the ocean floor with 12-km-thick overlying sediments. No tectonic features have been mapped in the near vicinity of its epicenter, the closest being the 85 degrees E ridge, located similar to 100 km to its west, hitherto regarded as seismically inactive. The few earthquakes that have occurred here in the past are clustered around its southern or northern limits, and a few are located midway, at around 10 degrees N. The 2014 earthquake, sourced close to the northern cluster, seems to be associated with a northwest-southeast-oriented fracture, located on the eastern flanks of the 85 degrees E ridge. If this causal association is possible, we believe that reactivation of fossil hotspot trails could be considered as another mechanism for oceanic intraplate seismicity.

Ocean atmosphere coupling over monsoon regions

In monsoon regions, the seasonal migration of the intertropical convergence zone (ITCZ) is manifested as a seasonal reversal of winds. Most of the summer monsoon rainfall over India occurs owing to synoptic and large-scale convection associated with the continental ITCZ (Fig. 1). We have investigated the interaction between these large-scale convective systems and the ocean over which they are generated1â3, concentrating on the relationship between organized convection over the Indian Ocean and sea surface temperature (SST). We report here that on a monthly basis the degree of cloudiness correlates well with SST for the relatively colder oceans, but when SST is maintained above 28 °C it ceases to be an important factor in determining the variability of cloudiness. Over the major regions of convection east of 70°E, which are warm year after year, the observed cloudiness cannot be correlated with variations in SST.

New species of Ptychozoon (Sauria: Gekkonidae) from the Nicobar Archipelago, Indian Ocean

A new species of Ptychozoon is described from the central portion of the Nicobar Archipelago, Bay of Bengal, India. It has been formerly referred to P. kuhli, a species widely distributed in Sundaland. Ptychozoon nicobarensis sp. nov. reaches an SVL of 100.3 mm, and is diagnosable from congeneric species in showing the following combination of characters: dorsum with a tan vertebral stripe, lacking dark transverse bars; supranasals in contact; cutaneous expansions on sides of head; absence of predigital notch in preantebrachial cutaneous expansion; imbricate parachute support scales; four irregular rows of low, rounded enlarged scales on dorsum; 20-29 scales across widest portion of tail terminus; three indistinct chevrons on dorsum; 7-11 pairs of preanal pores; femoral pores absent; tail with an expanded terminal flap and weak lobe fusion at proximal border of tail terminus. The curious distribution of the new species, centred around the central Nicobars is speculated to be the result of competition with and/or predation by large gekkonid species, to the north (Gekko verreauxi) and south (G. smithii) of the group of islands occupied by the new Ptychozoon from the central Nicobars.

Normal-mode sound propagation in an ocean with sinusoidal surface waves

The normal-mode solution to the problem of acoustic wave propagation in an isovelocity ocean with a wavy surface is considered. The surface wave amplitude is assumed to be small compared to the acoustic wavelength, and the method of multiple scales is employed to study the interaction between normal-mode acoustic waves and the surface waves. A nonresonant interaction causes small fluctuations of the amplitude and phase of the acoustic wave at a rate dependent on the frequency of the surface wave. Backscatter occurs if the wavenumber of the surface wave is larger than that of the acoustic wave. The interaction becomes resonant if appropriate phase-matching conditions are satisfied. In this case, two acoustic normal modes get coupled, resulting in a large-scale periodic exchange of energy from one mode to another.

The possible influence of curvature and rotation on ocean currents

Recent laboratory investigations have shown that rotation and (streamwise) curvature can have spectacular effects on momentum transport in turbulent shear flows. A simple model that takes account of these effects (based on an analogy with buoyant flows) utilises counterparts of the Richardson number Rg and the Monin-Oboukhov length. Estimates of Rg for meanders in ocean currents like the Gulf Stream show it to be of order 1 or more, while laboratory investigations reveal strong effects even at |Rg|∼0·1. These considerations lead to the conclusion that at a cyclonic bend in the Gulf Stream, a highly unstable flow in the outer half of the jet rides over a highly stable flow in the inner half. It is conjectured that the discrepancies noticed between observation and the various theories of Gulf Stream meanders, and such phenomena as the observed detachment of eddies from the Gulf Stream, may be due to the effects of curvature and rotation on turbulent transport.

Indian Ocean sea surface salinity variations in a coupled model

The variability of the sea surface salinity (SSS) in the Indian Ocean is studied using a 100-year control simulation of the Community Climate System Model (CCSM 2.0). The monsoon-driven seasonal SSS pattern in the Indian Ocean, marked by low salinity in the east and high salinity in the west, is captured by the model. The model overestimates runoff int the Bay of Bengal due to higher rainfall over the Himalayan-Tibetan regions which drain into the Bay of Bengal through Ganga-Brahmaputra rivers. The outflow of low-salinity water from the Bay of Bengal is to strong in the model. Consequently, the model Indian Ocean SSS is about 1 less than that seen in the climatology. The seasonal Indian Ocean salt balance obtained from the model is consistent with the analysis from climatological data sets. During summer, the large freshwater input into the Bay of Bengal and its redistribution decide the spatial pattern of salinity tendency. During winter, horizontal advection is the dominant contributor to the tendency term. The interannual variability of the SSS in the Indian Ocean is about five times larger than that in coupled model simulations of the North Atlantic Ocean. Regions of large interannual standard deviations are located near river mouths in the Bay of Bengal and in the eastern equatorial Indian Ocean. Both freshwater input into the ocean and advection of this anomalous flux are responsible for the generation of these anomalies. The model simulates 20 significant Indian Ocean Dipole (IOD) events and during IOD years large salinity anomalies appear in the equatorial Indian Ocean. The anomalies exist as two zonal bands: negative salinity anomalies to the north of the equator and positive to the south. The SSS anomalies for the years in which IOD is not present and for ENSO years are much weaker than during IOD years. Significant interannual SSS anomalies appear in the Indian Ocean only during IOD years.

Tree-ring based hydroclimate reconstruction over northern Vietnam from Fokienia hodginsii: eighteenth century mega-drought and tropical Pacific influence

We present here the first statistically calibrated and verified tree-ring reconstruction of climate from continental Southeast Asia.The reconstructed variable is March-May (MAM) Palmer Drought Severity Index (PDSI) based on ring widths from 22 trees (42 radial cores) of rare and long-lived conifer, Fokienia hodginsii (Po Mu as locally called) from northern Vietnam. This is the first published tree ring chronology from Vietnam as well as the first for this species. Spanning 535 years, this is the longest cross-dated tree-ring series yet produced from continental Southeast Asia. Response analysis revealed that the annual growth of Fokienia at this site was mostly governed by soil moisture in the pre-monsoon season. The reconstruction passed the calibration-verification tests commonly used in dendroclimatology, and revealed two prominent periods of drought in the mid-eighteenth and late-nineteenth enturies. The former lasted nearly 30 years and was concurrent with a similar drought over northwestern Thailand inferred from teak rings, suggesting a ``mega-drought'' extending across Indochina in the eighteenth century. Both of our reconstructed droughts are consistent with the periods of warm sea surface temperature (SST)anomalies in the tropical Pacific. Spatial correlation analyses with global SST indicated that ENSO-like anomalies might play a role in modulating droughts over the region, with El Nio (warm) phases resulting in reduced rainfall. However, significant correlation was also seen with SST over the Indian Ocean and the north Pacific,suggesting that ENSO is not the only factor affecting the climate of the area. Spectral analyses revealed significant peaks in the range of 53.9-78.8 years as well as in the ENSO-variability range of 2.0 to 3.2 years.