Comparison for the effects of different components of temperature variability on mortality: A multi-country time-series study.

BACKGROUND: Temperature variability (TV) is associated with increased mortality risk. However, it is still unknown whether intra-day or inter-day TV has different effects. OBJECTIVES: We aimed to assess the association of intra-day TV and inter-day TV with all-cause, cardiovascular, and respiratory mortality. METHODS: We collected data on total, cardiovascular, and respiratory mortality and meteorology from 758 locations in 47 countries or regions from 1972 to 2020. We defined inter-day TV as the standard deviation (SD) of daily mean temperatures across the lag interval, and intra-day TV as the average SD of minimum and maximum temperatures on each day. In the first stage, inter-day and intra-day TVs were modelled simultaneously in the quasi-Poisson time-series model for each location. In the second stage, a multi-level analysis was used to pool the location-specific estimates. RESULTS: Overall, the mortality risk due to each interquartile range [IQR] increase was higher for intra-day TV than for inter-day TV. The risk increased by 0.59% (95% confidence interval [CI]: 0.53, 0.65) for all-cause mortality, 0.64% (95% CI: 0.56, 0.73) for cardiovascular mortality, and 0.65% (95% CI: 0.49, 0.80) for respiratory mortality per IQR increase in intra-day TV0-7 (0.9 °C). An IQR increase in inter-day TV0-7 (1.6 °C) was associated with 0.22% (95% CI: 0.18, 0.26) increase in all-cause mortality, 0.44% (95% CI: 0.37, 0.50) increase in cardiovascular mortality, and 0.31% (95% CI: 0.21, 0.41) increase in respiratory mortality. The proportion of all-cause deaths attributable to intra-day TV0-7 and inter-day TV0-7 was 1.45% and 0.35%, respectively. The mortality risks varied by lag interval, climate area, season, and climate type. CONCLUSIONS: Our results indicated that intra-day TV may explain the main part of the mortality risk related to TV and suggested that comprehensive evaluations should be proposed in more countries to help protect human health.

Constructing synthetic populations in the age of big data.

BACKGROUND: To develop public health intervention models using micro-simulations, extensive personal information about inhabitants is needed, such as socio-demographic, economic and health figures. Confidentiality is an essential characteristic of such data, while the data should reflect realistic scenarios. Collection of such data is possible only in secured environments and not directly available for open-source micro-simulation models. The aim of this paper is to illustrate a method of construction of synthetic data by predicting individual features through models based on confidential data on health and socio-economic determinants of the entire Dutch population. METHODS: Administrative records and health registry data were linked to socio-economic characteristics and self-reported lifestyle factors. For the entire Dutch population (n = 16,778,708), all socio-demographic information except lifestyle factors was available. Lifestyle factors were available from the 2012 Dutch Health Monitor (n = 370,835). Regression model was used to sequentially predict individual features. RESULTS: The synthetic population resembles the original confidential population. Features predicted in the first stages of the sequential procedure are virtually similar to those in the original population, while those predicted in later stages of the sequential procedure carry the accumulation of limitations furthered by data quality and previously modelled features. CONCLUSIONS: By combining socio-demographic, economic, health and lifestyle related data at individual level on a large scale, our method provides us with a powerful tool to construct a synthetic population of good quality and with no confidentiality issues.

Global, regional, and national burden of mortality associated with short-term temperature variability from 2000-19: a three-stage modelling study.

