The Alfvén Wave Equation for Magnetospheric Plasmas

It is shown that besides the continuous spectrum which damps away as inverse power of time, the coupled Alfvén wave equation, which gives coupling between a shear Alfvén wave and a surface wave, can also admit a well behaved harmonic solution in the closed form for a set of initial conditions. This solution, though valid for finite time intervals, points out that the Alfvén surface waves can have a band of frequency (instead of a monochromatic frequency for a nonsheared magnetic field) within which the local field line resonance frequency can lie, and thus can excite magnetic pulsations with latitude-dependent frequency. By considering magnetic fields not only varying in magnitude but also in direction, it is shown that the time interval for the validity of the harmonic solution depend upon the angle between the magnetic field directions on either side of the magnetopause. For small values of the angle the time interval can become appreciably large.

Improved retrieval of aerosol size-resolved properties from moderate resolution imaging spectroradiometer over India: Role of aerosol model and surface reflectance

We have compared the total as well as fine mode aerosol optical depth (tau and tau(fine)) retrieved by Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua (2001-2005) with the equivalent parameters derived by Aerosol Robotic Network (AERONET) at Kanpur (26.45 degrees N, 80.35 degrees E), northern India. MODIS Collection 005 (C005)-derived tau(0.55) was found to be in good agreement with the AERONET measurements. The tau(fine) and eta (tau(fine)/tau) were, however, biased low significantly in most matched cases. A new set of retrieval with the use of absorbing aerosol model (SSA similar to 0.87) with increased visible surface reflectance provided improved tau and tau(fine) at Kanpur. The new derivation of eta also compares well qualitatively with an independent set of in situ measurements of accumulation mass fraction over much of the southern India. This suggests that though MODIS land algorithm has limited information to derive size properties of aerosols over land, more accurate parameterization of aerosol and surface properties within the existing C005 algorithm may improve the accuracy of size-resolved aerosol optical properties. The results presented in this paper indicate that there is a need to reconsider the surface parameterization and assumed aerosol properties in MODIS C005 algorithm over the Indian region in order to retrieve more accurate aerosol optical and size properties, which are essential to quantify the impact of human-made aerosols on climate.

Convergence feedback and unstable low frequency oscillations in a simple coupled ocean‐atmosphere model

The role of convergence feedback on the stability of a coupled ocean‐atmosphere system is studied using model III of Hirst (1986). It is shown that the unstable coupled mode found by Hirst is greatly modified by the convergence feedback. If the convergence feedback strength exceeds a critical value, several new unstable intraseasonal modes are also introduced. These modes have very weak dependence on the wave number. These results may explain the behaviour of some coupled models and to some extent provide a mechanism for the observed aperiodicity of the El‐Nino and Southern Oscillation (ENSO) events.

Comment on “Southward magnetic field in the neutral sheet produced by wavy motions Propagating in the dawn‐dusk direction”

In a recent paper Nakagawa and Nishida [1989] have suggested that wavy motions of the neutral sheet can be generated by the Kelvin‐Helmholtz instability if the dawn‐dusk flow of only several tens of km/s is present. However, their mathematical analysis is based on the choice of particular magnetic field directions in the three regions consisting of north, south lobes and the neutral sheet. In an earlier paper Uberoi [1986] discussed the Kelvin‐Helmholtz instability of a similar structured plasma layer without any assumptions either on velocity field directions or on the magnetic field directions, thus pointing out the angle effect due to variation in magnetic field directions on the instability criterion. The relevance of these results to the problem of wavy motions of the neutral sheet are pointed out. In particular it is found that when the y‐component of the magnetic field in each lobe is taken into consideration the Kelvin‐Helmholtz instability can be exicted only when the dawn‐dusk flow is of several hundreds of km/s a order of ten higher than that arrived in the analysis by Nakagawa and Nishida [1989].

Surprising observation of large anthropogenic aerosol fraction over the ``near-pristine'' southern Bay of Bengal: Climate implications

The Bay of Bengal (BoB), a small oceanic region surrounded by landmasses with distinct natural and anthropogenic activities and under the influence of seasonally changing airmass types, is characterized by a rather complex and highly heterogeneous aerosol environment. Concurrent measurements of the physical, optical, and chemical (offline analysis) properties of BoB aerosols, made onboard extensive ship-cruises and aircraft sorties during Integrated Campaign for Aerosols, gases and Radiation Budget of March-April 2006, and satellite-retrieved aerosol optical depths and derived parameters, were synthesized following a synergistic approach to delineate the anthropogenic fraction to the composite aerosol parameters and its spatial variation. Quite interestingly and contrary to the general belief, our studies revealed that, despite of the very high aerosol loading (in the marine atmospheric boundary layer as well as in the vertical column) over the northern BoB and a steep decreasing gradient toward the southern latitudes, the anthropogenic fraction showed a steady increase from North to South (where no obvious anthropogenic source regions exist). Consequently, the direct radiative forcing at the top of the atmosphere due to anthropogenic aerosols remained nearly constant over the entire BoB with values in the range from -3.3 to -3.6 Wm(-2). This interesting finding, beyond doubts calls for a better understanding of the complex aerosol system over the BoB through more focused field campaigns.

