Abstract
Lead (Pb) contaminated water remains a persistent problem in U.S. cities, disproportionately harming low-income households and communities of color. Although federal initiatives like Environmental Protection Agency’s (EPA) 2024 Lead and Copper Rule Improvements (LCRI) mandate the replacement of all lead service lines by 2037, equitable implementation of these policies hinges on how cities use discretion in prioritizing and sequencing replacement projects. In this study, we use data obtained through right-to-know requests to address the question: How and to what extent did Pittsburgh’s lead service line replacement (LSLR) program demonstrate equitable implementation? Using quantitative and spatial techniques, we show that Pittsburgh census tracts with greater percentages of low-income individuals, higher percentages of children with elevated blood lead levels (EBLL), and higher percentages of non-white people, were initially not prioritized relative to more affluent areas of the city. These early years of the program (2016–2018) reveal the consequences of equity as an afterthought: the most vulnerable areas of the city waited for replacements while whiter, wealthier areas were served first. Yet, as we observe, these patterns meaningfully shift post-2019 suggesting that advocacy through bodies like the Citizen Lead Response Advisory Committee (CLRAC) and the Pittsburgh Water Equity Task Force may redirect programs toward their stated equity goals. We end with a call for greater scrutiny and study of LSLR programs nationally to ensure that groups most historically affected by lead exposure are prioritized from the beginning of the process in every locale.
INTRODUCTION
In 2023, the Pittsburgh Water and Sewer Authority (PWSA) received the inaugural Environmental Justice and Equity Award from the Association of Metropolitan Water Agencies (AMWA) in recognition for its efforts to replace thousands of lead service lines. 1 For PWSA, this marked a high point in a long and challenging saga. Just six years earlier, routine water testing in 2016 revealed elevated lead (Pb) levels in Pittsburgh’s drinking water. In response to public pressure and a lawsuit filed by the Pennsylvania Attorney General alleging violation of the state’s Safe Drinking Water Act, 2 PWSA sought to restore confidence by making full lead service line replacement (LSLR) a central strategy for reducing lead exposure. A key component of this program was its stated focus on equity: prioritizing replacements in neighborhoods with the most vulnerable populations, and strengthening community engagement. 3 In 2019, under sustained pressure from local water justice advocates, PWSA established the Citizens Lead Response Advisory Committee (CLRAC) to incorporate community input and oversight into the implementation of its LSLR program. In the same year, PWSA also joined the Water Equity Task Force, along with members of local advocacy organizations and the University of Pittsburgh.
In this study, we analyze the translation of equity from rhetoric into practice, asking: How and to what extent did Pittsburgh’s LSLR program demonstrate equitable implementation? Using data recovered through right-to-know requests, quantitative analysis, and spatiotemporal mapping, we show that census tracts with greater percentages of low-income individuals, higher percentages of children with EBLL (≥5 µg/dL), and higher percentages of non-white people, were not initially prioritized in the LSLR process compared to more affluent census tracts.
BACKGROUND
With an estimated 9.2 million lead service lines still in use across the United States, lead-contaminated water threatens the well-being of millions of people. 4 Young children are particularly vulnerable to the impacts of lead exposure, which can cause harm to their developing bodies, especially affecting the brain. 5 In adults, lead has been linked to health problems such as cardiovascular disease 6 and psychiatric disorders. 7 Moreover, at the population level, some research links lead pollution to increased rates of crime (the “lead-crime” hypothesis) drawing connections between lead and both social and biological risk factors. 8
The ubiquitous nature of lead in our post-industrial environment means that individuals are often exposed via multiple routes including paint, 9 soil, 10 and water infrastructure. 11 However, as with most forms of environmental contamination, the risks of lead exposure are not uniformly distributed. Decades of empirical research and lived experience show that Black, Indigenous, People of Color (BIPOC) and low-income communities disproportionately suffer from lead exposure.12,13,14 Much of this is due to racialized patterns of residential segregation, which has resulted in marginalized groups being more likely to live in or near contaminated environments such as military bases, 15 Superfund sites, 16 and petrochemical infrastructure. 17 Moreover, when environmental contamination is detected, communities of color remain less likely to have their concerns addressed, due to unequal enforcement of environmental protections, a form of environmental racism. 18 Much critical environmental justice scholarship highlights the intersections between slow violence (incremental harm like toxic exposure, that unfolds gradually and often invisibly, over long periods of time), structural violence (harm produced by the way institutions and structures are organized and operate) and the role of the state in reinforcing and reproducing embedded social inequities. 19
The problem of lead exposure exemplifies all three of these dynamics because its health impacts may be misdiagnosed or undetected, resulting in an accumulation of harm over time. Lead hazards also reflect structural violence: substandard housing, 20 underfunded infrastructure, 21 and unequal regulatory enforcement 22 expose low-income and marginalized communities to disproportionate risk, not through singular decisions but through long-standing institutional arrangements. 23 The state plays a central role in this process, via historical policy choices such as discriminatory housing practices and uneven investment in public works, and through ongoing failures to remediate contaminated environments, thereby perpetuating and legitimizing these inequities over time. Together, these theories help to explain why and how the burden of lead exposure continues to disproportionately harm low-income people of color in the United States.
Regulating Pb
In 1991, the U.S. Environmental Protection Agency issued the Lead and Copper Rule (LCCR), 24 which regulates levels of Pb and copper in drinking water. Since its introduction, the rule has been revised and strengthened several times. As of this writing, the most current version—the Lead and Copper Rule Improvements (LCRI) released in October 2024—mandates that all drinking water systems across the U.S. must identify and replace lead service lines by 2037. 25 In 2021, after years of public outcry and media focus on the issue of lead-contaminated water in Flint, MI and other cities, the Biden Administration announced its Lead Pipe and Paint Action Plan (LPPAP), a federal plan allocating billions of dollars to the removal and replacement of lead services lines across the United States. 26 As the Biden administration elevated environmental justice across federal policy, equity and environmental justice became defining elements of the discourse on LSLR programs. Such framing marked a significant shift that positioned LSL removal not only as a technical and public health imperative, but as a central plank of the Biden administration’s (2020–2024) broader environmental justice agenda.
LSLR and equity
What constitutes equity in the context of lead service line replacement programs? According to the LSLR Collaborative, a national coalition of 28 organizations representing public health, water utilities, environmental groups, labor, consumer advocacy, housing, and various levels of government, “an equitable program will recognize that not everyone has the same societal and economic advantages, and provide support, not equally across the population, but rather as appropriate according to an individual’s circumstance.” 27 The LSLR Collaborative defines an equitable program as one that recognizes historic inequities, addresses income disparities, prioritizes tenants, and engages communities genuinely. 28
Although existing scholarship has focused heavily on the unequal distribution of lead exposure and its effects, relatively little attention has been paid to the question of equity within the context of LSLR. One exception is the work of Baehler et al. (2021), 29 which evaluates an LSLR program in Washington, D.C., and demonstrates that lower-income households and households of color were systematically less likely to obtain full service line replacements (the safer intervention) and more likely to receive unsafe partial replacements. In this study, we turn our focus to the sequencing and prioritization dimensions of equity in LSLR programs. Per the LSLR Collaborative, because each community has unique patterns of lead service line occurrence and social disparities, replacement planning offers a chance to prioritize replacements in ways that address environmental justice concerns (e.g., prioritizing households with multiple sources of lead exposure). Thus, in a context where widespread environmental remediation is required, questions of equity and environmental justice are closely intertwined with sequencing (i.e., which households in which geographical areas receive their service line replacement first).
The case of Pittsburgh, PA
In 2016, Pittsburgh residents learned that samples of the city’s water tested positive for elevated levels of lead, an exceedance of the EPA action level (15 ppb) at that time. 30 Subsequent investigation revealed that lead levels had been steadily rising for decades, a result of deteriorating water infrastructure, poor system maintenance, and controversial cost-cutting measures—including a change in water treatment chemicals—which increased water line corrosion. 31 As a result of this discovery, PWSA began the process of LSLR in 2016. In addition to the initial discovery of high lead levels in 2016, the water authority was sued by the state of Pennsylvania, and later entered into a $500k settlement agreement, for failing to notify and warn residents that LSLR work could (and did) cause temporary spikes in water lead levels. 32 Against this backdrop, PWSA (now known as Pittsburgh Water) took public strides to integrate equity into its general operations and LSLR program. In 2019, the water authority formed the Community Lead Response Advisory Committee (CLRAC). That same year, PWSA also convened a “Water Equity Learning Team” which joined the national Water Equity Taskforce, 33 a network of cities that collaborate to develop more equitable water policies and practices. In 2021, the Water Equity Learning Team (comprised of a number of community groups) published a report that outlined a water equity roadmap for Pittsburgh. The report discussed income and racial disparities and recommended that, “PWSA follow a replacement prioritization plan suggested by the Community Lead Response Advisory Committee. The plan ranks at-risk populations as a priority for LSLR.” 34 Since the beginning of the LSLR process in 2016, and as of this writing, PWSA has replaced over 14,000 lead service lines 35 within its water service area, as shown in Figure 1 and Table 1.

