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.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Frequency Domain Based Solution for Certain Class of Wave Equations: An exhaustive study of Numerical Solutions

The paper discusses the frequency domain based solution for a certain class of wave equations such as: a second order partial differential equation in one variable with constant and varying coefficients (Cantilever beam) and a coupled second order partial differential equation in two variables with constant and varying coefficients (Timoshenko beam). The exact solution of the Cantilever beam with uniform and varying cross-section and the Timoshenko beam with uniform cross-section is available. However, the exact solution for Timoshenko beam with varying cross-section is not available. Laplace spectral methods are used to solve these problems exactly in frequency domain. The numerical solution in frequency domain is done by discretisation in space by approximating the unknown function using spectral functions like Chebyshev polynomials, Legendre polynomials and also Normal polynomials. Different numerical methods such as Galerkin Method, Petrov- Galerkin method, Method of moments and Collocation method or the Pseudo-spectral method in frequency domain are studied and compared with the available exact solution. An approximate solution is also obtained for the Timoshenko beam with varying cross-section using Laplace Spectral Element Method (LSEM). The group speeds are computed exactly for the Cantilever beam and Timoshenko beam with uniform cross-section and is compared with the group speeds obtained numerically. The shear mode and the bending modes of the Timoshenko beam with uniform cross-section are separated numerically by applying a modulated pulse as the shear force and the corresponding group speeds for varying taper parameter in are obtained numerically by varying the frequency of the input pulse. An approximate expression for calculating group speeds corresponding to the shear mode and the bending mode, and also the cut-off frequency is obtained. Finally, we show that the cut-off frequency disappears for large in, for epsilon > 0 and increases for large in, for epsilon < 0.

A mathematical analysis of nonlinear waves in a fluid filled visco-elastic tube

Our investigations in this paper are centred around the mathematical analysis of a ldquomodal waverdquo problem. We have considered the axisymmetric flow of an inviscid liquid in a thinwalled viscoelastic tube under certain simplifying assumptions. We have first derived the propagation space equations in the long wave limit and also given a general procedure to derive these equations for arbitrary wave length, when the flow is irrotational. We have used the method of operators of multiple scales to derive the nonlinear Schrödinger equation governing the modulation of periodic waves and we have elaborated on the ldquolong modulated wavesrdquo and the ldquomodulated long wavesrdquo. We have also examined the existence and stability of Stokes waves in this system. This is followed by a discussion of the progressive wave solutions of the long wave equations. One of the most important results of our paper is that the propagation space equations are no longer partial differential equations but they are in terms of pseudo-differential operators.Die vorliegenden Untersuchungen beziehen sich auf die mathematische Behandlung des ldquorModalwellenrdquo-Problems. Die achsensymmetrische Strömung einer nichtviskosen Flüssigkeit in einem dünnwandigen viskoelastischen Rohr, unter bestimmten vereinfachenden Annahmen, wird betrachtet. Zuerst werden die Gleichungen des Ausbreitungsraumes im Langwellenbereich abgeleitet und eine allgemeine Methode zur Herleitung dieser Gleichungen für beliebige Wellenlängen bei nichtrotierender Strömung angegeben. Eine Operatorenmethode mit multiplem Maßstab wird verwendet zur Herleitung der nichtlinearen Schrödinger-Gleichung für die Modulation der periodischen Wellen, und die ldquorlangmodulierten Wellenrdquo sowie die ldquormodulierten Langwellenrdquo werden aufgezeigt. Weiters wird die Existenz und die Stabilität der Stokes-Wellen im System untersucht. Anschließend werden die progressiven Wellenlösungen der Langwellengleichungen diskutiert. Eines der wichtigsten Ergebnisse dieser Arbeit ist, daß die Gleichungen des Ausbreitungsraumes keine partiellen Differentialgleichungen mehr sind, sondern Ausdrücke von Pseudo-Differentialoperatoren.

Theory of Non-linear Waves in Multi-dimensions: with Special Reference to Surface Water Waves

The surface water waves are "modal" waves in which the "physical space" (t, x, y, z) is the product of a propagation space (t, x, y) and a cross space, the z-axis in the vertical direction. We have derived a new set of equations for the long waves in shallow water in the propagation space. When the ratio of the amplitude of the disturbance to the depth of the water is small, these equations reduce to the equations derived by Whitham (1967) by the variational principle. Then we have derived a single equation in (t, x, y)-space which is a generalization of the fourth order Boussinesq equation for one-dimensional waves. In the neighbourhood of a wave froat, this equation reduces to the multidimensional generalization of the KdV equation derived by Shen & Keller (1973). We have also included a systematic discussion of the orders of the various non-dimensional parameters. This is followed by a presentation of a general theory of approximating a system of quasi-linear equations following one of the modes. When we apply this general method to the surface water wave equations in the propagation space, we get the Shen-Keller equation.