Exact nonlinear travelling hydromagnetic wave solutions

Exact travelling wave solutions for hydromagnetic waves in an exponentially stratified incompressible medium are obtained. With the help of two integrals it becomes possible to reduce the system of seven nonlinear PDE's to a second order nonlinear ODE which describes an one dimensional harmonic oscillator with a nonlinear friction term. This equation is studied in detail in the phase plane. The travelling waves are periodic only when they propagate either horizontally or vertically. The reduced second order nonlinear differential equation describing the travelling waves in inhomogeneous conducting media has rather ubiquitous nature in that it also appears in other geophysical systems such as internal waves, Rossby waves and topographic Rossby waves in the ocean.

Impact of physical processes on the chlorophyll distribution in the Bay of Bengal

Satellite-derived chlorophyll a concentration (chl a) maps show three regions with high chl a in the Bay of Bengal. First among these is close to the coast, particularly off river mouths, with high values coinciding with the season of peak discharge; second is in the southwestern bay during the northeast monsoon, which is forced by local Ekman pumping; and the third is to the east of Sri Lanka in response to the summer monsoon winds. Chlorophyll-rich water from the mouths of rivers flows either along the coast or in an offshore direction, up to several hundred kilometers, depending on the prevailing ocean current pattern. The Irrawady River plume flows toward offshore and then turns northwestward during October–December, but it flows along the coast into the Andaman Sea for the rest of the year. From the Ganga-Brahmaputra river mouth, chl a–rich water flows directly southward into the open bay during spring but along the Indian coast during summer and winter. Along the Indian coast, the flow of chl a–rich water is determined by the East India Coastal Current (EICC). Whenever the EICC meanders off the Indian coast, it leads to an offshore outbreak of chl a–rich water from the coastal region into open ocean. The EICC as well as open ocean circulation in the bay is made up of several eddies, and these eddies show relatively higher chl a. Eddies near the coast, however, can often have higher chl a because of advection from the coastal region rather than generation within the eddy itself. The bay experiences several cyclones in a year, most of them occurring during October–November. These cyclones cause a drop in the sea surface temperature, a dip in the sea level, and a local increase in chl a. The impact of a cyclone is weaker in the northern part of the bay because of stronger stratification compared to the southern parts.

A possible origin of viscosity in Keplerian accretion disks due to secondary perturbation: Turbulent transport without magnetic fields

The origin of hydrodynamic turbulence in rotating shear flow is a long standing puzzle. Resolving it is especially important in astrophysics when the flow's angular momentum profile is Keplerian which forms an accretion disk having negligible molecular viscosity. Hence, any viscosity in such systems must be due to turbulence, arguably governed by magnetorotational instability, especially when temperature T greater than or similar to 10(5). However, such disks around quiescent cataclysmic variables, protoplanetary and star-forming disks, and the outer regions of disks in active galactic nuclei are practically neutral in charge because of their low temperature, and thus are not expected to be coupled with magnetic fields enough to generate any transport due to the magnetorotational instability. This flow is similar to plane Couette flow including the Coriolis force, at least locally. What drives their turbulence and then transport, when such flows do not exhibit any unstable mode under linear hydrodynamic perturbation? We demonstrate that the three-dimensional secondary disturbance to the primarily perturbed flow that triggers elliptical instability may generate significant turbulent viscosity in the range 0.0001 less than or similar to nu(t) less than or similar to 0.1, which can explain transport in accretion flows.

Paleoseismological evidence of surface faulting along the northeastern Himalayan front, India: Timing, size, and spatial extent of great earthquakes

