Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

From quantum to classical dynamics: The relativistic O ( N ) model in the framework of the real-time functional renormalization group

We investigate the transition from unitary to dissipative dynamics in the relativistic O(N) vector model with the λ(φ2)2 interaction using the nonperturbative functional renormalization group in the real-time formalism. In thermal equilibrium, the theory is characterized by two scales, the interaction range for coherent scattering of particles and the mean free path determined by the rate of incoherent collisions with excitations in the thermal medium. Their competition determines the renormalization group flow and the effective dynamics of the model. Here we quantify the dynamic properties of the model in terms of the scale-dependent dynamic critical exponent z in the limit of large temperatures and in 2≤d≤4 spatial dimensions. We contrast our results to the behavior expected at vanishing temperature and address the question of the appropriate dynamic universality class for the given microscopic theory.

Spatial Distribution of Stem Cell-Like Keratinocytes in Dissected Compound Hair Follicles of the Dog

Hair cycle disturbances are common in dogs and comparable to some alopecic disorders in humans. A normal hair cycle is maintained by follicular stem cells which are predominately found in an area known as the bulge. Due to similar morphological characteristics of the bulge area in humans and dogs, the shared particularity of compound hair follicles as well as similarities in follicular biomarker expression, the dog is a promising model to study human hair cycle and stem cell disorders. To gain insight into the spatial distribution of follicular keratinocytes with stem cell potential in canine compound follicles, we microdissected hair follicles in anagen and telogen from skin samples of freshly euthanized dogs. The keratinocytes isolated from different locations were investigated for their colony forming efficiency, growth and differentiation potential as well as clonal growth. Our results indicate that i) compound and single hair follicles exhibit a comparable spatial distribution pattern with respect to cells with high growth potential and stem cell-like characteristics, ii) the lower isthmus (comprising the bulge) harbors most cells with high growth potential in both, the anagen and the telogen hair cycle stage, iii) unlike in other species, colonies with highest growth potential are rather small with an irregular perimeter and iv) the keratinocytes derived from the bulbar region exhibit characteristics of actively dividing transit amplifying cells. Our results now provide the basis to conduct comparative studies of normal dogs and those with hair cycle disorders with the possibility to extend relevant findings to human patients.

The Effects of Density, Spatial Pattern, and Competitive Symmetry on Size Variation in Simulated Plant Populations

Patterns of size inequality in crowded plant populations are often taken to be indicative of the degree of size asymmetry of competition, but recent research suggests that some of the patterns attributed to size‐asymmetric competition could be due to spatial structure. To investigate the theoretical relationships between plant density, spatial pattern, and competitive size asymmetry in determining size variation in crowded plant populations, we developed a spatially explicit, individual‐based plant competition model based on overlapping zones of influence. The zone of influence of each plant is modeled as a circle, growing in two dimensions, and is allometrically related to plant biomass. The area of the circle represents resources potentially available to the plant, and plants compete for resources in areas in which they overlap. The size asymmetry of competition is reflected in the rules for dividing up the overlapping areas. Theoretical plant populations were grown in random and in perfectly uniform spatial patterns at four densities under size‐asymmetric and size‐symmetric competition. Both spatial pattern and size asymmetry contributed to size variation, but their relative importance varied greatly over density and over time. Early in stand development, spatial pattern was more important than the symmetry of competition in determining the degree of size variation within the population, but after plants grew and competition intensified, the size asymmetry of competition became a much more important source of size variation. Size variability was slightly higher at higher densities when competition was symmetric and plants were distributed nonuniformly in space. In a uniform spatial pattern, size variation increased with density only when competition was size asymmetric. Our results suggest that when competition is size asymmetric and intense, it will be more important in generating size variation than is local variation in density. Our results and the available data are consistent with the hypothesis that high levels of size inequality commonly observed within crowded plant populations are largely due to size‐asymmetric competition, not to variation in local density.

A Mathematical Model for Simplifying Representations of Objects in a Geographic Information System

