Data assimilation using ensemble transform Kalman filter (ETKF) in ROMS model for Indian Ocean

Study of Oceans dynamics and forecast is crucial as it influences the regional climate and other marine activities. Forecasting oceanographic states like sea surface currents, Sea surface temperature (SST) and mixed layer depth at different time scales is extremely important for these activities. These forecasts are generated by various ocean general circulation models (OGCM). One such model is the Regional Ocean Modelling System (ROMS). Though ROMS can simulate several features of ocean, it cannot reproduce the thermocline of the ocean properly. Solution to this problem is to incorporates data assimilation (DA) in the model. DA system using Ensemble Transform Kalman Filter (ETKF) has been developed for ROMS model to improve the accuracy of the model forecast. To assimilate data temperature and salinity from ARGO data has been used as observation. Assimilated temperature and salinity without localization shows oscillations compared to the model run without assimilation for India Ocean. Same was also found for u and v-velocity fields. With localization we found that the state variables are diverging within the localization scale.

A note on new indices for the equatorial Indian Ocean oscillation

It is now well known that there is a strong association of the extremes of the Indian summer monsoon rainfall (ISMR) with the El Nio and southern oscillation (ENSO) and the Equatorial Indian Ocean Oscillation (EQUINOO), later being an east-west oscillation in convection anomaly over the equatorial Indian Ocean. So far, the index used for EQUINOO is EQWIN, which is based on the surface zonal wind over the central equatorial Indian Ocean. Since the most important attribute of EQUINOO is the oscillation in convection/precipitation, we believe that the indices based on convection or precipitation would be more appropriate. Continuous and reliable data on outgoing longwave radiation (OLR), and satellite derived precipitation are now available from 1979 onwards. Hence, in this paper, we introduce new indices for EQUINOO, based on the difference in the anomaly of OLR/precipitation between eastern and western parts of the equatorial Indian Ocean. We show that the strong association of extremes of the Indian summer monsoon with ENSO and EQUINOO is also seen when the new indices are used to represent EQUINOO.

Narrowband signal detection techniques in shallow ocean by acoustic vector sensor array

This paper presents the formulation and performance analysis of four techniques for detection of a narrowband acoustic source in a shallow range-independent ocean using an acoustic vector sensor (AVS) array. The array signal vector is not known due to the unknown location of the source. Hence all detectors are based on a generalized likelihood ratio test (GLRT) which involves estimation of the array signal vector. One non-parametric and three parametric (model-based) signal estimators are presented. It is shown that there is a strong correlation between the detector performance and the mean-square signal estimation error. Theoretical expressions for probability of false alarm and probability of detection are derived for all the detectors, and the theoretical predictions are compared with simulation results. It is shown that the detection performance of an AVS array with a certain number of sensors is equal to or slightly better than that of a conventional acoustic pressure sensor array with thrice as many sensors.

Data acquisition and processing at ocean bottom for a Tsunami warning system

The design and development of a Bottom Pressure Recorder for a Tsunami Early Warning System is described here. The special requirements that it should satisfy for the specific application of deployment at ocean bed and pressure monitoring of the water column above are dealt with. A high-resolution data digitization and low circuit power consumption are typical ones. The implementation details of the data sensing and acquisition part to meet these are also brought out. The data processing part typically encompasses a Tsunami detection algorithm that should detect an event of significance in the background of a variety of periodic and aperiodic noise signals. Such an algorithm and its simulation are presented. Further, the results of sea trials carried out on the system off the Chennai coast are presented. The high quality and fidelity of the data prove that the system design is robust despite its low cost and with suitable augmentations, is ready for a full-fledged deployment at ocean bed. (C) 2013 Elsevier Ltd. All rights reserved.

Development of a regional model for the North Indian Ocean

