Towards convection-resolving, global atmospheric simulations with the Model for Prediction Across Scales (MPAS) v3.1: an extreme scaling experiment, links to sourcecode of model and grides

The Model for Prediction Across Scales (MPAS) is a novel set of Earth system simulation components and consists of an atmospheric model, an ocean model and a land-ice model. Its distinct features are the use of unstructured Voronoi meshes and C-grid discretisation to address shortcomings of global models on regular grids and the use of limited area models nested in a forcing data set, with respect to parallel scalability, numerical accuracy and physical consistency. This concept allows one to include the feedback of regional land use information on weather and climate at local and global scales in a consistent way, which is impossible to achieve with traditional limited area modelling approaches. Here, we present an in-depth evaluation of MPAS with regards to technical aspects of performing model runs and scalability for three medium-size meshes on four different high-performance computing (HPC) sites with different architectures and compilers. We uncover model limitations and identify new aspects for the model optimisation that are introduced by the use of unstructured Voronoi meshes. We further demonstrate the model performance of MPAS in terms of its capability to reproduce the dynamics of the West African monsoon (WAM) and its associated precipitation in a pilot study. Constrained by available computational resources, we compare 11-month runs for two meshes with observations and a reference simulation from the Weather Research and Forecasting (WRF) model. We show that MPAS can reproduce the atmospheric dynamics on global and local scales in this experiment, but identify a precipitation excess for the West African region. Finally, we conduct extreme scaling tests on a global 3?km mesh with more than 65 million horizontal grid cells on up to half a million cores. We discuss necessary modifications of the model code to improve its parallel performance in general and specific to the HPC environment. We confirm good scaling (70?% parallel efficiency or better) of the MPAS model and provide numbers on the computational requirements for experiments with the 3?km mesh. In doing so, we show that global, convection-resolving atmospheric simulations with MPAS are within reach of current and next generations of high-end computing facilities.

Age model and chemistry of sediment core STK7 from Steisslinger See

Stable isotope records for carbon and oxygen in bulk carbonates, carbon in bulk organic matter, and for total and chromium-reducible sulfur in a lacustrine sediment core from Lake Steisslingen (Southwest Germany) show several distinct and abrupt shifts during the last 15,000 years. Variations in the isotopic composition of authigenic carbonates indicate two major phases in the lake history. In the pre-Holocene, the hydrological budget of the lake was apparently stable. Variations of delta18O values of authigenic carbonates were, therefore, dominantly controlled by temperature changes. A decrease in the delta18Ocarb values of about 2 per mil at the Allerød/Younger Dryas transition is interpreted as a drop in mean annual air temperatures of approximately 5°C. An abrupt temperature increase of a similar magnitude is inferred at the Younger Dryas/Preboreal boundary. Throughout most of the Holocene, the isotopic composition of authigenic carbonates was influenced by marked changes in the hydrological budget of the lake. A major positive excursion in the delta13Ccarb and delta18Ocarb values at the beginning of the Atlantic and a smaller one in the Preboreal were related to evaporation effects, which indicate that dry climatic conditions must have prevailed at that time. A simultaneous increase in delta13C values of bulk organic matter at the beginning of the Atlantic suggests a high level of productivity in the lake. As a consequence, aqueous sulfate became limited as indicated by variations in the delta34S values of total and chromium-reducible sedimentary sulfur. Therefore, we conclude that the beginning of the Atlantic was characterized not only by dry but also by warm climatic conditions, which triggered a higher productivity in the lake. In the Subatlantic sediments, large variations in carbon, oxygen, and sulfur isotope ratios were observed as a result of human activities, causing considerable perturbations in the biogeochemical element cycling of Lake Steisslingen. Results obtained by the study of the continuous 15 ka record of Lake Steisslingen document clearly that isotopic proxy data from lacustrine sediments can provide useful information on environmental and climatic changes of local, regional, and in the case of the Younger Dryas event, of even hemispherical significance.

Holocene precipitation change in different monsoon sub-regions (time-slices and transient data) simulated by different global climate models, with links to model results

