Results of an aqueous source term model for a radiological risk assessment of the Drigg LLW Site, U.K.

A radionuclide source term model has been developed which simulates the biogeochemical evolution of the Drigg low level waste (LLW) disposal site. The DRINK (DRIgg Near field Kinetic) model provides data regarding radionuclide concentrations in groundwater over a period of 100,000 years, which are used as input to assessment calculations for a groundwater pathway. The DRINK model also provides input to human intrusion and gaseous assessment calculations through simulation of the solid radionuclide inventory. These calculations are being used to support the Drigg post closure safety case. The DRINK model considers the coupled interaction of the effects of fluid flow, microbiology, corrosion, chemical reaction, sorption and radioactive decay. It represents the first direct use of a mechanistic reaction-transport model in risk assessment calculations.

Consistent Dirichlet Boundary Conditions for Numerical Solution of Moving Boundary Problems

We consider the imposition of Dirichlet boundary conditions in the finite element modelling of moving boundary problems in one and two dimensions for which the total mass is prescribed. A modification of the standard linear finite element test space allows the boundary conditions to be imposed strongly whilst simultaneously conserving a discrete mass. The validity of the technique is assessed for a specific moving mesh finite element method, although the approach is more general. Numerical comparisons are carried out for mass-conserving solutions of the porous medium equation with Dirichlet boundary conditions and for a moving boundary problem with a source term and time-varying mass.

Application of the Lighthill-Ford theory of spontaneous imbalance to clear-air turbulence forecasting

A new method of clear-air turbulence (CAT) forecasting based on the Lighthill–Ford theory of spontaneous imbalance and emission of inertia–gravity waves has been derived and applied on episodic and seasonal time scales. A scale analysis of this shallow-water theory for midlatitude synoptic-scale flows identifies advection of relative vorticity as the leading-order source term. Examination of leading- and second-order terms elucidates previous, more empirically inspired CAT forecast diagnostics. Application of the Lighthill–Ford theory to the Upper Mississippi and Ohio Valleys CAT outbreak of 9 March 2006 results in good agreement with pilot reports of turbulence. Application of Lighthill–Ford theory to CAT forecasting for the 3 November 2005–26 March 2006 period using 1-h forecasts of the Rapid Update Cycle (RUC) 2 1500 UTC model run leads to superior forecasts compared to the current operational version of the Graphical Turbulence Guidance (GTG1) algorithm, the most skillful operational CAT forecasting method in existence. The results suggest that major improvements in CAT forecasting could result if the methods presented herein become operational.

On the generation mechanisms of short-scale unbalanced modes in rotating two-layer flows with vertical shear

We report on the results of a laboratory investigation using a rotating two-layer annulus experiment, which exhibits both large-scale vortical modes and short-scale divergent modes. A sophisticated visualization method allows us to observe the flow at very high spatial and temporal resolution. The balanced long-wavelength modes appear only when the Froude number is supercritical (i.e. $F\,{>}\,F_\mathrm{critical}\,{\equiv}\, \upi^2/2$), and are therefore consistent with generation by a baroclinic instability. The unbalanced short-wavelength modes appear locally in every single baroclinically unstable flow, providing perhaps the first direct experimental evidence that all evolving vortical flows will tend to emit freely propagating inertia–gravity waves. The short-wavelength modes also appear in certain baroclinically stable flows. We infer the generation mechanisms of the short-scale waves, both for the baro-clinically unstable case in which they co-exist with a large-scale wave, and for the baroclinically stable case in which they exist alone. The two possible mechanisms considered are spontaneous adjustment of the large-scale flow, and Kelvin–Helmholtz shear instability. Short modes in the baroclinically stable regime are generated only when the Richardson number is subcritical (i.e. $\hbox{\it Ri}\,{<}\,\hbox{\it Ri}_\mathrm{critical}\,{\equiv}\, 1$), and are therefore consistent with generation by a Kelvin–Helmholtz instability. We calculate five indicators of short-wave generation in the baroclinically unstable regime, using data from a quasi-geostrophic numerical model of the annulus. There is excellent agreement between the spatial locations of short-wave emission observed in the laboratory, and regions in which the model Lighthill/Ford inertia–gravity wave source term is large. We infer that the short waves in the baroclinically unstable fluid are freely propagating inertia–gravity waves generated by spontaneous adjustment of the large-scale flow.

Performance assessment of a volcanic ash transport model mini-ensemble used for inverse modeling of the 2010 Eyjafjallajökull eruption

