What is an Aquifer Thermal Energy Storage (ATES) doublet system?

An ATES doublet consists of two hydraulic wells drilled into a permeable confined or semi-confined aquifer layer: a 'warm well' (typically 18°C to 24°C) and a 'cold well' (typically 6°C to 10°C). During summer, cold groundwater is extracted to cool the greenhouse, absorbs heat via heat exchangers, and is reinjected into the warm well. In winter, flow direction is reversed: warm water is extracted as the evaporator source for industrial heat pumps, chilled down, and returned to the cold well.

How is thermal breakthrough prevented between ATES doublet wells?

According to VDI 4640 Blatt 4 standards, thermal breakthrough occurs if the thermal plume radius R_th overlaps between wells. The minimum well separation distance L must satisfy L ≥ 2.5 × R_th, where R_th = sqrt((c_w × V_pump) / (π × H × c_aq)), with c_w as volumetric heat capacity of water, c_aq as volumetric heat capacity of the saturated aquifer, H as aquifer thickness, and V_pump as seasonal pumped water volume.

What COP and energy savings can commercial greenhouses achieve with ATES and industrial heat pumps?

Because groundwater provides a stable source temperature of 12°C to 20°C year-round (compared to freezing winter ambient air), industrial water-to-water heat pumps achieve a Seasonal Coefficient of Performance (SCOP) between 4.8 and 5.6. Coupled with free-cooling in summer (COP > 25), commercial greenhouses reduce annual fossil natural gas consumption by 70% to 90% and primary energy costs by over 55%.

Industrial Geothermal Whitepaper • Vol. 75

Aquifer Thermal Energy Storage (ATES) & Industrial Heat Pump Integration in Commercial Greenhouses

✍️ By Inwoo Hwang 📅 October 1, 2026 ⏱️ 14 Min Read 🏷️ Geothermal HVAC • VDI 4640 • Decarbonization

Executive Engineering Summary: Modern Controlled Environment Agriculture (CEA) facilities and semi-closed Venlo greenhouses face massive thermodynamic asymmetry: an enormous surplus of solar thermal radiation in summer rejected via roof ventilation, followed by intense fossil gas heating loads in winter ($140\text{--}200\text{ W/m}^2$). Aquifer Thermal Energy Storage (ATES) solves this mismatch by utilizing deep permeable groundwater strata as an inter-seasonal thermal battery. This monograph provides rigorous engineering equations for Dupuit-Thiem doublet well hydraulics, VDI 4640 Blatt 4 thermal plume radius sizing, industrial water-to-water heat pump thermodynamic coupling ($\text{SCOP} \ge 5.2$), and capital funding capture through Germany’s BEW Modul 2 ($45\%$ grant) and US IRA Section 48C ($30\text{--}50\%$ ITC).

1. Fundamentals of ATES Doublet Architecture in CEA

Conventional commercial greenhouses discharge excess sensible and latent solar heat directly to the atmosphere during summer through continuous ridge vents. Concurrently, high-pressure fogging or mechanical chillers consume significant electrical energy to maintain sub-$28^\circ\text{C}$ canopies. In winter, the same structure burns fossil natural gas ($25\text{--}45\text{ m}^3/\text{m}^2/\text{year}$) or district heating to offset skin transmission losses.

An Aquifer Thermal Energy Storage (ATES) system balances this seasonal imbalance by operating two hydraulic wells—a Warm Well and a Cold Well—drilled into a saturated confined sand or gravel aquifer layer ($30\text{--}150\text{ m}$ depth).

🔁 Seasonal Reversible Doublet Operational Principle

Summer Cooling Mode: Chilled water ($6\text{--}10^\circ\text{C}$) is extracted from the Cold Well. It passes through plate heat exchangers in greenhouse Air Treatment Units (ATUs) to provide direct "free cooling" (EER > 25). The warmed water ($18\text{--}24^\circ\text{C}$) is injected into the Warm Well, charging the subterranean thermal plume.
Winter Heating Mode: Water at $18\text{--}24^\circ\text{C}$ is extracted from the Warm Well and directed into the evaporator of an industrial water-to-water heat pump. The heat pump extracts thermal energy, boosting supply temperatures to $40\text{--}45^\circ\text{C}$ for under-bench rail heating, while chilled groundwater ($5\text{--}8^\circ\text{C}$) is returned to the Cold Well.

2. Doublet Well Hydraulics & Dupuit-Thiem Drawdown Formulations

