Impact of mass lumping on gravity and Rossby waves in 2D finite-element shallow-water models

The goal of this study is to evaluate the effect of mass lumping on the dispersion properties of four finite-element velocity/surface-elevation pairs that are used to approximate the linear shallow-water equations. For each pair, the dispersion relation, obtained using the mass lumping technique, is computed and analysed for both gravity and Rossby waves. The dispersion relations are compared with those obtained for the consistent schemes (without lumping) and the continuous case. The P0-P1, RT0 and P-P1 pairs are shown to preserve good dispersive properties when the mass matrix is lumped. Test problems to simulate fast gravity and slow Rossby waves are in good agreement with the analytical results.

The interannual variability of energy transports within and over the Atlantic Ocean in a coupled climate model

To gain a new perspective on the interaction of the Atlantic Ocean and the atmosphere, the relationship between the atmospheric and oceanic meridional energy transports is studied in a version of HadCM3, the U.K. Hadley Centre's coupled climate model. The correlation structure of the energy transports in the atmosphere and Atlantic Ocean as a function of latitude, and the cross correlation between the two systems are analyzed. The processes that give rise to the correlations are then elucidated using regression analyses. In northern midlatitudes, the interannual variability of the Atlantic Ocean energy transport is dominated by Ekman processes. Anticorrelated zonal winds in the subtropics and midlatitudes, particularly associated with the North Atlantic Oscillation (NAO), drive anticorrelated meridional Ekman transports. Variability in the atmospheric energy transport is associated with changes in the stationary waves, but is only weakly related to the NAO. Nevertheless, atmospheric driving of the oceanic Ekman transports is responsible for a bipolar pattern in the correlation between the atmosphere and Atlantic Ocean energy transports. In the Tropics, the interannual variability of the Atlantic Ocean energy transport is dominated by an adjustment of the tropical ocean to coastal upwelling induced along the Venezuelan coast by a strengthening of the easterly trade winds. Variability in the atmospheric energy transport is associated with a cross-equatorial meridional overturning circulation that is only weakly associated with variability in the trade winds along the Venezuelan coast. In consequence, there is only very limited correlation between the atmosphere and Atlantic Ocean energy transports in the Tropics of HadCM3

Impact damage characterisation of thermoplastic matrix composites using transmission transient thermography

In this work, IR thermography is used as a non-destructive tool for impact damage characterisation on thermoplastic E-glass/polypropylene composites for automotive applications. The aim of this experimentation was to compare impact resistance and to characterise damage patterns of different laminates, in order to provide indications for their use in components. Two E-glass/polypropylene composites, commingled ®Twintex (with three different weave structures: directional, balanced and 3-D) and random reinforced GMT, were in particular characterised. Directional and balanced Twintex were also coupled in a number of hybrid configurations with GMT to evaluate the possible use of GMT/Twintex hybrids in high-energy absorption components. The laminates were impacted using a falling weight tower, with impact energies ranging from 15 J to penetration. Using IR thermography during cooling down following a long pulse (3 s), impact damaged areas were characterised and the influence of weave structure on damage patterns was studied. IR thermography offered good accuracy for laminates with thickness not exceeding 3.5 mm: this appears to be a limit for the direct use of this method on components, where more refined signal treatment would probably be needed for impact damage characterisation.

The impact of European legislative and technology measures to reduce air pollutants on air quality, human health and climate

