Seasonal variability of rainfall recorded in growth bands of the Giant African Land Snail Lissachatina fulica (Bowdich) from India

Seasonal rainfall patterns in Bangalore, India, have been reconstructed using stable isotopic ratios in the growth bands of Giant African Land Snail shells. The present study was conducted at Bangalore, India which receives rain during the summer months. The oxygen isotopic record in the rainwater samples collected during different months covering the period of the summer monsoon of the year 2008 is compared with the isotopic ratio in the gastropod growth bands deposited simultaneously. The chronology of the shell growth band is independently established assuming the growth rate observed in a chamber experiment maintaining similar relative humidity and temperature conditions. A consistent pattern observed in the isotopic ratio in the gastropod growth bands and rainwater is demonstrated and provides a novel approach for precipitation reconstruction at seasonal and weekly time scales. This approach of using isotopic ratios in the gastropod growth bands for rainfall can serve as a substitute for filling gaps in rainfall data and for cases where no rain records are available. In addition, they can be used to determine the frequencies and magnitudes of dry spells from the past records. (C) 2013 Elsevier B.V. All rights reserved.

Effect of SST Variation on ITCZ in APE Simulations

The effect of meridional variation of sea surface temperature (SST) on tropical atmospheric circulation is analyzed using Aqua-planet Experiment (APE) simulations. The meridional SST gradient around the narrow SST peak in CONTROL simulation favours a strong and single equatorial Intertropical Convergence Zone (ITCZ, defined by the maximum of zonally averaged total precipitation) in all APE models. In contrast, flat equatorial SST peak (FLAT simulation) favours split/double ITCZs flanking the SST maximum, in the majority of the APE models. Although there is reasonable agreement for SST sensitivity of ITCZ among the APE models in CONTROL, there exists disparity among them in FLAT case. Similarly, while the total and convective precipitation responses are consistent among the models, the large-scale precipitation response shows considerable inter-model variations in FLAT case. The APE intercomparison indicates that the occurrence and positioning of the ITCZ are primarily related to boundary layer moisture convergence as a response to the meridional variation of SST. Furthermore, the meridional gradient of tropospheric temperature is found to be an important factor that can influence the positioning of ITCZ. FLAT SST distribution is found to be similar to the observed distribution over the Indian region during summer season. Models that yield double ITCZs in this case simulate an easterly jet over the equatorial region (similar to 15 degrees equatorward of the ITCZ). This is analogous to the Tropical Easterly Jet (TEJ), which is a unique feature observed over the Indian region during summer monsoon season, with its core at 12 degrees N, equatorward of the seasonal convergence zone centered along 25 degrees N. In these models, positive meridional temperature gradient and the associated easterly shear in the atmosphere strengthened by moisture convergence penetrate up to the upper troposphere, with which TEJ is in thermal wind balance.

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.

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.

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.

Breeding in bamboo: a novel anuran reproductive strategy discovered in Rhacophorid frogs of the Western Ghats, India

Amphibians exhibit extraordinarily diverse sets of reproductive strategies among vertebrates. Understanding life history strategies in an evolutionary framework is lacking for many amphibian species in the tropics. Here, we report a novel reproductive mode where adult frogs enter hollow internodes of bamboo via a small opening, deposit direct developing eggs, and provide parental care. This behaviour is observed in two species of the frog genus Raorchestes. The first description of this unique life history and details of nest site characteristics and embryo development are provided along with ecological comparisons. Evolution of novel reproductive modes and parental care are discussed in context of natural selection. Dearth of natural history information on amphibians in the Western Ghats and much of the South-East Asian region is highlighted with suggestions for further studies.(c) 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 114, 1-11.

Impact of diurnal forcing on intraseasonal sea surface temperature oscillations in the Bay of Bengal

The diurnal cycle is an important mode of sea surface temperature (SST) variability in tropical oceans, influencing air-sea interaction and climate variability. Upper ocean mixing mechanisms are significant at diurnal time scales controlling the intraseasonal variability (ISV) of SST. Sensitivity experiments using an Ocean General Circulation Model (OGCM) for the summer monsoon of the year 2007 show that incorporation of diurnal cycle in the model atmospheric forcings improves the SST simulation at both intraseasonal and shorter time scales in the Bay of Bengal (BoB). The increase in SST-ISV amplitudes with diurnal forcing is approximate to 0.05 degrees C in the southern bay while it is approximate to 0.02 degrees C in the northern bay. Increased intraseasonal warming with diurnal forcing results from the increase in mixed layer heat gain from insolation, due to shoaling of the daytime mixed layer. Amplified intraseasonal cooling is dominantly controlled by the strengthening of subsurface processes owing to the nocturnal deepening of mixed layer. In the southern bay, intraseasonal variability is mainly determined by the diurnal cycle in insolation, while in the northern bay, diurnal cycle in insolation and winds have comparable contributions. Temperature inversions (TI) develop in the northern bay in the absence of diurnal variability in wind stress. In the northern bay, SST-ISV amplification is not as large as that in the southern bay due to the weaker diurnal variability of mixed layer depth (MLD) limited by salinity stratification. Diurnal variability of model MLD is not sufficient to create large modifications in mixed layer heat budget and SST-ISV in the northern bay.

