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.
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.
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.
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).
| Case | Definition | Purpose |
|---|---|---|
| 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.
| Case | Reference | Purpose |
|---|---|---|
| Small / moderate | Statistical design events (e.g., ESP/PSYCHIC 90th & 99th percentile fluence) | Frequent-event operational dose |
| August 1972 | King spectrum; classic hard-spectrum design case | Historical worst-case class for BFO dose |
| October 1989 | Event series incl. Oct 19 ground-level event | High-fluence, extended-duration case |
| Oct–Nov 2003 | "Halloween storms" event series | Multi-event sequence; recovery-interval stressor |
| September 2017 | Sept 10 GLE; modern well-instrumented event | Best-measured validation case |
| Carrington-class | Published fluence estimates for an 1859-scale event | Beyond-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 | Range | Why |
|---|---|---|
| Protons | 1 MeV – 100 GeV | Spectrum peaks near ~1 GeV; the multi-GeV tail penetrates multi-meter water and drives pion/muon production |
| Helium | 1 MeV/n – 100 GeV/n | Second-largest primary fluence; strong secondary-neutron source |
| Heavy ions (Z=3–28) | 1 MeV/n – 50 GeV/n | Dominant dose-equivalent contributors via high-LET tracks and fragmentation chains |
| Electrons / positrons | 100 keV – 10 GeV | Minor boundary term; completeness of the free-space field |
| Secondary neutrons | Thermal (0.025 eV) – 10 GeV | Dominant interior dose channel; must be transported to thermal capture (2.22 MeV n,γ on hydrogen) |
| Secondary γ, π, μ, fragments | Production threshold – full cascade | Complete accounting of the internally generated field at depth |
| ID | Function | Requirement |
|---|---|---|
| WSRE-F-001 | Environment 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-002 | SPE 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-003 | Species & 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-004 | Geometry 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-005 | Degraded 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-006 | Rotation 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-007 | Water 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-008 | Human 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-009 | Output 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-010 | Code 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-011 | Convergence & 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-012 | Experimental 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-013 | Standards 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-014 | Thickness 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.
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.
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 class | What is reported |
|---|---|
| Chronic exposure | Annual interior effective dose per occupancy location, per GCR case, with confidence interval, stated against then-current career-limit guidance |
| Shield sizing | The 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 duration | The maximum residency duration each location supports under applicable limits — an output of the evaluation, not an input to it |
| Storm exposure | BFO and effective dose per SPE design case, stated against short-term (30-day / annual) limits, with and without operational sheltering assumptions |
| Relative performance | Dose equivalent versus an equal-mass aluminum reference shield across the full environment envelope |
| Model confidence | Inter-code agreement per output quantity, and experiment–simulation agreement at validated beam conditions, each with documented bounds |
Usability requirement: outputs must be legible to program managers and human-rating reviewers, not only radiation physicists.
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.
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.
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.