Evaluating the urban heat island for cities of Emilia-Romagna region through numerical simulations

Lo scopo di questo studio è la comprensione della dinamica dello strato limite urbano per città dell’Emilia Romagna tramite simulazioni numeriche. In particolare, l’attenzione è posta sull’ effetto isola di calore, ovvero sulla differenza di temperatura dell’aria in prossimità del suolo fra zone rurali e urbane dovuta all’urbanizzazione. Le simulazioni sono state effettuate con il modello alla mesoscala "Weather Research and Forecasting" (WRF), accoppiato con le parametrizzazioni urbane "Building Effect Parametrization" (BEP) e "Building Energy Model" (BEM), che agiscono a vari livelli verticali urbani. Il periodo di studio riguarda sei giorni caldi e senza copertura nuvolosa durante un periodo di heat wave dell’anno 2015. La copertura urbana è stata definita con il "World Urban Databes and Access Portal Tools" (WUDAPT), un metodo che permette di classificare le aree urbane in dieci "urban climate zones" (UCZ), attraverso l’uso combinato di immagini satellitari e "training areas" manualmente definite con il software Google Earth. Sono state svolte diverse simulazioni a domini innestati, con risoluzione per il dominio più piccolo di 500 m, centrato sulla città di Bologna. Le differenze fra le simulazioni riguardano la presenza o l’assenza delle strutture urbane, il metodo di innesto e tipo di vegetazione rurale. Inoltre, è stato valutato l’effetto dovuto alla presenza di pannelli fotovoltaici sopra i tetti di ogni edificio e le variazioni che i pannelli esercitano sullo strato limite urbano. Per verificare la bontà del modello, i dati provenienti dalle simulazioni sono stati confrontati con misure provenienti da 41 stazioni all’interno dell’area di studio. Le variabili confrontate sono: temperatura, umidità relativa, velocità e direzione del vento. Le simulazioni sono in accordo con i dati osservativi e riescono a riprodurre l’effetto isola di calore: la differenza di temperatura fra città e zone rurali circostanti è nulla durante il giorno; al contrario, durante la notte l’isola di calore è presente, e in media raggiunge il massimo valore di 4°C alle 1:00. La presenza dei pannelli fotovoltaici abbassa la temperatura a 2 metri dell’aria al massimo di 0.8°C durante la notte, e l’altezza dello strato limite urbano dell’ordine 200mrispetto al caso senza pannelli. I risultati mostrano come l’uso di pannelli fotovoltaici all’interno del contesto urbano ha molteplici benefici: infatti, i pannelli fotovoltaici riescono a ridurre la temperatura durante un periodo di heat wave, e allo stesso tempo possono parzialmente sopperire all’alto consumo energetico, con una conseguente riduzione del consumo di combustibili fossili.

LPJ-LMfire simulations for Europe for the Last Glacial Maximum and preindustrial control, link to netCDF files

This data sets contains LPJ-LMfire dynamic global vegetation model output covering Europe and the Mediterranean for the Last Glacial Maximum (LGM; 21 ka) and for a preindustrial control simulation (20th century detrended climate). The netCDF data files are time averages of the final 30 years of the model simulation. Each netCDF file contains four or five variables: fractional cover of 9 plant functional types (PFTs; cover), total fractional coverage of trees (treecover), population density of hunter-gatherers (foragerPD; only for the "people" simulations), fraction of the gridcell burned on 30-year average (burnedf), and vegetation net primary productivity (NPP). The model spatial resolution is 0.5-degrees For the LGM simulations, LPJ-LMfire was driven by the PMIP3 suite of eight GCMs for which LGM climate simulations were available. Also provided in this archive is the result of an LPJ-LMfire run that was forced by the average climate of all GCMs (the "GCM-mean" files), and the average of each of the individual LPJ-LMfire runs over the eight LGM scenarios individually (the "LPJ-mean" files). The model simulations are provided that include the influence of human presence on the landscape (the "people" files), and in a "world without humans" scenario (the "natural" files). Finally this archive contains the preindustrial reference simulation with and without human influence ("PI_reference_people" and "PI_reference_nat", respectively). There are therefore 22 netCDF files in this archive: 8 each of LGM simulations with and without people (total 16) and the "GCM mean" simulation (2 files) and the "LPJ mean" aggregate (2 files), and finally the two preindustrial "control" simulations ("PI"), with and without humans (2 files). In addition to the LPJ-LMfire model output (netCDF files), this archive also contains a table of arboreal pollen percent calculated from pollen samples dated to the LGM at sites throughout (lgmAP.txt), and a table containing the location of archaeological sites dated to the LGM (LGM_archaeological_site_locations.txt).