BACKGROUND: Increased mortality risk is associated with short-term temperature variability. However, to our knowledge, there has been no comprehensive assessment of the temperature variability-related mortality burden worldwide. In this study, using data from the MCC Collaborative Research Network, we first explored the association between temperature variability and mortality across 43 countries or regions. Then, to provide a more comprehensive picture of the global burden of mortality associated with temperature variability, global gridded temperature data with a resolution of 0·5°â€ˆ× 0·5° were used to assess the temperature variability-related mortality burden at the global, regional, and national levels. Furthermore, temporal trends in temperature variability-related mortality burden were also explored from 2000-19. METHODS: In this modelling study, we applied a three-stage meta-analytical approach to assess the global temperature variability-related mortality burden at a spatial resolution of 0·5°â€ˆ× 0·5° from 2000-19. Temperature variability was calculated as the SD of the average of the same and previous days' minimum and maximum temperatures. We first obtained location-specific temperature variability related-mortality associations based on a daily time series of 750 locations from the Multi-country Multi-city Collaborative Research Network. We subsequently constructed a multivariable meta-regression model with five predictors to estimate grid-specific temperature variability related-mortality associations across the globe. Finally, percentage excess in mortality and excess mortality rate were calculated to quantify the temperature variability-related mortality burden and to further explore its temporal trend over two decades. FINDINGS: An increasing trend in temperature variability was identified at the global level from 2000 to 2019. Globally, 1 753 392 deaths (95% CI 1 159 901-2 357 718) were associated with temperature variability per year, accounting for 3·4% (2·2-4·6) of all deaths. Most of Asia, Australia, and New Zealand were observed to have a higher percentage excess in mortality than the global mean. Globally, the percentage excess in mortality increased by about 4·6% (3·7-5·3) per decade. The largest increase occurred in Australia and New Zealand (7·3%, 95% CI 4·3-10·4), followed by Europe (4·4%, 2·2-5·6) and Africa (3·3, 1·9-4·6). INTERPRETATION: Globally, a substantial mortality burden was associated with temperature variability, showing geographical heterogeneity and a slightly increasing temporal trend. Our findings could assist in raising public awareness and improving the understanding of the health impacts of temperature variability. FUNDING: Australian Research Council, Australian National Health & Medical Research Council.

Fluctuating temperature modifies heat-mortality association around the globe.

Studies have investigated the effects of heat and temperature variability (TV) on mortality. However, few assessed whether TV modifies the heat-mortality association. Data on daily temperature and mortality in the warm season were collected from 717 locations across 36 countries. TV was calculated as the standard deviation of the average of the same and previous days' minimum and maximum temperatures. We used location-specific quasi-Poisson regression models with an interaction term between the cross-basis term for mean temperature and quartiles of TV to obtain heat-mortality associations under each quartile of TV, and then pooled estimates at the country, regional, and global levels. Results show the increased risk in heat-related mortality with increments in TV, accounting for 0.70% (95% confidence interval [CI]: -0.33 to 1.69), 1.34% (95% CI: -0.14 to 2.73), 1.99% (95% CI: 0.29-3.57), and 2.73% (95% CI: 0.76-4.50) of total deaths for Q1-Q4 (first quartile-fourth quartile) of TV. The modification effects of TV varied geographically. Central Europe had the highest attributable fractions (AFs), corresponding to 7.68% (95% CI: 5.25-9.89) of total deaths for Q4 of TV, while the lowest AFs were observed in North America, with the values for Q4 of 1.74% (95% CI: -0.09 to 3.39). TV had a significant modification effect on the heat-mortality association, causing a higher heat-related mortality burden with increments of TV. Implementing targeted strategies against heat exposure and fluctuant temperatures simultaneously would benefit public health.

Geographical Variations of the Minimum Mortality Temperature at a Global Scale: A Multicountry Study.

BACKGROUND: Minimum mortality temperature (MMT) is an important indicator to assess the temperature-mortality association, indicating long-term adaptation to local climate. Limited evidence about the geographical variability of the MMT is available at a global scale. METHODS: We collected data from 658 communities in 43 countries under different climates. We estimated temperature-mortality associations to derive the MMT for each community using Poisson regression with distributed lag nonlinear models. We investigated the variation in MMT by climatic zone using a mixed-effects meta-analysis and explored the association with climatic and socioeconomic indicators. RESULTS: The geographical distribution of MMTs varied considerably by country between 14.2 and 31.1 °C decreasing by latitude. For climatic zones, the MMTs increased from alpine (13.0 °C) to continental (19.3 °C), temperate (21.7 °C), arid (24.5 °C), and tropical (26.5 °C). The MMT percentiles (MMTPs) corresponding to the MMTs decreased from temperate (79.5th) to continental (75.4th), arid (68.0th), tropical (58.5th), and alpine (41.4th). The MMTs indreased by 0.8 °C for a 1 °C rise in a community's annual mean temperature, and by 1 °C for a 1 °C rise in its SD. While the MMTP decreased by 0.3 centile points for a 1 °C rise in a community's annual mean temperature and by 1.3 for a 1 °C rise in its SD. CONCLUSIONS: The geographical distribution of the MMTs and MMTPs is driven mainly by the mean annual temperature, which seems to be a valuable indicator of overall adaptation across populations. Our results suggest that populations have adapted to the average temperature, although there is still more room for adaptation.

