An exploration framework to identify and track movement of cloud systems

We describe a framework to explore and visualize the movement of cloud systems. Using techniques from computational topology and computer vision, our framework allows the user to study this movement at various scales in space and time. Such movements could have large temporal and spatial scales such as the Madden Julian Oscillation (MJO), which has a spatial scale ranging from 1000 km to 10000 km and time of oscillation of around 40 days. Embedded within these larger scale oscillations are a hierarchy of cloud clusters which could have smaller spatial and temporal scales such as the Nakazawa cloud clusters. These smaller cloud clusters, while being part of the equatorial MJO, sometimes move at speeds different from the larger scale and in a direction opposite to that of the MJO envelope. Hitherto, one could only speculate about such movements by selectively analysing data and a priori knowledge of such systems. Our framework automatically delineates such cloud clusters and does not depend on the prior experience of the user to define cloud clusters. Analysis using our framework also shows that most tropical systems such as cyclones also contain multi-scale interactions between clouds and cloud systems. We show the effectiveness of our framework to track organized cloud system during one such rainfall event which happened at Mumbai, India in July 2005 and for cyclone Aila which occurred in Bay of Bengal during May 2009.

Elastic thickness structure of the Andaman subduction zone: Implications for convergence of the Ninetyeast Ridge

We use the Bouguer coherence (Morlet isostatic response function) technique to compute the spatial variation of effective elastic thickness (T-e) of the Andaman subduction zone. The recovered T-e map resolves regional-scale features that correlate well with known surface structures of the subducting Indian plate and the overriding Burma plate. The major structure on the India plate, the Ninetyeast Ridge (NER), exhibits a weak mechanical strength, which is consistent with the expected signature of an oceanic ridge of hotspot origin. However, a markedly low strength (0< T-e <3 km) in that region, where the NER is close to the Andaman trench (north of 10 N), receives our main attention in this study. The subduction geometry derived from the Bouguer gravity forward modeling suggests that the NER has indented beneath the Andaman arc. We infer that the bending stresses of the viscous plate, which were reinforced within the subducting oceanic plate as a result of the partial subduction of the NER buoyant load, have reduced the lithospheric strength. The correlation, T-e < T-s (seismogenic thickness) reveals that the upper crust is actively deforming beneath the frontal arc Andaman region. The occurrence of normal-fault earthquakes in the frontal arc, low Te zone, is indicative of structural heterogeneities within the subducting plate. The fact that the NER along with its buoyant root is subducting under the Andaman region is inhibiting the subduction processes, as suggested by the changes in trench line, interrupted back-arc volcanism, variation in seismicity mechanism, slow subduction, etc. The low T-e and thinned crustal structure of the Andaman back-arc basin are attributed to a thermomechanically weakened lithosphere. The present study reveals that the ongoing back-arc spreading and strike-slip motion along the West Andaman Fault coupled with the ridge subduction exerts an important control on the frequency and magnitude of seismicity in the Andaman region. (C) 2013 Elsevier Ltd. All rights reserved.

Nitrogen deposition: how important is it for global terrestrial carbon uptake?

Global carbon budget studies indicate that the terrestrial ecosystems have remained a large sink for carbon despite widespread deforestation activities. CO2 fertilization, N deposition and re-growth of mid-latitude forests are believed to be key drivers for land carbon uptake. In this study, we assess the importance of N deposition by performing idealized near-equilibrium simulations using the Community Land Model 4.0 (CLM4). In our equilibrium simulations, only 12-17% of the deposited nitrogen is assimilated into the ecosystem and the corresponding carbon uptake can be inferred from a C : N ratio of 20 : 1. We calculate the sensitivity of the terrestrial biosphere for CO2 fertilization, climate warming and N deposition as changes in total ecosystem carbon for unit changes in global mean atmospheric CO2 concentration, global mean temperature and Tera grams of nitrogen deposition per year, respectively. Based on these sensitivities, it is estimated that about 242 PgC could have been taken up by land due to the CO2 fertilization effect and an additional 175 PgC taken up as a result of the increased N deposition since the pre-industrial period. Because of climate warming, the terrestrial ecosystem could have lost about 152 PgC during the same period. Therefore, since pre-industrial times terrestrial carbon losses due to warming may have been more or less compensated by effects of increased N deposition, whereas the effect of CO2 fertilization is approximately indicative of the current increase in terrestrial carbon stock. Our simulations also suggest that the sensitivity of carbon storage to increased N deposition decreases beyond current levels, indicating that climate warming effects on carbon storage may overwhelm N deposition effects in the future.

