Evaluating COVID-19 Risk to Essential Workers by Occupational Group: A Case Study in Massachusetts.

Occupational exposure to SARS-CoV-2 varies by profession, but "essential workers" are often considered in aggregate in COVID-19 models. This aggregation complicates efforts to understand risks to specific types of workers or industries and target interventions, specifically towards non-healthcare workers. We used census tract-resolution American Community Survey data to develop novel essential worker categories among the occupations designated as COVID-19 Essential Services in Massachusetts. Census tract-resolution COVID-19 cases and deaths were provided by the Massachusetts Department of Public Health. We evaluated the association between essential worker categories and cases and deaths over two phases of the pandemic from March 2020 to February 2021 using adjusted mixed-effects negative binomial regression, controlling for other sociodemographic risk factors. We observed elevated COVID-19 case incidence in census tracts in the highest tertile of workers in construction/transportation/buildings maintenance (Phase 1: IRR 1.32 [95% CI 1.22, 1.42]; Phase 2: IRR: 1.19 [1.13, 1.25]), production (Phase 1: IRR: 1.23 [1.15, 1.33]; Phase 2: 1.18 [1.12, 1.24]), and public-facing sales and services occupations (Phase 1: IRR: 1.14 [1.07, 1.21]; Phase 2: IRR: 1.10 [1.06, 1.15]). We found reduced case incidence associated with greater percentage of essential workers able to work from home (Phase 1: IRR: 0.85 [0.78, 0.94]; Phase 2: IRR: 0.83 [0.77, 0.88]). Similar trends exist in the associations between essential worker categories and deaths, though attenuated. Estimating industry-specific risk for essential workers is important in targeting interventions for COVID-19 and other diseases and our categories provide a reproducible and straightforward way to support such efforts.

Impact of meteorology on indoor air quality, energy use, and health in a typical mid-rise multi-family home in the eastern United States.

Heating and cooling requirement differences across climates not only have carbon emissions and energy efficiency implications but also impact indoor air quality (IAQ) and health. Energy and IAQ building simulation models help understand tradeoffs or co-benefits, but these have not been applied to evaluate climate zone or multi-family home differences. We modeled a four-story multi-family home in six U.S. climate zones and quantified energy, IAQ, and health outcomes with EnergyPlus, CONTAM, and a pediatric asthma systems science model. Pollutant sources included cooking and ambient. Outputs were daily PM2.5 and NO2 indoor concentrations, infiltration, energy for heating and cooling, and asthma exacerbations, which were compared across climate zones, apartment units, and resident behaviors. Daily ambient-sourced PM2.5 decreased and cooking-sourced PM2.5 increased with higher ambient temperatures. Infiltration air changes per hour were higher on the first versus the fourth floor and in colder climates. Window opening during cooking led to decreases in total pollutant concentrations (11%-18% for PM2.5 and 9%-15% for NO2 ), 3%-4% decreases in asthma exacerbations within climate zones, and minimal impacts on cooling, but led to increased heating demand (4%-8%). Our results demonstrate the influence of meteorology, multi-family building characteristics, and resident behavior on IAQ, energy, and health, focused on multi-zone methodology.

Time-varying associations between COVID-19 case incidence and community-level sociodemographic, occupational, environmental, and mobility risk factors in Massachusetts.