Mortality risk attributable to wildfire-related PM2·5 pollution: a global time series study in 749 locations.

BACKGROUND: Many regions of the world are now facing more frequent and unprecedentedly large wildfires. However, the association between wildfire-related PM2·5 and mortality has not been well characterised. We aimed to comprehensively assess the association between short-term exposure to wildfire-related PM2·5 and mortality across various regions of the world. METHODS: For this time series study, data on daily counts of deaths for all causes, cardiovascular causes, and respiratory causes were collected from 749 cities in 43 countries and regions during 2000-16. Daily concentrations of wildfire-related PM2·5 were estimated using the three-dimensional chemical transport model GEOS-Chem at a 0·25°â€ˆ× 0·25° resolution. The association between wildfire-related PM2·5 exposure and mortality was examined using a quasi-Poisson time series model in each city considering both the current-day and lag effects, and the effect estimates were then pooled using a random-effects meta-analysis. Based on these pooled effect estimates, the population attributable fraction and relative risk (RR) of annual mortality due to acute wildfire-related PM2·5 exposure was calculated. FINDINGS: 65·6 million all-cause deaths, 15·1 million cardiovascular deaths, and 6·8 million respiratory deaths were included in our analyses. The pooled RRs of mortality associated with each 10 µg/m3 increase in the 3-day moving average (lag 0-2 days) of wildfire-related PM2·5 exposure were 1·019 (95% CI 1·016-1·022) for all-cause mortality, 1·017 (1·012-1·021) for cardiovascular mortality, and 1·019 (1·013-1·025) for respiratory mortality. Overall, 0·62% (95% CI 0·48-0·75) of all-cause deaths, 0·55% (0·43-0·67) of cardiovascular deaths, and 0·64% (0·50-0·78) of respiratory deaths were annually attributable to the acute impacts of wildfire-related PM2·5 exposure during the study period. INTERPRETATION: Short-term exposure to wildfire-related PM2·5 was associated with increased risk of mortality. Urgent action is needed to reduce health risks from the increasing wildfires. FUNDING: Australian Research Council, Australian National Health & Medical Research Council.

Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study.

BACKGROUND: Exposure to cold or hot temperatures is associated with premature deaths. We aimed to evaluate the global, regional, and national mortality burden associated with non-optimal ambient temperatures. METHODS: In this modelling study, we collected time-series data on mortality and ambient temperatures from 750 locations in 43 countries and five meta-predictors at a grid size of 0·5°â€ˆ× 0·5° across the globe. A three-stage analysis strategy was used. First, the temperature-mortality association was fitted for each location by use of a time-series regression. Second, a multivariate meta-regression model was built between location-specific estimates and meta-predictors. Finally, the grid-specific temperature-mortality association between 2000 and 2019 was predicted by use of the fitted meta-regression and the grid-specific meta-predictors. Excess deaths due to non-optimal temperatures, the ratio between annual excess deaths and all deaths of a year (the excess death ratio), and the death rate per 100 000 residents were then calculated for each grid across the world. Grids were divided according to regional groupings of the UN Statistics Division. FINDINGS: Globally, 5 083 173 deaths (95% empirical CI [eCI] 4 087 967-5 965 520) were associated with non-optimal temperatures per year, accounting for 9·43% (95% eCI 7·58-11·07) of all deaths (8·52% [6·19-10·47] were cold-related and 0·91% [0·56-1·36] were heat-related). There were 74 temperature-related excess deaths per 100 000 residents (95% eCI 60-87). The mortality burden varied geographically. Of all excess deaths, 2 617 322 (51·49%) occurred in Asia. Eastern Europe had the highest heat-related excess death rate and Sub-Saharan Africa had the highest cold-related excess death rate. From 2000-03 to 2016-19, the global cold-related excess death ratio changed by -0·51 percentage points (95% eCI -0·61 to -0·42) and the global heat-related excess death ratio increased by 0·21 percentage points (0·13-0·31), leading to a net reduction in the overall ratio. The largest decline in overall excess death ratio occurred in South-eastern Asia, whereas excess death ratio fluctuated in Southern Asia and Europe. INTERPRETATION: Non-optimal temperatures are associated with a substantial mortality burden, which varies spatiotemporally. Our findings will benefit international, national, and local communities in developing preparedness and prevention strategies to reduce weather-related impacts immediately and under climate change scenarios. FUNDING: Australian Research Council and the Australian National Health and Medical Research Council.