Simulations of charged droplet collisions in shear flow

Acknowledgments This research has been enabled by the use of computing resources provided by WestGrid, the Shared Hierarchical Academic Research Computing Network (SHARCNET: www.sharcnet.ca), and Compute/Calcul Canada. O.S. thanks NSERC for an Alexander Graham Bell Canada Graduate Scholarship.

Méthodes de simulations moléculaires accélérées : application et développement

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Risk from cattle trampling to nests of an endangered passerine evaluated using artificial nest experiments and simulations

Grasslands are often grazed by cattle and many grassland birds nest on the ground, potentially exposing nests to trampling. We tested for trampling risk introduced by cattle to nests of endangered Florida Grasshopper Sparrows (Ammodramus savannarum floridanus) using experimentally paired grids of artificial nests (i.e., clay targets) similar in size to nests of Florida Grasshopper Sparrows and counted the number of clay targets that were broken in paired grazed and ungrazed enclosures. Clay targets in grazed grids were trampled 3.9% more often than their respective ungrazed grids, and measurements of cattle presence or density were correlated with the number of broken clay targets, suggesting that excluding cattle during breeding is an important management recommendation for the Florida Grasshopper Sparrow. Trampling rates within grazed enclosures were spatially homogeneous with respect to cattle infrastructure such as supplemental feeding troughs and fences, and forests and stocking density were poor predictors of trampling rates when excluding ungrazed grids. We used population viability analysis to compare quasi-extinction rates, intrinsic growth rates, and median abundance in grazed and ungrazed Florida Grasshopper Sparrow aggregations to further understand the biological significance of management aimed at reducing trampling rates during the breeding season. Simulations indicated that trampling from grazing increased quasi-extinction rates by 41% while reducing intrinsic growth rates by 0.048, and reducing median abundance by an average of 214 singing males after 50 years. Management should avoid grazing enclosures occupied by Florida Grasshopper Sparrows during the nesting season to minimize trampling rates. Our methods that combine trampling experiments with population viability analysis provide a framework for testing effects from trampling on other grassland ground-nesting birds, and can directly inform conservation and management of the Florida Grasshopper Sparrow.

Numerical Simulations for the Design of Novel Integrated Parallel Reception, Excitation, and Shimming (iPRES) Coil Arrays

Magnetic field inhomogeneity results in image artifacts including signal loss, image blurring and distortions, leading to decreased diagnostic accuracy. Conventional multi-coil (MC) shimming method employs both RF coils and shimming coils, whose mutual interference induces a tradeoff between RF signal-to-noise (SNR) ratio and shimming performance. To address this issue, RF coils were integrated with direct-current (DC) shim coils to shim field inhomogeneity while concurrently emitting and receiving RF signal without being blocked by the shim coils. The currents applied to the new coils, termed iPRES (integrated parallel reception, excitation and shimming), were optimized in the numerical simulation to improve the shimming performance. The objectives of this work is to offer a guideline for designing the optimal iPRES coil arrays to shim the abdomen.

In this thesis work, the main field () inhomogeneity was evaluated by root mean square error (RMSE). To investigate the shimming abilities of iPRES coil arrays, a set of the human abdomen MRI data was collected for the numerical simulations. Thereafter, different simplified iPRES(N) coil arrays were numerically modeled, including a 1-channel iPRES coil and 8-channel iPRES coil arrays. For 8-channel iPRES coil arrays, each RF coil was split into smaller DC loops in the x, y and z direction to provide extra shimming freedom. Additionally, the number of DC loops in a RF coil was increased from 1 to 5 to find the optimal divisions in z direction. Furthermore, switches were numerically implemented into iPRES coils to reduce the number of power supplies while still providing similar shimming performance with equivalent iPRES coil arrays.

