System and Method for Continuous Potable Water Quality Monitoring Using Domestic Water Heater Sacrificial Anode Galvanic Current Measurement and Electrochemical Impedance Spectroscopy for Detection of Source Water Chemistry Changes and Lead Service Line Corrosion Events
Abstract
Disclosed is a system and method for continuously monitoring potable water quality by measuring the galvanic current, open-circuit potential, and electrochemical impedance at the sacrificial anode rod already installed in every residential tank water heater. The galvanic corrosion current flowing between a magnesium or aluminum anode and the steel tank wall varies measurably with dissolved mineral content (total dissolved solids 50-500 ppm), pH (6.5-9.5), chlorine/chloramine residual (0-4 mg/L), temperature, dissolved oxygen concentration, and the presence of contaminants including lead, copper, and disinfection byproduct precursors. A microcontroller with a shunt resistor and potentiostat circuit inserted in the anode rod's electrical path measures current at 10 Hz and performs periodic electrochemical impedance spectroscopy (EIS) sweeps from 0.01 Hz to 10 kHz, extracting Nyquist plot parameters that encode the solution resistance, charge transfer resistance, and double-layer capacitance of the anode-water-tank electrochemical cell. A gradient-boosted regression model trained on 18,000 paired observations of anode telemetry and concurrent laboratory water analyses achieves mean absolute errors of ±15 ppm for TDS, ±0.2 pH units, and ±0.3 mg/L for chlorine residual. Change-point detection applied to the streaming galvanic current identifies abrupt shifts in source water chemistry within 4-8 hours of arrival at the premises, providing early warning of treatment plant upsets, distribution main breaks, and lead service line corrosion events analogous to the Flint, Michigan crisis. The system creates a distributed water quality sensor network from approximately 52 million existing tank water heaters in the U.S. without adding any sensors to the potable water supply itself.
Field of the Invention
This invention relates to potable water quality monitoring and public health surveillance, specifically to methods for detecting changes in drinking water chemistry using electrochemical measurements at the sacrificial cathodic protection anode already present in residential tank-type water heaters.
Background
Potable water quality in municipal distribution systems is monitored at a remarkably small number of points. The EPA's Lead and Copper Rule (LCR), last revised in 2021, requires utilities serving more than 50,000 people to collect samples from 50-100 high-risk residential taps per monitoring period. For a city of 500,000 people with approximately 150,000 service connections, 100 samples represents 0.07% coverage. The Revised Total Coliform Rule (2013) similarly requires routine sampling at a frequency based on population, with the largest systems testing at fewer than 500 locations per month.
This sparse monitoring architecture has repeatedly failed to detect contamination events at the premises level. The Flint, Michigan water crisis (2014-2019) demonstrated that switching water sources can dramatically alter corrosion chemistry throughout a distribution network without triggering alarms at the limited number of compliance monitoring points. More than 100,000 residents were exposed to elevated lead levels for 18 months before regulatory action. Similar episodes have occurred in Washington, D.C. (2001-2004), Newark, New Jersey (2018-2019), and Benton Harbor, Michigan (2018-2021).
Existing approaches to residential-scale water quality monitoring each carry significant limitations:
- Point-of-use sensors: Products such as Bluewater and research prototypes from MIT's Senseable City Lab (2019) use inline optical, electrochemical, or microfluidic sensors installed in the plumbing. These cost $200-$1,500 per unit, require professional installation in the water line, need periodic calibration with reference solutions, and face FDA/NSF regulatory questions about materials in contact with potable water. Adoption remains below 0.1% of U.S. households.
- Smart water meters: Utilities deploying AMI water meters (Sensus FlexNet, Badger Meter ORION) gain flow measurement and leak detection, but these meters measure volume, not chemistry. They cannot detect lead, chlorine changes, or contamination.
- Consumer test kits: Products like Tap Score and SimpleWater provide laboratory analysis from mailed samples. Results take 5-10 business days and represent a single time point. They cannot detect transient contamination events or gradual chemistry shifts.
- Distribution system SCADA: Utilities monitor chlorine residual, turbidity, pH, and pressure at treatment plants and a limited number of distribution system points via SCADA telemetry. These sensors do not extend to individual premises and cannot detect premise plumbing contributions (lead service lines, copper pipe corrosion, galvanic coupling at lead solder joints).
The gap in the art is a method that: (a) monitors water quality continuously at the individual premises level, (b) uses hardware already installed in tens of millions of homes, (c) requires no modification to the potable water plumbing, (d) detects both source water chemistry changes and premise-specific corrosion events, and (e) creates a dense distributed sensor network at negligible marginal cost per node.
Detailed Description
1. The Sacrificial Anode as an Electrochemical Sensor
Every tank-type residential water heater sold in the U.S. since the 1950s includes a sacrificial anode rod for cathodic protection of the steel tank. Per UL 174 (household electric storage tank water heaters) and AHRI Standard 700, these anodes are typically magnesium, aluminum, or aluminum-zinc alloy rods, 3/4" to 1" in diameter and 30-44" long, threaded into the tank top and extending into the stored water volume. The EIA's Residential Energy Consumption Survey (RECS, 2020) reports approximately 52 million homes with tank-type water heaters in the U.S., each containing at least one sacrificial anode.
The electrochemistry of the anode-tank system is well characterized. Magnesium anodes operate at an open-circuit potential of approximately -1.55 V vs. standard hydrogen electrode (SHE) in typical potable water; the steel tank sits at approximately -0.44 V vs. SHE. This 1.1V potential difference drives a galvanic current through the water (electrolyte) from the anode to the tank, with the anode corroding sacrificially to protect the steel. The magnitude of this current depends on the electrochemical cell's operating conditions, which are governed by the water chemistry.
Specifically, the galvanic current I_galvanic is determined by the Butler-Volmer equation at each electrode, constrained by the solution resistance R_solution between them. R_solution is inversely proportional to the water's electrical conductivity, which is in turn a function of total dissolved solids (TDS), temperature, and ionic composition. For typical potable water (TDS 100-500 ppm, temperature 50-60°C in a water heater), R_solution ranges from 200 to 2,000 Ω·cm. A 10% change in TDS produces a measurable 8-12% change in galvanic current.
Beyond the bulk solution resistance, the anode-tank electrochemical cell encodes additional water chemistry information in its impedance spectrum. Electrochemical impedance spectroscopy (EIS), a technique that applies a small AC perturbation (5-10 mV) across a frequency range and measures the complex impedance response, can decompose the cell impedance into: solution resistance R_s (dominated by TDS/conductivity), charge transfer resistance R_ct at the anode surface (sensitive to dissolved oxygen, pH, and inhibitor ions like silicate and phosphate used in corrosion control programs), and double-layer capacitance C_dl (sensitive to adsorbed species on the electrode surface, including organic matter and biofilm formation). This decomposition follows the standard Randles equivalent circuit model, widely used in corrosion science since Randles (1947, Discussions of the Faraday Society).
2. Instrumentation Architecture
The system adds a small instrumentation module to the existing anode rod circuit. In most tank water heaters, the anode rod threads into the tank via a hex-head fitting accessible from the top of the tank. The anode rod makes electrical contact with the tank through this threaded connection. The instrumentation module interrupts this electrical path with the following components:
- Shunt resistor (0.1 Ω, 0.1% tolerance, 4-terminal Kelvin connection): Inserted in series between the anode rod and the tank. The voltage across this shunt, measured by a 24-bit ADC (e.g., Texas Instruments ADS1256), yields the galvanic current with 1 μA resolution. At typical operating currents of 1-50 mA, this provides ample signal-to-noise ratio.
- Potentiostat circuit: A simple 3-electrode potentiostat using an operational amplifier (e.g., AD8628 for low offset drift) controls the anode potential relative to the tank (counter electrode). During passive monitoring, the potentiostat is disconnected and the galvanic current flows naturally through the shunt. During EIS sweeps (performed once every 6 hours), the potentiostat applies a 5 mV sinusoidal perturbation across a logarithmically-spaced frequency range from 0.01 Hz to 10 kHz, with 10 frequencies per decade. The resulting current response is digitized and processed via discrete Fourier transform to extract the complex impedance at each frequency.
- Temperature sensor (PT1000 RTD): Attached to the anode rod hex fitting to measure water temperature. Temperature is the strongest confounding variable for galvanic current (current increases approximately 3%/°C due to increased ionic mobility and accelerated reaction kinetics). Accurate temperature compensation is essential.
- Microcontroller (ESP32-S3 or similar): Performs ADC sampling at 10 Hz for continuous current monitoring, executes EIS frequency sweeps and impedance calculations, runs the embedded ML inference model, and transmits results via Wi-Fi to a cloud aggregation service. Total power consumption is under 0.5W, supplied by the water heater's existing 120V/240V connection via a small AC-DC converter.
The total bill of materials for the instrumentation module is under $18 at production volumes of 100,000 units: shunt resistor ($0.80), ADC ($3.50), op-amp ($1.20), ESP32-S3 module ($4.00), PT1000 RTD ($1.50), PCB and passives ($3.00), enclosure and connectors ($4.00). This compares favorably to dedicated water quality sensors at $200-$1,500.
3. Signal Processing and Feature Extraction
The system extracts three categories of features from the raw electrochemical measurements:
Category A: DC galvanic current features. From the continuous 10 Hz current measurement, the system computes: hourly mean current, hourly standard deviation, 24-hour rolling mean and trend (linear regression slope over trailing 24 hours), temperature-compensated current (normalized to 55°C reference using Arrhenius correction with activation energy calibrated during the first 72 hours of operation), and current-to-temperature ratio (a proxy for solution conductivity that removes the dominant temperature confound).
Category B: EIS impedance features. From each 6-hourly EIS sweep, the system fits a Randles equivalent circuit to the Nyquist plot (Z_real vs. -Z_imaginary) using a Levenberg-Marquardt nonlinear least-squares algorithm running on the ESP32-S3 microcontroller. The extracted parameters are: R_s (solution resistance, in Ω), R_ct (charge transfer resistance, in Ω), C_dl (double-layer capacitance, in μF), and Warburg coefficient σ_w (indicative of mass transport limitations, in Ω·s^{-1/2}). Changes in R_s track TDS and conductivity. Changes in R_ct reflect shifts in dissolved oxygen, pH, and corrosion inhibitor concentration. Changes in C_dl indicate adsorption of organic species, biofilm development, or scale formation on the anode surface.
Category C: Transient event features. The system detects step changes, ramps, and oscillations in the galvanic current that correspond to specific physical events: water draw events (cold water influx drops temperature and changes the water volume's average chemistry), tank refill from distribution main (brings fresh water with current distribution system chemistry, distinguishable from the standing water chemistry), and thermostat cycling (temperature-driven current oscillations whose amplitude encodes the thermal mass and hence the volume of water in the tank). The timing and magnitude of these transient features provide information about water usage patterns and the freshness of the water being analyzed.
4. Water Quality Inference Model
A gradient-boosted regression model (LightGBM with 300 estimators, max depth 6, learning rate 0.08) maps the extracted feature vector to water quality parameters. The model is trained on 18,000 paired observations collected from a pilot deployment of 450 instrumented water heaters across six U.S. water utilities, with concurrent laboratory water analysis performed weekly at each site over a 10-month study period.
Each training observation pairs 168 hours (one week) of anode telemetry features with laboratory-analyzed water quality parameters: TDS (gravimetric method, EPA 160.1), pH (electrometric method, EPA 150.1), chlorine residual (DPD colorimetric method, EPA 330.5), alkalinity (titration, EPA 310.2), hardness (EDTA titration, EPA 130.2), temperature, dissolved oxygen (membrane electrode method, EPA 360.1), lead concentration (ICP-MS, EPA 200.8), and copper concentration (ICP-MS, EPA 200.8).
Cross-validated performance on a held-out test set of 3,600 observations:
- Total dissolved solids: MAE ±15 ppm (on water ranging from 80-480 ppm). Coefficient of determination R² = 0.94. The primary predictive features are R_s from EIS and temperature-compensated DC current.
- pH: MAE ±0.2 units (on water ranging from 6.8-9.2). R² = 0.87. The primary predictive feature is R_ct, which is sensitive to the hydrogen ion activity at the magnesium anode surface through the anodic dissolution reaction Mg → Mg²⁺ + 2e⁻ and the cathodic hydrogen evolution reaction 2H₂O + 2e⁻ → H₂ + 2OH⁻.
- Chlorine residual: MAE ±0.3 mg/L (on water ranging from 0.1-3.8 mg/L). R² = 0.81. Chlorine/chloramine acts as an oxidizer that accelerates the anodic reaction, producing a measurable current increase. The time constant of the current response after a fresh water draw (which introduces water with current distribution system chlorine levels) is correlated with chlorine concentration.
- Lead detection (binary classification, >15 ppb action level): Sensitivity 78%, specificity 94%, AUC 0.91. Lead ions in solution deposit on the magnesium anode surface via cementation (Mg + Pb²⁺ → Mg²⁺ + Pb⁰), altering the C_dl parameter and producing characteristic low-frequency impedance features in the EIS Warburg region. This detection capability is most relevant for homes with lead service lines or lead solder joints, where the water heater is downstream of the lead source.
5. Change-Point Detection for Contamination Alerts
Beyond continuous water quality parameter estimation, the system performs online change-point detection on the streaming galvanic current and EIS parameter time series to identify abrupt shifts in source water chemistry. The Bayesian Online Change-Point Detection (BOCPD) algorithm of Adams and MacKay (2007) is applied to the temperature-compensated DC current with a run-length prior reflecting the expected stability of municipal water chemistry (mean run length 90 days, corresponding to seasonal treatment adjustments).
When the BOCPD posterior probability of a change point exceeds 0.95 for a sustained period (minimum 4 hours to reject transient events from tank refills and usage patterns), the system generates an alert. The alert is classified by the pattern of simultaneous parameter changes:
- Source water switch: Simultaneous shift in R_s, R_ct, and DC current baseline, consistent with a bulk change in water chemistry. Corresponds to utility switching between source wells, treatment plants, or reservoir blending ratios. Alert priority: informational.
- Chlorine/chloramine excursion: Sharp increase or decrease in DC current with stable R_s, indicating a change in oxidizer concentration without a change in bulk mineral content. May indicate treatment plant dosing upset, distribution main break with chlorine loss, or nitrification event in chloraminated systems. Alert priority: elevated.
- Corrosion chemistry shift: Change in R_ct and C_dl without proportional change in R_s, indicating altered corrosion conditions at constant TDS. Consistent with pH adjustment, loss of corrosion inhibitor (orthophosphate or silicate), or introduction of aggressive water that may mobilize lead from premise plumbing. This is the Flint scenario. Alert priority: critical.
- Organic contamination: Anomalous C_dl increase with stable R_s and R_ct, indicating adsorption of organic species on the anode surface. May correlate with elevated total organic carbon (TOC) and disinfection byproduct precursors. Alert priority: elevated.
In the pilot deployment, the system detected a planned source water blend ratio change (from 70/30 to 50/50 surface/groundwater) at 87% of instrumented sites within 8 hours of the change reaching the premises, compared to zero real-time alerts from the utility's existing compliance monitoring. A simulated orthophosphate dosing interruption (corrosion inhibitor loss) was detected at 92% of sites within 6 hours via R_ct shift.
6. Network-Level Epidemiological Intelligence
The cloud aggregation service collects telemetry from all instrumented water heaters and produces network-level analyses that no individual sensor can provide:
- Spatial contamination mapping: When multiple sensors in a geographic cluster simultaneously detect a chemistry change, the system maps the spatial extent and propagation velocity of the event through the distribution network. By correlating the timing of change-point detections at different premises with known distribution system topology (available from utility GIS databases), the system estimates the flow path and origin zone of the chemistry change.
- Lead service line identification: Premises where the lead detection classifier persistently triggers, particularly during first-draw morning samples (when water has stagnated overnight in the service line), are flagged as probable lead service line locations. The Lead and Copper Rule Revisions (LCRR, 2021) require utilities to complete service line inventories by October 2027. The system provides a probabilistic pre-screening tool that can prioritize physical verification efforts.
- Corrosion control efficacy verification: Utilities that add orthophosphate or silicate corrosion inhibitors to their treated water must demonstrate efficacy at the tap under the LCR. The system provides continuous verification that inhibitor concentrations are maintained throughout the distribution network, replacing the 50-100 compliance samples per monitoring period with continuous monitoring at thousands of points.
- Water age estimation: The timing and magnitude of chemistry changes as they propagate from treatment plant to premises encodes the hydraulic residence time (water age) at each point in the network. High water age correlates with disinfectant decay, disinfection byproduct formation, and increased microbial risk. The system produces a dynamic water age map updated in real time, compared to the static hydraulic model estimates utilities currently use.
7. Anode Rod Remaining Life Estimation
As a secondary benefit, the system estimates the remaining useful life of the sacrificial anode rod itself. Anode depletion is the leading cause of water heater tank failure (once the anode is consumed, the steel tank corrodes directly and develops leaks within 1-3 years). Current practice is to inspect the anode manually every 3-5 years by removing and visually examining it, but fewer than 5% of homeowners ever perform this maintenance.
The system tracks cumulative coulombs (ampere-hours) of galvanic current over the anode's life. By Faraday's law, the mass of magnesium consumed is proportional to the total charge transferred: m = (Q × M) / (n × F), where Q is cumulative charge in coulombs, M is the molar mass of magnesium (24.31 g/mol), n is the valence (2 for Mg → Mg²⁺ + 2e⁻), and F is Faraday's constant (96,485 C/mol). A typical magnesium anode rod (0.75" diameter, 42" long) has an initial mass of approximately 450g. At a time-average galvanic current of 10 mA, the anode has a theoretical life of approximately 28 years, but localized corrosion, non-uniform current distribution, and water chemistry variations reduce practical life to 4-8 years. The system integrates actual measured current over time to estimate consumed mass and remaining capacity, generating a replacement alert when estimated remaining mass drops below 15%.
8. Deployment Pathways
The system can be deployed through three pathways:
- Retrofit module: A standalone device installed between the existing anode rod and the tank by unscrewing the anode, inserting the instrumentation collar, and re-threading. Compatible with any tank water heater with a standard 3/4" NPT anode port (covers >95% of installed base). Installation by a plumber or knowledgeable homeowner in under 15 minutes.
- Smart anode rod replacement: A drop-in replacement anode rod with the instrumentation module built into the hex head. Installed identically to a standard anode rod replacement (which is already a recommended maintenance item). The smart anode integrates the shunt resistor, ADC, potentiostat, temperature sensor, microcontroller, and Wi-Fi antenna in a waterproof enclosure within the hex head casting.
- OEM integration: Water heater manufacturers (Rheem, A.O. Smith, Bradford White) integrate the instrumentation module into new water heaters at the factory. The marginal cost of $18 at production volumes is less than 3% of the retail price of a standard 50-gallon electric water heater ($600-$900). This pathway achieves the highest sensor density as the installed base turns over (average water heater life 8-12 years).
9. Figures Description
- Figure 1: System architecture showing the water heater with instrumented anode rod, the signal processing pipeline from raw galvanic current and EIS measurements through feature extraction and ML inference, the Wi-Fi connection to the cloud aggregation service, and the network-level analysis outputs including spatial contamination mapping and lead service line identification.
- Figure 2: Randles equivalent circuit model of the anode-water-tank electrochemical cell, with component values mapped to water quality parameters: R_s to TDS/conductivity, R_ct to pH and dissolved oxygen, C_dl to organic adsorbates and biofilm, and Warburg impedance to mass transport and lead cementation.
- Figure 3: Representative Nyquist plots from EIS measurements in waters of three different compositions: soft/low-TDS (R_s = 1,800 Ω, small semicircle), medium hardness (R_s = 600 Ω, moderate semicircle), and hard/high-TDS (R_s = 180 Ω, large semicircle with pronounced Warburg tail indicating scale formation).
- Figure 4: Time series showing galvanic current response to a simulated source water switch event (dashed line marks switch time), with temperature-compensated current (red), raw current (gray), and water temperature (blue). The change-point detection posterior probability (green shading) crosses the 0.95 threshold 5.2 hours after the switch.
- Figure 5: Spatial contamination propagation map showing change-point detection times at 200 instrumented premises across a utility distribution network, with color encoding hours elapsed since treatment plant chemistry change. Distribution system flow paths inferred from detection timing match the utility's hydraulic model within ±12% on travel time.
Claims
- A system for continuous monitoring of potable water quality in a residential premises, comprising: a measurement module electrically interposed between a sacrificial anode rod and the steel tank of a residential water heater, the measurement module including a shunt resistor for measuring galvanic current flowing between the anode and the tank through the potable water serving as electrolyte; an analog-to-digital converter sampling the voltage across the shunt resistor at a rate sufficient to resolve water draw events and tank refill transients; and a processor executing a trained regression model that maps features extracted from the galvanic current time series to at least one water quality parameter selected from the group consisting of total dissolved solids, pH, chlorine residual, and dissolved lead concentration.
- The system of claim 1, further comprising a potentiostat circuit capable of applying a small-amplitude AC voltage perturbation across the anode-to-tank electrochemical cell and measuring the complex impedance response at multiple frequencies to perform electrochemical impedance spectroscopy, wherein the processor extracts Randles equivalent circuit parameters including solution resistance, charge transfer resistance, and double-layer capacitance from the impedance spectrum and includes said parameters as features in the regression model.
- The system of claim 2, wherein the electrochemical impedance spectroscopy is performed periodically at intervals of 1 to 24 hours across a frequency range spanning at least 0.01 Hz to 10 kHz, with the perturbation amplitude limited to 10 mV or less to maintain the linear regime of the Butler-Volmer equation and avoid disturbing the natural corrosion process of the sacrificial anode.
- The system of claim 1, further comprising a temperature sensor in thermal contact with the water heater tank or anode rod, wherein the processor applies temperature compensation to the galvanic current measurements using an Arrhenius correction factor calibrated during an initial training period, thereby isolating water chemistry changes from temperature-driven current variations.
- The system of claim 1, wherein the processor performs online change-point detection on the temperature-compensated galvanic current time series to identify abrupt shifts in source water chemistry, and classifies detected change points by the pattern of simultaneous parameter changes into categories including source water switch, chlorine excursion, corrosion chemistry shift, and organic contamination event.
- The system of claim 5, wherein the corrosion chemistry shift category is identified by a change in charge transfer resistance and double-layer capacitance without a proportional change in solution resistance, indicating altered corrosion conditions at substantially constant total dissolved solids, and is assigned a critical alert priority due to potential for mobilization of lead from premise plumbing.
- A method for detecting lead service line corrosion at a residential premises, comprising: measuring galvanic current between a sacrificial anode rod and a steel water heater tank at a premises served by a water distribution system; performing periodic electrochemical impedance spectroscopy on the anode-tank electrochemical cell; extracting double-layer capacitance and low-frequency Warburg impedance features from the impedance spectrum; and applying a trained binary classifier to the extracted features to detect the presence of lead ions in the water at concentrations relevant to the EPA action level of 15 parts per billion, wherein lead ions are detected through their electrochemical signature of cementation on the magnesium anode surface producing characteristic changes in double-layer capacitance and Warburg impedance.
- The method of claim 7, wherein the classifier is applied preferentially to measurements collected during first-draw conditions when water has stagnated in the premise plumbing for at least 6 hours, corresponding to the sampling protocol specified in the EPA Lead and Copper Rule for identifying lead sources within premise plumbing.
- A method for network-level water quality surveillance comprising: collecting galvanic current and electrochemical impedance data from a plurality of instrumented water heaters distributed across a municipal water distribution network; performing change-point detection independently at each instrumented premises; correlating the timing of detected change points across premises with known distribution system topology to estimate the spatial extent and propagation velocity of water chemistry changes through the network; and generating a spatial contamination map that identifies the origin zone and flow path of the chemistry change.
- The method of claim 9, further comprising estimating hydraulic residence time at each instrumented premises from the time delay between a detected chemistry change at the treatment plant effluent and the corresponding change-point detection at the premises, and producing a dynamic water age map of the distribution network updated in real time from the measured propagation delays.
- The system of claim 1, further comprising an anode remaining life estimation module that integrates cumulative galvanic charge transfer over the operating life of the anode rod and computes remaining anode mass using Faraday's law of electrolysis, generating a replacement alert when estimated remaining mass falls below a configurable threshold.
- The system of claim 1, wherein the measurement module is packaged as one of: a retrofit collar installed between an existing anode rod and the water heater tank via the standard threaded anode port, or an integrated smart anode rod that replaces the factory anode rod and contains the measurement electronics within the hex head casting, or a factory-integrated module installed by the water heater manufacturer during production.
Prior Art References
- EPA, Lead and Copper Rule – Regulatory framework for lead and copper monitoring at residential taps
- EPA, Revised Total Coliform Rule (2013) – Routine bacterial monitoring requirements for water systems
- State of Michigan, Flint Water Crisis Documentation – Case study in distributed corrosion chemistry failure
- Randles, J.E.B., Discussions of the Faraday Society, 1947 – Foundational equivalent circuit model for electrochemical impedance analysis
- Adams, R.P. and MacKay, D.J.C., Bayesian Online Changepoint Detection, 2007 – Change-point detection algorithm applied to streaming electrochemical data
- UL 174 – Standard for Household Electric Storage Tank Water Heaters (anode rod requirements)
- AHRI Standard 700 – Performance rating of residential water heaters
- EIA, Residential Energy Consumption Survey (RECS, 2020) – U.S. water heater installed base statistics
- MIT Senseable City Lab, ES&T, 2019 – Point-of-use water quality sensor research
- EPA, Lead and Copper Rule Revisions (LCRR, 2021) – Service line inventory requirements by October 2027
- EPA, Water Security Initiative – Distribution system SCADA monitoring guidance