Particulate air pollution from different sources and mortality in 7.5 million adults - The Dutch Environmental Longitudinal Study (DUELS).

BACKGROUND: Long-term exposure to particulate air pollution has been associated with mortality in urban cohort studies. Few studies have investigated the association between emission contributions from different particle sources and mortality in large-scale population registries, including non-urban populations. OBJECTIVES: The aim of the study was to evaluate the associations between long-term exposure to particulate air pollution from different source categories and non-accidental mortality in the Netherlands based on existing national databases. METHODS: We used existing Dutch national databases on mortality, individual characteristics, residence history, neighbourhood characteristics and modelled air pollution concentrations from different sources and air pollution components: particulate matter PM10, primary particulate matter PM10 (PPM10), particulate matter PM2.5, primary particulate matter PM2.5 (PPM2.5), elemental carbon (EC), nitrogen dioxide (NO2) and secondary inorganic aerosol (SIA) in PM10 (SIA10) or in PM2.5 (SIA2.5). We established a cohort of 7.5 million individuals 30 years or older. We followed the cohort for eight years (2008-2015). We applied Cox proportional hazard regression models adjusting for potential individual and area-specific confounders. RESULTS: We found statistically significant associations between total and primary particulate matter (PM10 and PM2.5), elemental carbon and mortality. Adjustment for nitrogen dioxide did not change the associations. Secondary inorganic aerosol showed less consistent associations. All primary PM sources were associated with mortality, except agricultural emissions and, depending on the statistical model, industrial PM emissions. CONCLUSIONS: We could not identify one or more specific source categories of particulate air pollution as main determinants of the mortality effects found in this and in a previous study. This suggests that present policy measures should be focussed on the wider spectrum of air pollution sources instead of on specific sources.

Estimating the prevalence of 26 health-related indicators at neighbourhood level in the Netherlands using structured additive regression.

BACKGROUND: Local policy makers increasingly need information on health-related indicators at smaller geographic levels like districts or neighbourhoods. Although more large data sources have become available, direct estimates of the prevalence of a health-related indicator cannot be produced for neighbourhoods for which only small samples or no samples are available. Small area estimation provides a solution, but unit-level models for binary-valued outcomes that can handle both non-linear effects of the predictors and spatially correlated random effects in a unified framework are rarely encountered. METHODS: We used data on 26 binary-valued health-related indicators collected on 387,195 persons in the Netherlands. We associated the health-related indicators at the individual level with a set of 12 predictors obtained from national registry data. We formulated a structured additive regression model for small area estimation. The model captured potential non-linear relations between the predictors and the outcome through additive terms in a functional form using penalized splines and included a term that accounted for spatially correlated heterogeneity between neighbourhoods. The registry data were used to predict individual outcomes which in turn are aggregated into higher geographical levels, i.e. neighbourhoods. We validated our method by comparing the estimated prevalences with observed prevalences at the individual level and by comparing the estimated prevalences with direct estimates obtained by weighting methods at municipality level. RESULTS: We estimated the prevalence of the 26 health-related indicators for 415 municipalities, 2599 districts and 11,432 neighbourhoods in the Netherlands. We illustrate our method on overweight data and show that there are distinct geographic patterns in the overweight prevalence. Calibration plots show that the estimated prevalences agree very well with observed prevalences at the individual level. The estimated prevalences agree reasonably well with the direct estimates at the municipal level. CONCLUSIONS: Structured additive regression is a useful tool to provide small area estimates in a unified framework. We are able to produce valid nationwide small area estimates of 26 health-related indicators at neighbourhood level in the Netherlands. The results can be used for local policy makers to make appropriate health policy decisions.

Assessment of spatial and physical neighborhood characteristics that influence sound quality and herewith well-being and health.