The study of operations on representations of objects is well documented in the realm of spatial engineering. However, the mathematical structure and formal proof of these operational phenomena are not thoroughly explored. Other works have often focused on query-based models that seek to order classes and instances of objects in the form of semantic hierarchies or graphs. In some models, nodes of graphs represent objects and are connected by edges that represent different types of coarsening operators. This work, however, studies how the coarsening operator "simplification" can manipulate partitions of finite sets, independent from objects and their attributes. Partitions that are "simplified first have a collection of elements filtered (removed), and then the remaining partition is amalgamated (some sub-collections are unified). Simplification has many interesting mathematical properties. A finite composition of simplifications can also be accomplished with some single simplification. Also, if one partition is a simplification of the other, the simplified partition is defined to be less than the other partition according to the simp relation. This relation is shown to be a partial-order relation based on simplification. Collections of partitions can not only be proven to have a partial- order structure, but also have a lattice structure and are complete. In regard to a geographic information system (GIs), partitions related to subsets of attribute domains for objects are called views. Objects belong to different views based whether or not their attribute values lie in the underlying view domain. Given a particular view, objects with their attribute n-tuple codings contained in the view are part of the actualization set on views, and objects are labeled according to the particular subset of the view in which their coding lies. Though the scope of the work does not mainly focus on queries related directly to geographic objects, it provides verification for the existence of particular views in a system with this underlying structure. Given a finite attribute domain, one can say with mathematical certainty that different views of objects are partially ordered by simplification, and every collection of views has a greatest lower bound and least upper bound, which provides the validity for exploring queries in this regard.

Public Commons for Geospatial Data: A Conceptual Model

A wide variety of spatial data collection efforts are ongoing throughout local, state and federal agencies, private firms and non-profit organizations. Each effort is established for a different purpose but organizations and individuals often collect and maintain the same or similar information. The United States federal government has undertaken many initiatives such as the National Spatial Data Infrastructure, the National Map and Geospatial One-Stop to reduce duplicative spatial data collection and promote the coordinated use, sharing, and dissemination of spatial data nationwide. A key premise in most of these initiatives is that no national government will be able to gather and maintain more than a small percentage of the geographic data that users want and desire. Thus, national initiatives depend typically on the cooperation of those already gathering spatial data and those using GIs to meet specific needs to help construct and maintain these spatial data infrastructures and geo-libraries for their nations (Onsrud 2001). Some of the impediments to widespread spatial data sharing are well known from directly asking GIs data producers why they are not currently involved in creating datasets that are of common or compatible formats, documenting their datasets in a standardized metadata format or making their datasets more readily available to others through Data Clearinghouses or geo-libraries. The research described in this thesis addresses the impediments to wide-scale spatial data sharing faced by GIs data producers and explores a new conceptual data-sharing approach, the Public Commons for Geospatial Data, that supports user-friendly metadata creation, open access licenses, archival services and documentation of parent lineage of the contributors and value- adders of digital spatial data sets.

Temporal and spatial regulation of cell division by the nucleoid in Escherichia coli

Selection of division sites and coordination of cytokinesis with other cell cycle events are critical for every organism to proliferate. In E. coli, the nucleoid is proposed to exclude division from the site of the chromosome (nucleoid occlusion model). We studied the effect of the nucleoid on timing and placement of cell division. An early cell division protein, FtsZ, was used to follow development of the division septum. FtsZ forms a ring structure (Z ring) at potential division sites. The dynamics of Z ring was visualized in live cells by fusing FtsZ with a green fluorescent protein (GFP). Emanating FtsZ-GFP polymers from the constricted septum or aggregates in daughter cells were also observed, probably representing the FtsZ depolymerization and immature FtsZ nucleation processes. We next examined the nucleoid occlusion model. Mutants carrying abnormally positioned chromosomes were employed. In chromosomal partition mutants, replicated chromosomes cannot segregate. The Z ring was excluded from midcell to the edge of the nucleoid. This negative effect of nucleoids was further confirmed in replication deficient dnaA mutants, in which only a single chromosome is present in the cell center. These results suggest that the nucleoid, replicating or not, inhibits division in the area where the chromosome occupies. In addition, increasing the level of FtsZ does not overcome nucleoid inhibition. Interestingly in anucleate cells produced by both mutants, the Z ring was localized in the central part of the cell, which indicates that the nucleoid is not required for FtsZ assembly. Relaxation of chromosomes by reducing the gyrase activity or disruption of protein translation/translocation did not abolish the division inhibition capacity of the nucleoid. However, preventing transcription did compromise the nucleoid occlusion effect, leading to formation of multiple FtsZ rings above the nucleoid. In summary, we demonstrate that nucleoids negatively regulate the timing and position of division by inhibiting FtsZ assembly at unselected sites. Relief of this inhibition at midcell is coincident with the completion of DNA replication. On the other hand, FtsZ assembly does not require the nucleoid. ^