On the nature of meandering of the springtime western boundary current in the Bay of Bengal

We present evidence that the springtime western boundary current (WBC) in the Bay of Bengal is a continuous northward-flowing current from about 12 degrees N to 17 degrees N, which then separates from the coast at around 18 degrees N. We first revisit a hydrographic data set collected in 1987 from a potential vorticity perspective, and then analyze absolute dynamic height maps from satellite altimeters during the period 2000-2010. The altimetric maps suggest that the mean configuration of the WBC is that of an intense current with two anticyclonic eddies on the offshore side, which are part of the basin-wide anticyclonic circulation. The WBC consistently separates from the coast at around 18 degrees N in all years between 2000 and 2010. The path of the eastward-flowing mean stream after separation appears to be consistent with isolines of f/H and with Ertel's potential vorticity, based on an analysis of the hydrographic data from 1987.

Spatial variations of effective elastic thickness over the Ninetyeast Ridge and implications for its structure and tectonic evolution

We present new data on the strength of oceanic lithosphere along the Ninetyeast Ridge (NER) from two independent methods: spectral analysis (Bouguer coherence) using the fan wavelet transform technique, and spatial analysis (flexure inversion) with the convolution method. The two methods provide effective elastic thickness (T-e) patterns that broadly complement each other, and correlate well with known surface structures and regional-scale features. Furthermore, our study presents a new high resolution database on the Moho configuration, which obeys flexural isostasy, and exhibit regional correlations with the T-e variations. A continuous ridge structure with a much lower T-e value than that of normal oceanic lithosphere provides strong support for the hotspot theory. The derived T-e values vary over the northern (higher T-e similar to 10-20 km), central (anomalously low T-e similar to 0-5 km), and southern (low T-e similar to 5 km) segments of the NER. The lack of correlation of the T-e value with the progressive aging of the lithosphere implies differences in thermo-mechanical setting of the crust and underlying mantle in different parts of the NER, again indicating diversity in their evolution. The anomalously low T-e and deeper Moho (similar to 22 km) estimates of the central NER (between 0.5 degrees N and 17 degrees S) are attributed to the interaction of a hotspot with the Wharton spreading ridge that caused significant thermal rejuvenation and hence weakening of the lithosphere. The higher mechanical strength values in the northern NER (north of 0.5 degrees N) may support the idea of off-ridge emplacement and a relatively large plate motion at the time of volcanism. The low T-e and deeper Moho (similar to 22 km) estimates in the southern part (south of 17 degrees S) suggest that the lithosphere was weak and therefore younger at the time of volcanism, and this supports the idea that the southern NER was emplaced on the edge of the Indian plate. (C) 2013 Elsevier B.V. All rights reserved.

Missing (in-situ) snow cover data hampers climate change and runoff studies in the Greater Himalayas

The Himalayas are presently holding the largest ice masses outside the polar regions and thus (temporarily) store important freshwater resources. In contrast to the contemplation of glaciers, the role of runoff from snow cover has received comparably little attention in the past, although (i) its contribution is thought to be at least equally or even more important than that of ice melt in many Himalayan catchments and (ii) climate change is expected to have widespread and significant consequences on snowmelt runoff. Here, we show that change assessment of snowmelt runoff and its timing is not as straightforward as often postulated, mainly as larger partial pressure of H2O, CO2, CH4, and other greenhouse gases might increase net long-wave input for snowmelt quite significantly in a future atmosphere. In addition, changes in the short-wave energy balance such as the pollution of the snow cover through black carbon or the sensible or latent heat contribution to snowmelt are likely to alter future snowmelt and runoff characteristics as well. For the assessment of snow cover extent and depletion, but also for its monitoring over the extremely large areas of the Himalayas, remote sensing has been used in the past and is likely to become even more important in the future. However, for the calibration and validation of remotely-sensed data, and even-more so in light of possible changes in snow-cover energy balance, we strongly call for more in-situ measurements across the Himalayas, in particular for daily data on new snow and snow cover water equivalent, or the respective energy balance components. Moreover, data should be made accessible to the scientific community, so that the latter can more accurately estimate climate change impacts on Himalayan snow cover and possible consequences thereof on runoff. (C) 2013 Elsevier B.V. All rights reserved.