We have developed a one-way nested Indian Ocean regional model. The model combines the National Oceanic and Atmospheric Administration (NOAA) Geophysical Fluid Dynamics Laboratory's (GFDL) Modular Ocean Model (MOM4p1) at global climate model resolution (nominally one degree), and a regional Indian Ocean MOM4p1 configuration with 25 km horizontal resolution and 1 m vertical resolution near the surface. Inter-annual global simulations with Coordinated Ocean-Ice Reference Experiments (CORE-II) surface forcing over years 1992-2005 provide surface boundary conditions. We show that relative to the global simulation, (i) biases in upper ocean temperature, salinity and mixed layer depth are reduced, (ii) sea surface height and upper ocean circulation are closer to observations, and (iii) improvements in model simulation can be attributed to refined resolution, more realistic topography and inclusion of seasonal river runoff. Notably, the surface salinity bias is reduced to less than 0.1 psu over the Bay of Bengal using relatively weak restoring to observations, and the model simulates the strong, shallow halocline often observed in the North Bay of Bengal. There is marked improvement in subsurface salinity and temperature, as well as mixed layer depth in the Bay of Bengal. Major seasonal signatures in observed sea surface height anomaly in the tropical Indian Ocean, including the coastal waveguide around the Indian peninsula, are simulated with great fidelity. The use of realistic topography and seasonal river runoff brings the three dimensional structure of the East India Coastal Current and West India Coastal Current much closer to observations. As a result, the incursion of low salinity Bay of Bengal water into the southeastern Arabian Sea is more realistic. (C) 2013 Elsevier Ltd. All rights reserved.

Evaluation of GODAS Using RAMA Mooring Observations from the Indian Ocean

We present a comparison of the Global Ocean Data Assimilation System (GODAS) five-day ocean analyses against in situ daily data from Research Moored Array for African-Asian-Australian Monsoon Analysis and Prediction (RAMA) moorings at locations 90 degrees E, 12 degrees N; 90 degrees E, 8 degrees N; 90 degrees E, 0 degrees N and 90 degrees E, 1.5 degrees S in the equatorial Indian Ocean and the Bay of Bengal during 2002-2008. We find that the GODAS temperature analysis does not adequately capture a prominent signal of Indian Ocean dipole mode of 2006 seen in the mooring data, particularly at 90 degrees E 0 degrees N and 90 degrees E 1.5 degrees S in the eastern India Ocean. The analysis, using simple statistics such as bias and root-mean-square deviation, indicates that standard GODAS temperature has definite biases and significant differences with observations on both subseasonal and seasonal scales. Subsurface salinity has serious deficiencies as well, but this may not be surprising considering the poorly constrained fresh water forcing, and possible model deficiencies in subsurface vertical mixing. GODAS reanalysis needs improvement to make it more useful for study of climate variability and for creating ocean initial conditions for prediction.

The April 2012 Indian Ocean earthquakes: Seismotectonic context and implications for their mechanisms

The 11 April 2012 earthquakes (M-w 8.6 and M-w 8.2) were sourced within the Northern Wharton Basin in the northeastern part of the Indo-Australian diffuse plate boundary. This unusually active oceanic intraplate region has generated many large earthquakes in the past, most of which are believed to have occurred by strike-slip motion, triggered by the NW-SE oriented compressional stresses acting across the Indian and Australian plates. In the aftermath of the 2004 megathrust earthquake along the nearby Sunda Trench, increased seismicity in the Northern Wharton Basin is attributed to the stress transfer from the Sumatra-Andaman plate boundary. Models proposed for the April 2012 earthquakes differ somewhat in details but partly attribute their complex rupture to the reactivation of pre-existing structures. These structures include previously mapped N-S trending fracture zones within the Northern Wharton Basin and E-W lineations across the Ninetyeast Ridge. In this paper, we review the regional tectonics and past seismicity on the Indo-Australian Plate in order to understand the seismotectonic setting of the April 2012 Indian Ocean earthquakes. (c) 2014 Elsevier B.V. All rights reserved.

Slow Slip Acceleration beneath Andaman Islands Triggered by the 11 April 2012 Indian Ocean Earthquakes