The requirement to forecast volcanic ash concentrations was amplified as a response to the 2010 Eyjafjallajökull eruption when ash safety limits for aviation were introduced in the European area. The ability to provide accurate quantitative forecasts relies to a large extent on the source term which is the emissions of ash as a function of time and height. This study presents source term estimations of the ash emissions from the Eyjafjallajökull eruption derived with an inversion algorithm which constrains modeled ash emissions with satellite observations of volcanic ash. The algorithm is tested with input from two different dispersion models, run on three different meteorological input data sets. The results are robust to which dispersion model and meteorological data are used. Modeled ash concentrations are compared quantitatively to independent measurements from three different research aircraft and one surface measurement station. These comparisons show that the models perform reasonably well in simulating the ash concentrations, and simulations using the source term obtained from the inversion are in overall better agreement with the observations (rank correlation = 0.55, Figure of Merit in Time (FMT) = 25–46%) than simulations using simplified source terms (rank correlation = 0.21, FMT = 20–35%). The vertical structures of the modeled ash clouds mostly agree with lidar observations, and the modeled ash particle size distributions agree reasonably well with observed size distributions. There are occasionally large differences between simulations but the model mean usually outperforms any individual model. The results emphasize the benefits of using an ensemble-based forecast for improved quantification of uncertainties in future ash crises.

The heliospheric Hale cycle over the last 300 years and its implications for a “lost” late 18th century solar cycle

A Hale cycle, one complete magnetic cycle of the Sun, spans two complete Schwabe cycles (also referred to as sunspot and, more generally, solar cycles). The approximately 22-year Hale cycle is seen in magnetic polarities of both sunspots and polar fields, as well as in the intensity of galactic cosmic rays reaching Earth, with odd- and even-numbered solar cycles displaying qualitatively different waveforms. Correct numbering of solar cycles also underpins empirical cycle-to-cycle relations which are used as first-order tests of stellar dynamo models. There has been much debate about whether the unusually long solar cycle 4 (SC4), spanning- 1784–1799, was actually two shorter solar cycles combined as a result of poor data coverage in the original Wolf sunspot number record. Indeed, the group sunspot number does show a small increase around 1794–1799 and there is evidence of an increase in the mean latitude of sunspots at this time, suggesting the existence of a cycle ‘‘4b’’. In this study, we use cosmogenic radionuclide data and associated reconstructions of the heliospheric magnetic field (HMF) to show that the Hale cycle has persisted over the last 300 years and that data prior to 1800 are more consistent with cycle 4 being a single long cycle (the ‘‘no SC4b’’ scenario). We also investigate the effect of cycle 4b on the HMF using an open solar flux (OSF) continuity model, in which the OSF source term is related to sunspot number and the OSF loss term is determined by the heliospheric current sheet tilt, assumed to be a simple function of solar cycle phase. The results are surprising; Without SC4b, the HMF shows two distinct peaks in the 1784–1799 interval, while the addition of SC4b removes the secondary peak, as the OSF loss term acts in opposition to the later rise in sunspot number. The timing and magnitude of the main SC4 HMF peak is also significantly changed by the addition of SC4b. These results are compared with the cosmogenic isotope reconstructions of HMF and historical aurora records. These data marginally favour the existence of SC4b (the ‘‘SC4b’’ scenario), though the result is less certain than that based on the persistence of the Hale cycle. Thus while the current uncertainties in the observations preclude any definitive conclusions, the data favour the ‘‘no SC4b’’ scenario. Future improvements to cosmogenic isotope reconstructions of the HMF, through either improved modelling or additional ice cores from well-separated geographic locations, may enable questions of the existence of SC4b and the phase of Hale cycle prior to the Maunder minimum to be settled conclusively.

Investigation of the impact of anthropogenic heat flux within an urban land surface model and PILPS-urban

Results from the first international urban model comparison experiment (PILPS-Urban) suggested that models which neglected the anthropogenic heat flux within the surface energy balance performed at least as well as models that include the source term, but this could not be explained. The analyses undertaken show that the results from PILPS-Urban were masked by the signal from including vegetation, which was identified in PILPS-Urban as being important. Including the anthropogenic heat flux does give improved performance, but the benefit is small for the site studied given the relatively small magnitude of this flux relative to other terms in the surface energy balance. However, there is no further benefit from including temporal variations in the flux at this site. The importance is expected to increase at sites with a larger anthropogenic heat flux and greater temporal variations.