The optimizations demonstrate that the shimming ability of an iPRES coil array increases with number of DC loops per RF coil. Furthermore, the z direction divisions tend to be more effective in reducing field inhomogeneity than the x and y divisions. Moreover, the shimming performance of an iPRES coil array gradually reach to a saturation level when the number of DC loops per RF coil is large enough. Finally, when switches were numerically implemented in the iPRES(4) coil array, the number of power supplies can be reduced from 32 to 8 while keeping the shimming performance similar to iPRES(3) and better than iPRES(1). This thesis work offers a guidance for the designs of iPRES coil arrays.

Méthodes de simulations moléculaires accélérées : application et développement

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Mechanisms Governing the Eyewall Replacement Cycle in Numerical Simulations of Tropical Cyclones

Eyewall replacement cycle (ERC) is frequently observed during the evolution of intensifying Tropical Cyclones (TCs). Although intensely studied in recent years, the underlying mechanisms of ERC are still poorly understood, and the forecast of ERC remains a great challenge. To advance our understanding of ERC and provide insights in improvement of numerical forecast of ERC, a series of numerical simulations is performed to investigate ERCs in TC-like vortices on a f-plane. The simulated ERCs possess key features similar to those observed in real TCs including the formation of a secondary tangential wind maximum associated with the outer eyewall. The Sawyer-Eliassen equation and tangential momentum budget analyses are performed to diagnose the mechanisms underlying the secondary eyewall formation (SEF) and ERC. Our diagnoses reveal crucial roles of outer rainband heating in governing the formation and development of the secondary tangential wind maximum and demonstrate that the outer rainband convection must reach a critical strength relative to the eyewall before SEF and the subsequent ERC can occur. A positive feedback among low-level convection, acceleration of tangential winds in the boundary layer, and surface evaporation that leads to the development of ERC and a mechanism for the demise of inner eyewall that involves interaction between the transverse circulations induced by eyewall and outer rainband convection are proposed. The tangential momentum budget indicates that the net tendency of tangential wind is a small residual resultant from a large cancellation between tendencies induced by the resolved and sub-grid scale (SGS) processes. The large SGS contribution to the tangential wind budget explains different characteristics of ERC shown in previous numerical studies and poses a great challenge for a timely correct forecast of ERC. The sensitivity experiments show that ERCs are strongly subjected to model physics, vortex radial structure and background wind. The impact of model physics on ERC can be well understood with the interaction among eyewall/outer rainband heating, radilal inflow in the boundary layer, surface layer turbulent processes, and shallow convection in the moat. However, further investigations are needed to fully understand the exhibited sensitivities of ERC to vortex radial structure and background wind.

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.

Field data, numerical simulations and probability analyses to assess lava flow hazards at Mount Etna

Improving lava flow hazard assessment is one of the most important and challenging fields of volcanology, and has an immediate and practical impact on society. Here, we present a methodology for the quantitative assessment of lava flow hazards based on a combination of field data, numerical simulations and probability analyses. With the extensive data available on historic eruptions of Mt. Etna, going back over 2000 years, it has been possible to construct two hazard maps, one for flank and the other for summit eruptions, allowing a quantitative analysis of the most likely future courses of lava flows. The effective use of hazard maps of Etna may help in minimizing the damage from volcanic eruptions through correct land use in densely urbanized area with a population of almost one million people. Although this study was conducted on Mt. Etna, the approach used is designed to be applicable to other volcanic areas.

Cenozoic global ice volume and temperature simulations over the past 40 million years

Variations in global ice volume and temperature over the Cenozoic era have been investigated with a set of one-dimensional (1-D) ice-sheet models. Simulations include three ice sheets representing glaciation in the Northern Hemisphere, i.e. in Eurasia, North America and Greenland, and two separate ice sheets for Antarctic glaciation. The continental mean Northern Hemisphere surface-air temperature has been derived through an inverse procedure from observed benthic d18O records. These data have yielded a mutually consistent and continuous record of temperature, global ice volume and benthic d18O over the past 35 Ma. The simple 1-D model shows good agreement with a comprehensive 3-D ice-sheet model for the past 3 Ma. On average, differences are only 1.0°C for temperature and 6.2 m for sea level. Most notably, over the 35 Ma period, the reconstructed ice volume-temperature sensitivity shows a transition from a climate controlled by Southern Hemisphere ice sheets to one controlled by Northern Hemisphere ice sheets. Although the transient behaviour is important, equilibrium experiments show that the relationship between temperature and sea level is linear and symmetric, providing limited evidence for hysteresis. Furthermore, the results show a good comparison with other simulations of Antarctic ice volume and observed sea level.