Characterization of healing tissue in a tooth extraction socket in a rabbit model

The histology of healing in a tooth extraction socket has been described in many studies. The focus of research in bone biology and healing is now centered on molecular events that regulate repair of injured tissue. Rapid progress in cellular and molecular biology has resulted in identification of many signaling molecules (growth factors and cytokines) associated with formation and repair of skeletal tissues. Some of these include members of the transforming growth factor-β superfamily (including the bone morphogenetic proteins), fibroblast growth factors, platelet derived growth factors and insulin like growth factors. ^ Healing of a tooth extraction socket is a complex process involving tissue repair and regeneration. It involves chemotaxis of appropriate cells into the wound, transformation of undifferentiated mesenchymal cells to osteoprogenitor cells, proliferation and differentiation of committed bone forming cells, extracellular matrix synthesis, mineralization of osteoid, maturation and remodeling of bone. Current data suggests that these cellular events are precisely controlled and regulated by specific signaling molecules. A plethora of cytokines; have been identified and studied in the past two decades. Some of these like transforming growth factor beta (TGF-β), vascular endothelial growth factor (VEGF), platelet derived growth factor (PDGF) and fibroblast growth factors (FGFs) are well conserved proteins involved in the initial response to injury and repair in soft and hard tissue. ^ The purpose of this study was to characterize the spatial and temporal localization of TGF-βl, VEGF, PDGF-A, FGF-2 and BMP-2, and secretory IgA in a tooth extraction socket model, and evaluate correlation of spatial and temporal changes of these growth factors to histological events. The results of this study showed positive correlation of histological events to spatial and temporal localization of TGF-β1, BMP-2, FGF-2, PDGF-A, and VEGF in a rabbit tooth extraction model. ^

Validation of a Spatially Continuous EDEN Water-Surface Model for the Everglades, Florida

The Everglades Depth Estimation Network (EDEN) is an integrated network of realtime water-level monitoring, ground-elevation modeling, and water-surface modeling that provides scientists and managers with current (2000-present), online water-stage and water-depth information for the entire freshwater portion of the Greater Everglades. Continuous daily spatial interpolations of the EDEN network stage data are presented on grid with 400-square-meter spacing. EDEN offers a consistent and documented dataset that can be used by scientists and managers to: (1) guide large-scale field operations, (2) integrate hydrologic and ecological responses, and (3) support biological and ecological assessments that measure ecosystem responses to the implementation of the Comprehensive Everglades Restoration Plan (CERP) (U.S. Army Corps of Engineers, 1999). The target users are biologists and ecologists examining trophic level responses to hydrodynamic changes in the Everglades. The first objective of this report is to validate the spatially continuous EDEN water-surface model for the Everglades, Florida developed by Pearlstine et al. (2007) by using an independent field-measured data-set. The second objective is to demonstrate two applications of the EDEN water-surface model: to estimate site-specific ground elevation by using the validated EDEN water-surface model and observed water depth data; and to create water-depth hydrographs for tree islands. We found that there are no statistically significant differences between model-predicted and field-observed water-stage data in both southern Water Conservation Area (WCA) 3A and WCA 3B. Tree island elevations were derived by subtracting field water-depth measurements from the predicted EDEN water-surface. Water-depth hydrographs were then computed by subtracting tree island elevations from the EDEN water stage. Overall, the model is reliable by a root mean square error (RMSE) of 3.31 cm. By region, the RMSE is 2.49 cm and 7.77 cm in WCA 3A and 3B, respectively. This new landscape-scale hydrological model has wide applications for ongoing research and management efforts that are vital to restoration of the Florida Everglades. The accurate, high-resolution hydrological data, generated over broad spatial and temporal scales by the EDEN model, provides a previously missing key to understanding the habitat requirements and linkages among native and invasive populations, including fish, wildlife, wading birds, and plants. The EDEN model is a powerful tool that could be adapted for other ecosystem-scale restoration and management programs worldwide.

Regional Integration and Industrial Location in a Landlocked Spatial Economy