Reliable continuous operation demands that pumping and injection rates do not exceed aquifer transmissivity, preventing screen cavitation, sand fluidization, or surface blowouts. Under steady-state radial flow in a confined aquifer of thickness $H$ [m] and hydraulic conductivity $K$ [m/s], the hydraulic transmissivity $T$ is defined as:

T = K \cdot H \quad [\text{m}^2/\text{s}]

The steady-state hydraulic drawdown (or injection head buildup) $s_w$ [m] at the well screen of radius $r_w$ [m] with total pumping flow rate $Q$ [$\text{m}^3/\text{s}$] and radius of hydraulic influence $R_0$ [m] is governed by the Dupuit-Thiem equation:

s_w = \frac{Q}{2\pi T} \ln\left(\frac{R_0}{r_w}\right) + \Delta s_{\text{skin}} + \Delta s_{\text{turb}}

Where $\Delta s_{\text{skin}}$ represents wellbore skin damage factor and $\Delta s_{\text{turb}} = B \cdot Q + C \cdot Q^2$ accounts for turbulent non-Darcian head losses across the gravel pack and Johnson screen slots.

Critical Well Screen Velocity Threshold

To prevent mineral incrustation, iron hydroxide precipitation ($\text{Fe}^{2+} \to \text{Fe(OH)}_3$), and mechanical clogging, the entrance velocity $v_{\text{entrance}}$ of groundwater passing through open screen area $A_{\text{open}}$ must be rigorously restricted below Siccardt's critical threshold:

v_{\text{entrance}} = \frac{Q}{A_{\text{open}} \cdot L_{\text{screen}}} \le \frac{\sqrt{K}}{30} \quad \text{and strictly} \le 0.03\text{ m/s}

3. VDI 4640 Blatt 4 Thermal Plume Sizing & Spacing Criteria

The volume of subterranean rock and water affected by seasonal heat injection forms a radial thermal cylinder. According to German standard VDI 4640 Blatt 4 (Thermal Use of the Underground - Aquifer Thermal Energy Storage), the theoretical thermal radius $R_{\text{th}}$ [m] is determined by:

R_{\text{th}} = \sqrt{\frac{c_w \cdot V_{\text{season}}}{\pi \cdot H \cdot c_{\text{aq}}}}

Where:

$c_w$: Volumetric heat capacity of water ($\approx 4.18\text{ MJ}/(\text{m}^3\cdot\text{K})$).
$c_{\text{aq}} = (1 - n)c_{\text{rock}} + n \cdot c_w$: Effective volumetric heat capacity of the saturated porous medium ($\approx 2.4\text{--}2.8\text{ MJ}/(\text{m}^3\cdot\text{K})$ for sandy gravels with porosity $n = 0.25\text{--}0.35$).
$V_{\text{season}}$: Total cumulative groundwater volume injected over the summer or winter half-year [$\text{m}^3$].
$H$: Screened net aquifer thickness [m].

Thermal Recovery Efficiency ($\eta_{\text{th}}$) & Separation Distance

Natural ambient groundwater flow ($v_{\text{darcy}} = K \cdot i$) advects and disperses the stored thermal bubble downstream. The seasonal thermal recovery factor $\eta_{\text{th}}$ represents the fraction of injected enthalpy retrievable during the reverse cycle:

\eta_{\text{th}} = \frac{E_{\text{extracted}}}{E_{\text{injected}}} = \frac{\int_0^{\tau_{\text{ext}}} \rho_w c_w Q(t) [T(t) - T_{\text{ambient}}] dt}{\int_0^{\tau_{\text{inj}}} \rho_w c_w Q(t) [T_{\text{inj}} - T_{\text{ambient}}] dt} \ge 80\text{--}85\%

To prevent premature thermal breakthrough—wherein the cold plume migrates into the warm well or vice versa, degrading system COP—the minimum inter-well doublet separation distance $L_{\text{doublet}}$ must satisfy:

L_{\text{doublet}} \ge 2.5 \cdot R_{\text{th}} + \frac{v_{\text{darcy}}}{n} \cdot \Delta t_{\text{inter-season}}

For typical commercial 1 to 5-hectare installations with $V_{\text{season}} = 150,000\text{ m}^3$ and $H = 20\text{ m}$, $R_{\text{th}}$ ranges between $55\text{--}75\text{ m}$, mandating well spacings of $L_{\text{doublet}} \ge 140\text{--}190\text{ m}$.

4. High-Lift Water-to-Water Heat Pump Integration