Data insertion in volcanic ash cloud forecasting

During the eruption of Eyjafjallajökull in April and May 2010, the London Volcanic Ash  Advisory Centre demonstrated the importance of infrared (IR) satellite imagery for monitoring volcanic ash and validating the Met Office operational model, NAME. This model is used to forecast ash dispersion and forms much of the basis of the advice given to civil aviation. NAME requires a source term describing the properties of the eruption plume at the volcanic source. Elements of the source term are often highly uncertain and significant effort has therefore been invested into the use of satellite observations of ash clouds to constrain them. This paper presents a data insertion method, where satellite observations of downwind ash clouds are used to create effective ‘virtual sources’ far from the vent. Uncertainty in the model output is known to increase over the duration of a model run, as inaccuracies in the source term, meteorological data and the parameterizations of the   modelled processes accumulate. This new technique, where the dispersion model (DM) is ‘reinitialized’ part-­way through a run, could go some way to addressing this.

Long-term changes in the solar photosphere associated with changes in the coronal source flux

Using sunspot observations from Greenwich and Mount Wilson, we show that the latitudinal spread of sunspot groups has increased since 1874, in a manner that closely mirrors the long-term (similar to 100 year) changes in the coronal source flux, F-s, as inferred from geomagnetic activity. This latitude spread is shown to be well correlated with the flux emergence rate required by the model of the coronal source flux variation by Solanki er al. [2000]. The time constant for the decay of this open flux is found to be 3.6 +/-0.8 years. Using this value, and quantifying the photospheric flux emergence rate using the latitudinal spread of sunspot groups, the model reproduces the observed coronal source flux variation. The ratio of the 100-year drift to the solar cycle amplitude for the flux emergence rate is found to be half of the same ratio for F-s.

Long-term drift of the coronal source magnetic flux and the total solar irradiance

We test the method of Lockwood et al. [1999] for deriving the coronal source flux from the geomagnetic aa index and show it to be accurate to within 12% for annual means and 4.5% for averages over a sunspot cycle. Using data from four solar constant monitors during 1981-1995, we find a linear relationship between this magnetic flux and the total solar irradiance. From this correlation, we show that the 131% rise in the mean coronal source field over the interval 1901-1995 corresponds to a rise in the average total solar irradiance of {\Delta}I = 1.65 +/- 0.23 Wm^{-2}.

Foods derived from animals: the impact of animal nutrition on their nutritive value and ability to sustain long-term health

Foods derived from domestic animals are a significant source of nutrients in the UK diet. However, certain aspects of some animal-derived foods, notably levels of saturated fatty acids, have given rise to concerns that these foods may contribute to the risk of cardiovascular disease, the metabolic syndrome and other conditions. However, the composition of the many animal-derived foods is not constant and can often be enhanced by manipulating the nutrition of the animal. This paper reviews these possibilities with particular attention to lipids, and draws attention to the fact that milk in particular, contains a number of compounds which may, for example, exert anti-carcinogenic effects. It is clear that the role of animal nutrition in creating foods closer to the optimum composition for long-term human health will not only be more relevant in the future, but will be vital in attempts to improve the health of the human population.

Effect of long-term olive oil dietary intervention on postprandial triacylglycerol and factor VII metabolism

Although the beneficial effects of Mediterranean-type diets, which are rich in olive oil, a good source of monounsaturated fatty acids (MUFAs), are generally accepted, little is known about the effects of long-term dietary MUFA intake on postprandial lipoprotein metabolism and hemostasis. This study used a single-blind, randomized, crossover design to investigate the relative effects of a long-term dietary olive oil intervention and a control [saturated fatty acid (SFA)-enriched] diet on postprandial triacylglycerol metabolism and factor VII activity. The postprandial response to a standard test meal was investigated in 23 healthy men who adhered to both diets for 8 wk. cis-MUFAs were successfully substituted for SFAs in the MUFA diet without affecting total dietary fat or energy intakes. The long-term dietary MUFA intervention significantly reduced plasma and LDL-cholesterol concentrations (P = 0.01). Postprandial triacylglycerol concentrations were significantly greater in the early postprandial period after the MUFA diet (P = 0.003). Postprandial factor VII activation and the concentration of the factor VII antigen were significantly lower after the MUFA diet (P = 0.04 and P = 0 006, respectively). This study showed that isoenergetic substitution of MUFAs for SFAs reduces plasma cholesterol and reduces the degree of postprandial factor VII activation. The alterations in the postprandial triacylglycerol response suggest a greater rate of dietary fat absorption and postprandial triacylglycerol metabolism after a diet rich in MUFAs. This study presents new insights into the biochemical basis of the beneficial effects associated with long-term dietary MUFA consumption, which may explain the lower rates of coronary mortality in Mediterranean regions.