Wave dispersion in a cellular composite with modulated microstructure

In this paper we report a modeling technique and analysis of wave dispersion in a cellular composite laminate with spatially modulated microstructure, which can be modeled by parameterization and homogenization in an appropriate length scale. Higher order beam theory is applied and the system of wave equations are derived. Homogenization of these equations are carried out in the scale of wavelength and frequency of the individual wave modes. Smaller scale scattering below the order of cell size are filtered out in the present approach. The longitudinal dispersion relations for different values of a modulation parameter are analyzed which indicates the existence of stop and pass band patterns. Dispersion relations for flexural-shear case are also analyzed which indicates a tendency toward forming the stop and pass bands for increasing values of a shear stiffness modulation parameter. The effect the phase angle (θ) of the incident wave indicates the existence more number of alternative stop bands and pass bands for θ = 45°.

Dynamics of a dilute sheared inelastic fluid. I. Hydrodynamic modes and velocity correlation functions

The hydrodynamic modes and the velocity autocorrelation functions for a dilute sheared inelastic fluid are analyzed using an expansion in the parameter epsilon=(1-e)(1/2), where e is the coefficient of restitution. It is shown that the hydrodynamic modes for a sheared inelastic fluid are very different from those for an elastic fluid in the long-wave limit, since energy is not a conserved variable when the wavelength of perturbations is larger than the ``conduction length.'' In an inelastic fluid under shear, there are three coupled modes, the mass and the momenta in the plane of shear, which have a decay rate proportional to k(2/3) in the limit k -> 0, if the wave vector has a component along the flow direction. When the wave vector is aligned along the gradient-vorticity plane, we find that the scaling of the growth rate is similar to that for an elastic fluid. The Fourier transforms of the velocity autocorrelation functions are calculated for a steady shear flow correct to leading order in an expansion in epsilon. The time dependence of the autocorrelation function in the long-time limit is obtained by estimating the integral of the Fourier transform over wave number space. It is found that the autocorrelation functions for the velocity in the flow and gradient directions decay proportional to t(-5/2) in two dimensions and t(-15/4) in three dimensions. In the vorticity direction, the decay of the autocorrelation function is proportional to t(-3) in two dimensions and t(-7/2) in three dimensions.

A fuzzy dynamic flood routing model for natural channels

A fuzzy dynamic flood routing model (FDFRM) for natural channels is presented, wherein the flood wave can be approximated to a monoclinal wave. This study is based on modification of an earlier published work by the same authors, where the nature of the wave was of gravity type. Momentum equation of the dynamic wave model is replaced by a fuzzy rule based model, while retaining the continuity equation in its complete form. Hence, the FDFRM gets rid of the assumptions associated with the momentum equation. Also, it overcomes the necessity of calculating friction slope (S-f) in flood routing and hence the associated uncertainties are eliminated. The fuzzy rule based model is developed on an equation for wave velocity, which is obtained in terms of discontinuities in the gradient of flow parameters. The channel reach is divided into a number of approximately uniform sub-reaches. Training set required for development of the fuzzy rule based model for each sub-reach is obtained from discharge-area relationship at its mean section. For highly heterogeneous sub-reaches, optimized fuzzy rule based models are obtained by means of a neuro-fuzzy algorithm. For demonstration, the FDFRM is applied to flood routing problems in a fictitious channel with single uniform reach, in a fictitious channel with two uniform sub-reaches and also in a natural channel with a number of approximately uniform sub-reaches. It is observed that in cases of the fictitious channels, the FDFRM outputs match well with those of an implicit numerical model (INM), which solves the dynamic wave equations using an implicit numerical scheme. For the natural channel, the FDFRM Outputs are comparable to those of the HEC-RAS model.

Equatorial modes observed in atmospheric variables

Wavenumber-frequency spectral analysis of different atmospheric variables has been carried Out using 25 years of data. The area considered is the tropical belt 25 degrees S-25 degrees N. A combined FFT wavelet analysis method has been used for this purpose. Variables considered are outgoing long wave radiation (OLR), 850 hPa divergence, zonal and meridional winds at 850, 500 and 200 hPa levels, sea level pressure and 850 hPa geopotential height. It is shown that the spectra of different variables have some common properties, but each variable also has few features diffe:rent from the rest. While Kelvin mode is prominent in OLR, and zonal winds, it is not clearly observed in pressure and geopotential height fields; the latter two have a dominant wavenumber zero mode not seen in other variables except in meridional wind at 200 hPa and 850 hPa divergences. Different dominant modes in the tropics show significant variations on sub-seasonal time scales.

