Skip to main content

Author: Julia Shipman

Air Quality and Pace of Aging among Older Adults in the United States

Air Quality and Pace of Aging among Older Adults in the United States

Investigators:

Arun Balachandran, Daniel W Belsky

Funding:

NIA R61AG086854 (CACHE)

Data sources:

Measures:

  • Health Measures: The Pace of Aging health measures used in the study is obtained from three blood biomarkers (HbA1c, C-reactive protein, cystatin-C), three physical assessments (diastolic blood pressure, peak-flow lung-function testing, waist circumference), and three functional tests (gait speed, balance, grip strength) available in the US Health and Retirement Study.
  • Aging Measures: Pace of Aging Age for all participants in HRS 2006-2016 who has at least two follow-ups of biomarker data (N=13, 358).
  • Climate Measures: Annual data of 5 exposure during 2002-16.

Project Summary:

Air pollution is emerging as a central public health threat to aging populations. Living in more polluted areas is associated with increased risk for a wide range of aging related diseases. There is emerging evidence that air pollution may accelerate the aging process itself, shortening healthy lifespans in already aging populations. Efforts are underway to reduce pollution and its harms. Metrics to monitor the impact of those efforts on population health are needed. Passively accumulated data such as hospitalizations for asthma or heart attacks are limited because they capture only the tip of the iceberg of latent morbidity caused by pollution. More sensitive and comprehensive measures are needed. If air pollution really does hasten the aging process, new methods to quantify the pace of biological aging could provide the answer.

The central hypothesis of this pilot proposal is that air pollution accelerates the pace of biological aging. We propose a one-year study to test this hypothesis and generate proof-of-concept for a method to monitor population health impacts of air pollution and efforts to reduce/mitigate it. Successful completion of this pilot study will position us to apply for R01 grants to expand our project to global scale, to develop an interactive toolkit for researchers and policymakers can use to evaluate how changes in air pollution levels will impact population aging, and to investigate the role of air pollution in social gradients in biological aging in the US.

Our pilot will analyze data from the US Health and Retirement Study (HRS), an ongoing longitudinal study of adults aged 50 and older and their spouses in the United States. HRS is ongoing since 1992. Survey data are collected every two years. Since 2006, biomarker data are collected every four years. Refresher panels are recruited periodically to replace study members who have died. The HRS has so far collected data on around 40,000 individuals, with roughly 20,00 participating at any given assessment wave. Our analysis will focus on a sample of 13,000 adults for whom we have previously conducted analysis to phenotype Pace of Aging, a longitudinal measurement of the rate of decline in the integrity of multiple bodily systems2. We will link Pace of Aging, sociodemographic, and morbidity and mortality data with small-area air pollution data within the HRS Contextual Data Resource (HRS-CDR) hosted within the Michigan Center for Demography of Aging (MiCDA) accessed via their Virtual Desktop Infrastructure (VDI). HRS-CDR is a collection of analysis-ready datasets that link HRS participant-level data with small-area characteristics relevant to health and wellbeing. Air pollution exposure data consist of average annual concentrations of PM 2.5 (mg/m3) and O3 (mg/m3) at the census-tract level for the period 2002-2016.

On the creation of the weather variables:

The air pollution data are integrated at the census-tract level, within the server of the Health and Retirement Study at Michigan Centre on Demography of Aging (MiCDA). The Virtual Data Enclave (VDE) supported by MiCDA gives the PM 2.5 (mg/m3) and O3 (mg/m3) at the census-tract level for the period 2002-2016, and the researchers made use of this.

Outputs:

  • Poster
  • Future publications and code

Continue reading

Effects of Outdoor Wildfire PM2.5 on Alzheimer’s Disease and Related Dementia

Effects of Outdoor Wildfire PM2.5 on Alzheimer’s Disease and Related Dementia

Investigators:

Jennifer Stowell, Chad Milando, and Greg Wellenius

Funding:

NIA R61AG086854 (CACHE)

Data sources:

Measures:

  • Health Measures: AD/ADRD ED visits and hospitalizations will be identified using the International Classification of Diseases versions 9 and 10.
  • Aging Measures: Age of event is recorded for all patients, and the analysis is restricted to persons aged at least 40 years.
  • Climate Measures: Daily WFS-specific PM2.5, heat, and relative humidity will allow us to examine un-biased associations between wildfire smoke and AD/ADRD. Meteorology will include multiple measures of temperature (i.e. absolute, heat index, wet bulb globe, etc.).

