Dust aerosols over India and adjacent continents retrieved using METEOSAT infrared radiance Part II: quantification of wind dependence and estimation of radiative forcing

Long-range transport of continental dust makes these particles a significant constituent even at locations far from their sources. It is important to study the temporal variations in dust loading over desert regions and the role of meteorology, in order to assess its radiative impact. In this paper, infrared radiance (10.5-12.5 mu m), acquired by the METEOSAT-5 satellite (similar to 5-km resolution) during 1999 and 2003 was used to quantify wind dependence of dust aerosols and to estimate the radiative forcing. Our analysis shows that the frequency of occurrence of dust events was higher during 2003 compared to 1999. Since the dust production function depends mainly on the surface wind speed over regions which are dry and without vegetation, the role of surface wind on IDDI was examined in detail. It was found that an increase of IDDI with wind speed was nearly linear and the rate of increase in IDDI with surface wind was higher during 2003 compared to 1999. It was also observed that over the Indian desert, when wind speed was the highest during monsoon months (June to August), the dust production rate was lower because of higher soil moisture (due to monsoon rainfall). Over the Arabian deserts, when the wind speed is the highest during June to August, the dust production rate is also highest, as soil moisture is lowest during this season. Even though nothing can be said precisely on the reason why 2003 had a greater number of dust events, examination of monthly mean soil moisture at source regions indicates that the occurrence of high winds simultaneous with high soil moisture could be the reason for the decreased dust production efficiency in 1999. It appears that the deserts of Northwest India are more efficient dust sources compared to the deserts of Saudi Arabia and Northeast Africa (excluding Sahara). The radiative impact of dust over various source regions is estimated, and the regionally and annually averaged top of the atmosphere dust radiative forcing (short wave, clear-sky and over land) over the entire study region (0-35 degrees N; 30 degrees-100 degrees E) was in the range of -0.9 to +4.5 W m(-2). The corresponding values at the surface were in the range of -10 to -25 W m(-2). Our studies demonstrate that neglecting the diurnal variation of dust can cause errors in the estimation of long wave dust forcing by as much as 50 to 100%, and nighttime retrieval of dust can significantly reduce the uncertainties. A method to retrieve dust aerosols during nighttime is proposed. The regionally and annually averaged long wave dust radiative forcing was +3.4 +/- 1.6 W m(-2).

Coherent rainfall zones: Case study for Karnataka

Generally average rainfall over meteorological subdivisions is used for assessment of the variability of monsoon rainfall. It is shown here that variations of seasonal rainfall over the meteorological subdivisions of interior Karnataka are not coherent. A methodology for delineating coherent rainfall zones is developed in this paper and applied to derive such zones for the State of Karnataka.

Retrieval of dust aerosols during night: improved assessment of long wave dust radiative forcing over Afro-Asian regions

Several investigators in the past have used the radiance depression (with respect to clear-sky infrared radiance), resulting from the presence of mineral dust aerosols in the atmosphere, as an index of dust aerosol load in the atmosphere during local noon. Here, we have used a modified approach to retrieve dust index during night since assessment of diurnal average infrared dust forcing essentially requires information on dust aerosols during night. For this purpose, we used infrared radiance (10.5-12.5 mu m), acquired from the METEOSAT-5 satellite (similar to 5 km resolution). We found that the `dust index' algorithm, valid for daytime, will no longer hold during the night because dust is then hotter than the theoretical dust-free reference. Hence we followed a `minimum reference' approach instead of a conventional `maximum reference' approach. A detailed analysis suggests that the maximum dust load occurs during the daytime. Over the desert regions of India and Africa, maximum change in dust load is as much as a factor of four between day and night and factor of two variations are commonly observed. By realizing the consequent impact on long wave dust forcing, sensitivity studies were carried out, which indicate that utilizing day time data for estimating the diurnally averaged long-wave dust radiative forcing results in significant errors (as much as 50 to 70%). Annually and regionally averaged long wave dust radiative forcing (which account for the diurnal variation of dust) at the top of the atmosphere over Afro-Asian region is 2.6 +/- 1.8 W m(-2), which is 30 to 50% lower than those reported earlier. Our studies indicate that neglecting diurnal variation of dust while assessing its radiative impact leads to an overestimation of dust radiative forcing, which in turn result in underestimation of the radiative impact of anthropogenic aerosols.

