Coupled subdaily and multiweek cycles during the lava dome erupt of Soufriere Hills Volcano, Montserrat

Observations of volcanoes extruding andesitic lava to produce lava domes often reveal cyclic behaviour. At Soufriere Hills Volcano, Montserrat, cycles with sub-daily and multi-week periods have been recognised on many occasions. These two types of cycle have been modelled separately as stick-slip magma flow at the junction between a dyke and an overlying cylindrical conduit (Costa et al. 2012), and as the filling and discharge of magma through the elastic-walled dyke (Costa et al., 2007a) respectively. Here, we couple these two models to simulate the behaviour over a period of well-observed multi-week cycles, with accompanying sub-daily cycles, from 13 May to 21 September 1997. The coupled model captures well the asymmetrical first-order behaviour: the first 40% of the multi-week cycle consists of high rates of lava extrusion during short period/high amplitude sub-daily cycles as the dyke reservoir discharges itself. The remainder of the cycle involves increasing pressurization as more magma is stored, and extrusion rate falls, followed by a gradual increase in the period of the sub-daily cycles.

Solutions of the fully compressible semi-geostrophic system

The fully compressible semi-geostrophic system is widely used in the modelling of large-scale atmospheric flows. In this paper, we prove rigorously the existence of weak Lagrangian solutions of this system, formulated in the original physical coordinates. In addition, we provide an alternative proof of the earlier result on the existence of weak solutions of this system expressed in the so-called geostrophic, or dual, coordinates. The proofs are based on the optimal transport formulation of the problem and on recent general results concerning transport problems posed in the Wasserstein space of probability measures.

Stability of an ice sheet on an elastic bed