Project Summary:

We will examine the impact of exposure to WFS-specific PM2.5 on emergency department visits and hospitalizations for incident AD/ADRD or exacerbations of AD/. We will link population-weighted exposure, meteorology, and demographic variables to AD/ADRD events across the contiguous US for 2006-2023 using a large medical claims dataset. We will accomplish this using distributed lag nonlinear models (DLNM) and conditional Poisson regression. We will explore multiple lag lengths to account for delays in exposure effects. These analyses will be conducted using a case-control study design where each case is matched to non-case days within the same month, year, and on the same day of the week. This design inherently controls for all time invariant confounders, and all models will include terms for confounding variables such as temperature, relative humidity, and holidays. We will repeat our analyses stratifying on measures of individual and community-level social determinants of health (SDOH) using age, sex, and select ACS variables.

This research will help to increase our understanding of the environmental factors associated with AD/ADRD. Our results will provide actionable evidence for public health practitioners, clinicians, and policymakers in future efforts to mitigate the impacts of climate change on AD/ADRD. Our future research will build on these results and inform an R01 proposal to examine the potential synergistic impacts of multiple extreme weather events (i.e. wildfire, drought, heat, etc.) and mixtures of pollutants on AD/ADRD in US adults.

Outputs:

  • Poster presentations
  • Grant proposal
  • Code
  • Future publications

Continue reading

What am I Reading? A Life Course Approach to Brain Health in a Changing Climate

What am I reading? A Life Course Approach to Brain Health in a Changing Climate

If you’re interested in contributing a short What Am I Reading post, we’d love to hear from you! Email us at cache@colorado.edu

Written by Kelly Perry and Jenna Merenstein

Calling for an exposome-informed approach to brain health across the life course

Older adults are a key group at risk from climate-related threats, including extreme heat, poor air quality, and flooding (EPA). While these hazards were once primarily seen as risks to the heart and lungs, new evidence shows they also play a significant role in brain health, aging, and an increased risk of dementia (Jones et al 2025). We highlight two recent notable studies that explore the links between adverse environmental exposures and impaired brain health and aging, with a broad theme underscoring the need for more equitable, justice-centered, and exposome-informed community-level and occupational interventions to mitigating environment-related acceleration of neurocognitive aging across the life course.

What does a life course, exposome-informed, justice-centered approach to understanding and protecting brain health look like in practice? First, it demands investment in transdisciplinary research that bridges cognitive science, neuroscience, urban planning, and environmental science. Second, it requires designing policies that protect the brain from environmental harm at all stages of life—from improved maternal housing to urban greening strategies for older adults. Third, it calls for integrating neuroimaging-based measures of brain structure and brain function into exposome-focused studies to identify early biomarkers of environmental impact. Doing so would help identify targetable properties for future exposome-based intervention work and ultimately help slow or delay the negative effects of the exposome on neurocognitive aging.

A quick, non-exhaustive review of recent articles discussing adverse environmental exposures and impaired brain health, with implications for an exposome-informed life course approach

Jones and colleagues (2025) conducted an umbrella review and meta-analysis of longitudinal studies that investigated environmental risk factors for dementia, and identified that exposure to the following nine factors was associated with increased relative risk of all-cause dementia (dementia resulting from any combination of underlying causes) compared to those who were unexposed: fine particulate matter (PM; e.g., PM less than or equal to 2.5 µg/m3), particulate matter (PM less than or equal to 10 µg/m3), nitrogen dioxide (NO2), nitrogen oxides (NOx), carbon monoxide (CO), shift work, night shift work, chronic noise, and extremely low-frequency magnetic fields (ELF-MF). They also highlighted community-level factors that were associated with a lower relative risk for dementia, such as neighborhood greenness. Regarding specific types of dementia, Jones and colleagues found the following factors were associated with increased relative risk of dementia of the Alzheimer’s type: PM2.5, ELF-MF, sulfur dioxide (SO2), chronic noise, and pesticides. Similarly, PM2.5, PM10, and chronic noise were associated with increased relative risk of vascular dementia.