The M-w 8.6 and 8.2 strike-slip earthquakes that struck the northeast Indian Ocean on 11 April 2012 resulted in coseismic deformation both at near and distant sites. The slip distribution, deduced using seismic-wave analysis for the orthogonal faults that ruptured during these earthquakes, is sufficient to predict the coseismic displacements at the Global Positioning System (GPS) sites, such as NTUS, PALK, and CUSV, but fall short at four continuous sites in the Andaman Islands region. Slip modeling, for times prior to the events, suggests that the lower portion of the thrust fault beneath the Andaman Islands has been slipping at least at the rate of 40 cm/yr, in response to the 2004 Sumatra-Andaman coseismic stress change. Modeling of GPS displacements suggests that the en echelon and orthogonal fault ruptures of the 2012 intraplate oceanic earthquakes could have possibly accelerated the ongoing slow slip, along the lower portion of the thrust fault beneath the islands with a month-long slip of 4-10 cm. The misfit to the coseismic GPS displacements along the Andaman Islands could be improved with a better source model, assuming that no local process contributed to this anomaly.

Three-dimensional localization of multiple acoustic sources in shallow ocean with non-Gaussian noise

In this paper, a low-complexity algorithm SAGE-USL is presented for 3-dimensional (3-D) localization of multiple acoustic sources in a shallow ocean with non-Gaussian ambient noise, using a vertical and a horizontal linear array of sensors. In the proposed method, noise is modeled as a Gaussian mixture. Initial estimates of the unknown parameters (source coordinates, signal waveforms and noise parameters) are obtained by known/conventional methods, and a generalized expectation maximization algorithm is used to update the initial estimates iteratively. Simulation results indicate that convergence is reached in a small number of (<= 10) iterations. Initialization requires one 2-D search and one 1-D search, and the iterative updates require a sequence of 1-D searches. Therefore the computational complexity of the SAGE-USL algorithm is lower than that of conventional techniques such as 3-D MUSIC by several orders of magnitude. We also derive the Cramer-Rao Bound (CRB) for 3-D localization of multiple sources in a range-independent ocean. Simulation results are presented to show that the root-mean-square localization errors of SAGE-USL are close to the corresponding CRBs and significantly lower than those of 3-D MUSIC. (C) 2014 Elsevier Inc. All rights reserved.

Sheltered coastal environments as archives of paleo-tsunami deposits: Observations from the 2004 Indian Ocean tsunami

The 2004 earthquake left several traces of coseismic land deformation and tsunami deposits, both on the islands along the plate boundary and distant shores of the Indian Ocean rim countries. Researchers are now exploring these sites to develop a chronology of past events. Where the coastal regions are also inundated by storm surges, there is an additional challenge to discriminate between the deposits formed by these two processes. Paleo-tsunami research relies largely on finding deposits where preservation potential is high and storm surge origin can be excluded. During the past decade of our work along the Andaman and Nicobar Islands and the east coast of India, we have observed that the 2004 tsunami deposits are best preserved in lagoons, inland streams and also on elevated terraces. Chronological evidence for older events obtained from such sites is better correlated with those from Thailand, Sri Lanka and Indonesia, reiterating their usefulness in tsunami geology studies. (C) 2014 Elsevier Ltd. All rights reserved.

Impact of river runoff into the ocean on Indian summer monsoon

Rivers of the world discharge about 36000 km 3 of freshwater into the ocean every year. To investigate the impact of river discharge on climate, we have carried out two 100 year simulations using the Community Climate System Model (CCSM3), one including the river runoff into the ocean and the other excluding it. When the river discharge is shut off, global average sea surface temperature (SST) rises by about 0.5 degrees C and the Indian Summer Monsoon Rainfall (ISMR) increases by about 10% of the seasonal total with large increase in the eastern Bay of Bengal and along the west coast of India. In addition, the frequency of occurrence of La Nina-like cooling events in the equatorial Pacific increases and the correlation between ISMR and Pacific SST anomalies become stronger. The teleconnection between the SST anomalies in the Pacific and monsoon is effected via upper tropospheric meridional temperature gradient and the North African-Asian Jet axis.

Stable isotopic signature of Southern Ocean deep water CO2 ventilation