Next generation e-learning platform: virtual laboratory simulations

Bioscience subjects require a significant amount of training in laboratory techniques to produce highly skilled science graduates. Many techniques which are currently used in diagnostic, research and industrial laboratories require expensive equipment for single users; examples of which include next generation sequencing, quantitative PCR, mass spectrometry and other analytical techniques. The cost of the machines, reagents and limited access frequently preclude undergraduate students from using such cutting edge techniques. In addition to cost and availability, the time taken for analytical runs on equipment such as High Performance Liquid Chromatography (HPLC) does not necessarily fit with the limitations of timetabling. Understanding the theory underlying these techniques without the accompanying practical classes can be unexciting for students. One alternative from wet laboratory provision is to use virtual simulations of such practical which enable students to see the machines and interact with them to generate data. The Faculty of Science and Technology at the University of Westminster has provided all second and third year undergraduate students with iPads so that these students all have access to a mobile device to assist with learning. We have purchased licences from Labster to access a range of virtual laboratory simulations. These virtual laboratories are fully equipped and require student responses to multiple answer questions in order to progress through the experiment. In a pilot study to look at the feasibility of the Labster virtual laboratory simulations with the iPad devices; second year Biological Science students (n=36) worked through the Labster HPLC simulation on iPads. The virtual HPLC simulation enabled students to optimise the conditions for the separation of drugs. Answers to Multiple choice questions were necessary to progress through the simulation, these focussed on the underlying principles of the HPLC technique. Following the virtual laboratory simulation students went to a real HPLC in the analytical suite in order to separate of asprin, caffeine and paracetamol. In a survey 100% of students (n=36) in this cohort agreed that the Labster virtual simulation had helped them to understand HPLC. In free text responses one student commented that "The terminology is very clear and I enjoyed using Labster very much”. One member of staff commented that “there was a very good knowledge interaction with the virtual practical”.

Westminster Serious Games Platform (wmin-SGP) a tool for real-time authoring of roleplay simulations for learning

The use of serious games in education and their pedagogical benefit is being widely recognized. However, effective integration of serious games in education depends on addressing two big challenges: the successful incorporation of motivation and engagement that can lead to learning; and the highly specialised skills associated with customised development to meet the required pedagogical objectives. This paper presents the Westminster Serious Games Platform (wmin-SGP) an authoring tool that allows educators/domain experts without games design and development technical skills to create bespoke roleplay simulations in three dimensional scenes featuring fully embodied virtual humans capable of verbal and non-verbal interaction with users fit for specific educational objectives. The paper presents the wmin-SGP system architecture and it evaluates its effectiveness in fulfilling its purpose via the implementation of two roleplay simulations, one for Politics and one for Law. In addition, it presents the results of two types of evaluation that address how successfully the wmin-SGP combines usability principles and game core drives based on the Octalysis gamification framework that lead to motivating games experiences. The evaluation results shows that the wmin-SGP: provides an intuitive environment and tools that support users without advanced technical skills to create in real-time bespoke roleplay simulations in advanced graphical interfaces; satisfies most of the usability principles; and provides balanced simulations based on the Octalysis framework core drives. The paper concludes with a discussion of future extension of this real time authoring tool and directions for further development of the Octalysis framework to address learning.

An implementation architecture for scenario-based simulations

The past years have witnessed an increased use of applied games for developing and evaluating communication skills. These skills benefit from in-terpersonal interactions. Providing feedback to students practicing communica-tion skills is difficult in a traditional class setting with one teacher and many students. This logistic challenge may be partly overcome by providing training using a simulation in which a student practices with communication scenarios. A scenario is a description of a series of interactions, where at each step the player is faced with a choice. We have developed a scenario editor that enables teachers to develop scenarios for practicing communication skills. A teacher can develop a scenario without knowledge of the implementation. This paper presents the implementation architecture for such a scenario-based simulation.