Canning and colleagues (2025) followed people from midlife to older age (43–69 years) to study how long-term exposure to air pollution affects the brain. They looked at common pollutants such as PM2.5, PM10, and nitrogen oxides (NOx), and measured cognitive thinking skills and acquired brain scans (i.e., structural magnetic resonance imaging, or MRI). A strength of this study is its long follow-up period and inclusion of adults over 65—an age group that is growing quickly worldwide but often excluded from studies of environmental exposure and neurocognitive aging. The team found that people exposed to more air pollution in mid-to-late life had lower cognitive performance, slower thinking speed, and greater age-related decreases in brain volume. Compared to individuals with lower exposure to NO2, NOx, and PM10, higher exposure to these pollutants was linked to larger brain ventricles and smaller hippocampal volume—and these changes are tied to memory abilities and overall brain health.

Interestingly, the team did not observe links between environmental exposures and measures of verbal memory, or between white matter hyperintensities (a marker for cardiovascular damage). These null findings could suggest that environmental risks may impact only certain aspects of how we think about the underlying neurobiology, although additional research is needed to confirm this. The authors also emphasize that these findings should be viewed within the bigger picture—where environmental exposures interact with genetics, cognition, and brain changes starting before birth and accumulating over the life course.

Researchers are now leveraging the idea of the “exposome”—the sum of environmental exposures an individual experiences across their lifetime—to offer a more holistic lens to brain health. For example, Legaz and colleagues (2025) propose an exposome framework that combines both social factors (e.g., education quality) and physical factors (e.g., air pollution) and connects them to brain outcomes measured with tools such as MRI (as shown by Canning et al., 2025). Legaz and colleagues’ framework highlights how systemic inequities—structures that promulgate an unequal distribution of resources and opportunities in communities—shape the pathways of brain aging. For instance, they report that greater structural inequities are associated with adverse brain outcomes, such as lower brain volume and reduced connectivity among different brain regions, and these outcomes have negative impacts on our cognitive abilities. These impacts are magnified in older adults and those living with dementia. We refer the reader to Legaz et al. for a descriptive figure of this model.

The takeaway for an exposome-informed approach for brain health over the life course

Taking an “exposome” approach means putting environmental justice at the center of studying brain health over the life course. Marginalized communities—often low-income or racially maligned—face the highest exposure to environmental risks while having the least resources to cope. As Legaz and colleagues (2025) emphasize, addressing this imbalance is both ethically imperative and essential to people’s health and wellbeing. Solutions include expanding equitable access to restorative green spaces (Besser et al., 2023), enforcing clean air and water quality standards in underserved areas, and investing in community-led initiatives. Implementing such steps would reduce harmful exposures, build community-level resilience, and promote brain health equity for all.

As air pollution, toxic environmental exposures, and climate change converge with a rapid increase in the global population of older adults (GBD Lancet 2022), we are at a critical juncture: such convergence demands a fundamental rethinking of how we understand and protect brain health across the life course.

Continue reading

Daily Temperature Data Processing and Analysis: An Example for New York City

Daily Temperature Data Processing and Analysis: An Example for New York City

Link to code

Click here

Date: September 2025


Authors/Creators/ Team Members: 

Author: Selen Ozdogan

Team Members: Frank Heiland, Deborah Balk, Jennifer Brite, Peter Marcotullio

Specific purpose of code:

This code aims to provide a comprehensive guide to acquiring and cleaning daily temperature and precipitation data for New York City between 2015-2022 using two primary data sources: Global Historical Climatology Network daily (GHCNd) from the U.S. National Centers for Environmental Information and ERA5-Land Reanalysis from the European Union’s Copernicus Project.

From these data sources, the code assembles daily air temperature and precipitation. It also calculates wet bulb temperature and creates temperature exposure variables with varying temporal resolutions. The extent of this example is New York City (NYC). Aggregation of the input data is necessary to generate estimates for all NYC.  

The code is embedded in an R Markdown pdf file.

General Application:

This is a guide to obtaining climate data and creating different temperature measures and temporal exposure lags. With minor tweaks, the code could be used for other locations/time-periods and can be merged with any daily data set for data analysis. Note that our example here is from 2015-2022, but that time period can be extended (as we also did in the underlying research); a short time period is given in this R Markdown package to facilitate the demonstration.

