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.

Inequities in COVID-19 vaccine and booster coverage across Massachusetts ZIP codes after the emergence of Omicron: A population-based cross-sectional study.

BACKGROUND: Inequities in Coronavirus Disease 2019 (COVID-19) vaccine and booster coverage may contribute to future disparities in morbidity and mortality within and between Massachusetts (MA) communities. METHODS AND FINDINGS: We conducted a population-based cross-sectional study of primary series vaccination and booster coverage 18 months into the general population vaccine rollout. We obtained public-use data on residents vaccinated and boosted by ZIP code (and by age group: 5 to 19, 20 to 39, 40 to 64, 65+) from MA Department of Public Health, as of October 10, 2022. We constructed population denominators for postal ZIP codes by aggregating census tract population estimates from the 2015-2019 American Community Survey. We excluded nonresidential ZIP codes and the smallest ZIP codes containing 1% of the state's population. We mapped variation in ZIP code-level primary series vaccine and booster coverage and used regression models to evaluate the association of these measures with ZIP code-level socioeconomic and demographic characteristics. Because age is strongly associated with COVID-19 severity and vaccine access/uptake, we assessed whether observed socioeconomic and racial/ethnic inequities persisted after adjusting for age composition and plotted age-specific vaccine and booster coverage by deciles of ZIP code characteristics. We analyzed data on 418 ZIP codes. We observed wide geographic variation in primary series vaccination and booster rates, with marked inequities by ZIP code-level education, median household income, essential worker share, and racial/ethnic composition. In age-stratified analyses, primary series vaccine coverage was very high among the elderly. However, we found large inequities in vaccination rates among younger adults and children, and very large inequities in booster rates for all age groups. In multivariable regression models, each 10 percentage point increase in "percent college educated" was associated with a 5.1 (95% confidence interval (CI) 3.9 to 6.3, p < 0.001) percentage point increase in primary series vaccine coverage and a 5.4 (95% CI 4.5 to 6.4, p < 0.001) percentage point increase in booster coverage. Although ZIP codes with higher "percent Black/Latino/Indigenous" and higher "percent essential workers" had lower vaccine coverage (-0.8, 95% CI -1.3 to -0.3, p < 0.01; -5.5, 95% CI -7.3 to -3.8, p < 0.001), these associations became strongly positive after adjusting for age and education (1.9, 95% CI 1.0 to 2.8, p < 0.001; 4.8, 95% CI 2.6 to 7.1, p < 0.001), consistent with high demand for vaccines among Black/Latino/Indigenous and essential worker populations within age and education groups. Strong positive associations between "median household income" and vaccination were attenuated after adjusting for age. Limitations of the study include imprecision of the estimated population denominators, lack of individual-level sociodemographic data, and potential for residential ZIP code misreporting in vaccination data. CONCLUSIONS: Eighteen months into MA's general population vaccine rollout, there remained large inequities in COVID-19 primary series vaccine and booster coverage across MA ZIP codes, particularly among younger age groups. Disparities in vaccination coverage by racial/ethnic composition were statistically explained by differences in age and education levels, which may mediate the effects of structural racism on vaccine uptake. Efforts to increase booster coverage are needed to limit future socioeconomic and racial/ethnic disparities in COVID-19 morbidity and mortality.

Assessing the impact of aircraft arrival on ambient ultrafine particle number concentrations in near-airport communities in Boston, Massachusetts.

Aircraft emissions contribute to overall ambient air pollution, including ultrafine particle (UFP) concentrations. However, accurately ascertaining aviation contributions to UFP is challenging due to high spatiotemporal variability along with intermittent aviation emissions. The objective of this study was to evaluate the impact of arrival aircraft on particle number concentration (PNC), a proxy for UFP, across six study sites 3-17 km from a major arrival aircraft flight path into Boston Logan International Airport by utilizing real-time aircraft activity and meteorological data. Ambient PNC at all monitoring sites was similar at the median but had greater variation at the 95th and 99th percentiles with more than two-fold increases in PNC observed at sites closer to the airport. PNC was elevated during the hours with high aircraft activity with sites closest to the airport exhibiting stronger signals when downwind from the airport. Regression models indicated that the number of arrival aircraft per hour was associated with measured PNC at all six sites, with a maximum contribution of 50% of total PNC at a monitor 3 km from the airport during hours with arrival activity on the flight path of interest (26% across all hours). Our findings suggest strong but intermittent contributions from arrival aircraft to ambient PNC in communities near airports.