Modeling non-stationarity in intensity, duration and frequency of extreme rainfall over India

Significant changes are reported in extreme rainfall characteristics over India in recent studies though there are disagreements on the spatial uniformity and causes of trends. Based on recent theoretical advancements in the Extreme Value Theory (EVT), we analyze changes in extreme rainfall characteristics over India using a high-resolution daily gridded (1 degrees latitude x 1 degrees longitude) dataset. Intensity, duration and frequency of excess rain over a high threshold in the summer monsoon season are modeled by non-stationary distributions whose parameters vary with physical covariates like the El-Nino Southern Oscillation index (ENSO-index) which is an indicator of large-scale natural variability, global average temperature which is an indicator of human-induced global warming and local mean temperatures which possibly indicate more localized changes. Each non-stationary model considers one physical covariate and the best chosen statistical model at each rainfall grid gives the most significant physical driver for each extreme rainfall characteristic at that grid. Intensity, duration and frequency of extreme rainfall exhibit non-stationarity due to different drivers and no spatially uniform pattern is observed in the changes in them across the country. At most of the locations, duration of extreme rainfall spells is found to be stationary, while non-stationary associations between intensity and frequency and local changes in temperature are detected at a large number of locations. This study presents the first application of nonstationary statistical modeling of intensity, duration and frequency of extreme rainfall over India. The developed models are further used for rainfall frequency analysis to show changes in the 100-year extreme rainfall event. Our findings indicate the varying nature of each extreme rainfall characteristic and their drivers and emphasize the necessity of a comprehensive framework to assess resulting risks of precipitation induced flooding. (C) 2014 Elsevier B.V. All rights reserved.

Investigation of orographic features influencing spatial distribution of rainfall over the Western Ghats of India using satellite data

Using remotely sensed Tropical Rainfall Measuring Mission (TRMM) 3B42 rainfall and topographic data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Digital Elevation Model (DEM), the impact of oroghraphical aspects such as topography, spatial variability of elevation and altitude of apexes are examined to investigate capacious summer monsoon rainfall over the Western Ghats (WG) of India. TRMM 3B42 v7 rainfall data is validated with Indian Meteorological Department (IMD) gridded rainfall data at 0.5 degrees resolution over the WG. The analysis of spatial pattern of monsoon rainfall with orography of the WG ascertains that the grade of orographic precipitation depends mainly on topography of the mountain barrier followed by steepness of windward side slope and altitude of the mountain. Longer and broader, i.e. cascaded topography, elevated summits and gradually increasing slopes impel the enhancement in precipitation. Comparing topography of various states of the WG, it has been observed that windward side of Karnataka receives intense rainfall in the WG during summer monsoon. It has been observed that the rainfall is enhanced before the peak of the mountain and confined up to the height about 800m over the WG. In addition to this, the spatial distribution of heavy and very heavy rainfall events in the last 14 years has also been explored. Heavy and very heavy rain events on this hilly terrain are categorized with a threshold of precipitation (R) in the range 150>R>120mmday(-1) and exceeding 150mmday(-1) using probability distribution of TRMM 3B42 v7 rainfall. The areas which are prone to heavy precipitation are identified. The study would help policy makers to manage the hazard scenario and, to improve weather predictions on mountainous terrain of the WG.

Meridional gradients in aerosol vertical distribution over Indian Mainland: Observations and model simulations

Multi-year observations from the network of ground-based observatories (ARFINET), established under the project `Aerosol Radiative Forcing over India' (ARFI) of Indian Space Research Organization and space-borne lidar `Cloud Aerosol Lidar with Orthogonal Polarization' (CALIOP) along with simulations from the chemical transport model `Goddard Chemistry Aerosol Radiation and Transport' (GOCART), are used to characterize the vertical distribution of atmospheric aerosols over the Indian landmass and its spatial structure. While the vertical distribution of aerosol extinction showed higher values close to the surface followed by a gradual decrease at increasing altitudes, a strong meridional increase is observed in the vertical spread of aerosols across the Indian region in all seasons. It emerges that the strong thermal convections cause deepening of the atmospheric boundary layer, which although reduces the aerosol concentration at lower altitudes, enhances the concentration at higher elevations by pumping up more aerosols from below and also helping the lofted particles to reach higher levels in the atmosphere. Aerosol depolarization ratios derived from CALIPSO as well as the GOCART simulations indicate the dominance of mineral dust aerosols during spring and summer and anthropogenic aerosols in winter. During summer monsoon, though heavy rainfall associated with the Indian monsoon removes large amounts of aerosols, the prevailing southwesterly winds advect more marine aerosols over to landmass (from the adjoining oceans) leading to increase in aerosol loading at lower altitudes than in spring. During spring and summer months, aerosol loading is found to be significant, even at altitudes as high as 4 km, and this is proposed to have significant impacts on the regional climate systems such as Indian monsoon. (C) 2015 Elsevier Ltd. All rights reserved.

