Importance of resolution and model configuration when downscaling extreme precipitation

Dynamical downscaling is frequently used to investigate the dynamical variables of extra-tropical cyclones, for example, precipitation, using very high-resolution models nested within coarser resolution models to understand the processes that lead to intense precipitation. It is also used in climate change studies, using long timeseries to investigate trends in precipitation, or to look at the small-scale dynamical processes for specific case studies. This study investigates some of the problems associated with dynamical downscaling and looks at the optimum configuration to obtain the distribution and intensity of a precipitation field to match observations. This study uses the Met Office Unified Model run in limited area mode with grid spacings of 12, 4 and 1.5 km, driven by boundary conditions provided by the ECMWF Operational Analysis to produce high-resolution simulations for the Summer of 2007 UK flooding events. The numerical weather prediction model is initiated at varying times before the peak precipitation is observed to test the importance of the initialisation and boundary conditions, and how long the simulation can be run for. The results are compared to raingauge data as verification and show that the model intensities are most similar to observations when the model is initialised 12 hours before the peak precipitation is observed. It was also shown that using non-gridded datasets makes verification more difficult, with the density of observations also affecting the intensities observed. It is concluded that the simulations are able to produce realistic precipitation intensities when driven by the coarser resolution data.

Multistep π dimerization of tetrakis(n-decyl)heptathienoacene radical cations: a combined experimental and theoretical study

Radical cations of a heptathienoacene a,b-substituted with four n-decyl side groups (D4T7C+) form exceptionally stable p-dimer dications already at ambient temperature (Chem. Comm. 2011, 47, 12622). This extraordinary p-dimerization process is investigated here with a focus on the ultimate[D4T7C+]2 p-dimer dication and yet-unreported transitoryspecies formed during and after the oxidation. To this end, we use a joint experimental and theoretical approach that combines cyclic voltammetry, in situ spectrochemistry and spectroelectrochemistry, EPR spectroscopy, and DFT calculations. The impact of temperature, thienoacene concentration, and the nature and concentration of counteranions on the p-dimerization process is also investigated in detail. Two different transitory species were detected in the course of the one-electron oxidation: 1) a different transient conformation of the ultimate [D4T7C+]2 p-dimer dications, the stability of which is strongly affected by the applied experimental conditions, and 2) intermediate [D4T7]2C+ p-dimer radical cations formed prior to the fully oxidized [D4T7]2C+ p-dimer dications. Thus, this comprehensive work demonstrates the formation of peculiar supramolecular species of heptathienoacene radical cations, the stability, nature, and structure of which have been successfully analyzed. We therefore believe that this study leads to a deeper fundamental understanding of the mechanism of dimer formation between conjugated aromatic systems.

Thiocyanate complexes of uranium in multiple oxidation states: a combined structural, magnetic, spectroscopic, spectroelectrochemical, and theoretical study

A comprehensive study of the complexes A4[U(NCS)8] (A = Cs, Et4N, nBu4N) and A3[UO2(NCS)5] (A = Cs, Et4N) is described, with the crystal structures of [nBu4N]4[U(NCS)8]·2MeCN and Cs3[UO2(NCS)5]·O0.5 reported. The magnetic properties of square antiprismatic Cs4[U(NCS)8] and cubic [Et4N]4[U(NCS)8] have been probed by SQUID magnetometry. The geometry has an important impact on the low-temperature magnetic moments: at 2 K, μeff = 1.21 μB and 0.53 μB, respectively. Electronic absorption and photoluminescence spectra of the uranium(IV) compounds have been measured. The redox chemistry of [Et4N]4[U(NCS)8] has been explored using IR and UV–vis spectroelectrochemical methods. Reversible 1-electron oxidation of one of the coordinated thiocyanate ligands occurs at +0.22 V vs Fc/Fc+, followed by an irreversible oxidation to form dithiocyanogen (NCS)2 which upon back reduction regenerates thiocyanate anions coordinating to UO22+. NBO calculations agree with the experimental spectra, suggesting that the initial electron loss of [U(NCS)8]4– is delocalized over all NCS– ligands. Reduction of the uranyl(VI) complex [Et4N]3[UO2(NCS)5] to uranyl(V) is accompanied by immediate disproportionation and has only been studied by DFT methods. The bonding in [An(NCS)8]4– (An = Th, U) and [UO2(NCS)5]3– has been explored by a combination of DFT and QTAIM analysis, and the U–N bonds are predominantly ionic, with the uranyl(V) species more ionic that the uranyl(VI) ion. Additionally, the U(IV)–NCS ion is more ionic than what was found for U(IV)–Cl complexes.

