Reference Concept • Shield Performance Verification

Water Shield Radiation Evaluation (WSRE)

A requirements-style definition of the simulation and experimental evaluation campaign required to establish, with quantified confidence, what continuous water shielding delivers for long-duration orbital habitation — across the full range of credible shield thicknesses, with thickness treated as an evaluated parameter rather than an assumed answer.

What this page is

This page defines a neutral evaluation class for water-based radiation shielding at habitat scale. It specifies what must be modeled, measured, and reported before a thick-water shield can be credited with a protection level in any human-rating argument. It is not a mission design, not a product announcement, and not a claim about the answer.

Full GCR + SPE envelope Cross-code validated Phantom-level dosimetry Beam-anchored Uncertainty-quantified
Core idea: "How much water is enough?" has no rule-of-thumb answer. At multi-meter areal densities, dose reduction is real but sub-exponential, secondary neutrons and light fragments dominate the residual dose, and transport codes disagree most in exactly this regime. The honest output is not a defense of any particular thickness — it is the dose-versus-thickness curve, with confidence bounds, from which any architecture selects its operating point against its own occupancy requirements.

The thick-shield question

For short missions, shielding analysis can lean on thin-shield intuition: more mass, less dose, roughly exponentially. A habitat intended for decades of occupancy under career-dose limits (NASA‑STD‑3001 currently sets 600 mSv career effective dose) lives in a different regime:

Water is among the best bulk shielding materials per unit mass (high hydrogen content, favorable fragmentation), which is why it is the reference material here. The evaluation exists to convert that qualitative advantage into a quantified, defensible protection curve over the credible thickness range (roughly 1–5 m, ~100–500 g/cm²) — and to map how any point on that curve degrades under realistic geometry, penetrations, and off-nominal fill states.

Radiation environment envelope

The evaluation shall cover the complete free-space environment at the habitat's operating location (cislunar / deep space; no geomagnetic shielding credited, no trapped-belt population).

Galactic cosmic rays