Model inter-comparison between statistical and dynamic model assessments of the long-term stability of blanket peat in Great Britain (1940-2099)

We compared output from 3 dynamic process-based models (DMs: ECOSSE, MILLENNIA and the Durham Carbon Model) and 9 bioclimatic envelope models (BCEMs; including BBOG ensemble and PEATSTASH) ranging from simple threshold to semi-process-based models. Model simulations were run at 4 British peatland sites using historical climate data and climate projections under a medium (A1B) emissions scenario from the 11-RCM (regional climate model) ensemble underpinning UKCP09. The models showed that blanket peatlands are vulnerable to projected climate change; however, predictions varied between models as well as between sites. All BCEMs predicted a shift from presence to absence of a climate associated with blanket peat, where the sites with the lowest total annual precipitation were closest to the presence/absence threshold. DMs showed a more variable response. ECOSSE predicted a decline in net C sink and shift to net C source by the end of this century. The Durham Carbon Model predicted a smaller decline in the net C sink strength, but no shift to net C source. MILLENNIA predicted a slight overall increase in the net C sink. In contrast to the BCEM projections, the DMs predicted that the sites with coolest temperatures and greatest total annual precipitation showed the largest change in carbon sinks. In this model inter-comparison, the greatest variation in model output in response to climate change projections was not between the BCEMs and DMs but between the DMs themselves, because of different approaches to modelling soil organic matter pools and decomposition amongst other processes. The difference in the sign of the response has major implications for future climate feedbacks, climate policy and peatland management. Enhanced data collection, in particular monitoring peatland response to current change, would significantly improve model development and projections of future change.

Aerosol contribution to the rapid warming of near-term climate under RCP 2.6

The importance of aerosol emissions for near term climate projections is investigated by analysing simulations with the HadGEM2-ES model under two different emissions scenarios: RCP2.6 and RCP4.5. It is shown that the near term warming projected under RCP2.6 is greater than under RCP4.5, even though the greenhouse gas forcing is lower. Rapid and substantial reductions in sulphate aerosol emissions due to a reduction of coal burning in RCP2.6 lead to a reduction in the negative shortwave forcing due to aerosol direct and indirect effects. Indirect effects play an important role over the northern hemisphere oceans, especially the subtropical northeastern Pacific where an anomaly of 5-10\,Wm$^{-2}$ develops. The pattern of surface temperature change is consistent with the expected response to this surface radiation anomaly, whilst also exhibiting features that reflect redistribution of energy, and feedbacks, within the climate system. These results demonstrate the importance of aerosol emissions as a key source of uncertainty in near term projections of global and regional climate.

Global and regional temperature-change potentials for near-term climate forcers

We examine the climate effects of the emissions of near-term climate forcers (NTCFs) from 4 continental regions (East Asia, Europe, North America and South Asia) using radiative forcing from the task force on hemispheric transport of air pollution source-receptor global chemical transport model simulations. These simulations model the transport of 3 aerosol species (sulphate, particulate organic matter and black carbon) and 4 ozone precursors (methane, nitric oxides (NOx), volatile organic compounds and carbon monoxide). From the equilibrium radiative forcing results we calculate global climate metrics, global warming potentials (GWPs) and global temperature change potentials (GTPs) and show how these depend on emission region, and can vary as functions of time. For the aerosol species, the GWP(100) values are −37±12, −46±20, and 350±200 for SO2, POM and BC respectively for the direct effects only. The corresponding GTP(100) values are −5.2±2.4, −6.5±3.5, and 50±33. This analysis is further extended by examining the temperature-change impacts in 4 latitude bands. This shows that the latitudinal pattern of the temperature response to emissions of the NTCFs does not directly follow the pattern of the diagnosed radiative forcing. For instance temperatures in the Arctic latitudes are particularly sensitive to NTCF emissions in the northern mid-latitudes. At the 100-yr time horizon the ARTPs show NOx emissions can have a warming effect in the northern mid and high latitudes, but cooling in the tropics and Southern Hemisphere. The northern mid-latitude temperature response to northern mid-latitude emissions of most NTCFs is approximately twice as large as would be implied by the global average.