Map of PWSA Service Areas. PWSA, Pittsburgh Water and Sewer Authority.
Timeline of Major Events
RESEARCH DESIGN AND METHODOLOGY
For this study, we combine quantitative and spatial analysis to examine the extent to which PWSA prioritized equity considerations in its LSLR replacement program (2016–2024).
Quantitative analysis
This analysis draws on primary data covering the replacement of more than 12,000 LSLs undertaken by PWSA between 2016 and 2024, acquired via a right-to-know request filed with PWSA. These records included essential information such as the address of the replaced lines, the date of replacement, and whether the replacement occurred on the public or private side of the property. To assess equity, this LSLR data was integrated with census-tract-level demographic indicators from the U.S. Census Bureau’s 2024 American Community Survey (ACS). Key indicators of vulnerability included median household income 36 (to differentiate affluent from less-affluent census tracts), percentage of BIPOC (particularly relevant given the documented racial disparities in environmental exposure), and the percent of children in a census tract with Elevated Blood Lead Levels (EBLL), data made available by the Allegheny County Health Department. 37 This comprehensive data collection strategy aligns with the principles highlighted by the LSLR Collaborative, which emphasizes the use of demographic, economic, and public health data to identify communities disproportionately at risk and to ensure equitable program implementation.
Spatio-temporal mapping
Because environmental inequality is inherently spatial, 38 we employed spatio-temporal mapping techniques to analyze the data. 39 All analyses were conducted at the census tract level to balance geographic precision with statistical power. Replacement records were geocoded and aggregated to tracts, and we calculated the average replacement date (expressed as a fractional year) for each tract, as well as the total number of replacements. To account for tract size, we normalized replacement counts by total population, creating a “replacements per capita” measure. Our focus on the sequencing dimension of equity helps us understand environmental justice in remediation efforts, as emphasized by the LSLR Collaborative, which notes the priorities of agencies conducting replacements.
We began with multivariate ordinary least squares (OLS) regression to model the average replacement date as a function of median income, percent of BIPOC residents, and percent of children with an EBLL, while controlling for population. Recognizing that a simple average can obscure temporal dynamics, we also constructed a binary “early” indicator (tracts replaced before versus after the median date) and fitted a logistic regression to examine whether demographic characteristics predicted being in the early-replacement group.
To capture how the influence of each factor might have shifted over the course of the program, we implemented a time-varying logistic regression on a tract-month panel dataset. This model interacted each demographic variable with a smooth spline function of time. This allowed, for example, the effect of the percent BIPOC variable on the probability of a replacement occurring in a given month to change continuously. Moreover, this approach directly addresses the question of whether initial prioritization patterns later changed over time.
We complemented these regression-based approaches with distributional and nonparametric methods. Lorenz curves and concentration indices were constructed by ranking tracts by each demographic variable and plotting the cumulative share of replacements against the cumulative population share. The concentration index summarizes whether replacements were disproportionately concentrated among tracts with higher or lower values of a given characteristic. We computed these separately for the pre-2019 and 2019-onward periods (using January 2019 as the cutoff) to visualize how the distribution shifted after the committee’s formation and membership in the Water Equity Task Force.
Spearman’s rank correlation provided a simple, robust measure of the monotonic relationship between each demographic variable and the average replacement timing, complementing the regression results.
Finally, to assess whether spatial clustering might be influencing our inferences, we computed Moran’s I on the residuals of our regression models and on the raw average replacement timing. We also produced Local Indicators of Spatial Association (LISA) cluster maps to identify tracts where replacement timing was significantly higher or lower than neighboring tracts, after accounting for demographics. Together, these methods allow us to examine equity not only as a final outcome, but as a dynamic process, revealing whether and how prioritization changed over the life of the program.
RESULTS AND FINDINGS
We restrict our presentation to the most salient patterns that emerged from the analyses, focusing on evidence of inequitable sequencing and its evolution over time. Full model outputs are available upon request and result summaries are provided in the Appendix.
Aggregate patterns
Over the study period, all 85 census tracts within the PWSA service area received at least one LSLR. The average replacement date (expressed as a fractional year) ranged from 2016.1 to 2022.6. 40 Spearman rank correlations (Table 2) show that, over the entire program, tracts with higher median income tended to be replaced earlier (ρ = –0.27, p = 0.01), although tracts with higher percentages of BIPOC residents tended to be replaced later (ρ = 0.26, p = 0.02). 41
Spearman Rank Correlations with the Timing_Avg Variable
To understand how the influence of race shifted over time, we estimated a time-varying logistic regression on a tract-month panel (8,667 observations). At the beginning of the program (early 2016), tracts with higher percent BIPOC were significantly less likely to receive a replacement (β = –0.31, p = 0.007). A significant positive interaction with the second spline term (β = 0.39, p = 0.002) indicates that this disadvantage diminished over time and eventually reversed. By late 2022, higher percent BIPOC tracts had become more likely to be replaced than whiter tracts. In contrast, the positive effect of EBLL on replacement probability remained constant throughout (β = 0.35, p < 0.001), with no evidence of change over time.
In the Lorenz curves in Figure 2, note that the closer a Lorenz curve lies to the diagonal, the more proportional the distribution. The Robin Hood (RH) (Hoover) index indicates the share of replacements needing reallocation to achieve equality. Before 2019, %BIPOC showed a modest bias toward whiter neighborhoods (CI = –0.16, RH = 0.02), requiring just 2% of replacements to be redistributed. After 2019, the distribution became nearly neutral (CI = –0.02, RH = 0.07). Income was a far stronger predictor early on as the curve sat well below the diagonal (CI = 0.13, RH = 0.25), meaning 25% of replacements would need to move from higher to lower-income tracts. This inequality was almost entirely corrected after 2019 (RH = 0.04). For child lead risk (EBLL), both periods showed positive concentration (CI = 0.08 pre, 0.39 post), with the Hoover index rising from 0.20 to 0.31 showing a deliberate and growing prioritization of high-risk areas after the CLRAC. Thus, class-based inequity was initially larger than racial inequity, yet both came close to resolution after 2019. Following Tsai’s caution 42 we note that the small number of census tracts (pre-CLRAC n = 77) may bias the RH and CI downward, leading to an overly optimistic impression of early racial equity. The more robust, statistically significant time-varying logistic model therefore serves as our primary evidence for early racial disadvantage, although the inequality metrics should be treated as suggestive.