Embracing Integrability in Stellar and Planetary Dynamics

Hamiltonian systems in stellar and planetary dynamics are typically near integrable. For example, Solar System planets are almost in two-body orbits, and in simulations of the Galaxy, the orbits of stars seem regular. For such systems, sophisticated numerical methods can be developed through integrable approximations. Following this theme, we discuss three distinct problems. We start by considering numerical integration techniques for planetary systems. Perturbation methods (that utilize the integrability of the two-body motion) are preferred over conventional "blind" integration schemes. We introduce perturbation methods formulated with Cartesian variables. In our numerical comparisons, these are superior to their conventional counterparts, but, by definition, lack the energy-preserving properties of symplectic integrators. However, they are exceptionally well suited for relatively short-term integrations in which moderately high positional accuracy is required. The next exercise falls into the category of stability questions in solar systems. Traditionally, the interest has been on the orbital stability of planets, which have been quantified, e.g., by Liapunov exponents. We offer a complementary aspect by considering the protective effect that massive gas giants, like Jupiter, can offer to Earth-like planets inside the habitable zone of a planetary system. Our method produces a single quantity, called the escape rate, which characterizes the system of giant planets. We obtain some interesting results by computing escape rates for the Solar System. Galaxy modelling is our third and final topic. Because of the sheer number of stars (about 10^11 in Milky Way) galaxies are often modelled as smooth potentials hosting distributions of stars. Unfortunately, only a handful of suitable potentials are integrable (harmonic oscillator, isochrone and Stäckel potential). This severely limits the possibilities of finding an integrable approximation for an observed galaxy. A solution to this problem is torus construction; a method for numerically creating a foliation of invariant phase-space tori corresponding to a given target Hamiltonian. Canonically, the invariant tori are constructed by deforming the tori of some existing integrable toy Hamiltonian. Our contribution is to demonstrate how this can be accomplished by using a Stäckel toy Hamiltonian in ellipsoidal coordinates.

Radiative interactions in boundary layers

A detailed description of radiative interactions in laminar compressible boundary layers for moderate Mach numbers is presented by way of asymptotic analysis and supporting solutions. The radiation field is described by the differential approximation. While the asymptotic analysis is valid for large N (the ratio of photon mean free path to molecular mean free path) and arbitrary Boltzmann number, Bo (the ratio of convective heat flux to radiation heat flux), the solutions are obtained for Bo [double less-than sign] 1, the case of strong radiative interactions. The asymptotic analysis shows the existence of an optically thin boundary layer for large N and all Bo. For Bo [double less-than sign] 1, two outer regions are observed — one optically thin (at short distances from the leading edge) and the other optically thick (at large distances from the leading edge). An interesting feature not pointed out in the previous literature is the existence of a wall layer at large distances from the leading edge where convective heat flux can be ignored to the leading order of approximation. The radiation field in all cases can be very well approximated by a one-dimensional description. The solutions have been constructed using the ideas of matched asymptotic expansions by approximate analytical procedures and numerical methods. It is shown that, to the leading order of approximation, the radiation slip method yields exactly the same result as the more complicated matching procedure. Both the cases of linear and nonlinear radiation have been considered, the former being of interest in developing approximate methods which are subsequently generalized to handle the nonlinear problem. Detailed results are presented for both cases.

Radiative effects of aerosols at an urban location in southern India: Observations versus model