Investigations of aerosol black carbon from a semi-urban site in the Indo-Gangetic Plain region

Long-term (2009-2012) data from ground-based measurements of aerosol black carbon (BC) from a semi-urban site, Pantnagar (29.0 degrees N, 79.5 degrees E, 231 m amsl), in the Indo-Gangetic Plain (IGP) near the Himalayan foothills are analyzed to study the regional characterization. Large variations are seen in BC at both diurnal and seasonal scales, associated with the mesoscale and synoptic meteorological processes, and local/regional anthropogenic activities. BC diurnal variations show two peaks (morning and evening) arising from the combined effects of the atmospheric boundary layer (ABL) dynamics and local emissions. The diurnal amplitudes as well as the rates of diurnal evolution are the highest in winter season, followed by autumn, and the lowest in summer-monsoon. BC exhibits nearly an inverse relation with mixing layer depth in all seasons; being strongest in winter (R-2 = 0.89) and weakest (R-2 = 0.33) in monsoon (July-August). Unlike BC, co-located aerosol optical depths (AOD) and aerosol absorption are highest in spring over IGP, probably due to the presence of higher abundances of aerosols (including dust) above the ABL (in the free troposphere). AOD (500 nm) showed annual peak (>0.6) in May-June, dominated by coarse mode, while fine mode aerosols dominated in late autumn and early winter. Aerosols profiles from CALIPSO show highest values close to the surface in winter/autumn, similar to the feature seen in surface BC, whereas at altitudes > 2 km, the extinction is maximum in spring/summer. WRF-Chem model is used to simulate BC temporal variations and then compared with observed BC. The model captures most of the important features of the diurnal and seasonal variations but significantly underestimated the observed BC levels, suggesting improvements in diurnal and seasonal varying BC emissions apart from the boundary layer processes. (C) 2015 Elsevier Ltd. All rights reserved.

Influence of the upwelling events on the delta C-13 and delta O-18 of the benthic bivalve shells of the South Western Continental Margin of India

In this study, the influence of the spatial and temporal variability of upwelling intensity and the associated biological productivity observed during different phases of summer monsoon along the southwestern continental margin of India (SWCMI) on the delta C-13 and delta O-18 of the inorganic biogenic carbonate shells was investigated. Multispecies benthic bivalve shells (1-5 mm) separated from ten surface sediment samples of SWCMI (off 12 degrees N, 10 degrees N and 9 degrees N) collected during the onset (OSM) and peak (PSM) phase of the summer monsoon of 2009 were analysed for delta C-13 and delta O-18. Sea surface temperature along the study region indicates prominent upwelling in PSM than in OSM. A comparison of analytical and predicted values for delta O-18 in the bivalve shells confirmed their in situ origin during both the sampling periods. During PSM, the delta C-13 values in the benthic bivalve shells were more depleted in C-13 than during OSM which recorded lower values of delta C-13 in dissolved inorganic carbon of bottom waters expected in the study region in PSM due to the upwelled waters, high surface productivity and the associated high degradation of the organic matter in the subsurface and bottom waters. However, this depletion of delta C-13 was not observed in benthic bivalve shells obtained from 10 degrees N, since it is influenced by high export fluxes of carbon from the Cochin estuary since early monsoon months.

PICES Press, Vol. 13, No. 2, July 2005

PICES in transition: The 3rd inter-sessional Science Board and Governing Council meeting (pdf 0.3 MB) New and upcoming PICES publications (pdf 0.3 MB) PICES hosts GLOBEC Symposium on “Climate Variability and Sub-Arctic Marine Ecosystems” in Victoria (pdf 0.4 MB) First CREAMS/PICES Workshop on East Asian Seas Time-series (pdf 0.4 MB) PICES workshop discusses Oceanic Ecodynamics COmparison in the Subarctic Pacific (OECOS) - a project proposal (pdf 0.3 MB) The state of the western North Pacific in the second half of 2004 (pdf 0.4 MB) The Bering Sea: Current status and recent events (pdf 0.3 MB) Recent trends in waters of the subarctic NE Pacific (pdf 0.3 MB) CFAME workshop on “Developing a working plan for CCCC synthesis” (pdf 0.5 MB) What is winter? (pdf 0.5 MB) The first specimens of Humboldt squid in British Columbia (pdf 0.4 MB) Obituary - Dr. Daniel M. Ware (pdf 0.3 MB) PICES Calendar (pdf 0.3 MB)