How does or could this code allow researchers to assess research questions related  to aging or life course?:

The output from this code, daily climate data, could be merged with any daily (or more aggregated temporal frequency) data to study the impact of extreme weather events on aging populations, so long as the underlying spatial resolution of the climate data and population data (from either administrative, census or survey data) are spatially and temporally compatible.

Data sets used: 

  • Climate, weather, disaster or environment data:

    Global Historical Climatology Network daily (GHCNd) – point location format.

    ERA5-Land Reanalysis data – grid format

  • All data are publicly available

Links to data:

  1. Global Historical Climatology Network daily (GHCNd)
  2. Climate Data Store

Coding Language:  R, Python

Tools and Packages used:

R: tidyverse, lubridate, magrittr, here, sf, raster, exactextractr, openxlsx, fixest, slider

Python: os, cdsapi, time, Path

Output(s): Dataset

Spatial extent: New York City (roughly 300 sq. miles or 778 sq. km.)

Temporal extent: 2015-2022

Comments: Replication package for the Demography article will be available here. 

Published papers that use this code:

Forthcoming paper “Extreme Weather and Mortality of Vulnerable Urban Populations:  An Examination of Temperature and Unclaimed Deaths in New York City”, in Demography (2026).

Link to PAA Poster

Related Content: 

Demonstration Project: Impact of Extreme Weather on Hard-to-Capture, Vulnerable Populations: Evidence from Hart Island — New York’s Public Burial Ground

Seminar: Measuring Extreme Temperatures and Thermal Comfort in Aging and Demographic Reseach 

Continue reading

Linking SVRGIS with FEMA Disaster Declarations and Census/ACS

Linking SVRGIS with FEMA Disaster Declarations and Census/ACS

Link to code

Click here

Date: September 2026


Author/Creator: Amy Read

Specific purpose of code: This code provides examples of how to pull, filter, and merge data from the NOAA/NWS Severe Weather GIS Database (SVRGIS), OpenFEMA Disaster Declarations Summaries, and data from the Decennial Census and American Community Survey (ACS). Walkthroughs are provided for A) merging FEMA disaster declarations to SVRGIS tornado paths data based on incident date and location, B) performing a spatial join between event paths and Census geographies (intersection of line and polygon) to identify geographic areas that were exposed to tornadoes, hail, and/or wind during the user-specified timeframe.

General Application: This template can be extended to access and merge other OpenFEMA datasets based on incident (such as Public Assistance or Individual Assistance summary data), and any other data on Census geographic boundaries that is of interest to the researcher. This code also allows the user to identify FEMA disaster declarations for tornado events at smaller geographic levels (tracts, block groups, etc.).

How does or could this code allow researchers to assess research questions related  to aging or life course?: This code could be used with any of the ACS/Census data subset by age group. Since this code focuses on the spatial join between tornado/wind/hail event paths and Census geographies, any other demographic/health datasets tracked by state, county, tract, block group, etc. could be merged by FIPS code into these data the same way one would combine them with Census data alone.

Data sets used: 

  • Population, socioeconomic, or health data: Decennial Census, ACS
  • Climate, weather, disaster or environment data: SVRGIS (Tornadoes, Wind, Hail) and FEMA Disaster Declarations Summaries

Are all the data publicly available or are some restricted-access? All are publicly available.  

Links to data:

Coding Language:  R 

Tools and Packages used: tidycensus, rfema, sf, tidyverse

Output(s): Merged dataset saved to .Rds format

Spatial extent: Contiguous United States

Temporal extent: Example focuses on 2000-2010 but explains how to filter/extend beyond that. SVRGIS data is available from 1950 for tornadoes and from 1955 for hail and wind. FEMA disaster declarations are available from 1953. 

Comments: This is a revised, streamlined, and more generalized version of the code used for the manuscript below. That code is also available on the author’s GitHub.

Published papers that use this code: Read, A. (2025). Repeated disaster and the economic valuation of place: Temporal dynamics of tornado effects on housing prices in the United States, 1980–2010. Population and Environment, 47(3), 29. https://doi.org/10.1007/s11111-025-00502-w