The radiative impact of aerosols is one of the largest sources of uncertainty in estimating anthropogenic climate perturbations. Here we have used independent ground-based radiometer measurements made simultaneously with comprehensive measurements of aerosol microphysical and optical properties at a highly populated urban site, Bangalore (13.02 degrees N, 77.6 degrees E) in southern India during a dedicated campaign during winter of 2004 and summer and pre-monsoon season of 2005. We have also used longer term measurements carried out at this site to present general features of aerosols over this region. The aerosol radiative impact assessments were made from direct measurements of ground reaching irradiance as well as by incorporating measured aerosol properties into a radiative transfer model. Large discrepancies were observed between measured and modeled (using radiative transfer models, which employed measured aerosol properties) radiative impacts. It appears that the presence of elevated aerosol layers and (or) inappropriate description of aerosol state of mixing are (is) responsible for the discrepancies. On a monthly scale reduction of surface irradiance due to the presence of aerosols (estimated using radiative flux measurements) varies from 30 to 65 W m(-2). The lowest values in surface radiative impact were observed during June when there is large reduction in aerosol as a consequence of monsoon rainfall. Large increase in aerosol-induced surface radiative impact was observed from winter to summer. Our investigations re-iterate the inadequacy of aerosol measurements at the surface alone and importance of representing column properties (using vertical profiles) accurately in order to assess aerosol-induced climate changes accurately. (C) 2010 Elsevier Ltd. All rights reserved.

Analytical and numerical solutions of inward spherical solidification of a superheated melt with radiative-convective heat transfer and density jump at freezing front

Analytical and numerical solutions of a general problem related to the radially symmetric inward spherical solidification of a superheated melt have been studied in this paper. In the radiation-convection type boundary conditions, the heat transfer coefficient has been taken as time dependent which could be infinite, at time,t=0. This is necessary, for the initiation of instantaneous solidification of superheated melt, over its surface. The analytical solution consists of employing suitable fictitious initial temperatures and fictitious extensions of the original region occupied by the melt. The numerical solution consists of finite difference scheme in which the grid points move with the freezing front. The numerical scheme can handle with ease the density changes in the solid and liquid states and the shrinkage or expansions of volumes due to density changes. In the numerical results, obtained for the moving boundary and temperatures, the effects of several parameters such as latent heat, Boltzmann constant, density ratios, heat transfer coefficients, etc. have been shown. The correctness of numerical results has also been checked by satisfying the integral heat balance at every timestep.

Coherent rainfall zones of the Indian region

We have delineated rainfall zones for the Indian region that are coherent with respect to the variations of the summer monsoon rainfall. Within each zone, the time series of the summer monsoon rainfall at every pair of stations are significantly positively correlated, and the mean interseries correlation for each zone is high. The interseries correlation data set is analysed in order to delineate the rainfall zones, using an objective method specifically developed for the purpose. Each of the zonal averages are shown to be representative of the zone as a whole. We suggest that this regionalization is appropriate for study of the variation of the summer monsoon rainfall over the Indian region on interannual and larger scales.

Meridional propagation of large-scale monsoon convective zones

Observational studies indicate that the convective activity of the monsoon systems undergo intraseasonal variations with multi-week time scales. The zone of maximum monsoon convection exhibits substantial transient behavior with successive propagating from the North Indian Ocean to the heated continent. Over South Asia the zone achieves its maximum intensity. These propagations may extend over 3000 km in latitude and perhaps twice the distance in longitude and remain as coherent entities for periods greater than 2-3 weeks. Attempts to explain this phenomena using simple ocean-atmosphere models of the monsoon system had concluded that the interactive ground hydrology so modifies the total heating of the atmosphere that a steady state solution is not possible, thus promoting lateral propagation. That is, the ground hydrology forces the total heating of the atmosphere and the vertical velocity to be slightly out of phase, causing a migration of the convection towards the region of maximum heating. Whereas the lateral scale of the variations produced by the Webster (1983) model were essentially correct, they occurred at twice the frequency of the observed events and were formed near the coastal margin, rather than over the ocean. Webster's (1983) model used to pose the theories was deficient in a number of aspects. Particularly, both the ground moisture content and the thermal inertia of the model were severely underestimated. At the same time, the sea surface temperatures produced by the model between the equator and the model's land-sea boundary were far too cool. Both the atmosphere and the ocean model were modified to include a better hydrological cycle and ocean structure. The convective events produced by the modified model possessed the observed frequency and were generated well south of the coastline. The improved simulation of monsoon variability allowed the hydrological cycle feedback to be generalized. It was found that monsoon variability was constrained to lie within the bounds of a positive gradient of a convective intensity potential (I). The function depends primarily on the surface temperature, the availability of moisture and the stability of the lower atmosphere which varies very slowly on the time scale of months. The oscillations of the monsoon perturb the mean convective intensity potential causing local enhancements of the gradient. These perturbations are caused by the hydrological feedbacks, discussed above, or by the modification of the air-sea fluxes caused by variations of the low level wind during convective events. The final result is the slow northward propagation of convection within an even slower convective regime. The ECMWF analyses show very similar behavior of the convective intensity potential. Although it is considered premature to use the model to conduct simulations of the African monsoon system, the ECMWF analysis indicates similar behavior in the convective intensity potential suggesting, at least, that the same processes control the low frequency structure of the African monsoon. The implications of the hypotheses on numerical weather prediction of monsoon phenomenon are discussed.

