Estimating the effects of ambient temperature on mortality : methodological challenges and proposed solutions

The health impacts of exposure to ambient temperature have been drawing increasing attention from the environmental health research community, government, society, industries, and the public. Case-crossover and time series models are most commonly used to examine the effects of ambient temperature on mortality. However, some key methodological issues remain to be addressed. For example, few studies have used spatiotemporal models to assess the effects of spatial temperatures on mortality. Few studies have used a case-crossover design to examine the delayed (distributed lag) and non-linear relationship between temperature and mortality. Also, little evidence is available on the effects of temperature changes on mortality, and on differences in heat-related mortality over time. This thesis aimed to address the following research questions: 1. How to combine case-crossover design and distributed lag non-linear models? 2. Is there any significant difference in effect estimates between time series and spatiotemporal models? 3. How to assess the effects of temperature changes between neighbouring days on mortality? 4. Is there any change in temperature effects on mortality over time? To combine the case-crossover design and distributed lag non-linear model, datasets including deaths, and weather conditions (minimum temperature, mean temperature, maximum temperature, and relative humidity), and air pollution were acquired from Tianjin China, for the years 2005 to 2007. I demonstrated how to combine the case-crossover design with a distributed lag non-linear model. This allows the case-crossover design to estimate the non-linear and delayed effects of temperature whilst controlling for seasonality. There was consistent U-shaped relationship between temperature and mortality. Cold effects were delayed by 3 days, and persisted for 10 days. Hot effects were acute and lasted for three days, and were followed by mortality displacement for non-accidental, cardiopulmonary, and cardiovascular deaths. Mean temperature was a better predictor of mortality (based on model fit) than maximum or minimum temperature. It is still unclear whether spatiotemporal models using spatial temperature exposure produce better estimates of mortality risk compared with time series models that use a single site’s temperature or averaged temperature from a network of sites. Daily mortality data were obtained from 163 locations across Brisbane city, Australia from 2000 to 2004. Ordinary kriging was used to interpolate spatial temperatures across the city based on 19 monitoring sites. A spatiotemporal model was used to examine the impact of spatial temperature on mortality. A time series model was used to assess the effects of single site’s temperature, and averaged temperature from 3 monitoring sites on mortality. Squared Pearson scaled residuals were used to check the model fit. The results of this study show that even though spatiotemporal models gave a better model fit than time series models, spatiotemporal and time series models gave similar effect estimates. Time series analyses using temperature recorded from a single monitoring site or average temperature of multiple sites were equally good at estimating the association between temperature and mortality as compared with a spatiotemporal model. A time series Poisson regression model was used to estimate the association between temperature change and mortality in summer in Brisbane, Australia during 1996–2004 and Los Angeles, United States during 1987–2000. Temperature change was calculated by the current day's mean temperature minus the previous day's mean. In Brisbane, a drop of more than 3 �C in temperature between days was associated with relative risks (RRs) of 1.16 (95% confidence interval (CI): 1.02, 1.31) for non-external mortality (NEM), 1.19 (95% CI: 1.00, 1.41) for NEM in females, and 1.44 (95% CI: 1.10, 1.89) for NEM aged 65.74 years. An increase of more than 3 �C was associated with RRs of 1.35 (95% CI: 1.03, 1.77) for cardiovascular mortality and 1.67 (95% CI: 1.15, 2.43) for people aged < 65 years. In Los Angeles, only a drop of more than 3 �C was significantly associated with RRs of 1.13 (95% CI: 1.05, 1.22) for total NEM, 1.25 (95% CI: 1.13, 1.39) for cardiovascular mortality, and 1.25 (95% CI: 1.14, 1.39) for people aged . 75 years. In both cities, there were joint effects of temperature change and mean temperature on NEM. A change in temperature of more than 3 �C, whether positive or negative, has an adverse impact on mortality even after controlling for mean temperature. I examined the variation in the effects of high temperatures on elderly mortality (age . 75 years) by year, city and region for 83 large US cities between 1987 and 2000. High temperature days were defined as two or more consecutive days with temperatures above the 90th percentile for each city during each warm season (May 1 to September 30). The mortality risk for high temperatures was decomposed into: a "main effect" due to high temperatures using a distributed lag non-linear function, and an "added effect" due to consecutive high temperature days. I pooled yearly effects across regions and overall effects at both regional and national levels. The effects of high temperature (both main and added effects) on elderly mortality varied greatly by year, city and region. The years with higher heat-related mortality were often followed by those with relatively lower mortality. Understanding this variability in the effects of high temperatures is important for the development of heat-warning systems. In conclusion, this thesis makes contribution in several aspects. Case-crossover design was combined with distribute lag non-linear model to assess the effects of temperature on mortality in Tianjin. This makes the case-crossover design flexibly estimate the non-linear and delayed effects of temperature. Both extreme cold and high temperatures increased the risk of mortality in Tianjin. Time series model using single site’s temperature or averaged temperature from some sites can be used to examine the effects of temperature on mortality. Temperature change (no matter significant temperature drop or great temperature increase) increases the risk of mortality. The high temperature effect on mortality is highly variable from year to year.