Continue reading

Weathering the Impact: ENSO-Driven Disasters, Power Disruption, and Health Outcomes in Medicare Populations

Weathering the Impact: ENSO-Driven Disasters, Power Disruption, and Health Outcomes in Medicare Populations

 

Investigators:

Sara Curran, Jeff Stanaway, Luanne Thompson, June Yang, Emmanuela Gakidou

Data sources:

ENSO indices from the National Oceanic and Atmospheric Administration Climate Prediction Center from 1980-2025.

SHELDUS: Special hazard events and losses database for the United States 1970 – 2023. Available by County and day.

Center for Medicare & Medicaid Services claim and beneficiary data from 1990-2023. Available by County and day.

HHS emPOWER: Count of medicare beneficiaries who rely on electricity-dependent durable medical and assistive equipment and devices from 2016-2025. Available by County and month.

DOE-417 archive data: Electric Emergency Incident and Disturbance Report from 2000 to 2023. Available by County and Day.

Spatial Coverage: United States

Temporal Coverage: 2016 – 2023.

Measures:

  • Climate Measurements:

    • ENSO Indices (NOAA): Seasonal measurements of sea surface temperature anomalies used to classify El Niño and La Niña events.
    • SHELDUS: Provides detailed records of natural disaster events (including type, date, and location), allowing alignment of local disasters with ENSO periods. It also includes loss data (e.g., property damage, crop loss, injuries) to quantify event severity.

    Power Infrastructure Measurements:

    • HHS emPOWER Dataset: Monthly, county-level counts of Medicare beneficiaries who rely on electricity-dependent durable medical equipment (DME), indicating vulnerability to power disruptions.
    • DOE OE-417 Archive: Incident-level reports on electrical disturbances, including type of event, date/time, geographic area, and population affected.

    Health Outcomes:

    • CMS Claims and Quality Measures: Data on Medicare beneficiaries covering:
      • Hospital readmission rates
      • Emergency department (ED) and inpatient utilization
      • Home health services
      • Minimum Data Set (MDS) indicators from nursing homes
      • Preventable hospitalizations (e.g., ambulatory care-sensitive conditions)

    Demographic and Socioeconomic Context:

    • County-Level Medicare Population Characteristics: Including race/ethnicity, age, sex, and dual eligibility for Medicare and Medicaid, enabling analysis of differential impacts across vulnerable groups.

Project Summary:

Disruptions driven by the El Niño-Southern Oscillation (ENSO), such as flooding, drought, and extreme temperatures, have long been implicated in public health threats across the United States. For example, during negative phases of ENSO (La Nina events) there is typically an increase in the number of hurricanes that can result in major floods. This also occurs during El Nino events in the Southwest US.  Major floods have been shown to elevate hospitalization rates among older adults for skin conditions, neurological illnesses, musculoskeletal disorders, and injuries in the weeks following exposure (Aggarwal et al. 2025). In addition, other impacts include health care disruptions due to displacement from homes and communities or through damage to health infrastructure. These might be particularly important for those requiring frequent care for chronic illnesses.

However,  mechanisms linking ENSO-related disasters with health outcomes remains poorly understood. For instance, the role of power infrastructure failure, such as prolonged outages disrupting electricity-dependent medical care is a plausible pathway for detrimental impacts on healthcare particularly for the elderly who are less mobile than the general population..

Our U.S.-based project investigates how ENSO-related disasters influence health outcomes among older adults, focusing especially on power disruptions as a critical explanatory pathway. By aligning the data sources mentioned above, we aim to trace if and how power outages amplify the health impact of ENSO-driven disasters on older adults.

We will employ a combination of panel regression models, Difference-in-Difference analysis, and causal mediation analysis to estimate the direct and indirect effects of ENSO-driven disasters, to quantify the role of infrastructure failure and identify populations at higher risk.

Comments

Papers on impacts of El Nino in the U.S. 

Outputs:

Peer-reviewed publications, grant proposals, conference presentations

References:

Aggarwal, Sarika, Jie K. Hu, Jonathan A. Sullivan, Robbie M. Parks, and Rachel C. Nethery. 2025. “Severe Flooding and Cause-Specific Hospitalisation among Older Adults in the USA: A Retrospective Matched Cohort Analysis.” The Lancet Planetary Health 9(7):101268. doi:10.1016/S2542-5196(25)00132-9.

