The Met Office unified model global atmosphere 4.0 and JULES global land 4.0 configurations

We describe Global Atmosphere 4.0 (GA4.0) and Global Land 4.0 (GL4.0): configurations of the Met Office Unified Model and JULES (Joint UK Land Environment Simulator) community land surface model developed for use in global and regional climate research and weather prediction activities. GA4.0 and GL4.0 are based on the previous GA3.0 and GL3.0 configurations, with the inclusion of developments made by the Met Office and its collaborators during its annual development cycle. This paper provides a comprehensive technical and scientific description of GA4.0 and GL4.0 as well as details of how these differ from their predecessors. We also present the results of some initial evaluations of their performance. Overall, performance is comparable with that of GA3.0/GL3.0; the updated configurations include improvements to the science of several parametrisation schemes, however, and will form a baseline for further ongoing development.

Workaround motivation model (WAMM): an adaptation of the Theory of Interpersonal Behaviour

The study of workarounds (WA) has increased in importance due to their impact on patient safety and efficiency. However, there are no adequate theories to explain the motivation to create and use a workaround in a healthcare sitting. Although theories of technology acceptance help to understand the reasons to accept or reject technology, they fail to explain drivers for alternatives. Also workarounds involve creators and performers that have different motivations. Models such as Theory of Planned Behaviour (TPB) or Theory of Reasoned Action (TRA) can help to explain the role of workaround users, but lack explanation of workaround creators’ dynamics. Our aim is to develop a theoretical foundation to explain workaround motivation behaviour models with norms that relate to sanctions to provide an integrated Workaround Motivation Model; WAMM. The development of WAMM model is explained in this paper based on workaround cases as part of further research to establish the model.

Evaluation of boundary-layer type in a weather forecast model utilising long-term Doppler lidar observations

Many studies evaluating model boundary-layer schemes focus either on near-surface parameters or on short-term observational campaigns. This reflects the observational datasets that are widely available for use in model evaluation. In this paper we show how surface and long-term Doppler lidar observations, combined in a way to match model representation of the boundary layer as closely as possible, can be used to evaluate the skill of boundary-layer forecasts. We use a 2-year observational dataset from a rural site in the UK to evaluate a climatology of boundary layer type forecast by the UK Met Office Unified Model. In addition, we demonstrate the use of a binary skill score (Symmetric Extremal Dependence Index) to investigate the dependence of forecast skill on season, horizontal resolution and forecast leadtime. A clear diurnal and seasonal cycle can be seen in the climatology of both the model and observations, with the main discrepancies being the model overpredicting cumulus capped and decoupled stratocumulus capped boundary-layers and underpredicting well mixed boundary-layers. Using the SEDI skill score the model is most skillful at predicting the surface stability. The skill of the model in predicting cumulus capped and stratocumulus capped stable boundary layer forecasts is low but greater than a 24 hr persistence forecast. In contrast, the prediction of decoupled boundary-layers and boundary-layers with multiple cloud layers is lower than persistence. This process based evaluation approach has the potential to be applied to other boundary-layer parameterisation schemes with similar decision structures.

Evaluating adjusted forcing and model spread for historical and future scenarios in the CMIP5 generation of climate models

We utilize energy budget diagnostics from the Coupled Model Intercomparison Project phase 5 (CMIP5) to evaluate the models' climate forcing since preindustrial times employing an established regression technique. The climate forcing evaluated this way, termed the adjusted forcing (AF), includes a rapid adjustment term associated with cloud changes and other tropospheric and land-surface changes. We estimate a 2010 total anthropogenic and natural AF from CMIP5 models of 1.9 ± 0.9 W m−2 (5–95% range). The projected AF of the Representative Concentration Pathway simulations are lower than their expected radiative forcing (RF) in 2095 but agree well with efficacy weighted forcings from integrated assessment models. The smaller AF, compared to RF, is likely due to cloud adjustment. Multimodel time series of temperature change and AF from 1850 to 2100 have large intermodel spreads throughout the period. The intermodel spread of temperature change is principally driven by forcing differences in the present day and climate feedback differences in 2095, although forcing differences are still important for model spread at 2095. We find no significant relationship between the equilibrium climate sensitivity (ECS) of a model and its 2003 AF, in contrast to that found in older models where higher ECS models generally had less forcing. Given the large present-day model spread, there is no indication of any tendency by modelling groups to adjust their aerosol forcing in order to produce observed trends. Instead, some CMIP5 models have a relatively large positive forcing and overestimate the observed temperature change.