Trends in aerosol optical depth over Indian region: Potential causes and impact indicators

The first regional synthesis of long-term (back to similar to 25 years at some stations) primary data (from direct measurement) on aerosol optical depth from the ARFINET (network of aerosol observatories established under the Aerosol Radiative Forcing over India (ARFI) project of Indian Space Research Organization over Indian subcontinent) have revealed a statistically significant increasing trend with a significant seasonal variability. Examining the current values of turbidity coefficients with those reported similar to 50 years ago reveals the phenomenal nature of the increase in aerosol loading. Seasonally, the rate of increase is consistently high during the dry months (December to March) over the entire region whereas the trends are rather inconsistent and weak during the premonsoon (April to May) and summer monsoon period (June to September). The trends in the spectral variation of aerosol optical depth (AOD) reveal the significance of anthropogenic activities on the increasing trend in AOD. Examining these with climate variables such as seasonal and regional rainfall, it is seen that the dry season depicts a decreasing trend in the total number of rainy days over the Indian region. The insignificant trend in AOD observed over the Indo-Gangetic Plain, a regional hot spot of aerosols, during the premonsoon and summer monsoon season is mainly attributed to the competing effects of dust transport and wet removal of aerosols by the monsoon rain. Contributions of different aerosol chemical species to the total dust, simulated using Goddard Chemistry Aerosol Radiation and Transport model over the ARFINET stations, showed an increasing trend for all the anthropogenic components and a decreasing trend for dust, consistent with the inference deduced from trend in Angstrom exponent.

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.

Glacial mass balance changes in the Karakoram and Himalaya based on CMIP5 multi-model climate projections