The effect of model configuration on modelled hillslope–riparian interactions

The transfer of hillslope water to and through the riparian zone forms a research area of importance in hydrological investigations. Numerical modelling schemes offer a way to visualise and quantify first-order controls on catchment runoff response and mixing. We use a two-dimensional Finite Element model to assess the link between model setup decisions (e.g. zero-flux boundary definitions, soil algorithm choice) and the consequential hydrological process behaviour. A detailed understanding of the consequences of model configuration is required in order to produce reliable estimates of state variables. We demonstrate that model configuration decisions can determine effectively the presence or absence of particular hillslope flow processes and, the magnitude and direction of flux at the hillslope–riparian interface. If these consequences are not fully explored for any given scheme and application, the resulting process inference may well be misleading.

A multipoint study of a substorm occurring on 7 December, 1992, and its theoretical implications

On 7 December 1992, a moderate substorm was observed by a variety of satellites and ground-based instruments. Ionospheric flows were monitored near dusk by the Goose Bay HF radar and near midnight by the EISCAT radar. The observed flows are compared here with magnetometer observations by the IMAGE array in Scandinavia and the two Greenland chains, the auroral distribution observed by Freja and the substorm cycle observations by the SABRE radar, the SAMNET magnetometer array and LANL geosynchronous satellites. Data from Galileo Earth-encounter II are used to estimate the IMF B-z component. The data presented show that the substorm onset electrojet at midnight was confined to closed field lines equatorward of the preexisting convection reversal boundaries observed in the dusk and midnight regions. No evidence of substantial closure of open flux was detected following this substorm onset. Indeed the convection reversal boundary on the duskside continued to expand equatorward after onset due to the continued presence of strong southward IMF, such that growth and expansion phase features were simultaneously present. Clear indications of closure of open flux were not observed until a subsequent substorm intensification 25 min after the initial onset. After this time, the substorm auroral bulge in the nightside hours propagated well poleward of the pre-existing convection reversal boundary, and strong flow perturbations were observed by the Goose Bay radar, indicative of flows driven by reconnection in the tail.

Managing change in complex projects: configuration management, asset information and big data

As we enter an era of ‘big data’, asset information is becoming a deliverable of complex projects. Prior research suggests digital technologies enable rapid, flexible forms of project organizing. This research analyses practices of managing change in Airbus, CERN and Crossrail, through desk-based review, interviews, visits and a cross-case workshop. These organizations deliver complex projects, rely on digital technologies to manage large data-sets; and use configuration management, a systems engineering approach with mid-20th century origins, to establish and maintain integrity. In them, configuration management has become more, rather than less, important. Asset information is structured, with change managed through digital systems, using relatively hierarchical, asynchronous and sequential processes. The paper contributes by uncovering limits to flexibility in complex projects where integrity is important. Challenges of managing change are discussed, considering the evolving nature of configuration management; potential use of analytics on complex projects; and implications for research and practice.

A coupled cloud physics–radiation parameterization of the bulk optical properties of cirrus and its impact on the Met Office unified model global atmosphere 5.0 configuration

A new coupled cloud physics–radiation parameterization of the bulk optical properties of ice clouds is presented. The parameterization is consistent with assumptions in the cloud physics scheme regarding particle size distributions (PSDs) and mass–dimensional relationships. The parameterization is based on a weighted ice crystal habit mixture model, and its bulk optical properties are parameterized as simple functions of wavelength and ice water content (IWC). This approach directly couples IWC to the bulk optical properties, negating the need for diagnosed variables, such as the ice crystal effective dimension. The parameterization is implemented into the Met Office Unified Model Global Atmosphere 5.0 (GA5) configuration. The GA5 configuration is used to simulate the annual 20-yr shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere (TOA), as well as the temperature structure of the atmosphere, under various microphysical assumptions. The coupled parameterization is directly compared against the current operational radiation parameterization, while maintaining the same cloud physics assumptions. In this experiment, the impacts of the two parameterizations on the SW and LW radiative effects at TOA are also investigated and compared against observations. The 20-yr simulations are compared against the latest observations of the atmospheric temperature and radiative fluxes at TOA. The comparisons demonstrate that the choice of PSD and the assumed ice crystal shape distribution are as important as each other. Moreover, the consistent radiation parameterization removes a long-standing tropical troposphere cold temperature bias but slightly warms the southern midlatitudes by about 0.5 K.