The recently proposed global monsoon hypothesis interprets monsoon systems as part of one global-scale atmospheric overturning circulation, implying a connection between the regional monsoon systems and an in-phase behaviour of all northern hemispheric monsoons on annual timescales (Trenberth et al., 2000). Whether this concept can be applied to past climates and variability on longer timescales is still under debate, because the monsoon systems exhibit different regional characteristics such as different seasonality (i.e. onset, peak, and withdrawal). To investigate the interconnection of different monsoon systems during the pre-industrial Holocene, five transient global climate model simulations have been analysed with respect to the rainfall trend and variability in different sub-domains of the Afro-Asian monsoon region. Our analysis suggests that on millennial timescales with varying orbital forcing, the monsoons do not behave as a tightly connected global system. According to the models, the Indian and North African monsoons are coupled, showing similar rainfall trend and moderate correlation in rainfall variability in all models. The East Asian monsoon changes independently during the Holocene. The dissimilarities in the seasonality of the monsoon sub-systems lead to a stronger response of the North African and Indian monsoon systems to the Holocene insolation forcing than of the East Asian monsoon and affect the seasonal distribution of Holocene rainfall variations. Within the Indian and North African monsoon domain, precipitation solely changes during the summer months, showing a decreasing Holocene precipitation trend. In the East Asian monsoon region, the precipitation signal is determined by an increasing precipitation trend during spring and a decreasing precipitation change during summer, partly balancing each other. A synthesis of reconstructions and the model results do not reveal an impact of the different seasonality on the timing of the Holocene rainfall optimum in the different sub-monsoon systems. They rather indicate locally inhomogeneous rainfall changes and show, that single palaeo-records should not be used to characterise the rainfall change and monsoon evolution for entire monsoon sub-systems.

Face recognition using a unified 3D morphable model

We address the problem of 3D-assisted 2D face recognition in scenarios when the input image is subject to degradations or exhibits intra-personal variations not captured by the 3D model. The proposed solution involves a novel approach to learn a subspace spanned by perturbations caused by the missing modes of variation and image degradations, using 3D face data reconstructed from 2D images rather than 3D capture. This is accomplished by modelling the difference in the texture map of the 3D aligned input and reference images. A training set of these texture maps then defines a perturbation space which can be represented using PCA bases. Assuming that the image perturbation subspace is orthogonal to the 3D face model space, then these additive components can be recovered from an unseen input image, resulting in an improved fit of the 3D face model. The linearity of the model leads to efficient fitting. Experiments show that our method achieves very competitive face recognition performance on Multi-PIE and AR databases. We also present baseline face recognition results on a new data set exhibiting combined pose and illumination variations as well as occlusion.

A framework for experimental determination of localised vertical pedestrian forces on full-scale structures using wireless attitude and heading reference systems

A major weakness among loading models for pedestrians walking on flexible structures proposed in recent years is the various uncorroborated assumptions made in their development. This applies to spatio-temporal characteristics of pedestrian loading and the nature of multi-object interactions. To alleviate this problem, a framework for the determination of localised pedestrian forces on full-scale structures is presented using a wireless attitude and heading reference systems (AHRS). An AHRS comprises a triad of tri-axial accelerometers, gyroscopes and magnetometers managed by a dedicated data processing unit, allowing motion in three-dimensional space to be reconstructed. A pedestrian loading model based on a single point inertial measurement from an AHRS is derived and shown to perform well against benchmark data collected on an instrumented treadmill. Unlike other models, the current model does not take any predefined form nor does it require any extrapolations as to the timing and amplitude of pedestrian loading. In order to assess correctly the influence of the moving pedestrian on behaviour of a structure, an algorithm for tracking the point of application of pedestrian force is developed based on data from a single AHRS attached to a foot. A set of controlled walking tests with a single pedestrian is conducted on a real footbridge for validation purposes. A remarkably good match between the measured and simulated bridge response is found, indeed confirming applicability of the proposed framework.

Generation and Characterisation of a Canine EGFP-HMGA2 Prostate Cancer In Vitro Model

