Effects of temperature on mortality in Chiang Mai city, Thailand: a time series study
© Guo et al.; licensee BioMed Central Ltd. 2012
Received: 20 June 2011
Accepted: 21 May 2012
Published: 9 July 2012
The association between temperature and mortality has been examined mainly in North America and Europe. However, less evidence is available in developing countries, especially in Thailand. In this study, we examined the relationship between temperature and mortality in Chiang Mai city, Thailand, during 1999–2008.
A time series model was used to examine the effects of temperature on cause-specific mortality (non-external, cardiopulmonary, cardiovascular, and respiratory) and age-specific non-external mortality (<=64, 65–74, 75–84, and > =85 years), while controlling for relative humidity, air pollution, day of the week, season and long-term trend. We used a distributed lag non-linear model to examine the delayed effects of temperature on mortality up to 21 days.
We found non-linear effects of temperature on all mortality types and age groups. Both hot and cold temperatures resulted in immediate increase in all mortality types and age groups. Generally, the hot effects on all mortality types and age groups were short-term, while the cold effects lasted longer. The relative risk of non-external mortality associated with cold temperature (19.35°C, 1st percentile of temperature) relative to 24.7°C (25th percentile of temperature) was 1.29 (95% confidence interval (CI): 1.16, 1.44) for lags 0–21. The relative risk of non-external mortality associated with high temperature (31.7°C, 99th percentile of temperature) relative to 28°C (75th percentile of temperature) was 1.11 (95% CI: 1.00, 1.24) for lags 0–21.
This study indicates that exposure to both hot and cold temperatures were related to increased mortality. Both cold and hot effects occurred immediately but cold effects lasted longer than hot effects. This study provides useful data for policy makers to better prepare local responses to manage the impact of hot and cold temperatures on population health.
KeywordsMortality Cardiovascular Respiratory Temperature Time series analysis
Many studies have demonstrated that both hot and cold temperatures had adverse effects on mortality. For example, elevated ambient temperature and heat-waves were associated with excess deaths in 86 US cities . High temperatures had significant impacts on deaths from all causes, chronic bronchitis, pneumonia, ischemic heart disease, and cerebrovascular disease in England and Wales . McMichael et al. found a U-shaped temperature-mortality relationship in developing countries, with strong evidence of increased deaths on hot and cold days . Chung et al. found an increase in mortality associated with elevated average temperature in Seoul, Beijing, Tokyo and Taipei in Asia .
Current knowledge of the health effects of temperature on mortality is mainly from developed countries [5, 6], particularly from regions with temperate climates. Few studies were conducted in developing countries, particularly in tropical regions . Developing countries are more sensitive to climate change, as they have poor public health infrastructure and more vulnerable populations .
Thailand is one of the countries experiencing a significant increase of temperature over recent years. According to the report of Meteorological Department of Thailand , ambient temperature significantly increased from 1975 to 2005 by approximately 1°C. Model projections show that the temperature will continue to increase in this century. Meanwhile, respiratory and cardiovascular diseases have become a serious public health issue in Thailand. They are among the top five leading causes of death and disability in Thailand . In particular, Chiang Mai city has the highest mortality and morbidity rates for respiratory disease in Thailand. The incidence and mortality of cardiovascular disease have increased slightly every year .
Few studies have examined the impacts of temperatures on cardiovascular, respiratory, and age-specific mortality in Chiang Mai city, Thailand. It is necessary to assess the relationship between temperature and cardiovascular and respiratory diseases in order to better manage this increasing health problem. In this study, we examined the effects of temperature on different types of mortality (including non-external, cardiopulmonary, cardiovascular, and respiratory) and age-specific non-external mortality in Chiang Mai city during 1999–2008.
Materials and method
Chiang Mai city is located in the north of Thailand (latitude 18°47'N and longitude 098°59'E). There were 1,603,220 population with 790,107 (49.28%) men and 813.113 (50.72%) women in 2008 . The major businesses in Chiang Mai city are trade, services, and agriculture. The topography of Chiang Mai city is mountain rimmed basin with poor ventilation. Chiang Mai city has grown rapidly in the past decade. It is the area with severe air pollution. Generally, the major ambient pollutants are particulate matter and ozone .
We collected daily mortality data from Chiang Mai city during 1999 to 2008. The mortality data were extracted from a computerized file of death certificates provided by the Bureau of Policy and Strategy, Ministry of Public Health, Thailand. These cases were residents in Chiang Mai city. The causes of death were coded according to the 10th Revision of the International Classification of Disease (ICD10). The mortality data were stratified into four cause-specific categories: non-external mortality (A00–R99), cardiopulmonary mortality (I00–I99 and J00–J99), cardiovascular mortality (I00–I99), and respiratory mortality (J00–J99). For cardiopulmonary and respiratory mortality, influenza (J10–J11) was excluded, as it was a potential confounder in the data analyses. We also stratified non-external mortality into four age groups (<=64, 65–74, 75–84, and > =85 years).
Daily data on minimum temperature (minimum value between 00:00 and 9:00 am) and maximum temperature (maximum value between 12:00 and 9: 00 pm) (degrees Celsius) and relative humidity were obtained from the Department of Meteorology at Chiang Mai station. Daily data on average particular matter ≤ 10 μm (PM10) (μg/m3) and daily average ozone (O3) (ppb) were collected from Yuparat School monitoring station and Provincial Hall monitoring station. The daily averaged concentrations of air pollution from the two stations were used in our analysis.
In this study, we used a Poisson regression model combined with a distributed lag non-linear model (DLNM) to examine the impact of temperature on mortality [12, 13]. We used a quasi-Poisson function that allows for over-dispersion in daily deaths. We controlled for PM10 (moving average of lags 0–7), O3 (moving average of lags 0–7) and relative humidity (moving average of lags 0–7) using a natural cubic spline with 3 degrees of freedom (df), as these variables are potential confounders of the association between temperature and mortality [14–16]. We controlled for day of the week as a category variable. We controlled for season and long-term trend using a natural cubic spline with 7 df per year for time, as the estimated effects of temperature were then stabilized.
Many studies have shown that the extreme temperature can not only affect current day’s mortality but also influence the several following days’ mortality (lag effect) [17–19]. Also, the relationship between temperature and mortality is non-linear . Therefore, we used a DLNM to examine the non-linear and delayed effects of temperature on cause-specific and age-specific mortality. The DLNM is developed on the basis of “cross-basis” function, which allows simultaneously estimating the non-linear effect of temperature at each lag and the non-linear effects across lags. The DLNM can show the relationship between temperature and mortality at each temperature point and lag. The DLNM can calculate cumulative effect in the existence of delayed contributions. .
We used a “natural cubic spline-natural cubic spline” DLNM to examine the non-linear temperature effect using a 5 df natural cubic spline and the lagged effect using a 4 df natural cubic spline. A maximum lag of 21 days was used to model the effect of temperature on mortality, because previous studies suggest that the effects of cold temperature last about ten days and the effects of high temperature are usually acute and short term [17, 19, 21]. We placed spline knots at equal spaces in the temperature range, and placed spline knots at equal intervals in the log scale of lags using the default setting of DLNM. The median value of temperature was used as the reference value to calculate the relative risks.
To examine the hot effect on cause-specific and age-specific mortality, we evaluated the relative risk of cause-specific and age-specific mortality associated with high temperature (31.7°C, 99th percentile of temperature) relative to the 75th percentile of temperature (28°C). To examine the cold effect on cause-specific and age-specific mortality, we evaluated the relative risk of cause-specific and age-specific mortality associated with cold temperature (19.35°C, 1st percentile of temperature) relative to the 25th percentile of temperature (24.7°C).
We evaluated the model fit using Akaike's Information Criterion for quasi-Poisson (Q-AIC). Our initial results showed that mean temperature (average of maximum and minimum temperatures) was a better predictor than maximum and minimum temperatures. In addition, mean temperature represents the exposure throughout whole day and night and can be straightforwardly interpreted for decision making purpose . Therefore, mean temperature was used in this study.
To check the main findings of this study, sensitivity analyses were performed by changing df (6–15 per year) for time to control for season. We changed df (4–7) for humidity, PM10, and O3, and varied the maximum lags from 22 to 30 days for the DLNM.
All the analyses were performed using R software (version 2.12.1). The “dlnm” package was used to perform DLNM .
Summary statistics for mean temperature, relative humidity, air pollutants, and mortality in Chiang Mai, Thailand, during 1999–2008
Relative humidity (%)
Non-external mortality < =64 years
> = 85 years
The cumulative effects of cold temperature on cause-specific mortality and age-specific non-external mortality, with 1 st percentile of temperature (19.35°C) relative to 25 th percentile of temperature (24.7°C)
Relative risk (95% CI)
> = 85 years
The cumulative effects of high temperature on cause-specific mortality and age-specific non-external mortality, with 99 th percentile of temperature (31.7°C) relative to 75 th percentile of temperature (28°C)
Relative risk (95% CI)
> = 85 years
We changed df (6–15 per year) for time to control for season, which gave similar results (data not shown). We changed df (3–7) for humidity, PM10, and O3, and the estimated effects of temperature were not substantially changed. In addition, we changed the maximum lag from 22 to 30 days, which gave similar results (data not shown). The models used in this study seemed to have adequately captured the main effects of temperature on mortality.
In this study, we examined the effects of temperature on cause-specific and age-specific mortality in Chiang Mai city, Thailand, during 1999–2008. We found that the temperature-mortality relationships were non-linear for all mortality types and age groups. Both extreme cold and hot temperatures had negative impacts on mortality. Both cold and hot effects occurred acutely for all mortality types. Cold effects lasted longer than hot effects.
In general, the relationship between temperature and cause-specific mortality is consistent with previous studies worldwide [3, 24]. This study has also explored the lag effects of temperature on case-specific mortality using DLNM up to 21 days. We found that both hot and cold temperatures were associated with increased non-external deaths. The strongest temperature-mortality relationship was at lag 0. Lag structures in the relationship between temperature and mortality have been reported in several studies [17, 25–27]. Previous studies also reported a similar lag effect of temperature on mortality [17, 19, 21, 28]. Many studies showed that hot temperatures had short term effects on mortality and morbidity [17, 25–27], while the cold effect were delayed and lasted for several days [6, 29]. The acute cold effects in this study might be due to the tropic climatic pattern in Chiang Mai city as people are not accustomed to cold weather.
Even though this study found that both high and low temperatures resulted in immediate increases in mortality, the effects of hot and cold temperatures were different. There were acute and short-term effects of hot temperature on mortality, while the effects of cold temperature persisted for several days. In addition, we found hot effects had mortality displacement for non-external and cardiopulmonary mortality, and people aged > = 85 years. This finding is consistent with several previous studies [17, 19, 30–32]. Hence, it suggests that even in a tropic climate, hot temperature has short term acute effect and cold temperature has relatively long, lagged effects on mortality. Studies in Europe and USA have confirmed that cold temperatures associated with risk in mortality, and that cold effects on mortality in warmer cities were higher than cold cities . People living in rooms without heating, and those wearing fewer clothers had higher cold related mortality .
Many studies have examined the relationship between temperature and mortality in Australia, Europe, and USA [35–39], but few are from Thailand [3, 40]. Studies showed that effects of temperature on mortality varied by population and region [37, 41–44]. The relationship between temperature and non-external mortality in this study was similar as a previous study conducted in Chiang Mai city . In Chiang Mai city, both extreme cold and high temperatures had significant effects on non-external mortality. Studies showed that cities with median or lower income (e.g., Bangkok, Mexico City, São Paulo, Delhi, Santiago, and Cape Town) had significant hot and cold effects on non-external mortality .
The findings in this study were biologically plausible. The human body regulates the heat exchange between the body and ambient temperature by physiological control. When the body is exposed to extreme high temperatures, the thermoregulation system may fail. The reason may be related to dehydration, salt depletion and increased surface blood circulation . High temperature is also related to elevated blood viscosity, cholesterol levels and sweating thresholds, which may also be the cause of heat-related mortality . Extreme cold temperatures are related to the increase in the heart rate, peripheral vasoconstriction, blood pressure, blood cholesterol levels, plasma fibrinogen concentrations, and platelet viscosity [47, 48].
This study has two major strengths. Firstly, this study investigated the effects of temperature on non-external, cardiopulmonary, cardiovascular and respiratory mortality, and age-specific non-external mortality using an advanced statistical approach (DLNM) in Chiang Mai city, Thailand. The advanced statistical approach can flexibly examine effects of temperature on mortality. For example, it can smooth temperature and lag at the same time. Secondly, we used ten years’ data which had high quality (no missing data for mortality). We also adjusted for a range of confounders including relative humidity and air pollutants. However, this study also has some limitations. This is an ecological study in which exposure misclassification might occur as detailed spatial information was unavailable in this study. We only used data from Chiang Mai city to examine the effects of temperature on mortality, so the findings may be difficult to generalize to other areas. We did not control for barometric pressure, wind direction, and wind speed, as these data were not available.
This study examined the effects of temperature on cause-specific and age-specific mortality in Chiang Mai city, Thailand. The relationships between temperature and cause-specific and age-specific mortality were non-linear. Both cold and hot temperatures were associated with increased risk in mortality. Both cold and hot effects occurred immediately. Cold effects lasted longer than hot effects. Cold related mortality risk in oldest people (aged > =85 years) was higher than other age groups. These results demonstrate that temperature is an important environmental hazard in Chiang Mai city.
Akaike's Information Criterion
Degree of freedom
Distributed lag non-linear model
International classification of cisease, tenth revision
Particulate matter ≤10 μm
We thank the Bureau of Policy and Strategy, Ministry of Public Health, and Department of Meteorology, Thailand for providing data. We thank Dr. Lyle Turner for his review.
GY is supported by the QUT Postgraduate Research Award (QUTPRA); TS is supported by a NHMRC Research Fellowship (#553043).
- Glover M: Mortality during periods of excessive temperature. Public Health Rep. 1938, 53: 1122-10.2307/4582590.View Article
- Langford I, Bentham G: The potential effects of climate change on winter mortality in England and Wales. International Journal of Biometeorology. 1995, 38: 141-147. 10.1007/BF01208491.View Article
- McMichael AJ, Wilkinson P, Kovats RS, Pattenden S, Hajat S, Armstrong B, Vajanapoom N, Niciu EM, Mahomed H, Kingkeow C: International study of temperature, heat and urban mortality: the 'ISOTHURM' project. Int J Epidemiol. 2008, 37: 1121-1131. 10.1093/ije/dyn086.View Article
- Chung J, Honda Y, Hong Y, Pan X, Guo Y, Kim H: Ambient temperature and mortality:An international study in four capital cities of East Asia. Science of the Total Environment. 2009, 408: 390-396. 10.1016/j.scitotenv.2009.09.009.View Article
- Guo Y, Barnett AG, Yu W, Pan X, Ye X, Huang C, Tong S: A large change in temperature between neighbouring days increases the risk of mortality. PLoS One. 2011, 6: e16511-10.1371/journal.pone.0016511.View Article
- Anderson G, Bell M: Weather-related mortality: how heat, cold, and heat waves affect mortality in the United States. Epidemiology. 2009, 20: 205-213. 10.1097/EDE.0b013e318190ee08.View Article
- Meteorological Department (Thailand): Thailand Annual Weather Summary. Meteorological Department (Thailand). 2009, [http://www.tmd.go.th/programs%5Cuploads%5CyearlySummary%5CAnnual2009_sum.pdf accessed 8/4/2012]
- World Health Organization:Country Health System Profile: Thailand. 2008, World Health Organization, Geneva, [http://www.searo.who.int/en/Section313/Section1525_10867.htm],
- Bureau of Policy and Strategy (Thailand):Public Health Statistic. Bureau of Policy and Strategy (Thailand). 2008, [http://bps.ops.moph.go.th/Healthinformation/index2.html accessed 8/4/2012],
- Ministry of Interior (Thailand): Chaing Mai Province Profile. In Book Chaing Mai Province Profile. 2008, Ministry of Interior (Thailand), Bangkok
- Pollution Control Department (PCD):State and Management of Air and Noise Pollution. 2003, Bangkok: Pollution Control Department (PCD), [http://www.pcd.go.th/indexeng.cfm],
- Armstrong B: Models for the relationship between ambient temperature and daily mortality. Epidemiology. 2006, 17: 624-631. 10.1097/01.ede.0000239732.50999.8f.View Article
- Gasparrini A, Armstrong B, Kenward MG: Distributed lag non-linear models. Stat Med. 2010, 29: 2224-2234. 10.1002/sim.3940.View Article
- Gouveia N, Fletcher T: Time series analysis of air pollution and mortality:effects by cause, age and socioeconomic status. Journal of Epidemiology Community Health. 2000, 54: 750-755. 10.1136/jech.54.10.750.View Article
- Bernard S, Samet J, Grambsch A, Ebi K, Romieu I: The potential impacts of climate variability and change on air pollution-related health effects in the united states. Environ Health Perspect. 2001, 109: 199-View Article
- Buadong D, Jinsaart W, Funatagawa I, Karita K, Yano E: Association between PM10 and O3 levels and hospital visits for cardiovascular Disease in Bangkok, Thailand. Journal Epidemiology. 2009, 19: 182-188. 10.2188/jea.JE20080047.View Article
- Braga A, Zanobetti A, Schwartz J: The effect of weather on respiratory and cardiovascular deaths in 12 U.S. cities. Environmental Health Perspective. 2002, 110: 859-863.View Article
- Zanobetti A: Generalized additive distributed lag models: quantifying mortality displacement. Biostatistic. 2000, 1: 279-292. 10.1093/biostatistics/1.3.279.View Article
- Guo Y, Barnett AG, Pan X, Yu W, Tong S: The impact of temperature on mortality in tianjin, china: a case-crossover design with a distributed lag nonlinear model. Environ Health Perspect. 2011, 119: 1719-1725. 10.1289/ehp.1103598.View Article
- Gasparrini A, Armstrong B, Kenward M: Distributed lag non-linear models. Statistics in Medicine. 2010, 29: 2224-2234. 10.1002/sim.3940.View Article
- Tong SL, Ren C, Becker N: Excess deaths during the 2004 heatwave in Brisbane, Australia. Int J Biometeorol. 2010, 54: 393-400. 10.1007/s00484-009-0290-8.View Article
- Yu W, Vaneckova P, Mengersen K, Pan X, Tong S: Is the association between temperature and mortality modified by age, gender and socio-economic status?. Sci Total Environ. 2010, 408: 3513-3518. 10.1016/j.scitotenv.2010.04.058.View Article
- Gasparrini A, Armstrong B: Distributed lag non-linear models in R: the package dlnm. [ftp://ftp.cn.debian.org/CRAN/web/packages/dlnm/vignettes/dlnmOverview.pdf]
- Pan W, Li L, Tsai M: Temperature extremes and mortality from coronary heart disease and cerebral infarction in elderly Chinese. the lancet. 1995, 345: 353-355. 10.1016/S0140-6736(95)90341-0.View Article
- Kunst A, Looman C, Mackenbatc J: Outdoor air temperature and mortality in the Netherlands: a time series analysis. American Journal of epidemiology. 1993, 137: 331-341.
- Hajat S, Haines A: Associations of cold temperatures with GP consultations for respiratory and cardiovascular disease amongst the elderly in London. International journal of Epidemiology. 2002, 31: 825-830. 10.1093/ije/31.4.825.View Article
- Basu R, Samet J: Relation between elevated ambient temperature and mortality: A review of the epidemiology evidence. Epidemiologic review. 2002, 82: 61-76.
- Hashizume M, Wagatsuma Y, Hayashi T, Saha S, Streatfield K, Yunus M: The effect of temperature on mortality in rural Bangladesh-a population-based time-series study. International Journal of Epidemiology. 2008, 38: 1689-1697.View Article
- Goodman PG, Dockery DW, Clancy L: Cause-specific mortality and the extended effects of particulate pollution and temperature exposure. Environmental Health Perspectives. 2004, 112: 179-185.View Article
- Deschenes O, Moretti E: Extreme weather events, mortality and migration: nber working paper. In Book Extreme Weather Events, Mortality and Migration: NBER Working Paper (Editor ed.^eds.). 2007, National Bureau of Economic Research, Cambridge, 1-45. 1–45, vol. 13227
- Hajat S, Armstrong B, Gouveia N, Wilkinson P: Mortality displacement of heat-related deaths: a comparison of Delhi, Sao Paulo, and London. Epidemiology. 2005, 16: 613-620. 10.1097/01.ede.0000164559.41092.2a.View Article
- Pattenden S, Nikiforov B, Armstrong BG: Mortality and temperature in Sofia and London. J Epidemiol Community Health. 2003, 57: 628-633. 10.1136/jech.57.8.628.View Article
- Analitis A, Katsouyanni K, Biggeri A, Baccini M, Forsberg B, Bisanti L, Kirchmayer U, Ballester F, Cadum E, Goodman PG: Effects of cold weather on mortality: results from 15 European cities within the PHEWE project. Am J Epidemiol. 2008, 168: 1397-1408. 10.1093/aje/kwn266.View Article
- Keatinge W, Donaldson G, Bucher K, Jendritsky G, Cordioli E, Martinelli M, Dardanoni L, Katsouyanni K, Kunst A, Mackenbach J: Cold exposure and winter mortality from ischaemic heart disease, cerebrovascular disease, respiratory disease, and all causes in warm and cold regions of Europe. Lancet. 1997, 349: 1341-1346.View Article
- Yu W, Mengersen K, Wang X, Ye X, Guo Y, Pan X, Tong S: Daily average temperature and mortality among the elderly: a meta-analysis and systematic review of epidemiological evidence. International Journal of Biometeorology. 2011, 56: 569-581. in pressView Article
- Yu W, Hu W, Mengersen K, Guo Y, Pan X, Connell D, Tong S: Time course of temperature effects on cardiovascular mortality in Brisbane, Australia. Heart. 2011, 97: 1089-1093. 10.1136/hrt.2010.217166.View Article
- Curriero FC, Heiner KS, Samet JM, Zeger SL, Strug L, Patz JA: Temperature and mortality in 11 cities of the eastern United States. Am J Epidemiol. 2002, 155: 80-87. 10.1093/aje/155.1.80.View Article
- Baccini M, Biggeri A, Accetta G, Kosatsky T, Katsouyanni K, Analitis A, Anderson HR, Bisanti L, D'Ippoliti D, Danova J: Heat effects on mortality in 15 European cities. Epidemiology. 2008, 19: 711-719. 10.1097/EDE.0b013e318176bfcd.View Article
- Rocklöv J, Forsberg B: The effect of temperature on mortality in Stockholm 1998—2003: A study of lag structures and heatwave effects. Scandinavian journal of public health. 2008, 36: 516-523. 10.1177/1403494807088458.View Article
- Hajat S, Kosatky T: Heat-related mortality: a review and exploration of heterogeneity. J Epidemiol Community Health. 2010, 64: 753-760. 10.1136/jech.2009.087999.View Article
- El-Zein A, Tewtel-Salem M, Nehme G: A time-series analysis of mortality and air temperature in Greater Beirut. Sci Total Environ. 2004, 330: 71-80. 10.1016/j.scitotenv.2004.02.027.View Article
- Revich B, Shaposhnikov D: Temperature-induced excess mortality in Moscow, Russia. Int J Biometeorol. 2008, 52: 367-374. 10.1007/s00484-007-0131-6.View Article
- Rocklov J, Forsberg B: The effect of temperature on mortality in Stockholm 1998–2003: a study of lag structures and heatwave effects. Scand J Public Health. 2008, 36: 516-523. 10.1177/1403494807088458.View Article
- Yu W, Mengersen K, Hu W, Guo Y, Pan X, Tong S: Assessing the relationship between global warming and mortality: Lag effects of temperature fluctuations by age and mortality categories. Environ Pollut. 2011, 159: 1789-1793. 10.1016/j.envpol.2011.03.039.View Article
- Bouchama A, Knochel JP: Heat stroke. N Engl J Med. 2002, 346: 1978-1988. 10.1056/NEJMra011089.View Article
- McGeehin M, Mirabelli M: The potential impacts of climate variability and change on temperature-related morbidity and mortality in the United States. Environmental Health Perspectives. 2001, 109: 185-189.View Article
- Carder M, McNamee R, Beverland I, Elton R, Cohen GR, Boyd J, Agius RM: The lagged effect of cold temperature and wind chill on cardiorespiratory mortality in Scotland. Occup Environ Med. 2005, 62: 702-710. 10.1136/oem.2004.016394.View Article
- Ballester F, Corella D, Perez-Hoyos S, Saez M, Hervas A: Mortality as a function of temperature, a study in Valencia, Spain, 1991–1993. International Journal of Epidemiology. 1997, 26: 551-561. 10.1093/ije/26.3.551.View Article
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.