An evaluation of the theoretical potential and practical opportunity for using greywater for domestic purposes in Ghana

Assessments concerning the effects of climate change, water resource availability and water deprivation in West Africa have not frequently considered the positive contribution to be derived from collecting and reusing water for domestic purposes. Where the originating water is taken from a clean water source and has been used the first time for washing or bathing, this water is commonly called “greywater”. Greywater is a prolific resource that is generated wherever people live. Treated greywater can be used for domestic cleaning, for flushing toilets where appropriate, for washing cars, sometimes for watering kitchen gardens, and for clothes washing prior to rinsing. Therefore, a large theoretical potential exists to increase total water resource availability if greywater were to be widely reused. Locally treated greywater reduces the distribution network requirement, lower construction effort and cost and, wherever possible, minimising the associated carbon footprint. Such locally treated greywater offers significant practical opportunities for increasing the total available water resources at a local level. The reuse of treated greywater is one important action that will help to mitigate the reducing availability of clean water supplies in some areas, and the expected mitigation required in future aligns well with WHO/UNICEF (2012) aspirations. The evaluation of potential opportunities for prioritising greywater systems to support water reuse takes into account the availability of water resources, water use indicators and published estimates in order to understand typical patterns of water demand. The approach supports knowledge acquisition regarding local conditions for enabling capacity building for greywater reuse, the understanding of systems that are most likely to encourage greywater reuse, and practices and future actions to stimulate greywater infrastructure planning, design and implementation. Although reuse might be considered to increase the uncertainty of achieving a specified quality of the water supply, robust methods and technologies are available for local treatment. Resource strategies for greywater reuse have the potential to consistently improve water efficiency and availability in water impoverished and water stressed regions of Ghana and West Africa. Untreated greywater is referred to as “greywater”; treated greywater is referred to as “treated greywater” in this paper.

A critical evaluation of the ability of the Spinning Enhanced Visible and Infrared Imager (SEVIRI) thermal infrared red-green-blue rendering to identify dust events: theoretical analysis

Using a combination of idealized radiative transfer simulations and a case study from the first field campaign of the Saharan Mineral Dust Experiment (SAMUM) in southern Morocco, this paper provides a systematic assessment of the limitations of the widely used Spinning Enhanced Visible and Infrared Imager (SEVIRI) red-green-blue (RGB) thermal infrared dust product. Both analyses indicate that the ability of the product to identify dust, via its characteristic pink coloring, is strongly dependent on the column water vapor, the lower tropospheric lapse rate, and dust altitude. In particular, when column water vapor exceeds ∼20–25 mm, dust presence, even for visible optical depths of the order 0.8, is effectively masked. Variability in dust optical properties also has a marked impact on the imagery, primarily as a result of variability in dust composition. There is a moderate sensitivity to the satellite viewing geometry, particularly in moist conditions. The underlying surface can act to confound the signal seen through variations in spectral emissivity, which are predominantly manifested in the 8.7μm SEVIRI channel. In addition, if a temperature inversion is present, typical of early morning conditions over the Sahara and Sahel, an increased dust loading can actually reduce the pink coloring of the RGB image compared to pristine conditions. Attempts to match specific SEVIRI observations to simulations using SAMUM measurements are challenging because of high uncertainties in surface skin temperature and emissivity. Recommendations concerning the use and interpretation of the SEVIRI RGB imagery are provided on the basis of these findings.

The Met Office Global Coupled model 2.0 (GC2) configuration

The latest coupled configuration of the Met Office Unified Model (Global Coupled configuration 2, GC2) is presented. This paper documents the model components which make up the configuration (although the scientific description of these components is detailed elsewhere) and provides a description of the coupling between the components. The performance of GC2 in terms of its systematic errors is assessed using a variety of diagnostic techniques. The configuration is intended to be used by the Met Office and collaborating institutes across a range of timescales, with the seasonal forecast system (GloSea5) and climate projection system (HadGEM) being the initial users. In this paper GC2 is compared against the model currently used operationally in those two systems. Overall GC2 is shown to be an improvement on the configurations used currently, particularly in terms of modes of variability (e.g. mid-latitude and tropical cyclone intensities, the Madden–Julian Oscillation and El Niño Southern Oscillation). A number of outstanding errors are identified with the most significant being a considerable warm bias over the Southern Ocean and a dry precipitation bias in the Indian and West African summer monsoons. Research to address these is ongoing.