A numerical model of the ionospheric signatures of time-varying magneticreconnection: I. ionospheric convection

This paper presents a numerical model for predicting the evolution of the pattern of ionospheric convection in response to general time-dependent magnetic reconnection at the dayside magnetopause and in the cross-tail current sheet of the geomagnetic tail. The model quantifies the concepts of ionospheric flow excitation by Cowley and Lockwood (1992), assuming a uniform spatial distribution of ionospheric conductivity. The model is demonstrated using an example in which travelling reconnection pulses commence near noon and then move across the dayside magnetopause towards both dawn and dusk. Two such pulses, 8 min apart, are used and each causes the reconnection to be active for 1 min at every MLT that they pass over. This example demonstrates how the convection response to a given change in the interplanetary magnetic field (via the reconnection rate) depends on the previous reconnection history. The causes of this effect are explained. The inherent assumptions and the potential applications of the model are discussed.

An AeroCom assessment of black carbon in Arctic snow and sea ice

Though many global aerosols models prognose surface deposition, only a few models have been used to directly simulate the radiative effect from black carbon (BC) deposition to snow and sea ice. Here, we apply aerosol deposition fields from 25 models contributing to two phases of the Aerosol Comparisons between Observations and Models (AeroCom) project to simulate and evaluate within-snow BC concentrations and radiative effect in the Arctic. We accomplish this by driving the offline land and sea ice components of the Community Earth System Model with different deposition fields and meteorological conditions from 2004 to 2009, during which an extensive field campaign of BC measurements in Arctic snow occurred. We find that models generally underestimate BC concentrations in snow in northern Russia and Norway, while overestimating BC amounts elsewhere in the Arctic. Although simulated BC distributions in snow are poorly correlated with measurements, mean values are reasonable. The multi-model mean (range) bias in BC concentrations, sampled over the same grid cells, snow depths, and months of measurements, are −4.4 (−13.2 to +10.7) ng g−1 for an earlier phase of AeroCom models (phase I), and +4.1 (−13.0 to +21.4) ng g−1 for a more recent phase of AeroCom models (phase II), compared to the observational mean of 19.2 ng g−1. Factors determining model BC concentrations in Arctic snow include Arctic BC emissions, transport of extra-Arctic aerosols, precipitation, deposition efficiency of aerosols within the Arctic, and meltwater removal of particles in snow. Sensitivity studies show that the model–measurement evaluation is only weakly affected by meltwater scavenging efficiency because most measurements were conducted in non-melting snow. The Arctic (60–90° N) atmospheric residence time for BC in phase II models ranges from 3.7 to 23.2 days, implying large inter-model variation in local BC deposition efficiency. Combined with the fact that most Arctic BC deposition originates from extra-Arctic emissions, these results suggest that aerosol removal processes are a leading source of variation in model performance. The multi-model mean (full range) of Arctic radiative effect from BC in snow is 0.15 (0.07–0.25) W m−2 and 0.18 (0.06–0.28) W m−2 in phase I and phase II models, respectively. After correcting for model biases relative to observed BC concentrations in different regions of the Arctic, we obtain a multi-model mean Arctic radiative effect of 0.17 W m−2 for the combined AeroCom ensembles. Finally, there is a high correlation between modeled BC concentrations sampled over the observational sites and the Arctic as a whole, indicating that the field campaign provided a reasonable sample of the Arctic.

The Met Office Unified Model Global Atmosphere 3.0/3.1 and JULES Global Land 3.0/3.1 configurations

We describe Global Atmosphere 3.0 (GA3.0): a configuration of the Met Office Unified Model (MetUM) developed for use across climate research and weather prediction activities. GA3.0 has been formulated by converging the development paths of the Met Office's weather and climate global atmospheric model components such that wherever possible, atmospheric processes are modelled or parametrized seamlessly across spatial resolutions and timescales. This unified development process will provide the Met Office and its collaborators with regular releases of a configuration that has been evaluated, and can hence be applied, over a variety of modelling régimes. We also describe Global Land 3.0 (GL3.0): a configuration of the JULES community land surface model developed for use with GA3.0.

Evaluating and optimizing oral formulations of live bacterial vaccines using a gastro-small intestine model