The similar to 2500 km long Himalayan arc has experienced three large to great earthquakes of M-w 7.8 to 8.4 during the past century, but none produced surface rupture. Paleoseismic studies have been conducted during the last decade to begin understanding the timing, size, rupture extent, return period, and mechanics of the faulting associated with the occurrence of large surface rupturing earthquakes along the similar to 2500 km long Himalayan Frontal Thrust (HFT) system of India and Nepal. The previous studies have been limited to about nine sites along the western two-thirds of the HFT extending through northwest India and along the southern border of Nepal. We present here the results of paleoseismic investigations at three additional sites further to the northeast along the HFT within the Indian states of West Bengal and Assam. The three sites reside between the meizoseismal areas of the 1934 Bihar-Nepal and 1950 Assam earthquakes. The two westernmost of the sites, near the village of Chalsa and near the Nameri Tiger Preserve, show that offsets during the last surface rupture event were at minimum of about 14 m and 12 m, respectively. Limits on the ages of surface rupture at Chalsa (site A) and Nameri (site B), though broad, allow the possibility that the two sites record the same great historical rupture reported in Nepal around A.D. 1100. The correlation between the two sites is supported by the observation that the large displacements as recorded at Chalsa and Nameri would most likely be associated with rupture lengths of hundreds of kilometers or more and are on the same order as reported for a surface rupture earthquake reported in Nepal around A.D. 1100. Assuming the offsets observed at Chalsa and Nameri occurred synchronously with reported offsets in Nepal, the rupture length of the event would approach 700 to 800 km. The easternmost site is located within Harmutty Tea Estate (site C) at the edges of the 1950 Assam earthquake meizoseismal area. Here the most recent event offset is relatively much smaller (<2.5 m), and radiocarbon dating shows it to have occurred after A.D. 1100 (after about A.D. 1270). The location of the site near the edge of the meizoseismal region of the 1950 Assam earthquake and the relatively lesser offset allows speculation that the displacement records the 1950 M-w 8.4 Assam earthquake. Scatter in radiocarbon ages on detrital charcoal has not resulted in a firm bracket on the timing of events observed in the trenches. Nonetheless, the observations collected here, when taken together, suggest that the largest of thrust earthquakes along the Himalayan arc have rupture lengths and displacements of similar scale to the largest that have occurred historically along the world's subduction zones.

Why is there a short-term increase in global precipitation in response to diminished CO2 forcing?

Recently, it was found that a reduction in atmospheric CO2 concentration leads to a temporary increase in global precipitation. We use the Hadley Center coupled atmosphere-ocean model, HadCM3L, to demonstrate that this precipitation increase is a consequence of precipitation sensitivity to changes in atmospheric CO2 concentrations through fast tropospheric adjustment processes. Slow ocean cooling explains the longer-term decrease in precipitation. Increased CO2 tends to suppress evaporation/precipitation whereas increased temperatures tend to increase evaporation/precipitation. When the enhanced CO2 forcing is removed, global precipitation increases temporarily, but this increase is not observed when a similar negative radiative forcing is applied as a reduction of solar intensity. Therefore, transient precipitation increase following a reduction in CO2-radiative forcing is a consequence of the specific character of CO2 forcing and is not a general feature associated with decreases in radiative forcing. Citation: Cao, L., G. Bala, and K. Caldeira (2011), Why is there a short-term increase in global precipitation in response to diminished CO2 forcing?, Geophys. Res. Lett., 38, L06703, doi:10.1029/2011GL046713.

An integrated model for optimal reservoir operation for irrigation of multiple crops

An integrated model is developed, based on seasonal inputs of reservoir inflow and rainfall in the irrigated area, to determine the optimal reservoir release policies and irrigation allocations to multiple crops. The model is conceptually made up of two modules, Module 1 is an intraseasonal allocation model to maximize the sum of relative yields of all crops, for a given state of the system, using linear programming (LP). The module takes into account reservoir storage continuity, soil moisture balance, and crop root growth with time. Module 2 is a seasonal allocation model to derive the steady state reservoir operating policy using stochastic dynamic programming (SDP). Reservoir storage, seasonal inflow, and seasonal rainfall are the state variables in the SDP. The objective in SDP is to maximize the expected sum of relative yields of all crops in a year. The results of module 1 and the transition probabilities of seasonal inflow and rainfall form the input for module 2. The use of seasonal inputs coupled with the LP-SDP solution strategy in the present formulation facilitates in relaxing the limitations of an earlier study, while affecting additional improvements. The model is applied to an existing reservoir in Karnataka State, India.

Hydrography and circulation in the western Bay of Bengal during the northeast monsoon

