Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

Geomagnetic Storms Over India: A Close Look at Historic Space Weather at Mumbai

The swirling colors of aurorae, familiar to many in polar communities, can occasionally be seen at middle latitudes in locations such as southern Canada and central Europe. But in rare instances, aurorae can even be seen in the tropics. On 6 February 1872, news of the sighting of one such aurora was carried by the Times of India newspaper. The aurora occurred on 4 February 1872 and, as noted, was also observed over the Middle East.

Universality of scaling laws in correlation between velocity and shear stress in turbulent boundary layers

In this study, we analyse simultaneous measurements (at 50 Hz) of velocity at several heights and shear stress at the surface made during the Utah field campaign for the presence of ranges of scales, where distinct scale-to-scale interactions between velocity and shear stress can be identified. We find that our results are similar to those obtained in a previous study [Venugopal et al., 2003] (contrary to the claim in V2003, that the scaling relations might be dependent on Reynolds number) where wind tunnel measurements of velocity and shear stress were analysed. We use a wavelet-based scale-to-scale cross-correlation to detect three ranges of scales of interaction between velocity and shear stress, namely, (a) inertial subrange, where the correlation is negligible; (b) energy production range, where the correlation follows a logarithmic law; and (c) for scales larger than the boundary layer height, the correlation reaches a plateau.

Influences of the springtime northern Indian biomass burning over the central Himalayas

The influences of the springtime northern Indian biomass burning are shown for the first time over the central Himalayas by using three years (2007-2009) of surface and space based observations along with a radiative transfer model. Near-surface ozone, black carbon (BC), spectral aerosol optical depths (AODs) and the meteorological parameters are measured at a high altitude site Nainital (29.37 degrees N, 79.45 degrees E, 1958 m amsl) located in the central Himalayas. The satellite observations include the MODIS derived fire counts and AOD (0.55 mu m), and OMI derived tropospheric column NO(2), ultraviolet aerosol index and single scattering albedo. MODIS fire counts and BC observations are used to identify the fire-impacted periods (372 h during 2007-2009) and hence the induced enhancements in surface BC, AOD (0.5 mu m) and ozone are estimated to be 1802 ng m(-3) (similar to 145%), 0.3 (similar to 150%) and 19 ppbv (similar to 34%) respectively. Large enhancements (53-100%) are also seen in the satellite derived parameters over a 2 degrees x 2 degrees region around Nainital. The present analysis highlights the northern Indian biomass burning induced cooling at the surface (-27 W m(-2)) and top of the atmosphere (-8 W m(-2)) in the lesser polluted high altitude regions of the central Himalayas. This cooling leads to an additional atmospheric warming of 19 W m(-2) and increases the lower atmospheric heating rate by 0.8 K day(-1). These biomass burning induced changes over the central Himalayan atmosphere during spring may also lead to enhanced short-wave absorption above clouds and might have an impact on the monsoonal rainfall.

Silver Nanoparticle Synthesis: Novel Route for Laser Triggering of Polyelectrolyte Capsules

We have demonstrated the synthesis of light-sensitive polyelectrolyte capsules (PECs) by utilizing a novel polyol reduction method and investigated its applicability as photosensitive drug delivery vehicle. The nanostructured capsules were prepared via layer by layer (LbL) assembly of poly(allylamine hydrochloride) (PAN) and dextran sulfate (DS) on silica particles followed by in-situ synthesis of silver nanoparticles (NPs). Capsules without silver NPs were permeable to low molecular weight (A(w), 479 g/mol) rhodamine but impermeable to higher molecular weight fluorescence labeled dextran (FITC-dextran). However, capsules synthesized with silver NPs showed porous morphology and were permeable to higher molecular weight (M(w) 70 kDa) FITC-dextran also. These capsules were loaded with FITC-dextran using thermal encapsulation method by exploiting temperature induced shrinking of the capsules. During heat treatment the porous morphology of the capsules transformed into smooth pore free structure which prevents the movement of dextran into bulk during the loading process. When these loaded capsules are exposed to laser pulses, the capsule wall ruptured, resulting in the release of the loaded drug/dye. The rupture of the capsules was dependent on particle size, laser pulse energy and exposure time. The release was linear with time when pulse energy of 400 mu J was used and burst release was observed when pulse energy increased to 600 mu J.