Spatiotemporal model or time series model for assessing city-wide temperature effects on mortality?

Most studies examining the temperature–mortality association in a city used temperatures from one site or the average from a network of sites. This may cause measurement error as temperature varies across a city due to effects such as urban heat islands. We examined whether spatiotemporal models using spatially resolved temperatures produced different associations between temperature and mortality compared with time series models that used non-spatial temperatures. We obtained daily mortality data in 163 areas across Brisbane city, Australia from 2000 to 2004. We used ordinary kriging to interpolate spatial temperature variation across the city based on 19 monitoring sites. We used a spatiotemporal model to examine the impact of spatially resolved temperatures on mortality. Also, we used a time series model to examine non-spatial temperatures using a single site and the average temperature from three sites. We used squared Pearson scaled residuals to compare model fit. We found that kriged temperatures were consistent with observed temperatures. Spatiotemporal models using kriged temperature data yielded slightly better model fit than time series models using a single site or the average of three sites' data. Despite this better fit, spatiotemporal and time series models produced similar associations between temperature and mortality. In conclusion, time series models using non-spatial temperatures were equally good at estimating the city-wide association between temperature and mortality as spatiotemporal models.

A meshfree-based local Galerkin method with condensation of degree of freedom

Condensation technique of degree of freedom is firstly proposed to improve the computational efficiency of meshfree method with Galerkin weak form. In present method, scattered nodes without connectivity are divided into several subsets by cells with arbitrary shape. The local discrete equations are established over each cell by using moving kriging interpolation, in which the nodes that located in the cell are used for approximation. Then, the condensation technique can be introduced into the local discrete equations by transferring equations of inner nodes to equations of boundary nodes based on cell. In the scheme of present method, the calculation of each cell is carried out by meshfree method with Galerkin weak form, and local search is implemented in interpolation. Numerical examples show that the present method has high computational efficiency and convergence, and good accuracy is also obtained.

Spatiotemporal patterns of Aedes aegyptipopulations in Cairns, Australia : assessing drivers of dengue transmission

OBJECTIVES To identify the meteorological drivers of dengue vector density and determine high- and low-risk transmission zones for dengue prevention and control in Cairns, Australia. METHODS Weekly adult female Ae. aegypti data were obtained from 79 double sticky ovitraps (SOs) located in Cairns for the period September 2007-May 2012. Maximum temperature, total rainfall and average relative humidity data were obtained from the Australian Bureau of Meteorology for the study period. Time series-distributed lag nonlinear models were used to assess the relationship between meteorological variables and vector density. Spatial autocorrelation was assessed via semivariography, and ordinary kriging was undertaken to predict vector density in Cairns. RESULTS Ae. aegypti density was associated with temperature and rainfall. However, these relationships differed between short (0-6 weeks) and long (0-30 weeks) lag periods. Semivariograms showed that vector distributions were spatially autocorrelated in September 2007-May 2008 and January 2009-May 2009, and vector density maps identified high transmission zones in the most populated parts of Cairns city, as well as Machans Beach. CONCLUSION Spatiotemporal patterns of Ae. aegypti in Cairns are complex, showing spatial autocorrelation and associations with temperature and rainfall. Sticky ovitraps should be placed no more than 1.2 km apart to ensure entomological coverage and efficient use of resources. Vector density maps provide evidence for the targeting of prevention and control activities. Further research is needed to explore the possibility of developing an early warning system of dengue based on meteorological and environmental factors.