The architectural transcription factor HMGA2 is abundantly expressed during embryonic development. In several malignant neoplasias including prostate cancer, high re-expression of HMGA2 is correlated with malignancy and poor prognosis. The let-7 miRNA family is described to regulate HMGA2 negatively. The balance of let-7 and HMGA2 is discussed to play a major role in tumour aetiology. To further analyse the role of HMGA2 in prostate cancer a stable and highly reproducible in vitro model system is precondition. Herein we established a canine CT1258-EGFP-HMGA2 prostate cancer cell line stably overexpressing HMGA2 linked to EGFP and in addition the reference cell line CT1258-EGFP expressing solely EGFP to exclude EGFP-induced effects. Both recombinant cell lines were characterised by fluorescence microscopy, flow cytometry and immunocytochemistry. The proliferative effect of ectopically overexpressed HMGA2 was determined via BrdU assays. Comparative karyotyping of the derived and the initial CT1258 cell lines was performed to analyse chromosome consistency. The impact of the ectopic HMGA2 expression on its regulator let-7a was analysed by quantitative real-time PCR. Fluorescence microscopy and immunocytochemistry detected successful expression of the EGFP-HMGA2 fusion protein exclusively accumulating in the nucleus. Gene expression analyses confirmed HMGA2 overexpression in CT1258-EGFP-HMGA2 in comparison to CT1258-EGFP and native cells. Significantly higher let-7a expression levels were found in CT1258-EGFP-HMGA2 and CT1258-EGFP. The BrdU assays detected an increased proliferation of CT1258-HMGA2-EGFP cells compared to CT1258-EGFP and native CT1258. The cytogenetic analyses of CT1258-EGFP and CT1258-EGFP-HMGA2 resulted in a comparable hyperdiploid karyotype as described for native CT1258 cells. To further investigate the impact of recombinant overexpressed HMGA2 on CT1258 cells, other selected targets described to underlie HMGA2 regulation were screened in addition. The new fluorescent CT1258-EGFP-HMGA2 cell line is a stable tool enabling in vitro and in vivo analyses of the HMGA2-mediated effects on cells and the development and pathogenesis of prostate cancer.

Quantification of the effects of architectural traits on dry mass production and light interception of tomato canopy under different temperature regimes using a dynamic functional-structural plant model

There is increasing interest in evaluating the environmental effects on crop architectural traits and yield improvement. However, crop models describing the dynamic changes in canopy structure with environmental conditions and the complex interactions between canopy structure, light interception, and dry mass production are only gradually emerging. Using tomato (Solanum lycopersicum L.) as a model crop, a dynamic functional-structural plant model (FSPM) was constructed, parameterized, and evaluated to analyse the effects of temperature on architectural traits, which strongly influence canopy light interception and shoot dry mass. The FSPM predicted the organ growth, organ size, and shoot dry mass over time with high accuracy (>85%). Analyses of this FSPM showed that, in comparison with the reference canopy, shoot dry mass may be affected by leaf angle by as much as 20%, leaf curvature by up to 7%, the leaf length: width ratio by up to 5%, internode length by up to 9%, and curvature ratios and leaf arrangement by up to 6%. Tomato canopies at low temperature had higher canopy density and were more clumped due to higher leaf area and shorter internodes. Interestingly, dry mass production and light interception of the clumped canopy were more sensitive to changes in architectural traits. The complex interactions between architectural traits, canopy light interception, dry mass production, and environmental conditions can be studied by the dynamic FSPM, which may serve as a tool for designing a canopy structure which is 'ideal' in a given environment.

Estimation of whole body muscle, adipose/fat mass, validation in health and during weight loss. Development of prediction equations using MRI as reference method