Theoretical insight into diagnosing observation error correlations using observation-minus-background and observation-minus-analysis statistics

To improve the quantity and impact of observations used in data assimilation it is necessary to take into account the full, potentially correlated, observation error statistics. A number of methods for estimating correlated observation errors exist, but a popular method is a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. The accuracy of the results it yields is unknown as the diagnostic is sensitive to the difference between the exact background and exact observation error covariances and those that are chosen for use within the assimilation. It has often been stated in the literature that the results using this diagnostic are only valid when the background and observation error correlation length scales are well separated. Here we develop new theory relating to the diagnostic. For observations on a 1D periodic domain we are able to the show the effect of changes in the assumed error statistics used in the assimilation on the estimated observation error covariance matrix. We also provide bounds for the estimated observation error variance and eigenvalues of the estimated observation error correlation matrix. We demonstrate that it is still possible to obtain useful results from the diagnostic when the background and observation error length scales are similar. In general, our results suggest that when correlated observation errors are treated as uncorrelated in the assimilation, the diagnostic will underestimate the correlation length scale. We support our theoretical results with simple illustrative examples. These results have potential use for interpreting the derived covariances estimated using an operational system.

Diruthenium complexes with bridging diethynyl polyaromatic ligands: synthesis, spectroelectrochemistry, and theoretical calculations

This work describes syntheses and electrochemical, spectroscopic, and bonding properties in a new series of dinuclear ruthenium(II) complexes bridged by polyaromatic (biphenyl, fluorene, phenanthrene, and pyrene) alkynyl ligands. Longitudinal expansion of the π-conjugated polyaromatic core of the bridging ligands caused a reduced potential difference between the anodic steps and reinforced their bridge-localized nature, as evidenced by UV/vis/near-IR and IR spectroelectrochemical data combined with DFT and TDDFT calculations. Importantly, the intricate multiple IR ν(CC) absorption bands for the singly oxidized states imply a thermal population of a range of conformers (rotamers) with distinct electronic character. This behavior was demonstrated with more accurate DFT calculations of selected nontruncated 1e− oxidized complexes in three different conformations. The combined experimental and theoretical data reveal that thermally populated rotamers featuring various mutual orientations of the ligated metal termini and the bridging diethynyl polyaromatic moieties have a significant impact on the electronic absorption and ν(CC) wavenumbers of the singly oxidized systems.

Ice-sheet configuration in the CMIP5/PMIP3 Last Glacial Maximum experiments

We describe the creation of a data set describing changes related to the presence of ice sheets, including ice-sheet extent and height, ice-shelf extent, and the distribution and elevation of ice-free land at the Last Glacial Maximum (LGM), which were used in LGM experiments conducted as part of the fifth phase of the Coupled Modelling Intercomparison Project (CMIP5) and the third phase of the Palaeoclimate Modelling Intercomparison Project (PMIP3). The CMIP5/PMIP3 data sets were created from reconstructions made by three different groups, which were all obtained using a model-inversion approach but differ in the assumptions used in the modelling and in the type of data used as constraints. The ice-sheet extent in the Northern Hemisphere (NH) does not vary substantially between the three individual data sources. The difference in the topography of the NH ice sheets is also moderate, and smaller than the differences between these reconstructions (and the resultant composite reconstruction) and ice-sheet reconstructions used in previous generations of PMIP. Only two of the individual reconstructions provide information for Antarctica. The discrepancy between these two reconstructions is larger than the difference for the NH ice sheets, although still less than the difference between the composite reconstruction and previous PMIP ice-sheet reconstructions. Although largely confined to the ice-covered regions, differences between the climate response to the individual LGM reconstructions extend over the North Atlantic Ocean and Northern Hemisphere continents, partly through atmospheric stationary waves. Differences between the climate response to the CMIP5/PMIP3 composite and any individual ice-sheet reconstruction are smaller than those between the CMIP5/PMIP3 composite and the ice sheet used in the last phase of PMIP (PMIP2).