Critical windows of susceptibility for the effects of prenatal exposure to heat and heat variability on gestational growth.

BACKGROUND: Studies have shown that prenatal heat exposure may impact fetal growth, but few studies have examined the critical windows of susceptibility. As extreme heat events and within season temperature variability is expected to increase in frequency, it is important to understand how this may impact gestational growth. OBJECTIVES: We investigated associations between various measures of weekly prenatal heat exposure (mean and standard deviation (SD) of temperature and heat index (HI), derived using temperature in °C and dew point) and term birthweight or odds of being born small for gestational age (SGA) to identify critical windows of susceptibility. METHODS: We analyzed data from mother-child dyads (n = 4442) in the Boston-based Children's HealthWatch cohort. Birthweights were collected from survey data and electronic health records. Daily temperature and HI values were obtained from 800 m gridded spatial climate datasets aggregated by the PRISM Climate Group. Distributed lag-nonlinear models were used to assess the effect of the four weekly heat metrics on measures of gestational growth (birthweight, SGA, and birthweight z-scores). Analyses were stratified by child sex and maternal homelessness status during pregnancy. RESULTS: HI variability was significantly associated with decreased term birthweight during gestational weeks 10-29 and with SGA for weeks 9-26. Cumulative effects for these time periods were -287.4 g (95% CI: -474.1 g, -100.8 g for birthweight and 4.7 (95% CI: 1.6, 14.1) for SGA. Temperature variability was also significantly associated with decreased birthweight between weeks 15 and 26. The effects for mean heat measures on term birthweight and SGA were not significant for any gestational week. Stratification by sex revealed a significant effect on term birthweight in females between weeks 23-28 and in males between weeks 9-26. Strongest effects of HI variability on term birthweight were found in children of mothers who experienced homelessness during pregnancy. Weekly HI variability was the heat metric most strongly associated with measures of gestational growth. The effects observed were largest in males and those who experienced homelessness during pregnancy. DISCUSSION: Given the impact of heat variability on birthweight and risk of SGA, it is important for future heat warnings to incorporate measure of heat index and temperature variability.

Fine particulate matter and polycystic ovarian morphology.

BACKGROUND: Polycystic ovary morphology (PCOM) is an ultrasonographic finding that can be present in women with ovulatory disorder and oligomenorrhea due to hypothalamic, pituitary, and ovarian dysfunction. While air pollution has emerged as a possible disrupter of hormone homeostasis, limited research has been conducted on the association between air pollution and PCOM. METHODS: We conducted a longitudinal cohort study using electronic medical records data of 5,492 women with normal ovaries at the first ultrasound that underwent a repeated pelvic ultrasound examination during the study period (2004-2016) at Boston Medical Center. Machine learning text algorithms classified PCOM by ultrasound. We used geocoded home address to determine the ambient annual average PM2.5 exposures and categorized into tertiles of exposure. We used Cox Proportional Hazards models on complete data (n = 3,994), adjusting for covariates, and additionally stratified by race/ethnicity and body mass index (BMI). RESULTS: Cumulative exposure to PM2.5 during the study ranged from 4.9 to 17.5 µg/m3 (mean = 10.0 µg/m3). On average, women were 31 years old and 58% were Black/African American. Hazard ratios and 95% confidence intervals (CI) comparing the second and third PM2.5 exposure tertile vs. the reference tertile were 1.12 (0.88, 1.43) and 0.89 (0.62, 1.28), respectively. No appreciable differences were observed across race/ethnicity. Among women with BMI ≥ 30 kg/m2, we observed weak inverse associations with PCOM for the second (HR: 0.93, 95% CI: 0.66, 1.33) and third tertiles (HR: 0.89, 95% CI: 0.50, 1.57). CONCLUSIONS: In this study of reproductive-aged women, we observed little association between PM2.5 concentrations and PCOM incidence. No dose response relationships were observed nor were estimates appreciably different across race/ethnicity within this clinically sourced cohort.