Regional integration proposals often require agreements between countries that differ in geographic size, resource endowments, transportation assets, technologies, and product quality. In this asymmetric setting, questions arise about the potential for mutual gains and the distribution of benefits among industries and workers in each country. This paper examines how regional integration between a small landlocked country and a large neighboring country--with a unique port facility that both nations must use to export goods--affects the wage and location decisions of firms, the allocation of labor, the welfare of each country's workers and firms, and aggregate measures of economic welfare in each country and the region. A simulated spatial labor market model is used to explore the economic effects of various stages of regional integration. Beginning with autarky as a benchmark case, we consider two forms of regional integration: partial mobility (mobile labor with geographically restricted firms); and full mobility (mobile labor and firms) with convergence of production technologies and product quality.

Spatial analysis of water body distance and West Nile Virus human cases in Houston, TX

West Nile Virus (WNV) is an arboviral disease that has affected hundreds of residents in Harris County, Texas since its introduction in 2002. Persistent infection, lingering sequelae and other long-term symptoms of patients reaffirm the need for prevention of this important vector-borne disease. This study aimed to determine if living within 400m of a water body increases one’s odds of infection with WNV. Additionally, we wanted to determine if one’s proximity to a particular water type or water body source increased one’s odds of infection with WNV.^ 145 cases’ addresses were abstracted from the initial interview and consent records from a cohort of patients (Epidemiology of Arboviral Encephalitis in Houston study, HSC-SPH-03-039). After applying inclusion criteria, 140 cases were identified for analysis. 140 controls were selected for analysis using a population proportionate to size model and US Census Bureau data. MapMarker USA v14 was used to geocode the cases’ addresses. Both cases’ and controls’ coordinates were uploaded onto a Harris County water shapefile in MapInfo Professional v9.5.1. Distance in meters to the closest water source, closest water source type, and closest water source name were recorded.^ Analysis of Variance (p=0.329, R2 = 0.0034) indicated no association between water body distance and risk of WNV disease. Living near a creek (x2 = 11.79, p < 0.001), or the combined group of creek and gully (x 2 = 14.02, p < 0.001) were found to be strongly associated with infection of WNV. Living near Cypress Creek and its feeders (x2 = 15.2, p < 0.001) was found to be strongly associated with WNV infection. We found that creek and gully habitats, particularly Cypress Creek, were preferential for the local disease transmitting Culex quinquefasciatus and reservoir avian population.^

A Bayesian analysis for spatial processes with application to mapping

In geographical epidemiology, maps of disease rates and disease risk provide a spatial perspective for researching disease etiology. For rare diseases or when the population base is small, the rate and risk estimates may be unstable. Empirical Bayesian (EB) methods have been used to spatially smooth the estimates by permitting an area estimate to "borrow strength" from its neighbors. Such EB methods include the use of a Gamma model, of a James-Stein estimator, and of a conditional autoregressive (CAR) process. A fully Bayesian analysis of the CAR process is proposed. One advantage of this fully Bayesian analysis is that it can be implemented simply by using repeated sampling from the posterior densities. Use of a Markov chain Monte Carlo technique such as Gibbs sampler was not necessary. Direct resampling from the posterior densities provides exact small sample inferences instead of the approximate asymptotic analyses of maximum likelihood methods (Clayton & Kaldor, 1987). Further, the proposed CAR model provides for covariates to be included in the model. A simulation demonstrates the effect of sample size on the fully Bayesian analysis of the CAR process. The methods are applied to lip cancer data from Scotland, and the results are compared. ^

Spatial analysis of cardiovascular disease mortality and air pollution in Harris County, Texas

The association between fine particulate matter air pollution (PM2.5) and cardiovascular disease (CVD) mortality was spatially analyzed for Harris County, Texas, at the census tract level. The objective was to assess how increased PM2.5 exposure related to CVD mortality in this area while controlling for race, income, education, and age. An estimated exposure raster was created for Harris County using Kriging to estimate the PM2.5 exposure at the census tract level. The PM2.5 exposure and the CVD mortality rates were analyzed in an Ordinary Least Squares (OLS) regression model and the residuals were subsequently assessed for spatial autocorrelation. Race, median household income, and age were all found to be significant (p<0.05) predictors in the model. This study found that for every one μg/m3 increase in PM2.5 exposure, holding age and education variables constant, an increase of 16.57 CVD deaths per 100,000 would be predicted for increased minimum exposure values and an increase of 14.47 CVD deaths per 100,000 would be predicted for increased maximum exposure values. This finding supports previous studies associating PM2.5 exposure with CVD mortality. This study further identified the areas of greatest PM2.5 exposure in Harris County as being the geographical locations of populations with the highest risk of CVD (i.e., predominantly older, low-income populations with a predominance of African Americans). The magnitude of the effect of PM2.5 exposure on CVD mortality rates in the study region indicates a need for further community-level studies in Harris County, and suggests that reducing excess PM2.5 exposure would reduce CVD mortality.^