Earthquakes of the Nepal Himalaya : towards a physical model of the seismic cycle

Home to hundreds of millions of souls and land of excessiveness, the Himalaya is also the locus of a unique seismicity whose scope and peculiarities still remain to this day somewhat mysterious. Having claimed the lives of kings, or turned ancient timeworn cities into heaps of rubbles and ruins, earthquakes eerily inhabit Nepalese folk tales with the fatalistic message that nothing lasts forever. From a scientific point of view as much as from a human perspective, solving the mysteries of Himalayan seismicity thus represents a challenge of prime importance. Documenting geodetic strain across the Nepal Himalaya with various GPS and leveling data, we show that unlike other subduction zones that exhibit a heterogeneous and patchy coupling pattern along strike, the last hundred kilometers of the Main Himalayan Thrust fault, or MHT, appear to be uniformly locked, devoid of any of the “creeping barriers” that traditionally ward off the propagation of large events. The approximately 20 mm/yr of reckoned convergence across the Himalaya matching previously established estimates of the secular deformation at the front of the arc, the slip accumulated at depth has to somehow elastically propagate all the way to the surface at some point. And yet, neither large events from the past nor currently recorded microseismicity nearly compensate for the massive moment deficit that quietly builds up under the giant mountains. Along with this large unbalanced moment deficit, the uncommonly homogeneous coupling pattern on the MHT raises the question of whether or not the locked portion of the MHT can rupture all at once in a giant earthquake. Univocally answering this question appears contingent on the still elusive estimate of the magnitude of the largest possible earthquake in the Himalaya, and requires tight constraints on local fault properties. What makes the Himalaya enigmatic also makes it the potential source of an incredible wealth of information, and we exploit some of the oddities of Himalayan seismicity in an effort to improve the understanding of earthquake physics and cipher out the properties of the MHT. Thanks to the Himalaya, the Indo-Gangetic plain is deluged each year under a tremendous amount of water during the annual summer monsoon that collects and bears down on the Indian plate enough to pull it away from the Eurasian plate slightly, temporarily relieving a small portion of the stress mounting on the MHT. As the rainwater evaporates in the dry winter season, the plate rebounds and tension is increased back on the fault. Interestingly, the mild waggle of stress induced by the monsoon rains is about the same size as that from solid-Earth tides which gently tug at the planets solid layers, but whereas changes in earthquake frequency correspond with the annually occurring monsoon, there is no such correlation with Earth tides, which oscillate back-and-forth twice a day. We therefore investigate the general response of the creeping and seismogenic parts of MHT to periodic stresses in order to link these observations to physical parameters. First, the response of the creeping part of the MHT is analyzed with a simple spring-and-slider system bearing rate-strengthening rheology, and we show that at the transition with the locked zone, where the friction becomes near velocity neutral, the response of the slip rate may be amplified at some periods, which values are analytically related to the physical parameters of the problem. Such predictions therefore hold the potential of constraining fault properties on the MHT, but still await observational counterparts to be applied, as nothing indicates that the variations of seismicity rate on the locked part of the MHT are the direct expressions of variations of the slip rate on its creeping part, and no variations of the slip rate have been singled out from the GPS measurements to this day. When shifting to the locked seismogenic part of the MHT, spring-and-slider models with rate-weakening rheology are insufficient to explain the contrasted responses of the seismicity to the periodic loads that tides and monsoon both place on the MHT. Instead, we resort to numerical simulations using the Boundary Integral CYCLes of Earthquakes algorithm and examine the response of a 2D finite fault embedded with a rate-weakening patch to harmonic stress perturbations of various periods. We show that such simulations are able to reproduce results consistent with a gradual amplification of sensitivity as the perturbing period get larger, up to a critical period corresponding to the characteristic time of evolution of the seismicity in response to a step-like perturbation of stress. This increase of sensitivity was not reproduced by simple 1D-spring-slider systems, probably because of the complexity of the nucleation process, reproduced only by 2D-fault models. When the nucleation zone is close to its critical unstable size, its growth becomes highly sensitive to any external perturbations and the timings of produced events may therefore find themselves highly affected. A fully analytical framework has yet to be developed and further work is needed to fully describe the behavior of the fault in terms of physical parameters, which will likely provide the keys to deduce constitutive properties of the MHT from seismological observations.