CaseDefinitionPurpose
Solar minimum Full GCR spectrum at deep solar minimum (e.g., 2009/2019‑class modulation), per the current NASA reference model (Badhwar–O'Neill 2020 or successor; ISO 15390 as cross-check) Worst-case chronic exposure; sizing case for annual and career dose
Solar maximum Full GCR spectrum at strong solar maximum modulation Lower bound on chronic GCR; pairs with SPE-dominated risk
Mission average Multi-decade average over full solar cycles Career-dose accounting for permanent habitation

All GCR cases are isotropic over 4π. Species: protons, helium, and all heavy ions Z = 3–28 individually (with C, N, O, Ne, Mg, Si, Fe explicitly resolved — iron is the single largest dose-equivalent contributor at depth), plus electrons and positrons as a boundary term.

Solar particle events

CaseReferencePurpose
Small / moderateStatistical design events (e.g., ESP/PSYCHIC 90th & 99th percentile fluence)Frequent-event operational dose
August 1972King spectrum; classic hard-spectrum design caseHistorical worst-case class for BFO dose
October 1989Event series incl. Oct 19 ground-level eventHigh-fluence, extended-duration case
Oct–Nov 2003"Halloween storms" event seriesMulti-event sequence; recovery-interval stressor
September 2017Sept 10 GLE; modern well-instrumented eventBest-measured validation case
Carrington-classPublished fluence estimates for an 1859-scale eventBeyond-design-basis margin assessment

Each SPE case shall be evaluated as event-integrated fluence; the worst-case (1972-class and Carrington-class) shall additionally be evaluated with a time-dependent profile to produce peak hourly dose rates relevant to operational response and 30-day limits.

Species and energy coverage

SpeciesRangeWhy
Protons1 MeV – 100 GeVSpectrum peaks near ~1 GeV; the multi-GeV tail penetrates multi-meter water and drives pion/muon production
Helium1 MeV/n – 100 GeV/nSecond-largest primary fluence; strong secondary-neutron source
Heavy ions (Z=3–28)1 MeV/n – 50 GeV/nDominant dose-equivalent contributors via high-LET tracks and fragmentation chains
Electrons / positrons100 keV – 10 GeVMinor boundary term; completeness of the free-space field
Secondary neutronsThermal (0.025 eV) – 10 GeVDominant interior dose channel; must be transported to thermal capture (2.22 MeV n,γ on hydrogen)
Secondary γ, π, μ, fragmentsProduction threshold – full cascadeComplete accounting of the internally generated field at depth

Evaluation objectives

Primary • Produce interior dose quantities (absorbed dose, dose equivalent, effective dose) as a function of shield areal density, at all defined occupancy locations, under every environment case
• Resolve the interior field by species and LET spectrum, not dose totals alone
• Quantify shield transmission and secondary production (neutron and gamma source terms)
• Attach an explicit uncertainty budget to every reported quantity
Recommended (program-dependent) • Map 3D dose contours through the full habitat volume
• Evaluate degraded and construction-phase shield states
• Rank design variants (thickness, additives, layering) by dose-per-mass efficiency
• Produce operational products: storm-shelter margins, peak-rate alarms, occupancy guidance
Success criterion: outputs support a quantitative human-rating and habitability argument — annual and career doses with confidence intervals, traceable to reference environments, cross-validated codes, and beam data — not qualitative "water shields well" statements.

Functional requirements (requirements-style)

IDFunctionRequirement
WSRE-F-001Environment coverage The evaluation shall model the full GCR envelope (solar minimum, solar maximum, mission average) using current accepted NASA reference environments, applied isotropically over 4π.
BON2020-classIsotropic
WSRE-F-002SPE design cases The evaluation shall model the defined SPE set (statistical design events; Aug 1972; Oct 1989; Oct–Nov 2003; Sep 2017; Carrington-class), event-integrated, with time-dependent profiles for the bounding cases.
Historical worst caseBeyond-design-basis
WSRE-F-003Species & energy All primary species and energy ranges in the coverage table shall be transported, and all secondary species (n, γ, π, μ, light fragments) tracked to absorption or thermalization. Truncating the neutron cascade above thermal energies is non-compliant.
Full cascadeThermal capture
WSRE-F-004Geometry fidelity Transport shall use the as-designed 3D habitat geometry at each evaluated shield thickness, including a manufacturing tolerance band about that thickness, bulkheads, pipes, structural supports, compartmentation, docking ports, windows, and utility penetrations. A uniform spherical-shell idealization is a scoping tool only and shall not be credited as verification.
As-designed CADPenetrations
WSRE-F-005Degraded states The evaluation shall include off-nominal shield states: entrained air/void distributions, partially drained maintenance sections, construction-phase partial fill, and docked-vehicle shadowing/streaming configurations.
VoidsPartial drainConstruction phase
WSRE-F-006Rotation averaging Interior dose quantities shall be computed as time-averages over habitat rotation in the isotropic external field, with directional (hub-axis vs. ring-plane) asymmetries and partial shadowing explicitly resolved before averaging.
Spin-averagedDirectional
WSRE-F-007Water properties Baseline transport shall assume liquid water at operating density and temperature; sensitivity cases shall cover density/temperature variation, thermal gradients, borated water, and frozen water.
Baseline + sensitivities
WSRE-F-008Human dosimetry Organ doses shall be computed using ICRP reference adult male and female voxel (or mesh) phantoms at defined occupancy locations, including blood-forming organs, lens, skin, CNS, and reproductive organs, in representative postures (standing, sleeping).
ICRP 110/145Organ dose
WSRE-F-009Output metrics Required outputs: absorbed dose, dose equivalent (ICRP Q(LET)), effective dose, LET spectra, species-resolved fluence, secondary n/γ production, 3D dose maps, annual dose, career dose, storm-event dose, peak hourly rate, and shield transmission fraction — each with uncertainty.
Species-resolvedWith uncertainty
WSRE-F-010Code cross-comparison Results shall be produced by at least two independent Monte Carlo transport codes (from GEANT4, PHITS, FLUKA, MCNP6-class tools) plus deterministic scoping (HZETRN-class), with documented physics-model selections and inter-code deltas reported per output quantity.
≥2 MC codesDocumented physics lists
WSRE-F-011Convergence & uncertainty Statistical uncertainty at every scoring location shall meet defined targets (e.g., <5% relative on dose equivalent), and a total uncertainty budget shall combine statistical, cross-section, environment-model, and geometry contributions.
Uncertainty budgetDecision-ready
WSRE-F-012Experimental anchoring Simulation results shall be anchored to accelerator measurements on water targets: heavy ions (He, C, Si, Fe) at flight-relevant energies at an NSRL-class facility, and high-energy protons (hundreds of MeV to the multi-GeV GCR tail) at a proton facility, using thick-target staged areal densities, multiple incidence angles, and instrumentation including TEPCs, ion chambers, silicon telescopes, neutron spectrometry, and activation foils.
Beam validationThick targetMulti-facility
WSRE-F-013Standards traceability All dose quantities, quality factors, and limits shall trace to current standards (NASA‑STD‑3001, ICRP 103/123), and a verification matrix shall link each mission objective to its verifying analysis, test, or inspection.
NASA-STD-3001Traceability matrix
WSRE-F-014Thickness parameterization Interior dose quantities shall be produced as a function of shield areal density across the credible design range (~100–500 g/cm²), such that shield sizing is an output of the evaluation. No specific thickness shall be presumed sufficient in advance of results.
Dose-vs-thickness curveSizing as output

Note: quantitative thresholds (dose targets, agreement bounds, confidence levels) are architecture-dependent and are set by the program's human-rating requirements; representative values appear below.

Interior dosimetry locations

Shield inner wallBounding exposure surface; highest secondary-field intensity
Living quarters & sleep stationsHighest occupancy fraction; sleeping-posture organ doses
Medical bay & command centerContinuous-occupancy work positions; storm-shelter candidates
Agriculture volumesCrew transit occupancy plus biological-systems exposure
Ring centerlineBest-case reference point; defines the interior dose gradient
Worst-case occupancyLocations adjacent to penetrations, ports, and windows — the streaming-path check

Graduated evaluation tiers

The campaign proceeds in tiers, each retiring risk for the next:

The terrestrial analog is standard practice: reactor and accelerator shielding is never certified on a single unvalidated code. The same discipline applies here; the physics is harder, not easier.

Sensitivity studies

Shield & configuration • Graded and layered thickness distributions (e.g., locally augmented zones)
• Gap, void, and drain-down distributions
• Docked-vehicle attachments and shadowing
• Structural materials and internal equipment mass
• Polyethylene and boron additions; borated water
Environment & operations • Solar-cycle phase and modulation uncertainty
• Rotation rate and orientation
• Orbital location variants
• Shield aging and water-quality drift
• Frozen vs. liquid phase, thermal gradients

Acceptance criterion classes

The evaluation shall report against the following criterion classes. Numerical thresholds are deliberately not asserted here: setting a target dose in advance of the evaluation would presume its conclusion. Thresholds derive from the program's human-rating requirements and current NASA‑STD‑3001 / ICRP guidance at the time of evaluation.

Criterion classWhat is reported
Chronic exposureAnnual interior effective dose per occupancy location, per GCR case, with confidence interval, stated against then-current career-limit guidance
Shield sizingThe minimum areal density meeting a given occupancy requirement, read from the dose-vs-thickness curve with its confidence bounds — the sizing decision follows the data
Occupancy durationThe maximum residency duration each location supports under applicable limits — an output of the evaluation, not an input to it
Storm exposureBFO and effective dose per SPE design case, stated against short-term (30-day / annual) limits, with and without operational sheltering assumptions
Relative performanceDose equivalent versus an equal-mass aluminum reference shield across the full environment envelope
Model confidenceInter-code agreement per output quantity, and experiment–simulation agreement at validated beam conditions, each with documented bounds

Data products

Minimum • Simulation report per environment case, with uncertainty
• Geometry files and CAD assumptions record
• Material and cross-section libraries used
• Dose maps and occupancy-location dose tables
• Verification matrix and validation report
• Executive summary legible to non-specialists
Optional / program-dependent • 3D interior dose contours
• Design-variant ranking (dose per unit mass)
• Operational storm-response products
• Risk assessment and recommendations
• Public benchmark dataset for cross-team reuse

Usability requirement: outputs must be legible to program managers and human-rating reviewers, not only radiation physicists.

Ballpark cost (ROM)

Order-of-magnitude, planning-grade only: a Tier 1–2 simulation campaign of this scope is a low single-digit $M effort dominated by analyst time and compute; a Tier 3 accelerator validation campaign is likewise low single-digit $M, dominated by beam time, targets, and instrumentation. Either figure is small against the cost of committing habitat-scale structure to an unverified protection level — or of over-building shield mass the physics does not require.

Relationship to mission risk

Absence of WSRE-class evaluation before committing to a shield architecture constitutes a known program risk in both directions: an under-performing shield discovered after construction is unrecoverable, and an over-built shield wastes the single most expensive commodity in the architecture — delivered mass. Chronic GCR exposure is the binding constraint on permanent habitation; it deserves the same evidentiary standard as structure and life support.

WSRE-class evaluation reduces epistemic uncertainty; it does not eliminate radiobiological uncertainty in risk models, and it does not replace onboard dosimetry and operational monitoring.

Standardization potential

Standardizing WSRE-class geometry files, environment definitions, and output formats enables cross-comparable results between independent teams and cumulative refinement as codes and cross-section libraries improve. Over time this constitutes a shared verification baseline for water-shielded habitats generally — a benchmark class, not a one-off study.