Environmental noise and health studies seldom address the positive effect of environments with high acoustic quality. Sound quality, in turn, is influenced by a large number of factors, including the spatial-physical characteristics of a neighborhood. In general, these characteristics cannot be retrieved from existing databases. In this article, we describe the design of an audit instrument and demonstrate its value for gathering information about these characteristics of neighborhoods. The audit instrument used was derived from research in other fields than environmental health. The instrument was tested in 33 neighborhoods in the Dutch cities of Amsterdam, Rotterdam, and Arnhem. In these neighborhoods, more or less homogeneous subareas were identified that were subject of the audit. The results show that the audit approach is suitable to gather neighborhood data that are relevant for the sound quality of neighborhoods. Together with survey data, they provide information that could further the field of soundscape and health. Several suggestions for improvement of the audit instrument were made.

Comparing the effect of a leaflet and a movie in preventing tick bites and Lyme disease in The Netherlands.

BACKGROUND: Lyme disease (LD) has become the most common vector borne illness in the Northern hemisphere. Prevention relies predominantly on fostering protective behaviors (e.g., avoiding tick areas, using protective clothing and repellent, and doing routine tick checks post-exposure). The objective of this study was to evaluate the effectiveness (in terms of knowledge, perceived severity and susceptibility, self-efficacy, response efficacy, intention, and behavior over time) and appreciation of a leaflet and a movie as tools for informing the public in the Netherlands about ticks and LD protective behaviors. METHODS: Participants (1,677 at t1 and 361 extra at t2) were members of a representative Internet panel (adults aged 18 years and above). A four group randomized controlled design was used to test the effect of an information leaflet and a movie (two intervention groups), compared to a control group and a follow-up only control group. Data were collected over two periods: July 15-29, 2013 (t1) and at follow-up 4 weeks later, August 16-31, 2013 (t2). RESULTS: Post-intervention results show all respondents in all groups possess good general basic knowledge of ticks and LD. Respondents in both the leaflet and movie groups knew more than respondents in the control group, and had greater awareness of best practices after a tick bite. Intention to perform protective behavior in future was stronger among respondents in the intervention groups. While respondents generally appreciated both the movie and the leaflet, they found the movie ran too long. Follow-up revealed no lasting positive effects from either the leaflet or the movie. CONCLUSIONS: Our results suggest that both the movie and the leaflet are valued and effective intervention tools for improving knowledge about tick bites and strengthening self-efficacy and intentions to perform protective behavior against ticks and LD . Achieving lasting effects, however, calls for more action.

Air Pollution and Mortality in Seven Million Adults: The Dutch Environmental Longitudinal Study (DUELS).

BACKGROUND: Long-term exposure to air pollution has been associated with mortality in urban cohort studies. Few studies have investigated this association in large-scale population registries, including non-urban populations. OBJECTIVES: The aim of the study was to evaluate the associations between long-term exposure to air pollution and nonaccidental and cause-specific mortality in the Netherlands based on existing national databases. METHODS: We used existing Dutch national databases on mortality, individual characteristics, residence history, neighborhood characteristics, and national air pollution maps based on land use regression (LUR) techniques for particulates with an aerodynamic diameter ≤ 10 µm (PM10) and nitrogen dioxide (NO2). Using these databases, we established a cohort of 7.1 million individuals ≥ 30 years of age. We followed the cohort for 7 years (2004-2011). We applied Cox proportional hazard models adjusting for potential individual and area-specific confounders. RESULTS: After adjustment for individual and area-specific confounders, for each 10-µg/m3 increase, PM10 and NO2 were associated with nonaccidental mortality [hazard ratio (HR) = 1.08; 95% CI: 1.07, 1.09 and HR = 1.03; 95% CI: 1.02, 1.03, respectively], respiratory mortality (HR = 1.13; 95% CI: 1.10, 1.17 and HR = 1.02; 95% CI: 1.01, 1.03, respectively), and lung cancer mortality (HR = 1.26; 95% CI: 1.21, 1.30 and HR = 1.10 95% CI: 1.09, 1.11, respectively). Furthermore, PM10 was associated with circulatory disease mortality (HR = 1.06; 95% CI: 1.04, 1.08), but NO2 was not (HR = 1.00; 95% CI: 0.99, 1.01). PM10 associations were robust to adjustment for NO2; NO2 associations remained for nonaccidental mortality and lung cancer mortality after adjustment for PM10. CONCLUSIONS: Long-term exposure to PM10 and NO2 was associated with nonaccidental and cause-specific mortality in the Dutch population of ≥ 30 years of age.