Spatial distribution of orofacial cleft birth defects and uranium-radium concentrations in tap water in Atascosa, Bee, Brooks, Calhoun, Duval, Goliad, Hidalgo, Jim Hogg, Jim Wells, Karnes, Kleberg, Live Oak, McMullen, Nueces, San Patricio, Refugio, Starr, Victoria, Webb, and Zavala Counties, Texas

Background Past and recent evidence shows that radionuclides in drinking water may be a public health concern. Developmental thresholds for birth defects with respect to chronic low level domestic radiation exposures, such as through drinking water, have not been definitely recognized, and there is a strong need to address this deficiency in information. In this study we examined the geographic distribution of orofacial cleft birth defects in and around uranium mining district Counties in South Texas (Atascosa, Bee, Brooks, Calhoun, Duval, Goliad, Hidalgo, Jim Hogg, Jim Wells, Karnes, Kleberg, Live Oak, McMullen, Nueces, San Patricio, Refugio, Starr, Victoria, Webb, and Zavala), from 1999 to 2007. The probable association of cleft birth defect rates by ZIP codes classified according to uranium and radium concentrations in drinking water supplies was evaluated. Similar associations between orofacial cleft birth defects and radium/radon in drinking water were reported earlier by Cech and co-investigators in another of the Gulf Coast region (Harris County, Texas).50, 55 Since substantial uranium mining activity existed and still exists in South Texas, contamination of drinking water sources with radiation and its relation to birth defects is a ground for concern. ^ Methods Residential addresses of orofacial cleft birth defect cases, as well as live births within the twenty Counties during 1999-2007 were geocoded and mapped. Prevalence rates were calculated by ZIP codes and were mapped accordingly. Locations of drinking water supplies were also geocoded and mapped. ZIP codes were stratified as having high combined uranium (≥30μg/L) vs. low combined uranium (<30μg/L). Likewise, ZIP codes having the uranium isotope, Ra-226 in drinking water, were also stratified as having elevated radium (≥3 pCi/L) vs. low radium (<3 pCi/L). A linear regression was performed using STATA® generalized linear model (GLM) program to evaluate the probable association between cleft birth defect rates by ZIP codes and concentration of uranium and radium via domestic water supply. These rates were further adjusted for potentially confounding variables such as maternal age, education, occupation, and ethnicity. ^ Results This study showed higher rates of cleft births in ZIP codes classified as having high combined uranium versus ZIP codes having low combined uranium. The model was further improved by adding radium stratified as explained above. Adjustment for maternal age and ethnicity did not substantially affect the statistical significance of uranium or radium concentrations in household water supplies. ^ Conclusion Although this study lacks individual exposure levels, the findings suggest a significant association between elevated uranium and radium concentrations in tap water and high orofacial birth defect rates by ZIP codes. Future case-control studies that can measure individual exposure levels and adjust for contending risk factors could result in a better understanding of the exposure-disease association.^

Using spatial linear models with SAR and CAR structure to examine Texas lung cancer incidence rates

Scholars have found that socioeconomic status was one of the key factors that influenced early-stage lung cancer incidence rates in a variety of regions. This thesis examined the association between median household income and lung cancer incidence rates in Texas counties. A total of 254 individual counties in Texas with corresponding lung cancer incidence rates from 2004 to 2008 and median household incomes in 2006 were collected from the National Cancer Institute Surveillance System. A simple linear model and spatial linear models with two structures, Simultaneous Autoregressive Structure (SAR) and Conditional Autoregressive Structure (CAR), were used to link median household income and lung cancer incidence rates in Texas. The residuals of the spatial linear models were analyzed with Moran's I and Geary's C statistics, and the statistical results were used to detect similar lung cancer incidence rate clusters and disease patterns in Texas.^