Gastrointestinal (GI) models that mimic physiological conditions in vitro are important tools for developing and optimizing biopharmaceutical formulations. Oral administration of live attenuated bacterial vaccines (LBV) can safely and effectively promote mucosal immunity but new formulations are required that provide controlled release of optimal numbers of viable bacterial cells, which must survive gastrointestinal transit overcoming various antimicrobial barriers. Here, we use a gastro-small intestine gut model of human GI conditions to study the survival and release kinetics of two oral LBV formulations: the licensed typhoid fever vaccine Vivotif comprising enteric coated capsules; and an experimental formulation of the model vaccine Salmonella Typhimurium SL3261 dried directly onto cast enteric polymer films and laminated to form a polymer film laminate (PFL). Neither formulation released significant numbers of viable cells when tested in the complete gastro-small intestine model. The poor performance in delivering viable cells could be attributed to a combination of acid and bile toxicity plus incomplete release of cells for Vivotif capsules, and to bile toxicity alone for PFL. To achieve effective protection from intestinal bile in addition to effective acid resistance, bile adsorbent resins were incorporated into the PFL to produce a new formulation, termed BR-PFL. Efficient and complete release of 4.4x107 live cells per dose was achieved from BR-PFL at distal intestinal pH, with release kinetics controlled by the composition of the enteric polymer film, and no loss in viability observed in any stage of the GI model. Use of this in vitro GI model thereby allowed rational design of an oral LBV formulation to maximize viable cell release.

Schooling and the intergenerational transmission of values

We present a model of the evolution of identity via dynamic interaction between the choice of education and the transmission of values in a community from parents to children, when parents care about the preservation of their traditional community values, different from the values of the host society. We compare the educational and socioeconomic outcomes in different scenarios (melting pot versus multiculturalism). If schooling shifts children’s identity away from their parents’ values, parents may choose lower levels of education for their children, at the cost of reducing their future earnings. We show how this effect can be attenuated and reversed when the school or, indeed, the host society are willing to accommodate the values of the community and/or to adjust to these values; otherwise the community gradually becomes alienated. This approach may be applied to the analysis of temporal changes in values and attitudes in a community of immigrants, as well as ethnic, religious, or other minority groups.

An idealized study of coupled atmosphere-ocean 4D-Var in the presence of model error

Atmosphere only and ocean only variational data assimilation (DA) schemes are able to use window lengths that are optimal for the error growth rate, non-linearity and observation density of the respective systems. Typical window lengths are 6-12 hours for the atmosphere and 2-10 days for the ocean. However, in the implementation of coupled DA schemes it has been necessary to match the window length of the ocean to that of the atmosphere, which may potentially sacrifice the accuracy of the ocean analysis in order to provide a more balanced coupled state. This paper investigates how extending the window length in the presence of model error affects both the analysis of the coupled state and the initialized forecast when using coupled DA with differing degrees of coupling. Results are illustrated using an idealized single column model of the coupled atmosphere-ocean system. It is found that the analysis error from an uncoupled DA scheme can be smaller than that from a coupled analysis at the initial time, due to faster error growth in the coupled system. However, this does not necessarily lead to a more accurate forecast due to imbalances in the coupled state. Instead coupled DA is more able to update the initial state to reduce the impact of the model error on the accuracy of the forecast. The effect of model error is potentially most detrimental in the weakly coupled formulation due to the inconsistency between the coupled model used in the outer loop and uncoupled models used in the inner loop.

Use of a predictive model for food insecurity estimates in Brazil

In 2004 the National Household Survey (Pesquisa Nacional par Amostras de Domicilios - PNAD) estimated the prevalence of food and nutrition insecurity in Brazil. However, PNAD data cannot be disaggregated at the municipal level. The objective of this study was to build a statistical model to predict severe food insecurity for Brazilian municipalities based on the PNAD dataset. Exclusion criteria were: incomplete food security data (19.30%); informants younger than 18 years old (0.07%); collective households (0.05%); households headed by indigenous persons (0.19%). The modeling was carried out in three stages, beginning with the selection of variables related to food insecurity using univariate logistic regression. The variables chosen to construct the municipal estimates were selected from those included in PNAD as well as the 2000 Census. Multivariate logistic regression was then initiated, removing the non-significant variables with odds ratios adjusted by multiple logistic regression. The Wald Test was applied to check the significance of the coefficients in the logistic equation. The final model included the variables: per capita income; years of schooling; race and gender of the household head; urban or rural residence; access to public water supply; presence of children; total number of household inhabitants and state of residence. The adequacy of the model was tested using the Hosmer-Lemeshow test (p=0.561) and ROC curve (area=0.823). Tests indicated that the model has strong predictive power and can be used to determine household food insecurity in Brazilian municipalities, suggesting that similar predictive models may be useful tools in other Latin American countries.