Long-term Exposure to PM2.5 and Mortality for the Older Population: Effect Modification by Residential Greenness.

BACKGROUND: Although many studies demonstrated reduced mortality risk with higher greenness, few studies examined the modifying effect of greenness on air pollution-health associations. We evaluated residential greenness as an effect modifier of the association between long-term exposure to fine particles (PM2.5) and mortality. METHODS: We used data from all Medicare beneficiaries in North Carolina (NC) and Michigan (MI) (2001-2016). We estimated annual PM2.5 averages using ensemble prediction models. We estimated mortality risk per 1 µg/m3 increase using Cox proportional hazards modeling, controlling for demographics, Medicaid eligibility, and area-level covariates. We investigated health disparities by greenness using the Normalized Difference Vegetation Index with measures of urbanicity and socioeconomic status. RESULTS: PM2.5 was positively associated with mortality risk. Hazard ratios (HRs) were 1.12 (95% confidence interval (CI) = 1.12 to 1.13) for NC and 1.01 (95% CI = 1.00 to 1.01) for MI. HRs were higher for rural than urban areas. Within each category of urbanicity, HRs were generally higher in less green areas. For combined disparities, HRs were higher in low greenness or low SES areas, regardless of the other factor. HRs were lowest in high-greenness and high-SES areas for both states. CONCLUSIONS: In our study, those in low SES and high-greenness areas had lower associations between PM2.5 and mortality than those in low SES and low greenness areas. Multiple aspects of disparity factors and their interactions may affect health disparities from air pollution exposures. Findings should be considered in light of uncertainties, such as our use of modeled PM2.5 data, and warrant further investigation.

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.

Descriptive characterization of sound levels in an environmental justice city before and during a global pandemic.

Many environmental justice communities face elevated exposures to multiple stressors, given biases in urban and environmental policy and planning. This paper aims to evaluate sound level exposure in a densely populated environmental justice city in close proximity to major roadways, a nearby airport and high levels of industrial activity. In this study we collected various sound level metrics to evaluate the loudness and frequency composition of the acoustical environment in Chelsea, Massachusetts, USA. A total of 29 week-long sites were collected from October 2019 to June 2020, a time period that also included the influence of the COVID-19 pandemic, which drastically altered activity patterns and corresponding sound level exposures. We found that Chelsea is exposed to high levels of sound, both day and night (65 dB (A), and 80 dB and 90 dB for low frequency, and infrasound sound levels). A spectral analysis shows that 63 Hz was the dominant frequency. Distance to major roads and flight activity (both arrivals and departures) were most strongly correlated with all metrics, most notably with metrics describing contributing from lower frequencies. Overall, we found similar patterns during the COVID-19 pandemic but at levels up to 10 dB lower. Our results demonstrate the importance of noise exposure assessments in environmental justice communities and the importance of using additional metrics to describe communities inundated with significant air, road, and industrial sound levels. It also provides a snapshot of how much quieter communities can be with careful and intentional urban and environmental policy and planning.

Multimodal Recruitment to Study Ovulation and Menstruation Health: Internet-Based Survey Pilot Study.