The impact of future climate change on the glaciers in the Karakoram and Himalaya (KH) is investigated using CMIP5 multi-model temperature and precipitation projections, and a relationship between glacial accumulation-area ratio and mass balance developed for the region based on the last 30 to 40 years of observational data. We estimate that the current glacial mass balance (year 2000) for the entire KH region is -6.6 +/- 1 Gta(-1), which decreases about sixfold to -35 +/- 2 Gta(-1) by the 2080s under the high emission scenario of RCP8.5. However, under the low emission scenario of RCP2.6 the glacial mass loss only doubles to -12 +/- 2 Gta(-1) by the 2080s. We also find that 10.6 and 27 % of the glaciers could face `eventual disappearance' by the end of the century under RCP2.6 and RCP8.5 respectively, underscoring the threat to water resources under high emission scenarios.

Nonlinear gravity-wave interactions in stratified turbulence

To investigate the dynamics of gravity waves in stratified Boussinesq flows, a model is derived that consists of all three-gravity-wave-mode interactions (the GGG model), excluding interactions involving the vortical mode. The GGG model is a natural extension of weak turbulence theory that accounts for exact three-gravity-wave resonances. The model is examined numerically by means of random, large-scale, high-frequency forcing. An immediate observation is a robust growth of the so-called vertically sheared horizontal flow (VSHF). In addition, there is a forward transfer of energy and equilibration of the nonzero-frequency (sometimes called ``fast'') gravity-wave modes. These results show that gravity-wave-mode interactions by themselves are capable of systematic interscale energy transfer in a stratified fluid. Comparing numerical simulations of the GGG model and the full Boussinesq system, for the range of Froude numbers (Fr) considered (0.05 a parts per thousand currency sign Fr a parts per thousand currency sign 1), in both systems the VSHF is hardest to resolve. When adequately resolved, VSHF growth is more vigorous in the GGG model. Furthermore, a VSHF is observed to form in milder stratification scenarios in the GGG model than the full Boussinesq system. Finally, fully three-dimensional nonzero-frequency gravity-wave modes equilibrate in both systems and their scaling with vertical wavenumber follows similar power-laws. The slopes of the power-laws obtained depend on Fr and approach -2 (from above) at Fr = 0.05, which is the strongest stratification that can be properly resolved with our computational resources.

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.

Neoarchean continental growth through arc magmatism in the Nilgiri Block, southern India

The formation and growth of continental crust in the Archean have been evaluated through models of subduction-accretion and mantle plume. The Nilgiri Block in southern India exposes exhumed Neoarchean lower crust, uplifted to heights of 2500 m above sea level along the north western margin of the Peninsula. Major lithologies in this block include charnockite with or without garnet, anorthosite-gabbro suite, pyroxenite, amphibolite and hornblende-biotite gneiss (TTG). All these rock types are closely associated as an arc magmatic suite, with diffuse boundaries and coeval nature. The charnockite and hornblende-biotite gneisses (TTG) show SiO2 content varying from 64 to 73 wt.%. The hornblende-biotite gneisses (TTG) are high-Al type with Al2O3 >15 wt.% whereas the charnockites show Al2O3 <15 wt.%. The composition of charnockite is mainly magnesian and calcic to calc-alkaline. The mafic-ultramafic rocks show composition close to that of tholeiitic series. The low values of K(2)o (<3 wt.%), (K/Rb)/K2O (<500), Zr/Ti, and trace element ratios like (La/Yb)n/(Sr/Y), (Y/Nb), (Y + Nb)/Rb, (Y+Ta)/Rb, Yb/Ta indicate a volcanic arc signature for these rocks. The geochemical signature is consistent with arc magmatic rocks generated through oceanic plate subduction. The primitive mantle normalized trace element patterns of these rocks display enrichment in large ion lithophile elements (LILE) and comparable high field strength elements (HFSE) in charnockite and hornblende-biotite gneisses (TTG) consistent with subduction-related origin. Primitive mantle normalized REE pattern displays an enrichment in LREE in the chamockite and hornblende-biotite gneisses (TTG) as compared to a flat pattern for the mafic rocks. The chondrite normalized REE patterns of zircons of all the rock types reveal cores with high HREE formed at ca. 2700 Ma and rims with low HREE formed at 2500-2450 Ma. Log-transformed La/Th-Nb/Th-Sm/Th-Yb/Th discrimination diagram for the mafic and ultramafic rocks from Nilgiri displays a transition from mid-oceanic ridge basalt (MORB) to island arc basalt (IAB) suggesting a MORB source. The U-Pb zircon data from the charnockites, mafic granulites and hornblende-biotite gneisses (TTG) presented in our study show that the magma generation during subduction and accretion events in this block occurred at 2700-2500 Ma. Together with the recent report on Neoarchean supra-subduction zone ophiolite suite at its southern margin, the Nilgiri Block provides one of the best examples for continental growth through vertical stacking and lateral accretion in a subduction environment during the Neoarchean. (c) 2014 Elsevier B.V. All rights reserved.

Role of orography in modulating space-time scales of convection over South Asia

Propagation of convective systems in the meridional direction during boreal summer is responsible for active and break phases of monsoon over south Asia. This region is unique in the world in its characteristics of monsoon variability and is in close proximity of mountains like the Himalayas. Here, using an atmospheric general circulation model, we try to understand the role of orography in determining spatial and temporal scales of these convective systems. Absence of orography (noGlOrog) decreased the simulated seasonal mean precipitation over India by 23 % due to delay in onset by about a month vis-a-vis the full-mountain case. In noGlOrog, poleward propagations were absent during the delayed period prior to onset. Post-onset, both simulations had similar patterns of poleward propagations. The spatial and temporal scales of propagating clouds bands were determined using wavelet analysis. These scales were found to be different in full-mountain and no-mountain experiments in June-July. However, after the onset of monsoon in noGlOrog, these scales become similar to that with orography. Simulations with two different sets of convection schemes confirmed this result. Further analysis shows that the absence (presence) of meridional propagations during early (late) phase of summer monsoon in noGlOrog was associated with weaker (stronger) vertical shear of zonal wind over south Asia. Our study shows that orography plays a major role in determining the time of onset over the Indian region. However, after onset, basic characteristics of propagating convective systems and therefore the monthly precipitation over India, are less sensitive to the presence of orography and are modulated by moist convective processes.