A modified rat model of neonatal anoxia: Development and evaluation by pulseoximetry, arterial gasometry and Fos immunoreactivity

Neonatal anoxia is a worldwide clinical problem that has serious and lasting consequences. The diversity of models does not allow complete reproducibility, so a standardized model is needed. In this study, we developed a rat model of neonatal anoxia that utilizes a semi-hermetic system suitable for oxygen deprivation. The validity of this model was confirmed using pulse oximetry, arterial gasometry, observation of skin color and behavior and analysis of Fos immunoreactivity in brain regions that function in respiratory control. For these experiments, 87 male albino neonate rats (Rattus norvegicus, lineage Wistar) aged approximate 30 postnatal hours were divided into anoxia and control groups. The pups were kept in an euthanasia polycarbonate chamber at 36 +/- 1 degrees C, with continuous 100% nitrogen gas flow at 3 L/min and 101.7 kPa for 25 min. The peripheral arterial oxygen saturation of the anoxia group decreased 75% from its initial value. Decreased pH and partial pressure of oxygen and increased partial pressure of carbon dioxide were observed in this group, indicating metabolic acidosis, hypoxia and hypercapnia. respectively. Analysis of neuronal activation showed Fos immunoreactivity in the solitary tract nucleus, the lateral reticular nucleus and the area postrema, confirming that those conditions activated areas related to respiratory control in the nervous system. Therefore, the proposed model of neonatal anoxia allows standardization and precise control of the anoxic condition, which should be of great value in indentifying both the mechanisms underlying neonatal anoxia and novel therapeutic strategies to combat or prevent this widespread public health problem. (C) 2011 Elsevier B.V. All rights reserved.

A Time Efficient Optical Model for GATE Simulation of a LYSO Scintillation Matrix Used in PET Applications

A time efficient optical model is proposed for GATE simulation of a LYSO scintillation matrix coupled to a photomultiplier. The purpose is to avoid the excessively long computation time when activating the optical processes in GATE. The usefulness of the model is demonstrated by comparing the simulated and experimental energy spectra obtained with the dual planar head equipment for dosimetry with a positron emission tomograph ( DoPET). The procedure to apply the model is divided in two steps. Firstly, a simplified simulation of a single crystal element of DoPET is used to fit an analytic function that models the optical attenuation inside the crystal. In a second step, the model is employed to calculate the influence of this attenuation in the energy registered by the tomograph. The use of the proposed optical model is around three orders of magnitude faster than a GATE simulation with optical processes enabled. A good agreement was found between the experimental and simulated data using the optical model. The results indicate that optical interactions inside the crystal elements play an important role on the energy resolution and induce a considerable degradation of the spectra information acquired by DoPET. Finally, the same approach employed by the proposed optical model could be useful to simulate a scintillation matrix coupled to a photomultiplier using single or dual readout scheme.

Synthesis of NiO-MgO-ZrO(2) catalysts and their performance in reforming of model biogas

Catalysts containing NiO/MgO/ZrO(2) mixtures were synthesized by the polymerization method in a single step. They were characterized by X-ray diffraction (XRD), temperature programmed reduction (TPR) and physisorption of N(2) (BET) and then tested in the reforming of a model biogas (1.5CH4:1CO(2)) in the presence of air (1.5CH(4) + 1CO(2) + 0.25O(2)) at 750 degrees C for 6h. It was observed that the catalyst Ni20MZ performed better in catalytic processes than the well known catalysts, Ni/ZrO(2) and Ni/MgO, synthesized under the same conditions. The formation of solid solutions, MgO-ZrO(2) and NiO-MgO, increased the rate of conversion of reactants (CH(4) and CO(2)) into synthesis gas (H(2) + CO). The formation of oxygen vacancies (in samples containing ZrO(2) and MgO) seems to promote removal of the coke deposited on the nickel surface. The values of the H(2)/CO ratio were generally found to be slightly lower than stoichiometric, owing to the reverse water gas shift reaction occurring in parallel. (C) 2011 Elsevier B.V. All rights reserved.

Chinese Basic Pension Substitution Rate: A Monte Carlo Demonstration of the Individual Account Model

At the end of 2005, the State Council of China passed ”The Decision on adjusting the Individual Account of Basic Pension System”, which adjusted the individual account in the 1997 basic pension system. In this essay, we will analyze the adjustment above, and use Life Annuity Actuarial Theory to establish the basic pension substitution rate model. Monte Carlo simulation is also used to prove the rationality of the model. Some suggestions are put forward associated with the substitution rate according to the current policy.