BACKGROUND: Multimodal recruitment strategies are a novel way to increase diversity in research populations. However, these methods have not been previously applied to understanding the prevalence of menstrual disorders such as polycystic ovary syndrome. OBJECTIVE: The purpose of this study was to test the feasibility of recruiting a diverse cohort to complete a web-based survey on ovulation and menstruation health. METHODS: We conducted the Ovulation and Menstruation Health Pilot Study using a REDCap web-based survey platform. We recruited 200 women from a clinical population, a community fair, and the internet. RESULTS: We recruited 438 women over 29 weeks between September 2017 and March 2018. After consent and eligibility determination, 345 enrolled, 278 started (clinic: n=43; community fair: n=61; internet: n=174), and 247 completed (clinic: n=28; community fair: n=60; internet: n=159) the survey. Among all participants, the median age was 25.0 (SD 6.0) years, mean BMI was 26.1 kg/m2 (SD 6.6), 79.7% (216/271) had a college degree or higher, and 14.6% (37/254) reported a physician diagnosis of polycystic ovary syndrome. Race and ethnicity distributions were 64.7% (176/272) White, 11.8% (32/272) Black/African American, 7.7% (21/272) Latina/Hispanic, and 5.9% (16/272) Asian individuals; 9.9% (27/272) reported more than one race or ethnicity. The highest enrollment of Black/African American individuals was in clinic (17/42, 40.5%) compared to 1.6% (1/61) in the community fair and 8.3% (14/169) using the internet. Survey completion rates were highest among those who were recruited from the internet (159/174, 91.4%) and community fairs (60/61, 98.4%) compared to those recruited in clinic (28/43, 65.1%). CONCLUSIONS: Multimodal recruitment achieved target recruitment in a short time period and established a racially diverse cohort to study ovulation and menstruation health. There were greater enrollment and completion rates among those recruited via the internet and community fair.

Impacts of high temperature on adverse birth outcomes in Seoul, Korea: Disparities by individual- and community-level characteristics.

BACKGROUND: Few studies have examined temperature's effect on adverse birth outcomes and relevant effect modifiers. OBJECTIVES: We investigated associations between heat and adverse birth outcomes and how individual and community characteristics affect these associations for Seoul, Korea, 2004-2012. METHODS: We applied logistic regression to estimate associations between heat index during pregnancy, 4 weeks before delivery, and 1 week before delivery and risk of preterm birth and term low birth weight. We investigated effect modification by individual (infant's sex, mother's age, and mother's educational level) and community characteristics (socioeconomic status (SES) and percentage of green areas near residence at the gu level, which is similar to borough in Western countries). We also evaluated associations by combinations of individual- and community-level SES. RESULTS: Heat exposure during whole pregnancy was significantly associated with risk of preterm birth. An interquartile (IQR) increase (5.5 °C) in heat index during whole pregnancy was associated with an odds ratio (OR) of 1.033 (95% CI 1.005, 1.061) with NO2 adjustment, and 1.028 (95% CI 0.998, 1.059) with PM10 adjustment, for preterm birth. We also found significant associations with heat exposure during 4 weeks before delivery and 1 week before delivery on preterm birth. We did not observe significant associations with term low birth weight. Higher risk of heat on preterm birth was associated with some individual characteristics such as infants with younger or older mothers and lower community-level SES. For combinations of individual- and community-level SES, the highest and most significant estimated effect was found for infants with low educated mothers living in low SES communities, with suggestions of effects of both individual-and community-level SES. CONCLUSIONS: Our findings have implications for evaluating impacts of high temperatures on birth outcomes, estimating health impacts of climate change, and identifying which subpopulations and factors are most relevant for disparities in this association.

Urban environment as an independent predictor of insulin resistance in a South Asian population.