Barriers to plantation activities in different agro-ecological zones of Southern India

This paper analyses environmental and socio-economic barriers for plantation activities on local and regional level and investigates the potential for carbon finance to stimulate the increased rates of forest plantation on wasteland, i.e., degraded lands, in southern India. Building on multidisciplinary field work and results from the model GCOMAP, the aim is to (1) identify and characterize the barriers to plantation activities in four agro-ecological zones in the state of Karnataka and (2) investigate what would be required to overcome these barriers and enhance the plantation rate and productivity. The results show that a rehabilitation of the wasteland based on plantation activities is not only possible but also anticipated by the local population and would lead to positive environmental and socio-economic effects at a local level. However, in many cases, the establishment of plantation activities is hindered by a lack of financial resources, low land productivity and water scarcity. Based on the model used and the results from the field work, it can be concluded that certified emission reductions such as carbon credits or other compensatory systems may help to overcome the financial barrier; however, the price needs to be significantly increased if these measures are to have any large-scale impact.

Aerosol radiative forcing due to enhanced black carbon at an urban site in India

[1] During a comprehensive aerosol field campaign, simultaneous measurements were made of aerosol spectral optical depths, black carbon mass concentration (M-b), total (M-t) and size segregated aerosol mass concentrations over an urban continental location, Bangalore (13 degreesN, 77 degreesE, 960 m msl), in India. Large amounts of BC were observed; both in absolute terms and fraction of total mass (similar to11%) and submicron mass (similar to23%) implying a significantly low single scatter albedo. The aerosol visible optical depth (tau(p)) was in the range 0.24 to 0.45. Estimated surface forcing is as high as -23 W m(-2) and top of the atmosphere (TOA) forcing is +5 Wm(-2) during relatively cleaner periods (tau(p) similar to 0.24). The net atmospheric absorption translates to an atmospheric heating of similar to0.8 K day(-1) for cleaner periods and similar to1.5 K day(-1) for less cleaner periods (tau(p) similar to 0.45). Our observations raise several issues, which may have impacts to regional climate and monsoon.

Aerosol radiative forcing over land: effect of surface and cloud reflection

It is now clearly understood that atmospheric aerosols have a significant impact on climate due to their important role in modifying the incoming solar and outgoing infrared radiation. The question of whether aerosol cools (negative forcing) or warms (positive forcing) the planet depends on the relative dominance of absorbing aerosols. Recent investigations over the tropical Indian Ocean have shown that, irrespective of the comparatively small percentage contribution in optical depth (similar to11%), soot has an important role in the overall radiative forcing. However, when the amount of absorbing aerosols such as soot are significant, aerosol optical depth and chemical composition are not the only determinants of aerosol climate effects, but the altitude of the aerosol layer and the altitude and type of clouds are also important. In this paper, the aerosol forcing in the presence of clouds and the effect of different surface types (ocean, soil, vegetation, and different combinations of soil and vegetation) are examined based on model simulations, demonstrating that aerosol forcing changes sign from negative (cooling) to positive (warming) when reflection from below (either due to land or clouds) is high.