Spatially-explicit modelling of grassland classes - an improved method of integrating a climate-based classification model with interpolated climate surfaces

Spatially-explicit modelling of grassland classes is important to site-specific planning for improving grassland and environmental management over large areas. In this study, a climate-based grassland classification model, the Comprehensive and Sequential Classification System (CSCS) was integrated with spatially interpolated climate data to classify grassland in Gansu province, China. The study area is characterized by complex topographic features imposed by plateaus, high mountains, basins and deserts. To improve the quality of the interpolated climate data and the quality of the spatial classification over this complex topography, three linear regression methods, namely an analytic method based on multiple regression and residues (AMMRR), a modification of the AMMRR method through adding the effect of slope and aspect to the interpolation analysis (M-AMMRR) and a method which replaces the IDW approach for residue interpolation in M-AMMRR with an ordinary kriging approach (I-AMMRR), for interpolating climate variables were evaluated. The interpolation outcomes from the best interpolation method were then used in the CSCS model to classify the grassland in the study area. Climate variables interpolated included the annual cumulative temperature and annual total precipitation. The results indicated that the AMMRR and M-AMMRR methods generated acceptable climate surfaces but the best model fit and cross validation result were achieved by the I-AMMRR method. Twenty-six grassland classes were classified for the study area. The four grassland vegetation classes that covered more than half of the total study area were "cool temperate-arid temperate zonal semi-desert", "cool temperate-humid forest steppe and deciduous broad-leaved forest", "temperate-extra-arid temperate zonal desert", and "frigid per-humid rain tundra and alpine meadow". The vegetation classification map generated in this study provides spatial information on the locations and extents of the different grassland classes. This information can be used to facilitate government agencies' decision-making in land-use planning and environmental management, and for vegetation and biodiversity conservation. The information can also be used to assist land managers in the estimation of safe carrying capacities which will help to prevent overgrazing and land degradation.

Developing sustainable fish farming in the Western Pacific : a viewpoint on potential reasons for why many attempts failed

Capture fisheries and aquaculture have been a major source of food and providers of economic benefits to many communities around the world for a very long time. While the history of aquaculture or fish farming can be traced back for more than 2000 years in some corners of the globe, notably in China, Japan and the Mediterranean, this is not true everywhere, where in general, fish farming is a relatively new industry. Rapid human population growth and increasing urbanisation over the last 20 to 40 years has meant that while fish consumption has doubled globally, returns from capture fisheries have remained static or have declined due to overexploitation and rising pollution levels, with some fisheries either closing or becoming economically unviable. Data from studies suggest that this trend is unlikely to be reversed unless appropriate fisheries management allows depleted wild stocks to rebuild. This has occurred during a time when demand for fish products has grown, in part due to improved purchasing power in some developing countries and changing dietary habits where fish are now considered to have a positive impact on health. Based on the projected population growth over the next two decades, Food and Agricultural Organization (FAO) estimates that at least an additional 40 million tonnes of aquatic food will be required to maintain the current per capita consumption (FAO 2006).

A sub‒domain smoothed Galerkin method for solid mechanics problems

A sub‒domain smoothed Galerkin method is proposed to integrate the advantages of mesh‒free Galerkin method and FEM. Arbitrarily shaped sub‒domains are predefined in problems domain with mesh‒free nodes. In each sub‒domain, based on mesh‒free Galerkin weak formulation, the local discrete equation can be obtained by using the moving Kriging interpolation, which is similar to the discretization of the high‒order finite elements. Strain smoothing technique is subsequently applied to the nodal integration of sub‒domain by dividing the sub‒domain into several smoothing cells. Moreover, condensation of DOF can also be introduced into the local discrete equations to improve the computational efficiency. The global governing equations of present method are obtained on the basis of the scheme of FEM by assembling all local discrete equations of the sub‒domains. The mesh‒free properties of Galerkin method are retained in each sub‒domain. Several 2D elastic problems have been solved on the basis of this newly proposed method to validate its computational performance. These numerical examples proved that the newly proposed sub‒domain smoothed Galerkin method is a robust technique to solve solid mechanics problems based on its characteristics of high computational efficiency, good accuracy, and convergence.