BACKGROUND: Developing countries, such as India, are experiencing rapid urbanization, which may have a major impact on the environment: including worsening air and water quality, noise and the problems of waste disposal. We used health data from an ongoing cohort study based in southern India to examine the relationship between the urban environment and homeostasis model assessment of insulin resistance (HOMA-IR). METHODS: We utilized three metrics of urbanization: distance from urban center; population density in the India Census; and satellite-based land cover. Restricted to participants without diabetes (N = 6350); we built logistic regression models adjusted for traditional risk factors to test the association between urban environment and HOMA-IR. RESULTS: In adjusted models, residing within 0-20 km of the urban center was associated with an odds ratio for HOMA-IR of 1.79 (95% CI 1.39, 2.29) for females and 2.30 (95% CI 1.64, 3.22) for males compared to residing in the furthest 61-80 km distance group. Similar statistically significant results were identified using the other metrics. CONCLUSIONS: We identified associations between urban environment and HOMA-IR in a cohort of adults. These associations were robust using various metrics of urbanization and adjustment for individual predictors. Our results are of public health concern due to the global movement of large numbers of people from rural to urban areas and the already large burden of diabetes.

Urban vegetation and heat-related mortality in Seoul, Korea.

Urban areas are particularly vulnerable to heat-related health outcomes. Simultaneous trends of climate change and urbanization may increase the urban heat-related health burden. We investigated the effects of urban vegetation on heat-related mortality, and evaluated whether different levels of vegetation and individuals' characteristics affect the temperature-mortality associations within Seoul, Korea 2000-2009. We used Normalized Difference Vegetation Index (NDVI) to assess the urban vegetation within Seoul. We applied an overdispersed Poisson generalized linear model with interaction term between temperature and indicator of NDVI group (categorized in 3 levels) to assess the effect modification of the temperature-mortality association by urban vegetation. We conducted stratified analysis to explore whether associations are affected by individual characteristics of sex and age. The association between total mortality and a 1°C increase in temperature above the 90th percentile (25.1°C) (the "heat effect") was the highest for gus with low NDVI. The heat effect was a 4.1% (95% confidence interval (CI) 2.3, 5.9%), 3.0% (95% CI 0.2, 5.9%), and 2.2% (95% CI -0.5, 5.0%) increase in mortality risk for low, medium, and high NDVI group, respectively. Estimated risks showed similar effects by sex and age. Our findings suggest a higher mortality effect of high temperature in areas with lower vegetation in Seoul, Korea.

Climate stress and anxiety, environmental context, and civic engagement: A nationally representative study.

There is increasing recognition that people are experiencing stress and anxiety around climate change, and that this climate stress/anxiety may be associated with more pro-environmental behavior. However, less is known about whether people's own environmental exposures affect climate stress/anxiety or the relationship between climate stress/anxiety and civic engagement. Using three waves of survey data (2020-2022) from the nationally representative Tufts Equity in Health, Wealth, and Civic Engagement Study of US adults (n = 1071), we assessed relationships among environmental exposures (county-level air pollution, greenness, number of toxic release inventory sites, and heatwaves), self-reported climate stress/anxiety, and civic engagement measures (canvasing behavior, collaborating to solve community problems, personal efficacy to solve community problems, group efficacy to solve community problems, voting behavior). Most participants reported experiencing climate stress/anxiety (61%). In general, the environmental exposures we assessed were not significantly associated with climate stress/anxiety or civic engagement metrics, but climate stress/anxiety was positively associated with most of the civic engagement outcomes (canvassing, personal efficacy, group efficacy, voter preference). Our results support the growing literature that climate stress/anxiety may spur constructive civic action, though do not suggest a consistent relationship between adverse environmental exposures and either climate stress/anxiety or civic engagement. Future research and action addressing the climate crisis should promote climate justice by ensuring mental health support for those who experience climate stress anxiety and by promoting pro-environmental civic engagement efforts.

Social Vulnerability and Biological Aging in New York City: An Electronic Health Records-Based Study.