Lorenz curves for % BIPOC, Median Income, and % EBLL, pre and post CLRAC. Curves below the diagonal indicate greater share to tracts with higher characteristic values; CI and RH quantify inequality.
The time-varying coefficient trajectories (Fig. 3) mirror these indices. The coefficient for percent BIPOC began negative (β = –0.31, p = 0.007) and, driven by a significant positive interaction (β = 0.39, p = 0.002), steadily increased, crossing zero around 2022. In contrast, the coefficient for EBLL remained positive and flat throughout (β = 0.35, p < 0.001), with no significant change over time. Together, these complementary analyses show that early racial inequity was corrected only after sustained community oversight, although health risk was consistently prioritized. 43

Results of Time-Varying Logistic Regression (Spline vs. Demographics).
Accounting for replacement density
When we controlled for the number of replacements per capita in an OLS model of average timing, EBLL emerged as a significant predictor of later replacement (β = 0.39, p = 0.04), even after holding replacement intensity constant. This suggests that high-EBLL tracts were not simply receiving fewer replacements; the replacements they did receive occurred later in the program. Race and income were not significant in this model, indicating that their aggregate effects may be mediated by replacement density (i.e., whiter, higher-income census tracts had more replacements overall, driving their earlier average timing).
Spatial clustering
A LISA analysis identified isolated clusters of tracts with significantly early or late replacement timing (Fig. 4). The low-low (early) cluster includes tracts in neighborhoods like Stanton Heights, Perry North, and Upper Lawrenceville, all with median incomes above $80,000. The high-high (late) cluster includes several Bloomfield tracts and one Shadyside tract, with median incomes between $47,000 and $68,000. No significant global spatial autocorrelation was detected (Moran’s I = 0.048, p = 0.199); therefore, these local patterns should be interpreted with caution. Our primary conclusions are drawn from the time-varying logistic regression and concentration indices, which do not rely on spatial assumptions.