Fussell, Elizabeth, Sara R. Curran, Matthew D. Dunbar, Michael A. Babb, Luanne Thompson, and Jacqueline Meijer-Irons. 2017. “Weather-Related Hazards and Population Change: A Study of Hurricanes and Tropical Storms in the United States, 1980–2012.” The ANNALS of the American Academy of Political and Social Science 669(1):146–67. doi:10.1177/0002716216682942.

Salas, Renee N., Laura G. Burke, Jessica Phelan, Gregory A. Wellenius, E. John Orav, and Ashish K. Jha. 2024. “Impact of Extreme Weather Events on Healthcare Utilization and Mortality in the United States.” Nature Medicine 30(4):1118–26. doi:10.1038/s41591-024-02833-x.

Continue reading

What Am I Reading: Disasters and Aging in Place

What am I reading? Disasters and Aging in Place

Link to article

If you’re interested in contributing a short What Am I Reading post, we’d love to hear from you! Email us at cache@colorado.edu

Written by Jenna Tipaldo, CUNY School of Public Health and CUNY Institute for Demographic Research, jenna.tipaldo09@sphmail.cuny.edu

A new report from Winkler & Mockrin (2025) entitled “Aging and wildfire risk to communities” explores the exposure of older populations to wildfires. A main findings is that “Nearly all (87 percent) of the population growth in higher wildfire risk locations between 2010 and 2020 was among people over the age of 60, many of whom had been living in higher risk places for years and are growing older (i.e., aging in place).” Each is relevant when thinking about older adults and exposure, vulnerability, and resilience to disasters.  Wildfires are just one of type of disaster to consider: recent evidence also suggests that coastal zones – areas at risk of storms and other seaward hazards such as flooding and tsunamis – are also aging faster than areas farther inland (Bukvic et al., 2018; Hauer et al., 2020; Tagtachian and Balk, 2023). 

No relocation: Forsyth & Molinsky (2020) note that to some, aging in place signifies remaining in their home while to others it may mean moving but within the same community, such as downsizing. Based on recent shifts in age distribution in Census blocks with “moderate-to-high” wildfire risk, Winkler & Mockrin (2025) conclude that the increase in older adult populations in fire-prone regions is likely attributable to populations aging in place rather than in-migration (Figure 4). They also note important spatial variation and also uncertainty about the relative contributions of migration and death. They also note that aging in place seems to be the “primary mechanism” in higher risk rural areas (Winkler & Mockrin, 2025). 

Source: Winkler & Mockrin (2025)

Health and Health Care : Winkler & Mockrin (2025) summarize the various ways in which older adults can be at higher risk due to wildfires including 1) physical limitations that are barriers to preparation or response, 2) factors like social isolation which can impact access to information and resources, and 3) higher rates of chronic diseases which are risk factors for adverse health outcomes due to fires and smoke. Furthermore, disasters can be disruptive to healthcare, not only by damaging facilities and displacing people from their homes but also by disrupting care which relies on movement. Examples include when patients are unable to travel to hospitals or medical providers, or if healthcare workers can‘t get to a patient’s home due to inaccessible roads (Tarabochia‐Gast et al., 2022) or suspended public transit systems. Rural areas face additional challenges with longer travel times for healthcare access, especially with high levels of hospital closures (Miler et al., 2020; McCarthy et al., 2021). Such patterns negatively impact health care access, emergency medical response, and transport times (GAO, 2021; Kaufman et al., 2016). On average, rural residents must travel about 20 miles farther for typical health care services – in non-disaster times (GAO 2021). While those miles may seem trivial, in emergencies they can mean loss of access to care and treatment. 

Personal choice: Aging in place can be a personal choice in support of maintaining one’s agency and independence by staying in one’s own home and community (Forsyth & Molinsky, 2020). Even so, staying in one’s home can also result from lack of choice due to limited resources and/or few desirable and affordable options. Modifications are expensive too. Even older adults who are relatively better off can struggle to pay for downsizing or modifying a new dwelling for care needs (Forsyth & Molinsky, 2020).  

From research to policy 