Chronological age is not an accurate predictor of morbidity and mortality risk, as individuals' aging processes are diverse. Phenotypic age acceleration (PhenoAgeAccel) is a validated biological age measure incorporating chronological age and biomarkers from blood samples commonly used in clinical practice that can better reflect aging-related morbidity and mortality risk. The heterogeneity of age-related decline is not random, as environmental exposures can promote or impede healthy aging. Social Vulnerability Index (SVI) is a composite index accounting for different facets of the social, economic, and demographic environment grouped into four themes: socioeconomic status, household composition and disability, minority status and language, and housing and transportation. We aim to assess the concurrent and combined associations of the four SVI themes on PhenoAgeAccel and the differential effects on disadvantaged groups. We use electronic health records data from 31,913 patients from the Mount Sinai Health System (116,952 person-years) and calculate PhenoAge for years with available laboratory results (2011-2022). PhenoAge is calculated as a weighted linear combination of lab results and PhenoAgeAccel is the differential between PhenoAge and chronological age. A decile increase in the mixture of SVI dimensions was associated with an increase of 0.23 years (95% CI: 0.21, 0.25) in PhenoAgeAccel. The socioeconomic status dimension was the main driver of the association, accounting for 61% of the weight. Interaction models revealed a more substantial detrimental association for women and racial and ethnic minorities with differences in leading SVI themes. These findings suggest that neighborhood-level social vulnerability increases the biological age of its residents, increasing morbidity and mortality risks. Socioeconomic status has the larger detrimental role amongst the different facets of social environment.

The Exposure Peaks of Traffic-Related Ultrafine Particles Associated with Inflammatory Biomarkers and Blood Lipid Profiles.

In this article, we explored the effects of ultrafine particle (UFP) peak exposure on inflammatory biomarkers and blood lipids using two novel metrics-the intensity of peaks and the frequency of peaks. We used data previously collected by the Community Assessment of Freeway Exposure and Health project from participants in the Greater Boston Area. The UFP exposure data were time-activity-adjusted hourly average concentration, estimated using land use regression models based on mobile-monitored ambient concentrations. The outcome data included C-reactive protein, interleukin-6 (IL-6), tumor necrosis factor-alpha receptor 2 (TNF-RII), low-density lipoprotein (LDL), high-density lipoprotein (HDL), triglycerides and total cholesterol. For each health indicator, multivariate regression models were used to assess their associations with UFP peaks (N = 364-411). After adjusting for age, sex, body mass index, smoking status and education level, an increase in UFP peak exposure was significantly (p < 0.05) associated with an increase in TNF-RII and a decrease in HDL and triglycerides. Increases in UFP peaks were also significantly associated with increased IL-6 and decreased total cholesterol, while the same associations were not significant when annual average exposure was used. Our work suggests that analysis using peak exposure metrics could reveal more details about the effect of environmental exposures than the annual average metric.

Modeling health and well-being measures using ZIP code spatial neighborhood patterns.

Individual-level assessment of health and well-being permits analysis of community well-being and health risk evaluations across several dimensions of health. It also enables comparison and rankings of reported health and well-being for large geographical areas such as states, metropolitan areas, and counties. However, there is large variation in reported well-being within such large spatial units underscoring the importance of analyzing well-being at more granular levels, such as ZIP codes. In this paper, we address this problem by modeling well-being data to generate ZIP code tabulation area (ZCTA)-level rankings through spatially informed statistical modeling. We build regression models for individual-level overall well-being index and scores from five subscales (Physical, Financial, Social, Community, Purpose) using individual-level demographic characteristics as predictors while including a ZCTA-level spatial effect. The ZCTA neighborhood information is incorporated by using a graph Laplacian matrix; this enables estimation of the effect of a ZCTA on well-being using individual-level data from that ZCTA as well as by borrowing information from neighboring ZCTAs. We deploy our model on well-being data for the U.S. states of Massachusetts and Georgia. We find that our model can capture the effects of demographic features while also offering spatial effect estimates for all ZCTAs, including ones with no observations, under certain conditions. These spatial effect estimates provide community health and well-being rankings of ZCTAs, and our method can be deployed more generally to model other outcomes that are spatially dependent as well as data from other states or groups of states.

Polycystic ovary syndrome underdiagnosis patterns by individual-level and spatial social vulnerability measures.