European air quality legislation has reduced emissions of air pollutants across Europe since the 1970s, affecting air quality, human health and regional climate. We used a coupled composition-climate model to simulate the impacts of European air quality legislation and technology measures implemented between 1970 and 2010. We contrast simulations using two emission scenarios; one with actual emissions in 2010 and the other with emissions that would have occurred in 2010 in the absence of technological improvements and end-of-pipe treatment measures in the energy, industrial and road transport sectors. European emissions of sulphur dioxide, black carbon (BC) and organic carbon in 2010 are 53%, 59% and 32% lower respectively compared to emissions that would have occurred in 2010 in the absence of legislative and technology measures. These emission reductions decreased simulated European annual mean concentrations of fine particulate matter(PM2.5) by 35%, sulphate by 44%, BC by 56% and particulate organic matter by 23%. The reduction in PM2.5 concentrations is calculated to have prevented 80 000 (37 000–116 000, at 95% confidence intervals) premature deaths annually across the European Union, resulting in a perceived financial benefit to society of US$232 billion annually (1.4% of 2010 EU GDP). The reduction in aerosol concentrations due to legislative and technology measures caused a positive change in the aerosol radiative effect at the top of atmosphere, reduced atmospheric absorption and also increased the amount of solar radiation incident at the surface over Europe. We used an energy budget approximation to estimate that these changes in the radiative balance have increased European annual mean surface temperatures and precipitation by 0.45 ± 0.11 °C and by 13 ± 0.8 mm yr−1 respectively. Our results show that the implementation of European legislation and technological improvements to reduce the emission of air pollutants has improved air quality and human health over Europe, as well as having an unintended impact on the regional radiative balance and climate.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part I: Dome shape and internal structure

Lava domes comprise core, carapace, and clastic talus components. They can grow endogenously by inflation of a core and/or exogenously with the extrusion of shear bounded lobes and whaleback lobes at the surface. Internal structure is paramount in determining the extent to which lava dome growth evolves stably, or conversely the propensity for collapse. The more core lava that exists within a dome, in both relative and absolute terms, the more explosive energy is available, both for large pyroclastic flows following collapse and in particular for lateral blast events following very rapid removal of lateral support to the dome. Knowledge of the location of the core lava within the dome is also relevant for hazard assessment purposes. A spreading toe, or lobe of core lava, over a talus substrate may be both relatively unstable and likely to accelerate to more violent activity during the early phases of a retrogressive collapse. Soufrière Hills Volcano, Montserrat has been erupting since 1995 and has produced numerous lava domes that have undergone repeated collapse events. We consider one continuous dome growth period, from August 2005 to May 2006 that resulted in a dome collapse event on 20th May 2006. The collapse event lasted 3 h, removing the whole dome plus dome remnants from a previous growth period in an unusually violent and rapid collapse event. We use an axisymmetrical computational Finite Element Method model for the growth and evolution of a lava dome. Our model comprises evolving core, carapace and talus components based on axisymmetrical endogenous dome growth, which permits us to model the interface between talus and core. Despite explicitly only modelling axisymmetrical endogenous dome growth our core–talus model simulates many of the observed growth characteristics of the 2005–2006 SHV lava dome well. Further, it is possible for our simulations to replicate large-scale exogenous characteristics when a considerable volume of talus has accumulated around the lower flanks of the dome. Model results suggest that dome core can override talus within a growing dome, potentially generating a region of significant weakness and a potential locus for collapse initiation.

A potential energy surface for the ground state of formaldehyde, H2CO(1A1)

A model potential energy function for the ground state of H2CO has been derived which covers the whole space of the six internal coordinates. This potential reproduces the experimental energy, geometry and quadratic force field of formaldehyde, and dissociates correctly to all possible atom, diatom and triatom fragments. Thus there are good reasons for believing it to be close to the true potential energy surface except in regions where both hydrogen atoms are close to the oxygen. It leads to the prediction that there should be a metastable singlet hydroxycarbene HCOH which has a planar trans structure and an energy of 2•31 eV above that of equilibrium formaldehyde. The reaction path for dissociation into H2 + CO is predicted to pass through a low symmetry transition state with an activation energy of 4•8 eV. Both of these predictions are in good agreement with recently published ab initio calculations.

Electronic transitions and bonding properties in a series of five-coordinate “16-electron” complexes [Mn(CO)3(L2)]− (L2=chelating redox-active π-donor ligand)