A meshfree-based local Galerkin method with condensation of degree of freedom for elastic dynamic analysis

Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dynamic analysis. In the present method, scattered nodes without connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approximation. Then local discrete equations can be simplified by condensation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by assembling all local discrete equations and are solved by using the standard implicit Newmark’s time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is implemented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.

Assessment of gene expression of intracellular calcium channels, pumps and exchangers with epidermal growth factor-induced epithelial-mesenchymal transition in a breast cancer cell line

Background Epithelial-mesenchymal transition (EMT) is a process implicated in cancer metastasis that involves the conversion of epithelial cells to a more mesenchymal and invasive cell phenotype. In breast cancer cells EMT is associated with altered store-operated calcium influx and changes in calcium signalling mediated by activation of cell surface purinergic receptors. In this study, we investigated whether MDA-MB-468 breast cancer cells induced to undergo EMT exhibit changes in mRNA levels of calcium channels, pumps and exchangers located on intracellular calcium storing organelles, including the Golgi, mitochondria and endoplasmic reticulum (ER). Methods Epidermal growth factor (EGF) was used to induce EMT in MDA-MB-468 breast cancer cells. Serum-deprived cells were treated with EGF (50 ng/mL) for 12 h and gene expression was assessed using quantitative RT-PCR. Results and conclusions These data reveal no significant alterations in mRNA levels of the Golgi calcium pump secretory pathway calcium ATPases (SPCA1 and SPCA2), or the mitochondrial calcium uniporter (MCU) or Na+/Ca2+ exchanger (NCLX). However, EGF-induced EMT was associated with significant alterations in mRNA levels of specific ER calcium channels and pumps, including (sarco)-endoplasmic reticulum calcium ATPases (SERCAs), and inositol 1,4,5-trisphosphate receptor (IP3R) and ryanodine receptor (RYR) calcium channel isoforms. The most prominent change in gene expression between the epithelial and mesenchymal-like states was RYR2, which was enriched 45-fold in EGF-treated MDA-MB-468 cells. These findings indicate that EGF-induced EMT in breast cancer cells may be associated with major alterations in ER calcium homeostasis.

Induction of epithelial-mesenchymal transition (EMT) in breast cancer cells is calcium signal dependent

Signals from the tumor microenvironment trigger cancer cells to adopt an invasive phenotype through epithelial-mesenchymal transition (EMT). Relatively little is known regarding key signal transduction pathways that serve as cytosolic bridges between cell surface receptors and nuclear transcription factors to induce EMT. A better understanding of these early EMT events may identify potential targets for the control of metastasis. One rapid intracellular signaling pathway that has not yet been explored during EMT induction is calcium. Here we show that stimuli used to induce EMT produce a transient increase in cytosolic calcium levels in human breast cancer cells. Attenuation of the calcium signal by intracellular calcium chelation significantly reduced epidermal growth factor (EGF)- and hypoxia-induced EMT. Intracellular calcium chelation also inhibited EGF-induced activation of signal transducer and activator of transcription 3 (STAT3), while preserving other signal transduction pathways such as Akt and extracellular signal-regulated kinase 1/2 (ERK1/2) phosphorylation. To identify calcium-permeable channels that may regulate EMT induction in breast cancer cells, we performed a targeted siRNA-based screen. We found that transient receptor potential-melastatin-like 7 (TRPM7) channel expression regulated EGF-induced STAT3 phosphorylation and expression of the EMT marker vimentin. Although intracellular calcium chelation almost completely blocked the induction of many EMT markers, including vimentin, Twist and N-cadherin, the effect of TRPM7 silencing was specific for vimentin protein expression and STAT3 phosphorylation. These results indicate that TRPM7 is a partial regulator of EMT in breast cancer cells, and that other calcium-permeable ion channels are also involved in calcium-dependent EMT induction. In summary, this work establishes an important role for the intracellular calcium signal in the induction of EMT in human breast cancer cells. Manipulation of calcium-signaling pathways controlling EMT induction in cancer cells may therefore be an important therapeutic strategy for preventing metastases.