OBJECTIVE: To use electronic health records (EHR) data at Boston Medical Center (BMC) to identify individual-level and spatial predictors of missed diagnosis, among those who meet diagnostic criteria for PCOS. METHODS: The BMC Clinical Data Warehouse was used to source patients who presented between October 1, 2003 and September 30, 2015 for any of the following: androgen blood tests, hirsutism, evaluation of menstrual regularity, pelvic ultrasound for any reason, or PCOS. Algorithm PCOS cases were identified as those with International Classification of disease (ICD) codes for irregular menstruation and either an ICD code for hirsutism, elevated testosterone lab, or polycystic ovarian morphology as identified using natural language processing on pelvic ultrasounds. Logistic regression models were used to estimate odds ratios (ORs) of missed PCOS diagnosis by age, race/ethnicity, education, primary language, body mass index (BMI), insurance type and social vulnerability index (SVI) score. RESULTS: In the 2003-2015 BMC-EHR PCOS at-risk cohort (n=23,786), there were 1,199 physician-diagnosed PCOS cases and 730 algorithm PCOS cases. In logistic regression models controlling for age, year, education, and SVI scores, Black/African American patients were more likely to have missed a PCOS diagnosis (OR = 1.69 [95% CI, 1.28, 2.24]) compared to non-Hispanic White patients, and relying on Medicaid or charity for insurance was associated with an increased odds of missed diagnosis when compared to private insurance (OR = 1.90 [95% CI, 1.47, 2.46], OR = 1.90 [95% CI, 1.41, 2.56], respectively). Higher SVI scores were associated with increased odds of missed diagnosis in univariate models. CONCLUSIONS: We observed individual-level and spatial disparities within the PCOS diagnosis. Further research should explore drivers of disparities for earlier intervention.

Quantifying urban park use in the USA at scale: empirical estimates of realised park usage using smartphone location data.

BACKGROUND: A large body of evidence connects access to greenspace with substantial benefits to physical and mental health. In urban settings where access to greenspace can be limited, park access and use have been associated with higher levels of physical activity, improved physical health, and lower levels of markers of mental distress. Despite the potential health benefits of urban parks, little is known about how park usage varies across locations (between or within cities) or over time. METHODS: We estimated park usage among urban residents (identified as residents of urban census tracts) in 498 US cities from 2019 to 2021 from aggregated and anonymised opted-in smartphone location history data. We used descriptive statistics to quantify differences in park usage over time, between cities, and across census tracts within cities, and used generalised linear models to estimate the associations between park usage and census tract level descriptors. FINDINGS: In spring (March 1 to May 31) 2019, 18·9% of urban residents visited a park at least once per week, with average use higher in northwest and southwest USA, and lowest in the southeast. Park usage varied substantially both within and between cities; was unequally distributed across census tract-level markers of race, ethnicity, income, and social vulnerability; and was only moderately correlated with established markers of census tract greenspace. In spring 2019, a doubling of walking time to parks was associated with a 10·1% (95% CI 5·6-14·3) lower average weekly park usage, adjusting for city and social vulnerability index. The median decline in park usage from spring 2019 to spring 2020 was 38·0% (IQR 28·4-46·5), coincident with the onset of physical distancing policies across much of the country. We estimated that the COVID-19-related decline in park usage was more pronounced for those living further from a park and those living in areas of higher social vulnerability. INTERPRETATION: These estimates provide novel insights into the patterns and correlates of park use and could enable new studies of the health benefits of urban greenspace. In addition, the availability of an empirical park usage metric that varies over time could be a useful tool for assessing the effectiveness of policies intended to increase such activities. FUNDING: Google.

Ambient air pollution and daily mortality in ten cities of India: a causal modelling study.