To help support healthy aging in place, Winkler & Mockrin (2025) suggest that existing programs that support older adults could be expanded to include wildfire risk reduction. An example is the USDA’s Section 504 Home Repair program which supports older low-income homeowners. In addition, organizations such as the AARP provide useful material for aging in place such as a checklist for people who are prepping their home. Such resources should be expanded to include disaster risk as a consideration.  

 

References:  

  • Bukvic, A., Gohlke, J., Borate, A., and Suggs, J. 2018. “Aging in Flood-Prone Coastal Areas: Discerning the Health and Well-Being Risk for Older Residents.” International Journal of Environmental Research and Public Health 15(12):2900. https://doi.org/10.3390/ijerph15122900.   
  • Hauer, Mathew E., Elizabeth Fussell, Valerie Mueller, Maxine Burkett, Maia Call, Kali Abel, Robert McLeman, and David Wrathall. 2020. “Sea-Level Rise and Human Migration.” Nature Reviews Earth & Environment 1(1):28–39. https://doi.org/10.1038/s43017-019-0002-9 
  • Kaufman, B.G., Thomas, S.R., Randolph, R.K., et al. The rising rate of rural hospital closures. The Journal of Rural Health. 2016;32(1):35-43. https://doi.org/10.1111/jrh.12128  
  • McCarthy, S., Moore, D., Smedley, W. A., Crowley, B. M., Stephens, S. W., Griffin, R. L., Tanner, L. C., & Jansen, J. O. (2021). Impact of Rural Hospital Closures on Health-Care Access. Journal of Surgical Research, 258, 170–178. https://doi.org/10.1016/j.jss.2020.08.055 
  • Miller, K.E.M., James, H.J., Holmes, G.M., Van Houtven, C.H. The effect of rural hospital closures on emergency medical service response and transport times. Health Serv Res. 2020;55(2):288-300. https://doi.org/10.1111/1475-6773.13254  
  • Tagtachian, D. and Balk, D., 2023. Uneven vulnerability: characterizing population composition and change in the low elevation coastal zone in the United States with a climate justice lens, 1990–2020. Frontiers in Environmental Science, 11, p.1111856. 
  • Tarabochia‐Gast, A. T., Michanowicz, D. R., & Bernstein, A. S. (2022). Flood Risk to Hospitals on the United States Atlantic and Gulf Coasts From Hurricanes and Sea Level Rise. GeoHealth, 6(10), e2022GH000651. https://doi.org/10.1029/2022GH000651 
  • Winkler, R. L., & Mockrin, M. H. (2025). Aging and wildfire risk to communities (Report No. EIB-284). U.S. Department of Agriculture, Economic Research Service. https://doi.org/10.32747/2025.9015828.ers 

 

Continue reading

Severe Heat Days using the Universal Thermal Comfort Index

Severe Heat Days using the Universal Thermal Comfort Index

Link to code

Click here

Date: September 2025


Authors/Creators/ Team Members: Marcial Yangali (marcialy@colmex.mx)

Reviewers: Landy Sánchez and Emerson Baptista

Specific purpose of code: It demonstrates how to construct severe heat measures using the UTCI data (ERA5-HEAT ). It starts by showing data manipulation from raster (grid data) to a tabular dataset that obtains UTCI values for each municipality in Mexico. Then, it presents how to map and analyze such data. Finally, it shows how to count the number of days of severe heat (32°C UTCI and above). This script is part of the demonstration CACHE project “Heat, Disability in older adults and Care” from El Colegio de Mexico.

General Application: This develops familiarity with heat data and fundamental skills to construct an environmental dataset that can be easily integrate with demographic information.

How does or could this code allow researchers to assess research questions related  to aging or life course?: This code help to analyze exposure to severe heat days by age structure at a subnational level. Particularly, to understand how aging realtes to higher exposure to extreme weather. Historical UTCI measures could be employ to assess cumulative heat impacts across the life course.

Data sets used: 

  • Two publicy available datasets are used:

    • Universal Thermal Climate Index (UTCI). Copernicous ERA5-HEAT
    • 2020 Mexican Municipalities and States boundaries. INEGI (Census Bureau office)

Coding Language:  R 