The Bay of Bengal, a semienclosed tropical basin that comes under the influence of monsoonal wind and freshwater influx, is distinguished by a strongly stratified surface layer and a seasonally reversing circulation. We discuss characteristics of these features in the western Bay during the northeast monsoon, when the East India Coastal Current (EICC) flows southward, using hydrographic data collected during December 1991. Vertical profiles show uniform temperature and salinity in a homogeneous surface layer, on average, 25 m deep but shallower northward and coastward. The halocline, immediately below, is approximately 50 m thick; salinity changes by approximately 3 parts per thousand. About two thirds of the profiles show temperature inversions in this layer. Salinity below the halocline hardly changes, and stratification is predominantly due to temperature variation, The halocline is noticeably better developed and the surface homogeneous layer is thinner in a low-salinity plume that hugs the coastline along the entire east coast of India, The plume is, on average, 50 km wide, with isohalines sloping down toward the coast. Most prominent in the geostrophic velocity field is the equatorward EICC. Its transport north of about 13 degrees N, computed with 1000 dbar as the level of reference, varies between 2.6 and 7.1 x 10(6) m(3) s(-1); just south of this latitude, a northwestward flow from offshore recurves and merges with the coastal current. At the southern end of the region surveyed, the transport is 7.7 x 10(6) m(3) s(-1). Recent model studies lead us to conclude that the EICC during the northeast monsoon is driven by winds along the east coast of India and Ekman pumping in the interior bay. In the south, Ekman pumping over the southwestern bay is responsible for the northwestward flow that merges with the EICC.

Interannual variations of Indian summer monsoon in a GCM: External conditions versus internal feedbacks

The potential predictability of the Indian summer monsoon due to slowly varying sea surface temperature (SST) forcing is examined. Factors responsible for limiting the predictability are also investigated. Three multiyear simulations with the R30 version of the Geophysical Fluid Dynamics Laboratory's climate model are carried out for this purpose, The mean monsoon simulated by this model is realistic including the mean summer precipitation over the Indian continent. The interannual variability of the large-scale component of the monsoon such as the "monsoon shear index" and its teleconnection with Pacific SST is well simulated by the model in a 15-yr integration with observed SST as boundary condition. On regional scales, the skill in simulating the interannual variability of precipitation over the Indian continent by the model is rather modest and its simultaneous correlation with eastern Pacific SST is negative but poor as observed. The poor predictability of precipitation over the Indian region in the model is related to the fact that contribution to the interannual variability over this region due to slow SST variations [El Nino-Southern Oscillation (ENSO) related] is comparable to those due to regional-scale fluctuations unrelated to ENSO SST. The physical mechanism through which ENSO SST tend to produce reduction in precipitation over the Indian continent is also elucidated. A measure of internal variability of the model summer monsoon is obtained from a 20-yr integration of the same model with fixed annual cycle SST as boundary conditions but with predicted soil moisture and snow cover. A comparison of summer monsoon indexes between this run and the observed SST run shows that the internal oscillations can account for a large fraction of the simulated monsoon variability. The regional-scale oscillations in the observed SST run seems to arise from these internal oscillations. It is discovered that most of the interannual internal variability is due to an internal quasi-biennial oscillation (QBO) of the model atmosphere. Such a QBO is also found in the author's third 18-yr simulation in which fixed annual cycle of SST as well as soil moisture and snow cover are prescribed. This shows that the model QBO is not due to land-surface-atmosphere interaction. It is proposed that the model QBO arises due to an interaction between nonlinear intraseasonal oscillations and the annual cycle. Spatial structure of the QBO and its role in limiting the predictability of the Indian summer monsoon is discussed.

Predictive uncertainty of chaotic daily streamflow using ensemble wavelet networks approach

Perfect or even mediocre weather predictions over a long period are almost impossible because of the ultimate growth of a small initial error into a significant one. Even though the sensitivity of initial conditions limits the predictability in chaotic systems, an ensemble of prediction from different possible initial conditions and also a prediction algorithm capable of resolving the fine structure of the chaotic attractor can reduce the prediction uncertainty to some extent. All of the traditional chaotic prediction methods in hydrology are based on single optimum initial condition local models which can model the sudden divergence of the trajectories with different local functions. Conceptually, global models are ineffective in modeling the highly unstable structure of the chaotic attractor. This paper focuses on an ensemble prediction approach by reconstructing the phase space using different combinations of chaotic parameters, i.e., embedding dimension and delay time to quantify the uncertainty in initial conditions. The ensemble approach is implemented through a local learning wavelet network model with a global feed-forward neural network structure for the phase space prediction of chaotic streamflow series. Quantification of uncertainties in future predictions are done by creating an ensemble of predictions with wavelet network using a range of plausible embedding dimensions and delay times. The ensemble approach is proved to be 50% more efficient than the single prediction for both local approximation and wavelet network approaches. The wavelet network approach has proved to be 30%-50% more superior to the local approximation approach. Compared to the traditional local approximation approach with single initial condition, the total predictive uncertainty in the streamflow is reduced when modeled with ensemble wavelet networks for different lead times. Localization property of wavelets, utilizing different dilation and translation parameters, helps in capturing most of the statistical properties of the observed data. The need for taking into account all plausible initial conditions and also bringing together the characteristics of both local and global approaches to model the unstable yet ordered chaotic attractor of a hydrologic series is clearly demonstrated.