BACKGROUND: The evidence for acute effects of air pollution on mortality in India is scarce, despite the extreme concentrations of air pollution observed. This is the first multi-city study in India that examines the association between short-term exposure to PM2·5 and daily mortality using causal methods that highlight the importance of locally generated air pollution. METHODS: We applied a time-series analysis to ten cities in India between 2008 and 2019. We assessed city-wide daily PM2·5 concentrations using a novel hybrid nationwide spatiotemporal model and estimated city-specific effects of PM2·5 using a generalised additive Poisson regression model. City-specific results were then meta-analysed. We applied an instrumental variable causal approach (including planetary boundary layer height, wind speed, and atmospheric pressure) to evaluate the causal effect of locally generated air pollution on mortality. We obtained an integrated exposure-response curve through a multivariate meta-regression of the city-specific exposure-response curve and calculated the fraction of deaths attributable to air pollution concentrations exceeding the current WHO 24 h ambient PM2·5 guideline of 15 µg/m3. To explore the shape of the exposure-response curve at lower exposures, we further limited the analyses to days with concentrations lower than the current Indian standard (60 µg/m3). FINDINGS: We observed that a 10 µg/m3 increase in 2-day moving average of PM2·5 was associated with 1·4% (95% CI 0·7-2·2) higher daily mortality. In our causal instrumental variable analyses representing the effect of locally generated air pollution, we observed a stronger association with daily mortality (3·6% [2·1-5·0]) than our overall estimate. Our integrated exposure-response curve suggested steeper slopes at lower levels of exposure and an attenuation of the slope at high exposure levels. We observed two times higher risk of death per 10 µg/m3 increase when restricting our analyses to observations below the Indian air quality standard (2·7% [1·7-3·6]). Using the integrated exposure-response curve, we observed that 7·2% (4·2%-10·1%) of all daily deaths were attributed to PM2·5 concentrations higher than the WHO guidelines. INTERPRETATION: Short-term PM2·5 exposure was associated with a high risk of death in India, even at concentrations well below the current Indian PM2·5 standard. These associations were stronger for locally generated air pollutants quantified through causal modelling methods than conventional time-series analysis, further supporting a plausible causal link. FUNDING: Swedish Research Council for Sustainable Development.

Impact of heatwaves on all-cause mortality in India: A comprehensive multi-city study.

BACKGROUND: Heatwaves are expected to increase with climate change, posing a significant threat to population health. In India, with the world's largest population, heatwaves occur annually but have not been comprehensively studied. Accordingly, we evaluated the association between heatwaves and all-cause mortality and quantifying the attributable mortality fraction in India. METHODS: We obtained all-cause mortality counts for ten cities in India (2008-2019) and estimated daily mean temperatures from satellite data. Our main extreme heatwave was defined as two-consecutive days with an intensity above the 97th annual percentile. We estimated city-specific heatwave associations through generalised additive Poisson regression models, and meta-analysed the associations. We reported effects as the percentage change in daily mortality, with 95% confidence intervals (CI), comparing heatwave vs non-heatwave days. We further evaluated heatwaves using different percentiles (95th, 97th, 99th) for one, two, three and five-consecutive days. We also evaluated the influence of heatwave duration, intensity and timing in the summer season on heatwave mortality, and estimated the number of heatwave-related deaths. FINDINGS: Among âˆ¼ 3.6 million deaths, we observed that temperatures above 97th percentile for 2-consecutive days was associated with a 14.7 % (95 %CI, 10.3; 19.3) increase in daily mortality. Alternative heatwave definitions with higher percentiles and longer duration resulted in stronger relative risks. Furthermore, we observed stronger associations between heatwaves and mortality with higher heatwave intensity. We estimated that around 1116 deaths annually (95 %CI, 861; 1361) were attributed to heatwaves. Shorter and less intense definitions of heatwaves resulted in a higher estimated burden of heatwave-related deaths. CONCLUSIONS: We found strong evidence of heatwave impacts on daily mortality. Longer and more intense heatwaves were linked to an increased mortality risk, however, resulted in a lower burden of heatwave-related deaths. Both definitions and the burden associated with each heatwave definition should be incorporated into planning and decision-making processes for policymakers.