BACKGROUND: Associations between community-level risk factors and COVID-19 incidence have been used to identify vulnerable subpopulations and target interventions, but the variability of these associations over time remains largely unknown. We evaluated variability in the associations between community-level predictors and COVID-19 case incidence in 351 cities and towns in Massachusetts from March to October 2020. METHODS: Using publicly available sociodemographic, occupational, environmental, and mobility datasets, we developed mixed-effect, adjusted Poisson regression models to depict associations between these variables and town-level COVID-19 case incidence data across five distinct time periods from March to October 2020. We examined town-level demographic variables, including population proportions by race, ethnicity, and age, as well as factors related to occupation, housing density, economic vulnerability, air pollution (PM2.5), and institutional facilities. We calculated incidence rate ratios (IRR) associated with these predictors and compared these values across the multiple time periods to assess variability in the observed associations over time. RESULTS: Associations between key predictor variables and town-level incidence varied across the five time periods. We observed reductions over time in the association with percentage of Black residents (IRR = 1.12 [95%CI: 1.12-1.13]) in early spring, IRR = 1.01 [95%CI: 1.00-1.01] in early fall) and COVID-19 incidence. The association with number of long-term care facility beds per capita also decreased over time (IRR = 1.28 [95%CI: 1.26-1.31] in spring, IRR = 1.07 [95%CI: 1.05-1.09] in fall). Controlling for other factors, towns with higher percentages of essential workers experienced elevated incidences of COVID-19 throughout the pandemic (e.g., IRR = 1.30 [95%CI: 1.27-1.33] in spring, IRR = 1.20 [95%CI: 1.17-1.22] in fall). Towns with higher proportions of Latinx residents also had sustained elevated incidence over time (IRR = 1.19 [95%CI: 1.18-1.21] in spring, IRR = 1.14 [95%CI: 1.13-1.15] in fall). CONCLUSIONS: Town-level COVID-19 risk factors varied with time in this study. In Massachusetts, racial (but not ethnic) disparities in COVID-19 incidence may have decreased across the first 8 months of the pandemic, perhaps indicating greater success in risk mitigation in selected communities. Our approach can be used to evaluate effectiveness of public health interventions and target specific mitigation efforts on the community level.

The impact of energy retrofits on pediatric asthma exacerbation in a Boston multi-family housing complex: a systems science approach.

BACKGROUND: Pediatric asthma is currently the most prevalent chronic disease in the United States, with children in lower income families disproportionately affected. This increased health burden is partly due to lower-quality and insufficient maintenance of affordable housing. A movement towards 'green' retrofits that improve energy efficiency and increase ventilation in existing affordable housing offers an opportunity to provide cost-effective interventions that can address these health disparities. METHODS: We combine indoor air quality modeling with a previously developed discrete event model for pediatric asthma exacerbation to simulate the effects of different types of energy retrofits implemented at an affordable housing site in Boston, MA. RESULTS: Simulation results show that retrofits lead to overall better health outcomes and healthcare cost savings if reduced air exchange due to energy-saving air tightening is compensated by mechanical ventilation. Especially when exposed to indoor tobacco smoke and intensive gas-stove cooking such retrofit would lead to an average annual cost saving of over USD 200, while without mechanical ventilation the same children would have experienced an increase of almost USD 200/year in health care utilization cost. CONCLUSION: The combination of indoor air quality modeling and discrete event modeling applied in this paper can allow for the inclusion of health impacts in cost-benefit analyses of proposed affordable housing energy retrofits.

Continental-scale quantification of landscape values using social media data.

Individuals, communities, and societies ascribe a diverse array of values to landscapes. These values are shaped by the aesthetic, cultural, and recreational benefits and services provided by those landscapes. However, across the globe, processes such as urbanization, agricultural intensification, and abandonment are threatening landscape integrity, altering the personally meaningful connections people have toward specific places. Existing methods used to study landscape values, such as social surveys, are poorly suited to capture dynamic landscape-scale processes across large geographic extents. Social media data, by comparison, can be used to indirectly measure and identify valuable features of landscapes at a regional, continental, and perhaps even worldwide scale. We evaluate the usefulness of different social media platforms-Panoramio, Flickr, and Instagram-and quantify landscape values at a continental scale. We find Panoramio, Flickr, and Instagram data can be used to quantify landscape values, with features of Instagram being especially suitable due to its relatively large population of users and its functional ability of allowing users to attach personally meaningful comments and hashtags to their uploaded images. Although Panoramio, Flickr, and Instagram have different user profiles, our analysis revealed similar patterns of landscape values across Europe across the three platforms. We also found variables describing accessibility, population density, income, mountainous terrain, or proximity to water explained a significant portion of observed variation across data from the different platforms. Social media data can be used to extend our understanding of how and where individuals ascribe value to landscapes across diverse social, political, and ecological boundaries.