Mechanisms of summer intraseasonal sea surface temperature oscillations in the Bay of Bengal

Intraseasonal variations (ISV) of sea surface temperature (SST) in the Bay of Bengal (BoB) is highest in its northwestern part. An Indian Ocean model forced by QuikSCAT winds and climatological river discharge (QR run) reproduces ISV of SST, albeit with weaker magnitude. Air-sea fluxes, in the presence of a shallow mixed layer, efficiently effect intraseasonal SST fluctuations. Warming during intraseasonal events is smaller (<1°C) for June - July period and larger (1.5° to 2°C) during September, the latter due to a thinner mixed layer. To examine the effect of salinity on ISV, the model was run by artificially increasing the salinity (NORR run) and by decreasing it (MAHA10 run). In NORR, both rainfall and river discharge were switched off and in MAHA10 the discharge by river Mahanadi was increased tenfold. The spatial pattern of ISV as well as its periodicity was similar in QR, NORR and MAHA10. The ISV was stronger in NORR and weaker in MAHA10, compared to QR. In NORR, both intraseasonal warming and cooling were higher than in QR, the former due to reduced air-sea heat loss as the mean SST was lower, and the latter due to enhanced subsurface processes resulting from weaker stratification. In MAHA10, both warming and cooling were lower than in QR, the former due to higher air-sea heat loss owing to higher mean SST, and the latter due to weak subsurface processes resulting from stronger stratification. These model experiments suggest that salinity effects are crucial in determining amplitudes of intraseasonal SST variations in the BoB.

Signatures of Indian ocean Dipole and El Nino-Southern Oscillation events in sea level variations in the Bay of Bengal

We investigate the impact of the Indian Ocean Dipole (IOD) and El Nino and the Southern Oscillation (ENSO) on sea level variations in the North Indian Ocean during 1957-2008. Using tide-gauge and altimeter data, we show that IOD and ENSO leave characteristic signatures in the sea level anomalies (SLAs) in the Bay of Bengal. During a positive IOD event, negative SLAs are observed during April-December, with the SLAs decreasing continuously to a peak during September-November. During El Nino, negative SLAs are observed twice (April-December and November-July), with a relaxation between the two peaks. SLA signatures during negative IOD and La Nina events are much weaker. We use a linear, continuously stratified model of the Indian Ocean to simulate their sea level patterns of IOD and ENSO events. We then separate solutions into parts that correspond to specific processes: coastal alongshore winds, remote forcing from the equator via reflected Rossby waves, and direct forcing by interior winds within the bay. During pure IOD events, the SLAs are forced both from the equator and by direct wind forcing. During ENSO events, they are primarily equatorially forced, with only a minor contribution from direct wind forcing. Using a lead/lag covariance analysis between the Nino-3.4 SST index and Indian Ocean wind stress, we derive a composite wind field for a typical El Nino event: the resulting solution has two negative SLA peaks. The IOD and ENSO signatures are not evident off the west coast of India.

On the basin-scale detection and attribution of human-induced climate change in monsoon precipitation and streamflow