Vave propagation and bandgaps in a parametrically modulated composite laminate

Wave propagation and its frequency bandgaps in a parametrically modulated composite laminate are reported in this paper. The modulated properties under considerations are due to periodic microstructure, for example honeycomb core sandwich composite, which can be parameterized and homogenized in a suitable scale. Wave equations are derived by assuming a third-order shear deformation theory. Homogenization of the wave equations is carried out in the scale of wavelength. In-plane wave and flexural-shear wave dispersions are obtained for a range of values of a stiffness modulation coefficient (alpha). A clear pattern of stop-bands is observed for alpha >= 4. To validate the band-gap phenomena, we take recourse to time domain response obtained from finite element simulation. As predicted by the proposed analytical technique, a distinct correlation between the chosen frequency band and the simulated wave arrival time and amplitude reduction is found. This promises practical applications of the proposed analytical technique to designing parametrically modulated composite laminate for wave suppression. (C) 2009 Elsevier B.V. All rights reserved.

The vibration spectrum of rutile

The long-wave lattice dynamics of rutile has been studied using a rigid ion model. The vibration frequencies for the zero wavevector have been calculated using the expressions for the frequencies of the normal modes derived group theoretically. The observed Raman and infrared frequencies have been explained.

Raman spectra of alkali halides

In this article, we present a comparative study of the Raman spectra of alkali halides in relation to the lattice dynamics ofBorn andRaman. It is shown that the experimentally observed limit of the second-order spectra in almost all the cases can be explained well by the Lyddane-Sachs-Teller relation. It is also seen, while, an explanation of the second-order Raman spectrum of a crystal of diamond or zinc blende structure requires the frequencies from the critical points,W, Gamma, X andL inBorn's analysis, the frequencies fromGamma, X andL alone are sufficient and necessary for an interpretation of the same onRaman's model. Some similarities in the determination of the long wave properties of crystals like elastic constants and limiting frequencies of the lattice vibrations in the symmetry directions in both the models are pointed out.

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).

Diagnosis of dominant forcing factors for large-scale vertical velocities during active and break phases of the monsoon

An attempt to diagnose the dominant forcings which drive the large-scale vertical velocities over the monsoon region has been made by computing the forcings like diabatic heating fields,etc. and the large-scale vertical velocities driven by these forcings for the contrasting periods of active and break monsoon situations; in order to understand the rainfall variability associated with them. Computation of diabatic heating fields show us that among different components of diabatic heating it is the convective heating that dominates at mid-tropospheric levels during an active monsoon period; whereas it is the sensible heating at the surface that is important during a break period. From vertical velocity calculations we infer that the prime differences in the large-scale vertical velocities seen throughout the depth of the atmosphere are due to the differences in the orders of convective heating; the maximum rate of latent heating being more than 10 degrees Kelvin per day during an active monsoon period; whereas during a break monsoon period it is of the order of 2 degrees Kelvin per day at mid-tropospheric levels. At low levels of the atmosphere, computations show that there is large-scale ascent occurring over a large spatial region, driven only by the dynamic forcing associated with vorticity and temperature advection during an active monsoon period. However, during a break monsoon period such large-scale spatial organization in rising motion is not seen. It is speculated that these differences in the low-level large-scale ascent might be causing differences in convective heating because the weaker the low level ascent, the lesser the convective instability which produces deep cumulus clouds and hence lesser the associated latent heat release. The forcings due to other components of diabatic heating, namely, the sensible heating and long wave radiative cooling do not influence the large-scale vertical velocities significantly.

Mechanisms of variability and predictability of the tropical coupled ocean-atmosphere system

A conceptual model is proposed to explain the observed aperiodicity in the short term climate fluctuations of the tropical coupled ocean-atmosphere system. This is based on the evidence presented here that the tropical coupled ocean-atmosphere system sustains a low frequency inter-annual mode and a host of higher frequency intra-seasonal unstable modes. At long wavelengths, the low frequency mode is dominant while at short wavelengths, the high frequency modes are dominant resulting in the co-existence of a long wave low frequency mode with some short wave intra-seasonal modes in the tropical coupled system. It is argued that due to its long wavelength, the low frequency mode would behave like a linear oscillator while the higher frequency short wave modes would be nonlinear. The conceptual model envisages that an interaction between the low frequency linear oscillator and the high frequency nonlinear oscillations results in the observed aperiodicity of the tropical coupled system. This is illustrated by representing the higher frequency intra-seasonal oscillations by a nonlinear low order model which is then coupled to a linear oscillator with a periodicity of four years. The physical mechanism resulting in the aperiodicity in the low frequency oscillations and implications of these results on the predictability of the coupled system are discussed.