Background: Body composition is affected by diseases, and affects responses to medical treatments, dosage of medicines, etc., while an abnormal body composition contributes to the causation of many chronic diseases. While we have reliable biochemical tests for certain nutritional parameters of body composition, such as iron or iodine status, and we have harnessed nuclear physics to estimate the body’s content of trace elements, the very basic quantification of body fat content and muscle mass remains highly problematic. Both body fat and muscle mass are vitally important, as they have opposing influences on chronic disease, but they have seldom been estimated as part of population health surveillance. Instead, most national surveys have merely reported BMI and waist, or sometimes the waist/hip ratio; these indices are convenient but do not have any specific biological meaning. Anthropometry offers a practical and inexpensive method for muscle and fat estimation in clinical and epidemiological settings; however, its use is imperfect due to many limitations, such as a shortage of reference data, misuse of terminology, unclear assumptions, and the absence of properly validated anthropometric equations. To date, anthropometric methods are not sensitive enough to detect muscle and fat loss. Aims: The aim of this thesis is to estimate Adipose/fat and muscle mass in health disease and during weight loss through; 1. evaluating and critiquing the literature, to identify the best-published prediction equations for adipose/fat and muscle mass estimation; 2. to derive and validate adipose tissue and muscle mass prediction equations; and 3.to evaluate the prediction equations along with anthropometric indices and the best equations retrieved from the literature in health, metabolic illness and during weight loss. Methods: a Systematic review using Cochrane Review method was used for reviewing muscle mass estimation papers that used MRI as the reference method. Fat mass estimation papers were critically reviewed. Mixed ethnic, age and body mass data that underwent whole body magnetic resonance imaging to quantify adipose tissue and muscle mass (dependent variable) and anthropometry (independent variable) were used in the derivation/validation analysis. Multiple regression and Bland-Altman plot were applied to evaluate the prediction equations. To determine how well the equations identify metabolic illness, English and Scottish health surveys were studied. Statistical analysis using multiple regression and binary logistic regression were applied to assess model fit and associations. Also, populations were divided into quintiles and relative risk was analysed. Finally, the prediction equations were evaluated by applying them to a pilot study of 10 subjects who underwent whole-body MRI, anthropometric measurements and muscle strength before and after weight loss to determine how well the equations identify adipose/fat mass and muscle mass change. Results: The estimation of fat mass has serious problems. Despite advances in technology and science, prediction equations for the estimation of fat mass depend on limited historical reference data and remain dependent upon assumptions that have not yet been properly validated for different population groups. Muscle mass does not have the same conceptual problems; however, its measurement is still problematic and reference data are scarce. The derivation and validation analysis in this thesis was satisfactory, compared to prediction equations in the literature they were similar or even better. Applying the prediction equations in metabolic illness and during weight loss presented an understanding on how well the equations identify metabolic illness showing significant associations with diabetes, hypertension, HbA1c and blood pressure. And moderate to high correlations with MRI-measured adipose tissue and muscle mass before and after weight loss. Conclusion: Adipose tissue mass and to an extent muscle mass can now be estimated for many purposes as population or groups means. However, these equations must not be used for assessing fatness and categorising individuals. Further exploration in different populations and health surveys would be valuable.

Experimental and diagnostic protocol for the physical component of the CMIP6 Ocean Model Intercomparison Project (OMIP)

The Ocean Model Intercomparison Project (OMIP) aims to provide a framework for evaluating, understanding, and improving the ocean and sea-ice components of global climate and earth system models contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses these aims in two complementary manners: (A) by providing an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing, (B) by providing a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) offering details for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows that of the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II have become the standard method to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP (Scenario MIP), as well as the ocean-sea ice OMIP simulations. The bulk of this paper offers scientific rationale for saving these diagnostics.

OMIP contribution to CMIP6: experimental and diagnostic protocol for the physical component of the Ocean Model Intercomparison Project

The Ocean Model Intercomparison Project (OMIP) is an endorsed project in the Coupled Model Intercomparison Project Phase 6 (CMIP6). OMIP addresses CMIP6 science questions, investigating the origins and consequences of systematic model biases. It does so by providing a framework for evaluating (including assessment of systematic biases), understanding, and improving ocean, sea-ice, tracer, and biogeochemical components of climate and earth system models contributing to CMIP6. Among the WCRP Grand Challenges in climate science (GCs), OMIP primarily contributes to the regional sea level change and near-term (climate/decadal) prediction GCs. OMIP provides (a) an experimental protocol for global ocean/sea-ice models run with a prescribed atmospheric forcing; and (b) a protocol for ocean diagnostics to be saved as part of CMIP6. We focus here on the physical component of OMIP, with a companion paper (Orr et al., 2016) detailing methods for the inert chemistry and interactive biogeochemistry. The physical portion of the OMIP experimental protocol follows the interannual Coordinated Ocean-ice Reference Experiments (CORE-II). Since 2009, CORE-I (Normal Year Forcing) and CORE-II (Interannual Forcing) have become the standard methods to evaluate global ocean/sea-ice simulations and to examine mechanisms for forced ocean climate variability. The OMIP diagnostic protocol is relevant for any ocean model component of CMIP6, including the DECK (Diagnostic, Evaluation and Characterization of Klima experiments), historical simulations, FAFMIP (Flux Anomaly Forced MIP), C4MIP (Coupled Carbon Cycle Climate MIP), DAMIP (Detection and Attribution MIP), DCPP (Decadal Climate Prediction Project), ScenarioMIP, HighResMIP (High Resolution MIP), as well as the ocean/sea-ice OMIP simulations.