Detecting and quantifying the presence of human-induced climate change in regional hydrology is important for studying the impacts of such changes on the water resources systems as well as for reliable future projections and policy making for adaptation. In this article a formal fingerprint-based detection and attribution analysis has been attempted to study the changes in the observed monsoon precipitation and streamflow in the rain-fed Mahanadi River Basin in India, considering the variability across different climate models. This is achieved through the use of observations, several climate model runs, a principal component analysis and regression based statistical downscaling technique, and a Genetic Programming based rainfall-runoff model. It is found that the decreases in observed hydrological variables across the second half of the 20th century lie outside the range that is expected from natural internal variability of climate alone at 95% statistical confidence level, for most of the climate models considered. For several climate models, such changes are consistent with those expected from anthropogenic emissions of greenhouse gases. However, unequivocal attribution to human-induced climate change cannot be claimed across all the climate models and uncertainties in our detection procedure, arising out of various sources including the use of models, cannot be ruled out. Changes in solar irradiance and volcanic activities are considered as other plausible natural external causes of climate change. Time evolution of the anthropogenic climate change ``signal'' in the hydrological observations, above the natural internal climate variability ``noise'' shows that the detection of the signal is achieved earlier in streamflow as compared to precipitation for most of the climate models, suggesting larger impacts of human-induced climate change on streamflow than precipitation at the river basin scale.

Computing Reeb Graphs as a Union of Contour Trees

The Reeb graph of a scalar function tracks the evolution of the topology of its level sets. This paper describes a fast algorithm to compute the Reeb graph of a piecewise-linear (PL) function defined over manifolds and non-manifolds. The key idea in the proposed approach is to maximally leverage the efficient contour tree algorithm to compute the Reeb graph. The algorithm proceeds by dividing the input into a set of subvolumes that have loop-free Reeb graphs using the join tree of the scalar function and computes the Reeb graph by combining the contour trees of all the subvolumes. Since the key ingredient of this method is a series of union-find operations, the algorithm is fast in practice. Experimental results demonstrate that it outperforms current generic algorithms by a factor of up to two orders of magnitude, and has a performance on par with algorithms that are catered to restricted classes of input. The algorithm also extends to handle large data that do not fit in memory.

Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics

The way in which basal tractions, associated with mantle convection, couples with the lithosphere is a fundamental problem in geodynamics. A successful lithosphere-mantle coupling model for the Earth will satisfy observations of plate motions, intraplate stresses, and the plate boundary zone deformation. We solve the depth integrated three-dimensional force balance equations in a global finite element model that takes into account effects of both topography and shallow lithosphere structure as well as tractions originating from deeper mantle convection. The contribution from topography and lithosphere structure is estimated by calculating gravitational potential energy differences. The basal tractions are derived from a fully dynamic flow model with both radial and lateral viscosity variations. We simultaneously fit stresses and plate motions in order to delineate a best-fit lithosphere-mantle coupling model. We use both the World Stress Map and the Global Strain Rate Model to constrain the models. We find that a strongly coupled model with a stiff lithosphere and 3-4 orders of lateral viscosity variations in the lithosphere are best able to match the observational constraints. Our predicted deviatoric stresses, which are dominated by contribution from mantle tractions, range between 20-70 MPa. The best-fitting coupled models predict strain rates that are consistent with observations. That is, the intraplate areas are nearly rigid whereas plate boundaries and some other continental deformation zones display high strain rates. Comparison of mantle tractions and surface velocities indicate that in most areas tractions are driving, although in a few regions, including western North America, tractions are resistive. Citation: Ghosh, A., W. E. Holt, and L. M. Wen (2013), Predicting the lithospheric stress field and plate motions by joint modeling of lithosphere and mantle dynamics.

Buildup of aerosols over the Indian Region

Climate change has great significance globally in general and South Asia in particular. Here we have used data from a network of 35 aerosol observatories over the Indian region to generate the first time regional synthesis using primary data and estimate the aerosol trends. On an average, aerosol optical depth (AOD) was found increasing at a rate of 2.3% (of its value in 1985) per year and more rapidly (similar to 4%) during the last decade. If the trends continue so, AOD at several locations would nearly double and approach unity in the next few decades leading to an enhancement in aerosol-induced lower atmospheric warming by a factor of two. However, a regionally averaged scenario can be ascertained only in the coming years, when longer and denser data would become available. The regional and global climate implications of such trends in the forcing elements need to be better assessed using GCMs.