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.

Trellis coded modulation scheme for the fading relay channel

A decode and forward protocol based Trellis Coded Modulation (TCM) scheme for the half-duplex relay channel, in a Rayleigh fading environment, is presented. The proposed scheme can achieve any spectral efficiency greater than or equal to one bit per channel use (bpcu). A near-ML decoder for the suggested TCM scheme is proposed. It is shown that the high Signal to Noise Ratio (SNR) performance of this near-ML decoder approaches the performance of the optimal ML decoder. Based on the derived Pair-wise Error Probability (PEP) bounds, design criteria to maximize the diversity and coding gains are obtained. Simulation results show a large gain in SNR for the proposed TCM scheme over uncoded communication as well as the direct transmission without the relay.

Nonlinear dielectric response of BaBi4Ti4O15 ceramics

Nonlinear dielectric response of BaBi4Ti4O15 ceramics synthesized via the conventional solid-state reaction route has been monitored over a wide range of electric field strengths (E-0 = 0.5 - 5 kV/cm). Dielectric permittivity was found to increase linearly within the range of applied field. Rayleigh relations were employed to interpret the nonlinear dielectric response and the contribution of irreversible domain wall motion to the macroscopic permittivity was separated. The values of room temperature Rayleigh dielectric coefficient (alpha) and relative initial permittivity (epsilon'(init)) were found to be 2.28 +/- 0.02 cm/kV and 146.10 +/- 0.07, respectively. A reasonable agreement between the simulated and measured polarization-electric field (P-E) hysteresis loops was observed at an applied electric field of 5 kV/cm.

Power allocation for interference channels

We propose power allocation algorithms for increasing the sum rate of two and three user interference channels. The channels experience fast fading and there is an average power constraint on each transmitter. Our achievable strategies for two and three user interference channels are based on the classification of the interference into very strong, strong and weak interferences. We present numerical results of the power allocation algorithm for two user Gaussian interference channel with Rician fading with mean total power gain of the fade Omega = 3 and Rician factor kappa = 0.5 and compare the sum rate with that obtained from ergodic interference alignment with water-filling. We show that our power allocation algorithm increases the sum rate with a gain of 1.66dB at average transmit SNR of 5dB. For the three user Gaussian interference channel with Rayleigh fading with distribution CN(0, 0.5), we show that our power allocation algorithm improves the sum rate with a gain of 1.5dB at average transmit SNR of 5dB.

Can Superflares Occur on Our Sun?

Recent observations of Sun-like stars, similar to our Sun in their surface temperature (5600-6000 K) and slow rotation (rotational period > 10 d), using the Kepler satellite by Maehara et al. (2012, Nature, 485, 478) have revealed the existence of superflares (with energy of 10(33)-10(35) erg). From statistical analyses of these superflares, it was found that superflares with energy of 10(34) erg occur once in 800 yr, and superflares with 10(35) erg occur once in 5000 yr. In this paper, we examine whether superflares with energy of 10(33)-10(35) erg could occur on the present Sun through the use of simple order-of-magnitude estimates based on current ideas related to the mechanisms of the solar dynamo. If magnetic flux is generated by differential rotation at the base of the convection zone, as assumed in typical dynamo models, it is possible that the present Sun would generate a large sunspot with a total magnetic flux of similar to 2 x 10(23) Mx (= G cm(2)) within one solar cycle period, and lead to superflares with an energy of 10(34) erg. To store a total magnetic flux of similar to 10(24) Mx, necessary for generating 10(35) erg superflares, it would take similar to 40 yr. Hot Jupiters have often been argued to be a necessary ingredient for the generation of superflares, but we found that they do not play any essential role in the generation of magnetic flux in the star itself, if we consider only the magnetic interaction between the star and the hot Jupiter. This seems to be consistent with Maehara et al.'s finding of 148 superflare-generating solar-type stars that do not have a hot Jupiter-like companion. Altogether, our simple calculations, combined with Maehara et al.'s analysis of superflares on Sun-like stars, show that there is a possibility that superflares of 10(34) erg would occur once in 800 yr on our present Sun.

Universal flame propagation behavior in packed bed of biomass

This article aims at seeking the universal behavior of propagation rate variation with air superficial velocity (V-s) in a packed bed of a range of biomass particles in reverse downdraft mode while also resolving the differing and conflicting explanations in the literature. Toward this, measurements are made of exit gas composition, gas phase and condensed phase surface temperature (T-g and T-s), and reaction zone thickness for a number of biomass with a range of properties. Based on these data, two regimes are identified: gasificationvolatile oxidation accompanied by char reduction reactions up to 16 +/- 1cm/s of V-s and above this, and char oxidationsimultaneous char oxidation and gas phase combustion. In the gasification regime, the measured T-s is less than T-g; a surface heat balance incorporating a diffusion controlled model for flaming combustion gives and matches with the experimental results to within 5%. In the char oxidation regime, T-g and T-s are nearly equal and match with the equilibrium temperature at that equivalence ratio. Drawing from a recent study of the authors, the ash layer over the oxidizing char particle is shown to play a critical role in regulating the radiation heat transfer to fresh biomass in this regime and is shown to be crucial in explaining the observed propagation behavior. A simple model based on radiation-convection balance that tracks the temperature-time evolution of a fresh biomass particle is shown to support the universal behavior of the experimental data on reaction front propagation rate from earlier literature and the present work for biomass with ash content up to 10% and moisture fraction up to 10%. Upstream radiant heat transfer from the ash-laden hot char modulated by the air flow is shown to be the dominant feature of this model.

Dynamics of the North American continent

The forces that cause deformation of western North America have been debated for decades. Recent studies, primarily based on analysis of crustal stresses in the western United States, have suggested that the deformation of the region is mainly controlled by gravitational potential energy (GPE) variations and boundary loads, with basal tractions due to mantle flow playing a relatively minor role. We address these issues by modelling the deviatoric stress field over western North America from a 3-D finite element mantle circulation model with lateral viscosity variations. Our approach takes into account the contribution from both topography and shallow lithosphere structure (GPE) as well as that from deeper mantle flow in one single model, as opposed to separate lithosphere and circulation models, as has been done so far. In addition to predicting the deviatoric stresses we also jointly fit the constraints of geoid, dynamic topography and plate motion both globally and over North America, in order to ensure that the forces that arise in our models are dynamically consistent. We examine the sensitivity of the dynamic models to different lateral viscosity variations. We find that circulation models that include upper mantle slabs yield a better fit to observed plate velocities. Our results indicate that a model of GPE variations coupled with mantle convection gives the best fit to the observational constraints. We argue that although GPE variations control a large part of the deformation of the western United States, deeper mantle tractions also play a significant role. The average deviatoric stress magnitudes in the western United States range 30-40 MPa. The cratonic region exhibits higher coupling to mantle flow than the rest of the continent. We find that a relatively strong San Andreas fault gives a better fit to the observational constraints, especially that of plate velocity in western North America.

Magnetohydrodynamic stability of stochastically driven accretion flows

We investigate the evolution of magnetohydrodynamic (or hydromagnetic as coined by Chandrasekhar) perturbations in the presence of stochastic noise in rotating shear flows. The particular emphasis is the flows whose angular velocity decreases but specific angular momentum increases with increasing radial coordinate. Such flows, however, are Rayleigh stable but must be turbulent in order to explain astrophysical observed data and, hence, reveal a mismatch between the linear theory and observations and experiments. The mismatch seems to have been resolved, at least in certain regimes, in the presence of a weak magnetic field, revealing magnetorotational instability. The present work explores the effects of stochastic noise on such magnetohydrodynamic flows, in order to resolve the above mismatch generically for the hot flows. We essentially concentrate on a small section of such a flow which is nothing but a plane shear flow supplemented by the Coriolis effect, mimicking a small section of an astrophysical accretion disk around a compact object. It is found that such stochastically driven flows exhibit large temporal and spatial autocorrelations and cross-correlations of perturbation and, hence, large energy dissipations of perturbation, which generate instability. Interestingly, autocorrelations and cross-correlations appear independent of background angular velocity profiles, which are Rayleigh stable, indicating their universality. This work initiates our attempt to understand the evolution of three-dimensional hydromagnetic perturbations in rotating shear flows in the presence of stochastic noise.

Wet and dry spell characteristics of global tropical rainfall

In this study, we analyze satellite-based daily rainfall observations to compare and contrast the wet and dry spell characteristics of tropical rainfall. Defining a wet (dry) spell as the number of consecutive rainy (nonrainy) days, we find that the distributions of wet spells appear to exhibit universality in the following sense. While both ocean and land regions with high seasonal rainfall accumulation (humid regions; e. g., India, Amazon, Pacific Ocean) show a predominance of 2-4 day wet spells, those regions with low seasonal rainfall accumulation (arid regions; e. g., South Atlantic, South Australia) exhibit a wet spell duration distribution that is essentially exponential in nature, with a peak at 1 day. The behavior that we observed for wet spells is reversed for the dry spell characteristics. In other words, the main contribution to the dry part of the season, in terms of the number of nonrainy days, appears to come from 3-4 day dry spells in the arid regions, as opposed to 1 day dry spells in the humid regions. The total rainfall accumulated in each wet spell has also been analyzed, and we find that the major contribution to seasonal rainfall for arid regions comes from 1-5 day wet spells; however, for humid regions, this contribution comes from wet spells of duration as long as 30 days. We also explore the role of chance as well as the influence of organized convection in determining some of the observed features.

Wireless Network Coding for MIMO Two-Way Relaying

The design of modulation schemes for the physical layer network-coded two-way MIMO relaying scenario is considered, with the denoise-and-forward protocol which employs two phases: Multiple Access phase and Broadcast phase. It is shown that for MIMO two-way relaying, the minimum distance of the effective constellation at the relay becomes zero when all the rows of the channel fade coefficient matrix belong to a finite number of vector subspaces referred to as the singular fade subspaces. The singular fade subspaces can be classified into two kinds based on whether their harmful effects can be removed or not: (i) the removable and (ii) the non-removable singular fade subspaces. It is shown that network coding maps obtained by the completion of appropriate partially filled Latin Rectangles can remove the harmful effects of all the removable singular fade subspaces. For 2(lambda)-PSK signal set, the removable and non-removable singular fade subspaces are characterized and, it is shown that the number of non-removable singular fade subspaces is a small fraction of the total number of singular fade subspaces and this fraction tends to zero as the constellation size tends to infinity. The Latin Rectangles for the case when the end nodes use different number of antennas are shown to be obtainable from the Latin Squares for the case when they use the same number of antennas. Also, the network coding maps which remove all the removable singular singular fade subspaces are shown to be obtainable from a small set of Latin Squares. The removal of all the singular fade subspaces by properly choosing the network coding map, provides a gain of 5.5 dB over the conventional Exclusive-OR network coding, in a Rayleigh fading scenario with 2 antennas at the end nodes and one antenna at the relay node, for 4-PSK signal set.

Enhanced DMT-optimality criterion for STBC schemes for asymmetric MIMO systems

For any n(t) transmit, n(r) receive antenna (n(t) x n(r)) multiple-input multiple-output (MIMO) system in a quasi-static Rayleigh fading environment, it was shown by Elia et al. that linear space-time block code schemes (LSTBC schemes) that have the nonvanishing determinant (NVD) property are diversity-multiplexing gain tradeoff (DMT)-optimal for arbitrary values of n(r) if they have a code rate of n(t) complex dimensions per channel use. However, for asymmetric MIMO systems (where n(r) < n(t)), with the exception of a few LSTBC schemes, it is unknown whether general LSTBC schemes with NVD and a code rate of n(r) complex dimensions per channel use are DMT optimal. In this paper, an enhanced sufficient criterion for any STBC scheme to be DMT optimal is obtained, and using this criterion, it is established that any LSTBC scheme with NVD and a code rate of min {n(t), n(r)} complex dimensions per channel use is DMT optimal. This result settles the DMT optimality of several well-known, low-ML-decoding-complexity LSTBC schemes for certain asymmetric MIMO systems.

Single electron transfer-driven multi-dimensional signal read-out function of TCNQ as an ``off-the-shelf'' detector for cyanide

Herein we report the first applications of TCNQ as a rapid and highly sensitive off-the-shelf cyanide detector. As a proof-of-concept, we have applied a kinetically selective single-electron transfer (SET) from cyanide to deep-lying LUMO orbitals of TCNQ to generate a persistently stable radical anion (TCNQ(center dot-)), under ambient condition. In contrast to the known cyanide sensors that operate with limited signal outputs, TCNQ(center dot-) offers a unique multiple signaling platform. The signal readability is facilitated through multichannel absorption in the UV-vis-NIR region and scattering-based spectroscopic methods like Raman spectroscopy and hyper Rayleigh scattering techniques. Particularly notable is the application of the intense 840 nm NIR absorption band to detect cyanide. This can be useful for avoiding background interference in the UV-vis region predominant in biological samples. We also demonstrate the fabrication of a practical electronic device with TCNQ as a detector. The device generates multiorder enhancement in current with cyanide because of the formation of the conductive TCNQ(center dot-).

A note on new indices for the equatorial Indian Ocean oscillation

It is now well known that there is a strong association of the extremes of the Indian summer monsoon rainfall (ISMR) with the El Nio and southern oscillation (ENSO) and the Equatorial Indian Ocean Oscillation (EQUINOO), later being an east-west oscillation in convection anomaly over the equatorial Indian Ocean. So far, the index used for EQUINOO is EQWIN, which is based on the surface zonal wind over the central equatorial Indian Ocean. Since the most important attribute of EQUINOO is the oscillation in convection/precipitation, we believe that the indices based on convection or precipitation would be more appropriate. Continuous and reliable data on outgoing longwave radiation (OLR), and satellite derived precipitation are now available from 1979 onwards. Hence, in this paper, we introduce new indices for EQUINOO, based on the difference in the anomaly of OLR/precipitation between eastern and western parts of the equatorial Indian Ocean. We show that the strong association of extremes of the Indian summer monsoon with ENSO and EQUINOO is also seen when the new indices are used to represent EQUINOO.

Dynamics of a flexible splitter plate in the wake of a circular cylinder

Rigid splitter plates in the wake of bluff bodies are known to suppress the primary vortex shedding. In the present work, we experimentally study the problem of a flexible splitter plate in the wake of a circular cylinder. In this case, the splitter plate is free to continuously deform along its length due to the fluid forces acting on it; the flexural rigidity (EI) of the plate being an important parameter. Direct visualizations of the splitter plate motions, for very low values of flexural rigidity (EI), indicate periodic traveling wave type deformations of the splitter plate with maximum tip amplitudes of the order of I cylinder diameter. As the Reynolds number based on cylinder diameter is varied, two regimes of periodic splitter plate motions are found that are referred to as mode I and mode II, with a regime of aperiodic motions between them. The frequency of plate motions in both periodic modes is found to be close to the plane cylinder Strouhal number of about 0.2, while the average frequencies in the non-periodic regime are substantially lower. The measured normalized phase speed of the traveling wave for both periodic modes is also close to the convection speed of vortices in the plane cylinder wake. As the flexural rigidity of the plate (EI) is increased, the response of the plate was found to shift to the right when plotted with flow speed or Re. To better capture the effect of varying EI, we define and use a non-dimensional bending stiffness, K*, similar to the ones used in the flag flutter problem, K*=EI/(0.5 rho(UL3)-L-2), where U is the free-stream velocity and L is the splitter plate length. Amplitude data for different EI cases when plotted against this parameter appear to collapse on to a single curve for a given splitter plate length. Measurements of the splitter plate motions for varying splitter plate lengths indicate that plates that are substantially larger than the formation length of the plane cylinder wake have similar responses, while shorter plates show significant differences.

An equivalent undercooling model to account for flow effect on binary alloy dendrite growth rate

A theoretical analysis is carried out to observe the influence of important flow parameters such as Nusselt number and Sherwood number on the tip speed of an equiaxed dendrite growing in a convecting alloy melt. The effect of thermal and solutal transfer at the interface due to convection is equated to an undercooling of the melt, and an expression is derived for this equivalent undercooling in terms of the flow Nusselt number and Sherwood number. Results for the equivalent undercooling are compared with corresponding numerical values obtained by performing simulations based on the enthalpy method. This method represents a relatively simple procedure to analyze the effects of melt convection on the growth rate of dendrites. (C) 2013 Elsevier Ltd. All rights reserved.