Tools and Packages used: R packages for data manipulation (mainly dplyr but also janitor, stringr, tidyr and forcats), visualization (ggplot2, patchwork), and geospatial operations (sf).

Output(s): Analysis results, maps and graphs

Spatial extent:Mexico (33.0°N (North), -118.5°E (East), 14.0°N (South), and -86.0°E (West)

Temporal extent: May 2019

Key words: heat, severe weather, thermal index, Mexico

Continue reading

Aging and disability in the Mexican Population

Aging and disability in the Mexican Population

Link to code

Click here

Date: September 2025


Authors/Creators/ Team Members: Marcial Yangali (marcialy@colmex.mx) and Mariana Ramos

Reviewers: Emerson Baptista and Landy Sánchez

Specific purpose of code: To show how to construct disability measures using the international recommendation of the Washington Group. It also demonstrates how to evaluate age and sex composition of the population with disability and their territorial distribution.

General Application: This demonstrates basic manipulation of demographic data with R, particularly useful for those with limited familiarity with population measures.

How does or could this code allow researchers to assess research questions related  to aging or life course?: This code examines age structure in census data. Specifically, explores how the prevalence of disability condition increases with age.

Data sets used: 

  • All are publicy available datasets:

    • 2020 Mexican Housing and Population Data (IPUMS International)

Coding Language:  R 

Tools and Packages used: R packages for data manipulation (mainly dplyr but also stringr, ipumsr, tidyr, srvyr, gtsummary, purr, and labelled), visualization (ggplot2, patchwork, geofacet, scales and plotly), and geospatial operations (sf, biscale, ggtern).

Output(s): Analysis results, maps and graphs

Spatial extent: Mexico (33.0°N (North), -118.5°E (East), 14.0°N (South), and -86.0°E (West)

Temporal extent: 2022

Key words: disability, Mexico, age structure, sex composition, aging

Continue reading

Joining ACAG Annual Estimates of PM2.5 with Social Determinants of Health (SDOH) data

Joining ACAG Annual Estimates of PM2.5 with Social Determinants of Health (SDOH) data

Link to code

Click here

Date: July 2025


Authors/Creators/Team Members: Zoé Haskell-Craig and Priyanka deSouza

Specific purpose of code: This code aggregates gridded annual average PM2.5 concentration estimates produced by the Atmospheric Composition Analysis Group (ACAG) to the census tract level, producing a variable containing the average PM2.5 exposure for each tract. This is then combined with socioeconomic and demographic data available at the tract level from the social determinants of health (SDOH) database produced by the Agency for Healthcare Research and Quality (AHRQ).

General Application: This code takes advantage of the `tigris` package to aggregate high resolution (fine spatial scale) modelled estimates of PM2.5 pollution to the administrative boundaries at which demographic and SDOH data are available. As an example, here we demonstrate computing the annual average PM2.5 concentrations in 2020 for census tracts and combining this with SDOH data on race/ethnicity and income. With minor changes, this code can be used for other years and temporal resolutions (i.e. monthly estimates of PM2.5) and for other administrative units (ZCTAs, blockgroups, counties, etc).

How does or could this code allow researchers to assess research questions related to aging or life course?: While the dataset output from this code does not contain variables on age, the raw SDOH dataset contains census information on age which could be included with minor edits to the code. Also, the aggregation of PM2.5 exposure to the census tract (or other administrative units) allows researchers to combine this exposure with other information available from the census, such as income by age breakdowns.

Data sets used: 

Are all the data publicly available or are some restricted-access? Publicly available 

Coding Language:  R 

Tools and Packages used: tidyverse, readxl, sf, raster, ncdf4, exactextractr, tigris, viridis

Output(s): Dataset with census tracts as unit of analysis and a map of the average PM2.5 per census tract in 2020. 

Spatial extent: Continental US

Temporal extent: Single-year, 2020 (code can be modified to produce data for any year from 1998 – 2023, or monthly for any month in that period). 

Additional Comments: Journal article using this code is forthcoming. 

Published papers that use this code: Zoé Haskell-Craig, Kevin P. Josey, Patrick L. Kinney, and Priyanka deSouza. (2025). Equity in the Distribution of Regulatory PM2.5 Monitors. Environmental Science & Technology. DOI: 10.1021/acs.est.4c12915

Continue reading