Health beliefs affect the correct replacement of daily disposable contact lenses:predicting compliance with the Health Belief Model and the Theory of Planned Behaviour

Purpose: To assess the compliance of Daily Disposable Contact Lenses (DDCLs) wearers with replacing lenses at a manufacturer-recommended replacement frequency. To evaluate the ability of two different Health Behavioural Theories (HBT), The Health Belief Model (HBM) and The Theory of Planned Behaviour (TPB), in predicting compliance. Method: A multi-centre survey was conducted using a questionnaire completed anonymously by contact lens wearers during the purchase of DDCLs. Results: Three hundred and fifty-four questionnaires were returned. The survey comprised 58.5% females and 41.5% males (mean age 34. ±. 12. years). Twenty-three percent of respondents were non-compliant with manufacturer-recommended replacement frequency (re-using DDCLs at least once). The main reason for re-using DDCLs was "to save money" (35%). Predictions of compliance behaviour (past behaviour or future intentions) on the basis of the two HBT was investigated through logistic regression analysis: both TPB factors (subjective norms and perceived behavioural control) were significant (p. <. 0.01); HBM was less predictive with only the severity (past behaviour and future intentions) and perceived benefit (only for past behaviour) as significant factors (p. <. 0.05). Conclusions: Non-compliance with DDCLs replacement is widespread, affecting 1 out of 4 Italian wearers. Results from the TPB model show that the involvement of persons socially close to the wearers (subjective norms) and the improvement of the procedure of behavioural control of daily replacement (behavioural control) are of paramount importance in improving compliance. With reference to the HBM, it is important to warn DDCLs wearers of the severity of a contact-lens-related eye infection, and to underline the possibility of its prevention.

Development of the Simulation Model for the CoSES Laboratory Test Microgrid in Modelica

Evolution of the traditional consumer in a power system to a prosumer has posed many problems in the traditional uni-directional grid. This evolution in the grid model has made it important to study the behaviour of microgrids. This thesis deals with the laboratory microgrid setup at the Munich School of Engineering, built to assist researchers in studying microgrids. The model is built in Dymola which is a tool for the OpenModelica language. Models for the different components were derived, suiting the purpose of this study. The equivalent parameters were derived from data sheets and other simulation programs such as PSCAD. The parameters were entered into the model grid and tested at steady state, firstly. This yielded satisfactory results that were similar to the reference results from MATPOWER power flow. Furthermore, fault conditions at several buses were simulated to observe the behaviour of the grid under these conditions. Recommendations for further developing this model to include more detailed models for components, such as power electronic converters, were made at the end of the thesis.

Discharge characteristics and discharge simulation model of long air gaps containing a floating conductor

Long air gaps containing a floating conductor are common insulation types in power grids. During the transmission line live-line work, the process of lineman entering the transmission line air gap constitutes a live-line work combined air gap, which is a typical long air gap containing a floating conductor. This thesis investigates the discharge characteristics, the discharge mechanism and a discharge simulation model of long air gaps containing a floating conductor in order to address the engineering issues in live-line work. The innovative achievements of the thesis are as follows: (1) The effect of the gap distance, the floating electrode structure, the switching impulse wavefront time, the altitude, and the deviation of the floating conductor from the axis on the breakdown voltage was determined. (2) The physical process of the discharges in long air gaps containing a floating conductor was determined. The reason why the discharge characteristics of long air gaps containing a floating electrode with complex geometrics and sharp protrusions and long air gaps with a rod-shaped floating electrode are similar has been studied. The formation mechanism of the lowest breakdown voltage area of a long air gap containing a floating conductor is explained. (3) A simulation discharge model of long air gaps containing a floating conductor was established, which can describe the physical process and predict the breakdown voltage. The model can realize the accurate prediction of the breakdown voltage of typical long air gaps containing a floating conductor and live-line work combined air gaps in transmission lines. The findings of the study can provide theoretical reference and technical support for improving the safety of live-line work.