Modelling the chemical weathering fluxes at the watershed scale in the Tropics (Mule Hole, South India): Relative contribution of the smectite/kaolinite assemblage versus primary minerals

We investigate the chemical weathering processes and fluxes in a small experimental watershed (SEW) through a modelling approach. The study site is the Mule Hole SEW developed on a gneissic basement located in the climatic gradient of the Western Ghats, South India. The model couples a lumped hydrological model simulating the water budget at the watershed scale to the WITCH model estimating the dissolution/precipitation rates of minerals using laboratory kinetic laws. Forcing functions and parameters of the simulation are defined by the field data. The coupled model is calibrated with stream and groundwater compositions through the testing of a large range of smectite solubility and abundance in the soil horizons. We found that, despite the low abundance of smectite in the dominant soil type of the watershed (4 vol.%), their net dissolution provides 75% of the export of dissolved silica, while primary silicate mineral dissolution releases only 15% of this flux. Overall, smectites (modelled as montmorillonites) are not stable under the present day climatic conditions. Furthermore, the dissolution of trace carbonates in the saprolitic horizon provides 50% of the calcium export at the watershed scale. Modelling results show the contrasted behavior of the two main soil types of the watershed: red soils (88% of the surface) are provider of calcium, while black soils (smectite-rich and characterized by a lower drainage) consumes calcium through overall carbonate precipitation. Our model results stress the key role played by minor/accessory minerals and by the thermodynamic properties of smectite minerals, and by the drainage of the weathering profiles on the weathering budget of a tropical watershed. (C) 2010 Elsevier B.V. All rights reserved.

Sensitivity of the simulated precipitation to changes in convective relaxation time scale

The paper describes the sensitivity of the simulated precipitation to changes in convective relaxation time scale (TAU) of Zhang and McFarlane (ZM) cumulus parameterization, in NCAR-Community Atmosphere Model version 3 (CAM3). In the default configuration of the model, the prescribed value of TAU, a characteristic time scale with which convective available potential energy (CAPE) is removed at an exponential rate by convection, is assumed to be 1 h. However, some recent observational findings suggest that, it is larger by around one order of magnitude. In order to explore the sensitivity of the model simulation to TAU, two model frameworks have been used, namely, aqua-planet and actual-planet configurations. Numerical integrations have been carried out by using different values of TAU, and its effect on simulated precipitation has been analyzed. The aqua-planet simulations reveal that when TAU increases, rate of deep convective precipitation (DCP) decreases and this leads to an accumulation of convective instability in the atmosphere. Consequently, the moisture content in the lower-and mid-troposphere increases. On the other hand, the shallow convective precipitation (SCP) and large-scale precipitation (LSP) intensify, predominantly the SCP, and thus capping the accumulation of convective instability in the atmosphere. The total precipitation (TP) remains approximately constant, but the proportion of the three components changes significantly, which in turn alters the vertical distribution of total precipitation production. The vertical structure of moist heating changes from a vertically extended profile to a bottom heavy profile, with the increase of TAU. Altitude of the maximum vertical velocity shifts from upper troposphere to lower troposphere. Similar response was seen in the actual-planet simulations. With an increase in TAU from 1 h to 8 h, there was a significant improvement in the simulation of the seasonal mean precipitation. The fraction of deep convective precipitation was in much better agreement with satellite observations.

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (epsilon) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.

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.

Piezoelectrically actuated insect scale flapping wing

An energy method is used in order to derive the non-linear equations of motion of a smart flapping wing. Flapping wing is actuated from the root by a PZT unimorph in the piezofan configuration. Dynamic characteristics of the wing, having the same size as dragonfly Aeshna Multicolor, are analyzed using numerical simulations. It is shown that flapping angle variations of the smart flapping wing are similar to the actual dragonfly wing for a specific feasible voltage. An unsteady aerodynamic model based on modified strip theory is used to obtain the aerodynamic forces. It is found that the smart wing generates sufficient lift to support its own weight and carry a small payload. It is therefore a potential candidate for flapping wing of micro air vehicles.

Strong nonlocalization induced by small scale parameter on terahertz flexural wave dispersion characteristics of a monolayer graphene

This paper presents the strong nonlocal scale effect on the flexural wave propagation in a monolayer graphene sheet. The graphene is modeled as an isotropic plate of one atom thick. Nonlocal governing equation of motion is derived and wave propagation analysis is performed using spectral analysis. The present analysis shows that the flexural wave dispersion in graphene obtained by local and nonlocal elasticity theories is quite different. The nonlocal elasticity calculation shows that the wavenumber escapes to infinite at certain frequency and the corresponding wave velocity tends to zero at that frequency indicating localization and stationary behavior. This behavior is captured in the spectrum and dispersion curves. The cut-off frequency of flexural wave not only depend on the axial wavenumber but also on the nonlocal scaling parameter. The effect of axial wavenumber on the wave behavior in graphene is also discussed in the present manuscript. (C) 2010 Elsevier B.V. All rights reserved.

Length-scale-dependent ensemble-averaged conductance of a 1D disordered conductor: Conductance minimum

An exact numerical calculation of ensemble-averaged length-scale-dependent conductance for the one-dimensional Anderson model is shown to support an earlier conjecture for a conductance minimum. The numerical results can be understood in terms of the Thouless expression for the conductance and the Wigner level-spacing statistics.

On the out of phase appearance of large-scale diffuse magnetic field of the sun with respect to sunspots

We assume the large-scale diffuse magnetic field of the Sun to originate from the poloidal component of a dynamo operating at the base of the convection zone, whereas the sunspots are due to the toroidal component. The evolution of the poloidal component is studied to model the poleward migration of the diffuse field seen on the solar surface and the polar reversal at the time of sunspot maxima (Dikpati and Choudhuri 1994, 1995).