Influence of geospatial resolution on sociodemographic predictors of COVID-19 in Massachusetts.

PURPOSE: When studying health risks across a large geographic region such as a state or province, researchers often assume that finer-resolution data on health outcomes and risk factors will improve inferences by avoiding ecological bias and other issues associated with geographic aggregation. However, coarser-resolution data (e.g., at the town or county-level) are more commonly publicly available and packaged for easier access, allowing for rapid analyses. The advantages and limitations of using finer-resolution data, which may improve precision at the cost of time spent gaining access and processing data, have not been considered in detail to date. METHODS: We systematically examine the implications of conducting town-level mixed-effect regression analyses versus census-tract-level analyses to study sociodemographic predictors of COVID-19 in Massachusetts. In a series of negative binomial regressions, we vary the spatial resolution of the outcome, the resolution of variable selection, and the resolution of the random effect to allow for more direct comparison across models. RESULTS: We find stability in some estimates across scenarios, changes in magnitude, direction, and significance in others, and tighter confidence intervals on the census-tract level. Conclusions regarding sociodemographic predictors are robust when regions of high concentration remain consistent across town and census-tract resolutions. CONCLUSIONS: Inferences about high-risk populations may be misleading if derived from town- or county-resolution data, especially for covariates that capture small subgroups (e.g., small racial minority populations) or are geographically concentrated or skewed (e.g., % college students). Our analysis can help inform more rapid and efficient use of public health data by identifying when finer-resolution data are truly most informative, or when coarser-resolution data may be misleading.

Mapping the gaps between cooling benefits of urban greenspace and population heat vulnerability.

We provide a novel method to assess the heat mitigation impacts of greenspace though studying the mechanisms of ecosystems responsible for benefits and connecting them to heat exposure metrics. We demonstrate how the ecosystem services framework can be integrated into current practices of environmental health research using supply/demand state-of-the-art methods of ecological modeling of urban greenspace. We compared the supply of cooling ecosystem services in Boston measured through an indicator of high resolution evapotranspiration modeling, with the demand for benefits from cooling measured as a heat exposure risk score based on exposure, hazard and population characteristics. The resulting evapotranspiration indicator follows a pattern similar to conventional greenspace indicators based on vegetation abundance, except in warmer areas such as those with higher levels of impervious surface. We identified demand-supply mismatch areas across the city of Boston, some coinciding with affordable housing complexes and long term care facilities. This novel ES-framework provides cross-disciplinary methods to prioritize urban areas where greenspace interventions can have the most impact based on heat-related demand.

Community predictors of COVID-19 cases and deaths in Massachusetts: Evaluating changes over time using geospatially refined data.

BACKGROUND: The COVID-19 pandemic has highlighted the need for targeted local interventions given substantial heterogeneity within cities and counties. Publicly available case data are typically aggregated to the city or county level to protect patient privacy, but more granular data are necessary to identify and act upon community-level risk factors that can change over time. METHODS: Individual COVID-19 case and mortality data from Massachusetts were geocoded to residential addresses and aggregated into two time periods: "Phase 1" (March-June 2020) and "Phase 2" (September 2020 to February 2021). Institutional cases associated with long-term care facilities, prisons, or homeless shelters were identified using address data and modeled separately. Census tract sociodemographic and occupational predictors were drawn from the 2015-2019 American Community Survey. We used mixed-effects negative binomial regression to estimate incidence rate ratios (IRRs), accounting for town-level spatial autocorrelation. RESULTS: Case incidence was elevated in census tracts with higher proportions of Black and Latinx residents, with larger associations in Phase 1 than Phase 2. Case incidence associated with proportion of essential workers was similarly elevated in both Phases. Mortality IRRs had differing patterns from case IRRs, decreasing less substantially between Phases for Black and Latinx populations and increasing between Phases for proportion of essential workers. Mortality models excluding institutional cases yielded stronger associations for age, race/ethnicity, and essential worker status. CONCLUSIONS: Geocoded home address data can allow for nuanced analyses of community disease patterns, identification of high-risk subgroups, and exclusion of institutional cases to comprehensively reflect community risk.