Characterising the spatial dynamics of sympatric Aedes aegypti and Aedes albopictus populations in the Philippines

Entomological surveillance and control are essential to the management of dengue fever (DF). Hence, understanding the spatial and temporal patterns of DF vectors, Aedes (Stegomyia) aegypti (L.) and Ae. (Stegomyia) albopictus (Skuse), is paramount. In the Philippines, resources are limited and entomological surveillance and control are generally commenced during epidemics, when transmission is difficult to control. Recent improvements in spatial epidemiological tools and methods offer opportunities to explore more efficient DF surveillance and control solutions: however, there are few examples in the literature from resource-poor settings. The objectives of this study were to: (i) explore spatial patterns of Aedes populations and (ii) predict areas of high and low vector density to inform DF control in San Jose village, Muntinlupa city, Philippines. Fortnightly, adult female Aedes mosquitoes were collected from 50 double-sticky ovitraps (SOs) located in San Jose village for the period June-November 2011. Spatial clustering analysis was performed to identify high and low density clusters of Ae. aegypti and Ae. albopictus mosquitoes. Spatial autocorrelation was assessed by examination of semivariograms, and ordinary kriging was undertaken to create a smoothed surface of predicted vector density in the study area. Our results show that both Ae. aegypti and Ae. albopictus were present in San Jose village during the study period. However, one Aedes species was dominant in a given geographic area at a time, suggesting differing habitat preferences and interspecies competition between vectors. Density maps provide information to direct entomological control activities and advocate the development of geographically enhanced surveillance and control systems to improve DF management in the Philippines.

Remodeling of purinergic receptor-mediated Ca 2+ signaling as a consequence of EGF-induced epithelial-mesenchymal transition in breast cancer cells

Background The microenvironment plays a pivotal role in tumor cell proliferation, survival and migration. Invasive cancer cells face a new set of environmental challenges as they breach the basement membrane and colonize distant organs during the process of metastasis. Phenotypic switching, such as that which occurs during epithelial-mesenchymal transition (EMT), may be associated with a remodeling of cell surface receptors and thus altered responses to signals from the tumor microenvironment. Methodology/Principal Findings We assessed changes in intracellular Ca 2+ in cells loaded with Fluo-4 AM using a fluorometric imaging plate reader (FLIPR TETRA) and observed significant changes in the potency of ATP (EC 50 0.175 μM (-EGF) versus 1.731 μM (+EGF), P<0.05), and the nature of the ATP-induced Ca 2+ transient, corresponding with a 10-fold increase in the mesenchymal marker vimentin (P<0.05). We observed no change in the sensitivity to PAR2-mediated Ca 2+ signaling, indicating that these alterations are not simply a consequence of changes in global Ca 2+ homeostasis. To determine whether changes in ATP-mediated Ca 2+ signaling are preceded by alterations in the transcriptional profile of purinergic receptors, we analyzed the expression of a panel of P2X ionotropic and P2Y metabotropic purinergic receptors using real-time RT-PCR and found significant and specific alterations in the suite of ATP-activated purinergic receptors during EGF-induced EMT in breast cancer cells. Our studies are the first to show that P2X 5 ionotropic receptors are enriched in the mesenchymal phenotype and that silencing of P2X 5 leads to a significant reduction (25%, P<0.05) in EGF-induced vimentin protein expression. Conclusions The acquisition of a new suite of cell surface purinergic receptors is a feature of EGF-mediated EMT in MDA-MB-468 breast cancer cells. Such changes may impart advantageous phenotypic traits and represent a novel mechanism for the targeting of cancer metastasis.

Spatial prediction on a river network

This article develops methods for spatially predicting daily change of dissolved oxygen (Dochange) at both sampled locations (134 freshwater sites in 2002 and 2003) and other locations of interest throughout a river network in South East Queensland, Australia. In order to deal with the relative sparseness of the monitoring locations in comparison to the number of locations where one might want to make predictions, we make a classification of the river and stream locations. We then implement optimal spatial prediction (ordinary and constrained kriging) from geostatistics. Because of their directed-tree structure, rivers and streams offer special challenges. A complete approach to spatial prediction on a river network is given, with special attention paid to environmental exceedances. The methodology is used to produce a map of Dochange predictions for 2003. Dochange is one of the variables measured as part of the Ecosystem Health Monitoring Program conducted within the Moreton Bay Waterways and Catchments Partnership.