In this article we present for the first time accurate density functional theory (DFT) and time-dependent (TD) DFT data for a series of electronically unsaturated five-coordinate complexes [Mn(CO)(3)(L-2)](-), where L-2 stands for a chelating strong pi-donor ligand represented by catecholate, dithiolate, amidothiolate, reduced alpha-diimine (1,4-dialkyl-1,4-diazabutadiene (R-DAB), 2,2'-bipyridine) and reduced 2,2'-biphosphinine types. The single-crystal X-ray structure of the unusual compound [Na(BPY)][Mn(CO)(3)(BPY)]center dot Et2O and the electronic absorption spectrum of the anion [Mn(CO)(3)(BPY)](-) are new in the literature. The nature of the bidentate ligand determines the bonding in the complexes, which varies between two limiting forms: from completely pi-delocalized diamagnetic {(CO)(3)Mn-L-2}(-) for L-2 = alpha-diimine or biphosphinine, to largely valence-trapped {(CO)(3)Mn-1-L-2(2-)}(-) for L-2(2-) = catecholate, where the formal oxidation states of Mn and L-2 can be assigned. The variable degree of the pi-delocalization in the Mn(L-2) chelate ring is indicated by experimental resonance Raman spectra of [Mn(CO)(3)(L-2)](-) (L-2=3,5-di-tBu-catecholate and iPr-DAB), where accurate assignments of the diagnostically important Raman bands have been aided by vibrational analysis. The L-2 = catecholate type of complexes is known to react with Lewis bases (CO substitution, formation of six-coordinate adducts) while the strongly pi-delocalized complexes are inert. The five-coordinate complexes adopt usually a distorted square pyramidal geometry in the solid state, even though transitions to a trigonal bipyramid are also not rare. The experimental structural data and the corresponding DFT-computed values of bond lengths and angles are in a very good agreement. TD-DFT calculations of electronic absorption spectra of the studied Mn complexes and the strongly pi-delocalized reference compound [Fe(CO)(3)(Me-DAB)] have reproduced qualitatively well the experimental spectra. Analyses of the computed electronic transitions in the visible spectroscopic region show that the lowest-energy absorption band always contains a dominant (in some cases almost exclusive) contribution from a pi(HOMO) -> pi*(LUMO) transition within the MnL2 metallacycle. The character of this optical excitation depends strongly on the composition of the frontier orbitals, varying from a partial L-2 -> Mn charge transfer (LMCT) through a fully delocalized pi(MnL2) -> pi*(MnL2) situation to a mixed (CO)Mn -> L-2 charge transfer (LLCT/MLCT). The latter character is most apparent in the case of the reference complex [Fe(CO)(3)(Me-DAB)]. The higher-lying, usually strongly mixed electronic transitions in the visible absorption region originate in the three lower-lying occupied orbitals, HOMO - 1 to HOMO - 3, with significant metal-d contributions. Assignment of these optical excitations to electronic transitions of a specific type is difficult. A partial LLCT/MLCT character is encountered most frequently. The electronic absorption spectra become more complex when the chelating ligand L-2, such as 2,2'-bipyridine, features two or more closely spaced low-lying empty pi* orbitals.

Post-occupancy studies of an office environment: Energy performance and occupants’ satisfaction

It is generally accepted that the physical workplace environment affects employees’ satisfaction and, consequently, their perceived productivity and well-being. This study investigated whether employee “satisfaction” variables can predict perceived productivity, well-being and enjoyment at work, and if so, to what extent. The study also explored whether limiting employees’ control over their environment could save energy without compromising employees’ satisfaction and perceived productivity. Preoccupancy and post-occupancy evaluation studies were conducted, in terms of both energy consumption and employee perceptions, to make comparisons between a company’s old and current headquarters buildings, both located in the same area of London. The results showed that employees were more satisfied with their work environment at their new HQ, in general, than with that of their previous office. Also, employees’ self-reported productivity, well-being and enjoyment at work improved after the move. It was revealed that the combination of employees’ level of satisfaction with “interior use of space” and “physical conditions” was the best predictor of their perceived productivity, while satisfaction with “indoor facilities” was not a good predictor. In terms of energy performance, although the new HQ’s energy consumption per m2 was significantly less than that of the previous building, there was still a gap between the refurbishment design target and the actual performance of the building. The findings suggest that this gap could be due to a number of factors, including an ineffective use of interior space, and occupants’ behaviour.

Rapid assessment of historic, current and future habitat quality for biodiversity around UK Natura 2000 sites

Changes in landscape composition and structure may impact the conservation and management of protected areas. Species that depend on specific habitats are at risk of extinction when these habitats are degraded or lost. Designing robust methods to evaluate landscape composition will assist decision- and policy-making in emerging landscapes. This paper describes a rapid assessment methodology aimed at evaluating landcover quality for birds, plants, butterflies and bees around seven UK Natura 2000 sites. An expert panel assigned quality values to standard Coordination of Information on the Environment (CORINE) landcover classes for each taxonomic group. Quality was assessed based on historical (1950, 1990), current (2000) and future (2030) land-cover data, the last projected using three alternative scenarios: a growth applied strategy (GRAS), a business-as-might-beusual (BAMBU) scenario, and sustainable European development goal (SEDG) scenario. A quantitative quality index weighted the area of each land-cover parcel with a taxa-specific quality measure. Land parcels with high quality for all taxonomic groups were evaluated for temporal changes in area, size and adjacency. For all sites and taxonomic groups, the rate of deterioration of land-cover quality was greater between 1950 and 1990 than current rates or as modelled using the alternative future scenarios (2000– 2030). Model predictions indicated land-cover quality stabilized over time under the GRAS scenario, and was close to stable for the BAMBU scenario. The SEDG scenario suggested an ongoing loss of quality, though this was lower than the historical rate of c. 1% loss per decade. None of the future scenarios showed accelerated fragmentation, but rather increases in the area, adjacency and diversity of high quality land parcels in the landscape.

Impact of atmospheric forcing data on simulations of the Laptev Sea polynya dynamics using the sea-ice ocean model FESOM

The polynyas of the Laptev Sea are regions of particular interest due to the strong formation of Arctic sea-ice. In order to simulate the polynya dynamics and to quantify ice production, we apply the Finite Element Sea-Ice Ocean Model FESOM. In previous simulations FESOM has been forced with daily atmospheric NCEP (National Centers for Environmental Prediction) 1. For the periods 1 April to 9 May 2008 and 1 January to 8 February 2009 we examine the impact of different forcing data: daily and 6-hourly NCEP reanalyses 1 (1.875° x 1.875°), 6-hourly NCEP reanalyses 2 (1.875° x 1.875°), 6-hourly analyses from the GME (Global Model of the German Weather Service) (0.5° x 0.5°) and high-resolution hourly COSMO (Consortium for Small-Scale Modeling) data (5 km x 5 km). In all FESOM simulations, except for those with 6-hourly and daily NCEP 1 data, the openings and closings of polynyas are simulated in principle agreement with satellite products. Over the fast-ice area the wind fields of all atmospheric data are similar and close to in situ measurements. Over the polynya areas, however, there are strong differences between the forcing data with respect to air temperature and turbulent heat flux. These differences have a strong impact on sea-ice production rates. Depending on the forcing fields polynya ice production ranges from 1.4 km3 to 7.8 km3 during 1 April to 9 May 2011 and from 25.7 km3 to 66.2 km3 during 1 January to 8 February 2009. Therefore, atmospheric forcing data with high spatial and temporal resolution which account for the presence of the polynyas are needed to reduce the uncertainty in quantifying ice production in polynyas.

A cell by cell anisotropic adaptive mesh ALE scheme for the numerical solution of the Euler equations

In this paper a cell by cell anisotropic adaptive mesh technique is added to an existing staggered mesh Lagrange plus remap finite element ALE code for the solution of the Euler equations. The quadrilateral finite elements may be subdivided isotropically or anisotropically and a hierarchical data structure is employed. An efficient computational method is proposed, which only solves on the finest level of resolution that exists for each part of the domain with disjoint or hanging nodes being used at resolution transitions. The Lagrangian, equipotential mesh relaxation and advection (solution remapping) steps are generalised so that they may be applied on the dynamic mesh. It is shown that for a radial Sod problem and a two-dimensional Riemann problem the anisotropic adaptive mesh method runs over eight times faster.

A review of the theoretical basis for bulk mass flux convective parameterization

Most parameterizations for precipitating convection in use today are bulk schemes, in which an ensemble of cumulus elements with different properties is modelled as a single, representative entraining-detraining plume. We review the underpinning mathematical model for such parameterizations, in particular by comparing it with spectral models in which elements are not combined into the representative plume. The chief merit of a bulk model is that the representative plume can be described by an equation set with the same structure as that which describes each element in a spectral model. The equivalence relies on an ansatz for detrained condensate introduced by Yanai et al. (1973) and on a simplified microphysics. There are also conceptual differences in the closure of bulk and spectral parameterizations. In particular, we show that the convective quasi-equilibrium closure of Arakawa and Schubert (1974) for spectral parameterizations cannot be carried over to a bulk parameterization in a straightforward way. Quasi-equilibrium of the cloud work function assumes a timescale separation between a slow forcing process and a rapid convective response. But, for the natural bulk analogue to the cloud-work function (the dilute CAPE), the relevant forcing is characterised by a different timescale, and so its quasi-equilibrium entails a different physical constraint. Closures of bulk parameterization that use the non-entraining parcel value of CAPE do not suffer from this timescale issue. However, the Yanai et al. (1973) ansatz must be invoked as a necessary ingredient of those closures.

Sonic to ultrasonic Q of sandstones and limestones: Laboratory measurements at in situ pressures

Laboratory measurements of the attenuation and velocity dispersion of compressional and shear waves at appropriate frequencies, pressures, and temperatures can aid interpretation of seismic and well-log surveys as well as indicate absorption mechanisms in rocks. Construction and calibration of resonant-bar equipment was used to measure velocities and attenuations of standing shear and extensional waves in copper-jacketed right cylinders of rocks (30 cm in length, 2.54 cm in diameter) in the sonic frequency range and at differential pressures up to 65 MPa. We also measured ultrasonic velocities and attenuations of compressional and shear waves in 50-mm-diameter samples of the rocks at identical pressures. Extensional-mode velocities determined from the resonant bar are systematically too low, yielding unreliable Poisson's ratios. Poisson's ratios determined from the ultrasonic data are frequency corrected and used to calculate the sonic-frequency compressional-wave velocities and attenuations from the shear- and extensional-mode data. We calculate the bulk-modulus loss. The accuracies of attenuation data (expressed as 1000/Q, where Q is the quality factor) are +/- 1 for compressional and shear waves at ultrasonic frequency, +/- 1 for shear waves, and +/- 3 for compressional waves at sonic frequency. Example sonic-frequency data show that the energy absorption in a limestone is small (Q(P) greater than 200 and stress independent) and is primarily due to poroelasticity, whereas that in the two sandstones is variable in magnitude (Q(P) ranges from less than 50 to greater than 300, at reservoir pressures) and arises from a combination of poroelasticity and viscoelasticity. A graph of compressional-wave attenuation versus compressional-wave velocity at reservoir pressures differentiates high-permeability (> 100 mD, 9.87 X 10(-14) m(2)) brine-saturated sandstones from low-permeability (< 100 mD, 9.87 X 10 (14) m(2)) sandstones and shales.

Modelling the lava dome extruded at Soufrière Hills Volcano, Montserrat, August 2005–May 2006. Part II: Rockfall activity and talus deformation

During many lava dome-forming eruptions, persistent rockfalls and the concurrent development of a substantial talus apron around the foot of the dome are important aspects of the observed activity. An improved understanding of internal dome structure, including the shape and internal boundaries of the talus apron, is critical for determining when a lava dome is poised for a major collapse and how this collapse might ensue. We consider a period of lava dome growth at the Soufrière Hills Volcano, Montserrat, from August 2005 to May 2006, during which a 100 × 106 m3 lava dome developed that culminated in a major dome-collapse event on 20 May 2006. We use an axi-symmetrical Finite Element Method model to simulate the growth and evolution of the lava dome, including the development of the talus apron. We first test the generic behaviour of this continuum model, which has core lava and carapace/talus components. Our model describes the generation rate of talus, including its spatial and temporal variation, as well as its post-generation deformation, which is important for an improved understanding of the internal configuration and structure of the dome. We then use our model to simulate the 2005 to 2006 Soufrière Hills dome growth using measured dome volumes and extrusion rates to drive the model and generate the evolving configuration of the dome core and carapace/talus domains. The evolution of the model is compared with the observed rockfall seismicity using event counts and seismic energy parameters, which are used here as a measure of rockfall intensity and hence a first-order proxy for volumes. The range of model-derived volume increments of talus aggraded to the talus slope per recorded rockfall event, approximately 3 × 103–13 × 103 m3 per rockfall, is high with respect to estimates based on observed events. From this, it is inferred that some of the volumetric growth of the talus apron (perhaps up to 60–70%) might have occurred in the form of aseismic deformation of the talus, forced by an internal, laterally spreading core. Talus apron growth by this mechanism has not previously been identified, and this suggests that the core, hosting hot gas-rich lava, could have a greater lateral extent than previously considered.

Redox chemistry and electronic properties of 2,3,5,6-tetrakis(2-pyridyl)pyrazine-bridged diruthenium complexes controlled by N,C,N '-biscyclometalated ligands

To investigate the consequences of cyclometalation for electronic communication in dinuclear ruthenium complexes, a series of 2,3,5,6-tetrakis(2-pyridyl)pyrazine (tppz) bridged diruthenium complexes was prepared and studied. These complexes have a central tppz ligand bridging via nitrogen-to-ruthenium coordination bonds, while each ruthenium atom also binds either a monoanionic, N,C,N'-terdentate 2,6-bis(2'-pyridyl)phenyl (R-N boolean AND C boolean AND N) ligand or a 2,2':6',2 ''-terpyridine (tpy) ligand. The N,C,N'-, that is, biscyclometalation, instead of the latter N,N', N ''-bonding motif significantly changes the electronic properties of the resulting complexes. Starting from well-known [{Ru(tpy)}(2)(mu-tppz)](4+) (tpy = 2,2':2 '',6-terpyridine) ([3](4+)) as a model compound, the complexes [{Ru(R-N boolean AND C boolean AND N)}(mu-tppz){Ru(tpy)}](3+) (R-N boolean AND C(H)boolean AND N = 4-R-1,3-dipyridylbenzene, R = H ([4a](3+)), CO2Me ([4b](3+))), and [{Ru(R-N boolean AND C boolean AND N)}(2)(mu-tppz)](2+), (R = H ([5a](2+)), CO2Me ([5b](2+))) were prepared with one or two N,C,N'-cyclometalated terminal ligands. The oxidation and reduction potentials of cyclometalated [4](3+) and [5](2+) are shifted negatively compared to non-cyclometalated [3](4+), the oxidation processes being affected more significantly. Compared to [3](4+), the electronic spectra of [5](2+) display large bathochromic shifts of the main MLCT transitions in the visible spectral region with low-energy absorptions tailing down to the NIR region. One-electron oxidation of [3](4+) and [5](2+) gives rise to low-energy absorption bands. The comproportionation constants and NIR band shape correspond to delocalized Robin-Day class III compounds. Complexes [4a](3+) (R = H) and [4b](3+) (R = CO2Me) also exhibit strong electronic communication, and notwithstanding the large redox-asymmetry the visible metal-to-ligand charge-transfer absorption is assigned to originate from both metal centers. The potential of the first, ruthenium-based, reversible oxidation process is strongly negatively shifted. On the contrary, the second oxidation is irreversible and cyclometalated ligand-based. Upon one-electron oxidation, a weak and low-energy absorption arises.