Opens in a new tab
Home/Example Projects/Environmental Monitoring Study for a Stormwater Retrofit Program

ILLUSTRATIVE EXAMPLE

Environmental Monitoring Study for a
Stormwater Retrofit Program

Monthly nutrient series at about 12 stream stations, analyzed with segmented models, seasonal terms, censored values below the detection limit and block bootstrap intervals, reported in journal format.

Education & Social SciencesEnvironmental & Sustainability StudiesRInterrupted time series

This example does not describe a real client or project. Every engagement is confidential, during the project and after it ends.

1

RESEARCH QUESTION

What the Study Set Out to Ask

This environmental monitoring study example follows an environmental research team that is asked to evaluate a stormwater retrofit program. Bioretention basins and permeable pavement have been added across a set of urban catchments, and the program's sponsors want to know whether nutrient concentrations in the streams below them changed after construction was finished, compared with what the earlier pattern and the untreated streams would suggest.

The monitoring data were collected for routine regulatory reporting, not for this question, and the streams respond to weather, season and land use as well as to the retrofits. The question is framed as a change relative to the projected pre-retrofit trend, read against the comparison streams, and not as proof that the retrofits caused it. The plan fixes the completion month, any phase-in period and the outcome before the post-retrofit months are examined.

2

STUDY DESIGN AND DATA

Monthly Series at About 12 Stations

Monthly samples from about 12 stations, six downstream of retrofitted catchments and six downstream of untreated catchments matched on size and land cover, for 36 months before and 36 months after the month construction is complete. Each record holds the station, the month, the total phosphorus concentration, a flag for values reported only as below the detection limit, and rainfall in the weeks before sampling. The laboratory changed its analytical method partway through the record, and with it the reporting limit, so the limit is not the same for every record.

  • Design: treated and comparison stations observed monthly before and after a known completion month
  • Outcome: monthly total phosphorus concentration, with values below the reporting limit kept as censored
  • Software requested: R, with scripts and a record of package versions
3

METHOD AND ASSUMPTIONS

Why a Segmented Model, and What Must Hold

Each station's series is modeled with a pre-retrofit trend, a change in level at completion and a change in slope after it, with seasonal terms and a rainfall covariate so that weather and the annual cycle are not mistaken for an effect. Concentrations are analyzed on the log scale because they are positive and right-skewed. Successive months are correlated, which makes conventional intervals too narrow, so intervals come from a moving-block bootstrap that resamples runs of consecutive months and keeps the correlation intact, with the block length chosen from the residual autocorrelation. Station-level estimates are combined in a random-effects summary, so the interval for the average change reflects how much the stations differ.

Values below the detection limit are left-censored, and the reporting limit differs before and after the laboratory's method change, so setting every such value to half its limit would distort both the trend and the variance, and could create a step at the date of the change. The primary model uses a censored-data likelihood that accepts a different limit for each record, and half-limit substitution is run as a sensitivity analysis. The design assumes that the comparison stations show what the treated stations would have done without the retrofits, which the plan examines in the pre-retrofit months. The completion month is also checked against the date of the laboratory change, so the two cannot be confused.

4

OUTPUT AND INTERPRETATION

How the Results Are Read

The output has a plot for each station with the fitted segments and the projected pre-retrofit trend, the changes in level and in slope for the treated and the comparison stations, and the primary estimate: the difference between the two groups' pooled changes, with its 95% confidence interval. That contrast is also expressed as a percent difference from the projected concentration at set times after completion, such as 12 and 24 months, so the size of the change can be judged and not only its direction. The share of values below the limit at each station is tabulated next to the estimates.

Sensitivity analyses vary the phase-in period, replace the censored-data model with half-limit substitution, add an autoregressive error model and drop any station near a change in land use. A change in the treated streams that the comparison streams do not share is described as consistent with an effect of the retrofits, with other events and the length of the record named as limits. The report gives the sponsors and the research team the size of the change and its uncertainty; any decision to extend the retrofits is theirs.

5

REPORTING

What the Researcher Received

A written report in journal format, with the design, the handling of values below the detection limit and the dependence model written as methods-section text a reviewer can follow choice by choice. The monitoring question, the environmental interpretation and the conclusions remain with the research team.

Delivered package

Interpreted written report (journal format), with tables and figures
Station plots with fitted segments and the projected pre-retrofit trend
Estimates with 95% confidence intervals and the share of values below the limit by station
Sensitivity analyses, with the detection-limit decision explained
Annotated R scripts and a record of package versions for reproducibility

NEXT STEP

Have a Similar
Research Need?

Share your project and receive a personalized assessment from a statistically trained Research Coordinator, with a fixed fee, usually the same day.

No Commitment Required

Strictly Confidential

Fixed Project Fee

In-House PhD-Level Statisticians

Same-Day Response

Serving Researchers Since 2014