Lung cancer risk and past exposure to emissions from a large steel plant.

We studied the spatial distribution of cancer incidence rates around a large steel plant and its association with historical exposure. The study population was close to 600,000. The incidence data was collected for 1995-2006. From historical emission data the air pollution concentrations for polycyclic aromatic hydrocarbons (PAH) and metals were modelled. Data were analyzed using Bayesian hierarchical Poisson regression models. The standardized incidence ratio (SIR) for lung cancer was up to 40% higher than average in postcodes located in two municipalities adjacent to the industrial area. Increased incidence rates could partly be explained by differences in socioeconomic status (SES). In the highest exposure category (approximately 45,000 inhabitants) a statistically significant increased relative risk (RR) of 1.21 (1.01-1.43) was found after adjustment for SES. The elevated RRs were similar for men and women. Additional analyses in a subsample of the population with personal smoking data from a recent survey suggested that the observed association between lung cancer and plant emission, after adjustment for SES, could still be caused by residual confounding. Therefore, we cannot indisputably conclude that past emissions from the steel plant have contributed to the increased risk of lung cancer.

Trends in relative risk estimates for the association between air pollution and mortality in The Netherlands, 1992-2006.

BACKGROUND: Daily variations in the levels of air pollution are well known to be associated with daily variations in mortality counts. Given the large number of time-series studies, there is little need for simple replication of these results in additional locations. However, additional analyses of time-series data might be useful in elucidating remaining questions on the role of air pollution on mortality. OBJECTIVES: Because of ongoing issues related to causality, changing toxicity, the difficulty in isolating the independent effects of individual pollutants, the availability of new methods to detect effect thresholds, and questions about the extent to which effects are restricted to frail members of the population, additional analyses of time-series data might be helpful in addressing these issues. We show an example where additional time-series analyses can be helpful in elucidating specific questions in the field of air pollution epidemiology. METHODS: We analysed daily mortality and air pollution data using Poisson regression in generalised additive models. Air pollution data for the overall period 1992-2006 and for four different periods were analysed to assess the overall risk estimates for the whole period and to assess variability over time for the different effect estimates. RESULTS AND CONCLUSION: We found some statistically significant upward trends, but this was only the case for a few associations without a consistent pattern over the cause-specific deaths. Whether these findings are consistent over time or whether our findings are merely the result of statistical chance can only be elucidated by continuation of monitoring of the relative risks over time in the future. Although these results may indicate that both photochemical and particulate matter air pollution might have become more toxic, the lack of a clear pattern in the results makes these conclusions speculative.

The association between noise exposure and blood pressure and ischemic heart disease: a meta-analysis.

It has been suggested that noise exposure is associated with blood pressure changes and ischemic heart disease risk, but epidemiologic evidence is still limited. Furthermore, most reviews investigating these relations were not carried out in a systematic way, which makes them more prone to bias. We conducted a meta-analysis of 43 epidemiologic studies published between 1970 and 1999 that investigate the relation between noise exposure (both occupational and community) and blood pressure and/or ischemic heart disease (International Classification of Diseases, Ninth Revision, codes 410-414). We studied a wide range of effects, from blood pressure changes to a myocardial infarction. With respect to the association between noise exposure and blood pressure, small blood pressure differences were evident. Our meta-analysis showed a significant association for both occupational noise exposure and air traffic noise exposure and hypertension: We estimated relative risks per 5 dB(A) noise increase of 1.14 (1.01-1.29) and 1.26 (1.14-1.39), respectively. Air traffic noise exposure was positively associated with the consultation of a general practitioner or specialist, the use of cardiovascular medicines, and angina pectoris. In cross-sectional studies, road traffic noise exposure increases the risk of myocardial infarction and total ischemic heart disease. Although we can conclude that noise exposure can contribute to the prevalence of cardiovascular disease, the evidence for a relation between noise exposure and ischemic heart disease is still inconclusive because of the limitations in exposure characterization, adjustment for important confounders, and the occurrence of publication bias.