Analysis and modelling of the hydraulic conductivity in aquitards : application to the Galilee Basin and the Great Artesian Basin, Australia

In large sedimentary basins with layers of different rocks, the groundwater flow between aquifers depends on the hydraulic conductivity (K) of the separating low-permeable rocks, or aquitards. Three methods were developed to evaluate K in aquitards for areas with limited field data: • Coherence and harmonic analysis: estimates the regional-scale K based on water-level fluctuations in adjacent aquifers. • Cokriging and Bayes' rule: infers K from downhole geophysical logs. • Fluvial process model: reproduces the lithology architecture of sediment formations which can be converted to K. These proposed methods enable good estimates of K and better planning of further drillholes.

Recovering Functional Integrity in Pastures Dominated by an Invasive Grass Species

Exotic grasses have been introduced in countries worldwide for pasture improvement, soil stabilisation and ornamental purposes. Some of these introductions have proven successful, but many have not (Cook & Dias 2006). In Australia, the Commonwealth Plant Introduction Scheme was initiated in 1929, and over-time introduced more than 5000 species of grasses, legumes and other forage and browse plants (Cook & Dias 2006). Lonsdale (1994) suggested that, in tropical Australia, 13% of introductions have become a problem, with only 5% being considered useful for agriculture. Low (1997) suggested that 5 out of 18 of Australia's worst tropical environmental weeds were intentionally introduced as pasture grasses. The spread and dominance of invasive grass species that degrade the quality of pastures for production can impact significantly on the livelihoods of small proprietors. Although Livestock grazing contributes only a small percentage to the world's GDP (1.5%), maintaining the long-term stability of this industry is crucial because of the high social and environmental consequence of a collapse. One billion of the world's poor are dependent on livestock grazing for food and income with this industry occupying more than 25% of the world's land base (Steinfeld et al. 2006). The ling-term sustainability of livestock grazing is also crucial for the environment. A recent FAO report attributed livestock production as a major cause of five of the most serious environmental problems: global warming, land degredation, air and water pollution, and the loss of biodiversity (Steinfeld et al. 2006). For these reasons, finding more effective approaches that guide the sustainable management of pastures is urgently needed. In Australia more than 55% of land use is for livestock grazing by sheelp and/or cattle. This land use dominate in the semi-arid and arid regions where rainfall and soil conditions are marginal for production (Commonwealth of Australia 2004). Although the level of agriculture production by conglomerates is increasing, the majority of livestock grazing within Australia remains family owned and operated (Commonwealth of Australia 2004). The sustainability of production from a grazed pasture is dependent on its botanical composition (Kemp & Dowling 1991, Kemp et al. 1996). In a grazed pasture, the dominance of an invasive grass species can impact on the functional integrity of the ecosystem, including production and nutrient cycling; wwhich will in turn, affect the income of proprietors and the ability of the system to recover from disturbance and environmental change. In Australia, $0.3 billion is spent on weed control in livestock production, but despite this substantial investment $1.9 billion is still lost in yield as a result of weeds (Sinden et al. 2004). In this paper, we adaprt a framework proposed for the restoration of degraded rainforest communities (Lamb & Gilmour 2003, Lamb et al. 2005) to compare and contrast options for recovering function integrity (i.e. a diverse set of desirable plant species that maintain key ecological processes necessary for sustainable production and nutrient cycling) within pasture communities dominated by an invasive grass species. To do this, we uase a case-study of the invasion of Eragrostis curvula (Africal lovegrss; hereafter, Lovegrass), a serious concern in Australian agricultural communities (Parsons and Cuthbertson 1992). The spread and dominance of Lovegrass is a problem because its low palatability, low nutritional content and competitiveness affect the livelihood of graziers by reducing the diversity of other plant species. We conclude by suggesting modifications to this framework for pasture ecosystems to help increase the effiency of strategies to protect functional integrity and balance social/economic and biodiversity values.