Because groundwater extracted from the Warm Well enters the evaporator at $18\text{--}24^\circ\text{C}$ (substantially higher than typical ambient winter air at $-5\text{--}5^\circ\text{C}$ or surface water at $4^\circ\text{C}$), the Carnot temperature lift $\Delta T_{\text{lift}} = T_{\text{cond}} - T_{\text{evap}}$ is cut by half.

\text{COP}_{\text{Carnot}} = \frac{T_{\text{cond}}}{T_{\text{cond}} - T_{\text{evap}}} \quad \implies \quad \text{COP}_{\text{actual}} = \eta_{\text{Carnot}} \cdot \text{COP}_{\text{Carnot}} \ge 5.0\text{--}5.6

Where $\eta_{\text{Carnot}} \approx 0.52\text{--}0.58$ for industrial semi-hermetic twin-screw or magnetic-bearing centrifugal chillers using low-GWP refrigerants (such as HFO-1234ze, R513A, or natural Ammonia R717 / CO₂ R744 transcritical loops).

Heating Technology	Heat Source Temp (°C)	Supply Temp (°C)	Seasonal COP / Efficiency	Primary Energy CO₂ (kg/MWh)
Natural Gas Boiler (Condensing)	Chemical Combustion	70 / 50	94% (HHV)	202 kg CO₂/MWh
Air-Source Industrial HP (ASHP)	Ambient Air (-2°C avg)	45 / 38	2.8 ~ 3.4	90 ~ 110 kg CO₂/MWh
Closed-Loop Borehole (BHE)	Ground Brine (4°C ~ 8°C)	45 / 38	3.8 ~ 4.2	72 ~ 80 kg CO₂/MWh
ATES Doublet + Industrial W2W HP	Warm Aquifer (18°C ~ 22°C)	42 / 35	5.2 ~ 5.8 (Free Cool > 25)	35 ~ 48 kg CO₂/MWh

5. Decarbonization Subsidies: German BEW & US IRA Sizing

While ATES and industrial heat pumps deliver unparalleled operational expenditure (OPEX) reductions, capital expenditure (CAPEX) for deep hydrogeological exploratory drilling, double-screened stainless well installation, and megawatt-scale heat pumps is substantial (€350,000 to €850,000 per doublet pair). Fortunately, major federal decarbonization programs in 2026 drastically de-risk capital deployment:

🇩🇪 German BEW (Bundesförderung für effiziente Wärmenetze)

Modul 2 (Investitionsförderung): Provides up to 40% to 45% non-repayable direct grants for the installation of large-scale geothermal heat pumps, ATES doublet drilling, and district heating connections for commercial greenhouse networks.
KfW 295 / 432 Concessional Loans: Long-term low-interest debt financing with debt relief redemption bonuses up to 20% for projects eliminating more than 70% of fossil natural gas baseload.

🇺🇸 US Inflation Reduction Act (IRA Section 48 & 48C)

Section 48 Investment Tax Credit (ITC): 30% base tax credit for commercial geothermal heat pumps and subterranean thermal storage, stackable with +10% Domestic Content and +10% Energy Community bonuses for a total ITC up to 50%.
USDA REAP (Rural Energy for America Program): Direct grants up to $500,000 or 50% of total project cost for agricultural producers adopting renewable thermal energy systems.

Run Live Doublet Well & Financial Payback Simulations

Calculate exact Dupuit-Thiem drawdown ($s_w$), VDI 4640 thermal plume radius ($R_{\text{th}}$), inter-well spacing, seasonal heating/cooling COP, and BEW/IRA grant net NPV with our verified engineering engine.

🚀 Launch Tool 38: ATES Simulator & Payback Engine →

6. Field Engineering Checklist & Water Chemistry Protection

A successful ATES deployment requires strict adherence to geochemical and physical protection protocols:

Continuous Nitrogen Blanketing or Closed Degassing: Ensure 100% airtight groundwater loops. Exposure to atmospheric oxygen oxidizes dissolved ferrous iron ($\text{Fe}^{2+}$) into ferric hydroxide ($\text{Fe(OH)}_3$) flocculants, permanently clogging injection well screens within weeks.
Dual Backwash Redundancy: Incorporate automated high-rate backwash cycles (1.5× normal production flow rate) equipped with 50 µm hydrocyclone sand separators to flush fine silt and bio-matting.
Submersible VFD Well Pumps with Non-Return Valves: Modulate pumping flow directly with greenhouse heating/cooling demand to prevent cycling pressure transients (water hammer) across deep sandstone formations.
Dual Plate Heat Exchangers (Titanium / SMO 254): Isolate groundwater hydraulically from the internal greenhouse heating loops to safeguard boilers, ATU coils, and crops from mineral scaling and corrosion.