LISA Cluster map of average replacement timing by census tract.
DISCUSSION AND CONCLUSION
In November 2023, PWSA accepted an award for its commitment to equitable water policy. 44 As the analyses presented here reveal, however, the agency originally fell short of prioritizing the most vulnerable areas for the earliest period of intervention. Yet with community input through the CLRAC and the Water Equity Task Force, PWSA was able to correct course and modify its implementation approach.
Limitations
We acknowledge several data limitations. The EBLL data, obtained from the Allegheny County Health Department, reflect only children who were tested and may not capture the full population at risk; moreover, sample sizes at the census tract level can be small, necessitating caution in interpreting tract-level percentages. 45 American Community Survey 2024 (ACS) estimates are subject to sampling variability, particularly for smaller tracts, which may introduce measurement error in our demographic predictors. Although the OLS model’s overall fit is marginally non-significant, the individual coefficient for percentage of BIPOC residents remains statistically significant (p = 0.046), as individual variables can reveal valid relationships even in social science models where low r2 values are common due to unpredictable institutional behaviors and unobserved local factors in cross-sectional data. 46 We address this limitation by corroborating our findings across multiple lines of evidence, including Spearman rank correlations and time-varying logistic regression, both of which independently confirm our findings. Additionally, although replacement records were obtained through a formal right-to-know request and are presumed to be complete for PWSA’s reported work, we cannot rule out the possibility of missing or mis-coded entries. As PWSA was unable to produce records about the distribution of LSLs, we chose to normalize our data by population rather than by existing LSLs, which may be unevenly distributed within the city. Despite these limitations, the consistency of findings across multiple analytical approaches—time-varying logistic regression, concentration indices, spatial clustering, and non-parametric correlations—strengthens confidence in the overall patterns we report.
Conclusion
This study provides empirical evidence that PWSA’s early LSLR implementation (2016–2018) did not prioritize equity. Instead, during this period replacements were disproportionately concentrated in tracts with higher percentages of white residents and higher average median incomes. The time-varying logistic model demonstrates that the disadvantage for BIPOC tracts was statistically significant at the program’s outset and only diminished over time. Concentration indices confirm this pattern: before the CLRAC’s formation and PWSA’s membership in the Water Equity Task Force in 2019, replacements were concentrated in census tract areas with higher percentages of white residents (CI = –0.14); afterward, the distribution became racially neutral. Income showed a similar but weaker trajectory, shifting from favoring higher-income tracts pre-2019 to favoring lower-income tracts 2019-onward.
In Pittsburgh’s case, early inequity was only revealed through independent analysis and was not publicly documented by the utility itself, highlighting the need for transparency and accountability mechanisms. In the case of the ongoing clean water transition, community engagement is of paramount importance as repairs cannot be completed without the public’s cooperation and buy-in. This emphasis on engagement aligns with a robust body of environmental justice literature that identifies meaningful public participation as a central tenet, with scholars underscoring that equitable outcomes depend on inclusive decision-making processes that empower affected communities.47,48,49
We do not believe Pittsburgh’s LSLR program is an outlier. Rather, it mirrors challenges observed in other cities. Washington D.C.’s cost-sharing model, which linked lead line replacement access to household wealth, illustrates how neoliberal policy design can create disparities. However, as our time-varying analysis shows, the eventual correction of racial disparities from 2019 onward suggests that such failures can be addressed with sustained advocacy and oversight. Federal funding mechanisms anticipating partial homeowner contributions similarly risk cementing racial and class inequities, as seen in both cities. With changing federal priorities, such funding mechanisms themselves are in flux.
Some key takeaways for other municipalities considering LSLR include: First, equity should be embedded in prioritization criteria from a program’s inception, rather than applied later as a corrective measure. Second, transparency in prioritization decisions is essential. Public dashboards showing how sequencing aligns with stated equity goals can build trust and enable accountability. Third, equity and public health objectives need not conflict. Despite early disparities in race and income, PWSA’s eventual correction of the imbalances led to sustained prioritization of high-EBLL areas, thereby demonstrating that improving equity does not come at the expense of protecting children from lead exposure, rather, these goals can be pursued in tandem.
Looking ahead, many unanswered questions remain with regard to how equity has or has not shaped decision-making in the thousands of cities now undertaking LSLRs. To date, no comparative, systematic studies have examined equity outcomes across different cities or regions. One major barrier is the absence of standardized, publicly available data. Future research could critically examine this data gap itself, asking what it reveals about transparency, accountability, and the political will to address environmental justice concerns at scale. Another direction could be efforts directed toward the development of an equity-weighted prioritization tool to assist water systems in the task of strategically sequencing repairs in ways that deliver the greatest benefit to historically marginalized and high-risk communities. Embedding community co-design into the creation of weighting criteria and decision objectives will be essential to ensure that equity-focused algorithmic tools reflect lived experiences, enhance procedural justice, and support more accountable infrastructure planning.
AUTHORS’ CONTRIBUTIONS
Conceptualization: C.T. and A.V. Methodology: A.V. Formal analysis: A.V. Investigation: C.T. and A.V. Data curation: C.T. and A.V. Writing original draft: C.T. and A.V. Writing—review and editing: C.T. and A.V. Supervision: C.T. and A.V. Project Administration: C.T. and A.V.
Footnotes
AUTHOR DISCLOSURE STATEMENT
No competing financial interests exist.
FUNDING INFORMATION
No funding was received for this article.
Appendix
Data Dictionary
| Variable | Description | Data type | Notes/Source |
|---|---|---|---|
| a. master_demographics_final.csv | |||
| Tract_Number | Census tract number (without state/county prefix) | String | Identifier; may require conversion to full GEOID |
| Median_Income | Median household income in the census tract | Numeric (float) | Source: ACS 2024 5-year estimates |
| Total_Population | Total population of the census tract | Numeric (float) | Source: ACS 2024 5-year estimates |
| Pct_BIPOC | Percentage of population that is Black, Indigenous, or People of Color (i.e., non-White) | Numeric (float) | Computed from ACS race/ethnicity categories |
| Pct_EBLL | Percentage of children under 6 years tested with elevated blood lead levels (≥ 3.5 μg/dL) | Numeric (float) | Source: Allegheny County Health Department; based on tested children only |
| GEOID_11 | Full 11-digit census tract FIPS code (42 + county + tract) | String | Used for merging with replacement records |
| GEOID | Cleaned version of GEOID_11 (string, 11 digits, no decimals) | String | Derived field for consistent merging |
| b. pipe_replacements_with_tract.csv | |||
| Address | Street address of the property where replacement occurred | String | As provided by PWSA |
| Site Verification Public | Material type of the public side of the service line before replacement | String | Lead, Copper, Galvanized Iron, etc. |
| Site Verification Private | Material type of the private side of the service line before replacement | String | Lead, Copper, etc. |
| Replaced Material Public | Material used to replace the public side | String | PEX, Copper, etc. |
| Replaced Material Private | Material used to replace the private side | String | PEX, Copper, etc. |
| Replacement Date Public | Date of public-side replacement | Date (string) | May be missing if not replaced |
| Replacement Date Private | Date of private-side replacement | Date (string) | May be missing if not replaced |
| Latitude | Latitude coordinate of the address | Numeric (float) | Geocoded from address |
| Longitude | Longitude coordinate of the address | Numeric (float) | Geocoded from address |
| Neighborhood | Neighborhood name from the original geocoding | String | May be missing or inconsistent |
| GEOID | Full 11-digit census tract FIPS code assigned to the address | String | Assigned via spatial join to census tract boundaries |
| Date | Cleaned replacement date (using public date if available, else private) | Date (datetime) | Derived field |
| Timing_val | Fractional year representing replacement date (e.g., 2016.5 = July 2016) | Numeric (float) | Computed as year + (month-1)/12 |