The stability of stationary flow of a two-dimensional ice sheet is studied when the ice obeys a power flow law (Glen's flow law). The mass accumulation rate at the top is assumed to depend on elevation and span and the bed supporting the ice sheet consists of an elastic layer lying on a rigid surface. The normal perturbation of the free surface of the ice sheet is a singular eigenvalue problem. The singularity of the perturbation at the front of the ice sheet is considered using matched asymptotic expansions, and the eigenvalue problem is seen to reduce to that with fixed ice front. Numerical solution of the perturbation eigenvalue problem shows that the dependence of accumulation rate on elevation permits the existence of unstable solutions when the equilibrium line is higher than the bed at the ice divide. Alternatively, when the equilibrium line is lower than the bed, there are only stable solutions. Softening of the bed, expressed through a decrease of its elastic modulus, has a stabilising effect on the ice sheet.

Flow synthesis using gaseous ammonia in a Teflon AF-2400 tube-in-tube reactor: Paal–Knorr pyrrole formation and gas concentration measurement by inline flow titration

Using a simple and accessible Teflon AF-2400 based tube-intube reactor, a series of pyrroles were synthesised in flow using the Paal–Knorr reaction of 1,4-diketones with gaseous ammonia. An inline flow titration technique allowed measurement of the ammonia concentration and its relationship to residence time and temperature.

Continuous-flow processing of gaseous ammonia using a Teflon AF-2400 tube-in-tube reactor: synthesis of thioureas and in-line titrations

A simple tube-in-tube reactor based on the gas-permeable membrane Teflon AF-2400 was used in the continuous flow reaction of gaseous ammonia with isothiocyanates and one isocyanate. A colourimetric in-line titration technique is also reported as a simple method to quantify the amount of ammonia taken up by the solvent in the system.

Elastic net orthogonal forward regression

An efficient two-level model identification method aiming at maximising a model׳s generalisation capability is proposed for a large class of linear-in-the-parameters models from the observational data. A new elastic net orthogonal forward regression (ENOFR) algorithm is employed at the lower level to carry out simultaneous model selection and elastic net parameter estimation. The two regularisation parameters in the elastic net are optimised using a particle swarm optimisation (PSO) algorithm at the upper level by minimising the leave one out (LOO) mean square error (LOOMSE). There are two elements of original contributions. Firstly an elastic net cost function is defined and applied based on orthogonal decomposition, which facilitates the automatic model structure selection process with no need of using a predetermined error tolerance to terminate the forward selection process. Secondly it is shown that the LOOMSE based on the resultant ENOFR models can be analytically computed without actually splitting the data set, and the associate computation cost is small due to the ENOFR procedure. Consequently a fully automated procedure is achieved without resort to any other validation data set for iterative model evaluation. Illustrative examples are included to demonstrate the effectiveness of the new approaches.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Resilience of rice (Oryzaspp.) pollen germination and tube growth to temperature stress

Resilience of rice cropping systems to potential global climate change will partly depend on temperature tolerance of pollen germination (PG) and tube growth (PTG). Germination of pollen of high temperature susceptible Oryza glaberrima Steud. (cv. CG14) and O. sativa L. ssp. indica (cv. IR64) and high temperature tolerant O. sativa ssp. aus (cv. N22), was assessed on a 5.6-45.4°C temperature gradient system. Mean maximum PG was 85% at 27°C with 1488 μm PTG at 25°C. The hypothesis that in each pollen grain, minimum temperature requirements (Tn) and maximum temperature limits (Tx) for germination operate independently was accepted by comparing multiplicative and subtractive probability models. The maximum temperature limit for PG in 50% of grains (Tx(50)) was lowest (29.8°C) in IR64 compared with CG14 (34.3°C) and N22 (35.6°C). Standard deviation (sx) of Tx was also low in IR64 (2.3°C) suggesting that the mechanism of IR64's susceptibility to high temperatures may relate to PG. Optimum germination temperatures and thermal times for 1mm PTG were not linked to tolerating high temperatures at anthesis. However, the parameters Tx(50) and sx in the germination model define new pragmatic criteria for successful and resilient PG, preferable to the more traditional cardinal (maximum and minimum) temperatures.

Studying the contact point and interface moving in a sinusoidal tube with lattice Boltzmann method

The multicomponent nonideal gas lattice Boltzmann model by Shan and Chen (S-C) is used to study the immiscible displacement in a sinusoidal tube. The movement of interface and the contact point (contact line in three-dimension) is studied. Due to the roughness of the boundary, the contact point shows "stick-slip" mechanics. The "stick-slip" effect decreases as the speed of the interface increases. For fluids that are nonwetting, the interface is almost perpendicular to the boundaries at most time, although its shapes at different position of the tube are rather different. When the tube becomes narrow, the interface turns a complex curves rather than remains simple menisci. The velocity is found to vary considerably between the neighbor nodes close to the contact point, consistent with the experimental observation that the velocity is multi-values on the contact line. Finally, the effect of three boundary conditions is discussed. The average speed is found different for different boundary conditions. The simple bounce-back rule makes the contact point move fastest. Both the simple bounce-back and the no-slip bounce-back rules are more sensitive to the roughness of the boundary in comparison with the half-way bounce-back rule. The simulation results suggest that the S-C model may be a promising tool in simulating the displacement behaviour of two immiscible fluids in complex geometry.

Lattice Boltzmann simulation of viscous fluid systems with elastic boundaries

A lattice Boltzmann model able to simulate viscous fluid systems with elastic and movable boundaries is proposed. By introducing the virtual distribution function at the boundary, the Galilean invariance is recovered for the full system. As examples of application, the how in elastic vessels is simulated with the pressure-radius relationship similar to that of the pulmonary blood vessels. The numerical results for steady how are in good agreement with the analytical prediction, while the simulation results for pulsative how agree with the experimental observation of the aortic flows qualitatively. The approach has potential application in the study of the complex fluid systems such as the suspension system as well as the arterial blood flow.

Switchable amplification of vibrational circular dichroism as a probe of local chiral structure

A new method to detect the vibrational circular dichroism (VCD) of a localized part of a chiral molecular system is reported. A local VCD amplifier was implemented, and the distance dependence of the amplification was investigated in a series of peptides. The results indicate a characteristic distance of 2.0±0.3 bonds, which suggests that the amplification is a localized phenomenon. The amplifier can be covalently coupled to a specific part of a molecule, and can be switched ON and OFF electrochemically. By subtracting the VCD spectra obtained when the amplifier is in the ON and OFF states, the VCD of the local environment of the amplifier can be separated from the total VCD spectrum. Switchable local VCD amplification thus makes it possible to “zoom in” on a specific part of a chiral molecule.

Elucidating the structure of chiral molecules by using amplified vibrational circular dichroism: from theory to experimental realization

Recent experimental observations of enhanced vibrational circular dichroism (VCD) in molecular systems with low-lying electronically excited states suggest interesting new applications of VCD spectroscopy. The theory describing VCD enhancement through vibronic coupling schemes was derived by Nafie in 1983, but only recently experimental evidence of VCD amplification has demonstrated the extent to which this effect can be exploited as a structure elucidation tool to probe local structure. In this Concept paper, we give an overview of the physics behind vibrational circular dichroism, in particular the equations governing the VCD amplification effect, and review the latest experimental developments with a prospective view on the application of amplified VCD to locally probe biomolecular structure.

Dehydrodipeptide hydrogelators containing naproxen N‑Capped tryptophan: self-assembly, hydrogel characterization, and evaluation as potential drug nanocarriers

In this work, we introduce dipeptides containing tryptophan N-capped with the nonsteroidal anti-inflammatory drug naproxen and C-terminal dehydroamino acids, dehydrophenylalanine (ΔPhe), dehydroaminobutyric acid (ΔAbu), and dehydroalanine (ΔAla) as efficacious protease resistant hydrogelators. Optimized conditions for gel formation are reported. Transmission electron microscopy experiments revealed that the hydrogels consist of networks of micro/nanosized fibers formed by peptide self-assembly. Fluorescence and circular dichroism spectroscopy indicate that the self-assembly process is driven by stacking interactions of the aromatic groups. The naphthalene groups of the naproxen moieties are highly organized in the fibers through chiral stacking. Rheological experiments demonstrated that the most hydrophobic peptide (containing C-terminal ΔPhe) formed more elastic gels at lower critical gelation concentrations. This gel revealed irreversible breakup, while the C-terminal ΔAbu and ΔAla gels, although less elastic, exhibited structural recovery and partial healing of the elastic properties. A potential antitumor thieno[3,2-b]pyridine derivative was incorporated (noncovalently) into the gel formed by the hydrogelator containing C-terminal ΔPhe residue. Fluorescence and Förster resonance energy transfer measurements indicate that the drug is located in a hydrophobic environment, near/associated with the peptide fibers, establishing this type of hydrogel as a good drug-nanocarrier candidate.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.