The link between atmospheric CO2 level and ventilation state of the deep ocean is poorly understood due to the lack of coherent observations on the partitioning of carbon between atmosphere and ocean. In this Southern Ocean study, we have classified the Southern Ocean into different zones based on its hydrological features and have binned the variability in latitudinal air-CO2 concentration and its isotopic ratios. Together with air-CO2, we analysed the surface water for the isotopic ratios in dissolved inorganic carbon (DIC). Using the binary mixing approach on the isotopic ratio of atmospheric CO2 and its concentration, we identified the delta C-13 value of source CO2. The isotopic composition of source CO2 was around -9.22 +/- 0.26 parts per thousand for the year 2011 and 2012, while a composition of -13.49 +/- 4.07 parts per thousand was registered for the year 2013. We used the delta C-13 of DIC to predict the CO2 composition in air under equilibrium and compared our estimates with actual observations. We suggest that the degeneration of the DIC in presence of warm water in the region was the factor responsible for adding the CO2 to the atmosphere above. The place of observation coincides with the zone of high wind speed which promotes the process of CO2 exsolution from sea water. (C) 2015 Elsevier Ltd. All rights reserved.

Effect of Surface Energy Anisotropy on the Stability of Growth Fronts in Multiphase Alloys

Eutectic growth offers a variety of examples for pattern formation which are interesting both for theoreticians as well as experimentalists. One such example of patterns is ternary eutectic colonies which arise as a result of instabilities during growth of two solid phases. Here, in addition to the two major components being exchanged between the solid phases during eutectic growth, there is an impurity component which is rejected by both solid phases. During progress of solidification, there develops a boundary layer of the third impurity component ahead of the solidification front of the two solid phases. Similar to Mullins-Sekerka type instabilities, such a boundary layer tends to make the global solidification envelope unstable to morphological perturbations giving rise to two-phase cells. This phenomenon has been studied numerically in two dimensions for the conditions of directional solidification, by Plapp and Karma (Phys Rev E 66:061608, 2002) using phase-field simulations. While, in the work by Plapp and Karma (Phys Rev E 66:061608, 2002) all interfaces are isotropic, in our presentation, we extend the phase-field model by considering interfacial anisotropy in the solid-solid and solid-liquid interfaces and characterize the role of interfacial anisotropy on the stability of the growth front through phase-field simulations in two dimensions.

Prediction of Indian rainfall during the summer monsoon season on the basis of links with equatorial Pacific and Indian Ocean climate indices

Interannual variation of Indian summer monsoon rainfall (ISMR) is linked to El Nino-Southern oscillation (ENSO) as well as the Equatorial Indian Ocean oscillation (EQUINOO) with the link with the seasonal value of the ENSO index being stronger than that with the EQUINOO index. We show that the variation of a composite index determined through bivariate analysis, explains 54% of ISMR variance, suggesting a strong dependence of the skill of monsoon prediction on the skill of prediction of ENSO and EQUINOO. We explored the possibility of prediction of the Indian rainfall during the summer monsoon season on the basis of prior values of the indices. We find that such predictions are possible for July-September rainfall on the basis of June indices and for August-September rainfall based on the July indices. This will be a useful input for second and later stage forecasts made after the commencement of the monsoon season.

Isotopic disequilibrium in Globigerina bulloides and carbon isotope response to productivity increase in Southern Ocean

Oxygen and carbon isotope ratios in planktonic foraminifera Globigerina bulloides collected from tow samples along a transect from the equatorial Indian ocean to the Southern Ocean (45 degrees E and 80 degrees E and 10 degrees N to 53 degrees S) were analysed and compared with the equilibrium delta O-18 and delta C-13 values of calcite calculated using the temperature and isotopic composition of the water column. The results agree within similar to 0.25% for the region between 10 degrees N and 40 degrees S and 75-200 m water depth which is considered to be the habitat of Globigerina bulloides. Further south (from 40 degrees S to 55 degrees S), however, the measured delta O-18 and delta C-13 values are higher than the expected values by similar to 2% and similar to 1% respectively. These enrichments can be attributed to either a `vital effect' or a higher calcification rate. An interesting pattern of increase in the delta C-13(DIC) value of the surface water with latitude is observed between 35 degrees S and similar to 60 degrees S, with a peak at similar to 42 degrees S. This can be caused by increased organic matter production and associated removal. A simple model accounting for the increase in the delta C-13(DIC) values is proposed which fits well with the observed chlorophyll abundance as a function of latitude.