On the Large-Scale Diffuse Magnetic Field of the Sun

Although the sunspots migrate towards the equator, the large-scale weak diffuse magnetic fields of the Sun migrate poleward with the solar cycle, the polar field reversing at the time of the sunspot maxima. We apply the vector model of Dikpati and Choudhuri (1994, Paper I) to fit these observations. The dynamo layer at the base of the convection zone is taken to be the source of the diffuse field, which is then evolved in the convection zone subject to meridional circulation and turbulent diffusion. We find that the longitudinally averaged observational data can be fitted reasonably well both for positive and negative values of the alpha-effect by adjusting the subsurface meridional flow suitably. The model will be extended in a future paper to include the decay of active regions as an extra source of the diffuse field, which may be necessary to explain the probable phase lag between B-tau and B-phi at lower latitudes.

The actor-critic algorithm as multi-time-scale stochastic approximation

The actor-critic algorithm of Barto and others for simulation-based optimization of Markov decision processes is cast as a two time Scale stochastic approximation. Convergence analysis, approximation issues and an example are studied.

A Dual Scale Model for Macrosegregation in Alloy Solidification

A numerical approach for coupling the temperature and concentration fields using a micro/macro dual scale model for a solidification problem is presented. The dual scale modeling framework is implemented on a hybrid explicit-implicit solidification scheme. The advantage of this model lies in more accurate consideration of microsegregation occurring at micro-scale using a subgrid model. The model is applied to the case of solidification of a Pb-40% Sn alloy in a rectangular cavity. The present simulation results are compared with the corresponding experimental results reported in the literature, showing improvement in macrosegregation predictions. Subsequently, a comparison of macrosegregation prediction between the results of the present method with those of a parameter model is performed, showing similar trends.

On the large-scale diffuse magnetic field of the sun - II. The contribution of active regions

Dikpati and Choudhuri (1993, 1995) developed a model for the poleward migration of the weak diffuse magnetic field on the Sun's surface. This field was identified with the poloidal component produced by the solar dynamo operating at the base of the convection zone, and its evolution was studied by considering the effects of meridional circulation and turbulent diffusion. The earlier model is extended in this paper by incorporating the flux from, the decay of tilted active regions near the solar surface as an additional source of the poloidal field. This extended model can now explain various low-latitude features in the time-latitude diagram of the weak diffuse fields. These low-latitude features could not be accounted for in the earlier model, which was very successful in modeling the behavior at high latitudes. The time-latitude diagrams show that regions of a particular polarity often have 'tongues' of opposite polarity. Such tongues can be produced in the theoretical model by incorporating fluctuations in the source term arising out of the decaying active regions.

Comparison of the ultrafast to slow time scale dynamics of three liquid crystals in the isotropic phase

The dynamics of three liquid crystals, 4'(pentyloxy)-4-biphenylcarbonitrile (5-OCB), 4'-pentyl-4-biphenylcarbonitrile (5-CB), and 1-isothiocyanato-(4-propylcyclohexyl)benzene (3-CHBT), are investigated from very short time (similar to1 ps) to very long time (>100 ns) as a function of temperature using optical heterodyne detected optical Kerr effect experiments. For all three liquid crystals, the data decay exponentially only on the longest time scale (> several ns). The temperature dependence of the long time scale exponential decays is described well by the Landau-de Gennes theory of the randomization of pseudonematic domains that exist in the isotropic phase of liquid crystals near the isotropic to nematic phase transition. At short time, all three liquid crystals display power law decays. Over the full range of times, the data for all three liquid crystals are fit with a model function that contains a short time power law. The power law exponents for the three liquid crystals range between 0.63 and 0.76, and the power law exponents are temperature independent over a wide range of temperatures. Integration of the fitting function gives the empirical polarizability-polarizability (orientational) correlation function. A preliminary theoretical treatment of collective motions yields a correlation function that indicates that the data can decay as a power law at short times. The power law component of the decay reflects intradomain dynamics. (C) 2002 American Institute of Physics.

Mapping of large scale 158 mu m [CII] line emission: Orion A

We present the first results of an observational programme undertaken to map the fine structure line emission of singly ionized carbon ([ CII] 157 : 7409 mum) over extended regions using a Fabry Perot spectrometer newly installed at the focal plane of a 100 cm balloon- borne far- infrared telescope. This new combination of instruments has a velocity resolution of similar to 200 km s(-1) and an angular resolution of 1.'5. During the first flight, an area of 30' x 15' in Orion A was mapped. These observations extend over a larger area than previous observations, the map is fully sampled and the spectral scanning method used enables reliable estimation of the continuum emission at frequencies adjacent to the [ CII] line. The total [ CII] line luminosity, calculated by considering up to 20% of the maximum line intensity is 0.04% of the luminosity of the far- infrared continuum. We have compared the [ CII] intensity distribution with the velocity- integrated intensity distributions of (CO)-C-13(1- 0), CI(1- 0) and CO( 3- 2) from the literature. Comparison of the [ CII], [ CI] and the radio continuum intensity distributions indicates that the largescale [ CII] emission originates mainly from the neutral gas, except at the position of M 43, where no [ CI] emission corresponding to the [ CII] emission is seen. Substantial part of the [ CII] emission from here originates from the ionized gas. The observed line intensities and ratios have been analyzed using the PDR models by Kaufman et al. ( 1999) to derive the incident UV flux and volume density at a few selected positions. The models reproduce the observations reasonably well at most positions excepting the [ CII] peak ( which coincides with the position of theta(1) Ori C). Possible reason for the failure could be the simplifying assumption of a homogeneous plane parallel slab in place of a more complicated geometry.