!Hot springs can cause serious injury or death. Stay on marked paths, follow local signs, never cross barriers, and do not rely on this site for safe navigation. Safety information →

Methodologies · 1000 Springs Project

Sampling & Analysis Protocols

Field sampling techniques, laboratory processing procedures, sequencing protocols, and health & safety guidelines used across the 1000 Springs Project.

1,000geothermal features targeted
~70,000sequences read per sample
97%OTU similarity threshold
250 bptrimmed read length

Overview diagram

1000 Springs Project data flow, from field sampling through analysis and computing to the online database
Fig. 1. Data flow from field sampling through laboratory analysis and computing to the online database.

01: Feature Sampling

Sample site selection & water collection

Site selection

The research targets a broad range of physicochemical conditions in geothermally-influenced springs within the Taupō Volcanic Zone. Water columns are sampled from springs. Sediments, biofilms, mudpots, and heated soils are excluded.

The goal is to collect 1,000 samples, primarily one per individual feature. Repeated sampling occurs where conditions vary significantly, for temporal studies, or for quality control purposes.

Collection technique

A telescopic sampling pole fitted with sterile containers or custom-built samplers collects water across the full depth of the spring, ensuring homogenous and representative samples.

A field tablet and custom app record physical, chemical, and metadata in the field, eliminating double-handling of field notes and reducing transcription errors. Metadata uploads instantly to the database.

Sample water distribution

Table 1. Volume and container allocation per sample.

VolumeContainerPurpose
2,150 mLSterile 2 L polypropylene bottleMicrobial diversity & chemistry analysis
~330 mLRubber-necked bottleH₂S, bicarbonate & chloride analysis
500 mLSterile polypropylene bottlePhysical attributes (conductivity, dissolved oxygen)
~60 mLSyringeDissolved gas content
~10 gSediment containerSediment sample (where possible), earmarked for future research
Collecting water samples using the telescopic sampling pole.
Fig. 2. Collecting water samples using the telescopic sampling pole.
Field sampling documentation and mobile lab setup.
Fig. 3. Field sampling documentation and mobile lab setup.

Field metadata recorded

Feature type (spring, stream, etc.)
Area description
Sampling date and staff
Feature location & GPS coordinates
Spring size
Ebullition (bubbling / water upflow)
Spring colour
Images and bird's-eye-view sketches for future resampling

02: Field & Laboratory Processing

Measurements, filtration & chemistry

Field measurements

Processing occurs in a 4WD field van configured as a mobile laboratory to minimise temperature changes, dissolved oxygen changes, and microbial activity effects.

pHHanna field multiparameter meter
Oxidation-reduction potential (ORP / redox)Hanna field multiparameter meter
Conductivity (salinity proxy)Hanna field multiparameter meter
Total dissolved solidsHanna field multiparameter meter
TurbidityHanna field multiparameter meter
Dissolved oxygenHanna field multiparameter meter
Ferrous iron (Fe²⁺)Field spectrophotometer (immediate)

Microorganism collection

A battery-operated peristaltic pump pushes up to 2 litres through a 0.2 µm Sterivex column filter (Millipore). After extensive field testing, Sterivex filters demonstrated higher filtration volume capacity than tested alternatives.

Filtrate water is collected in 15 mL and 50 mL tubes for additional chemistry analyses.

All filtrate samples, filters, and soil samples are cooled to 4 °C and transported immediately to the Wairakei laboratory for storage until analysis.

Chemistry analysis locations

NZGALNZ Geothermal Analytical Laboratory, GNS Science
ERGExtremophile Research Group, GNS Science
UWaikatoUniversity of Waikato, School of Science
The 4WD mobile laboratory used for field processing.
Fig. 4. The 4WD mobile laboratory used for field processing.
Filtering samples through the Sterivex column.
Fig. 5. Filtering samples through the Sterivex column.

Sample processing summary

Table 2. Container, processing method, storage, and analysis location per parameter.

ContainerProcessingStorageParameterAnalytical methodLocation
2 L bottleFiltered−20 °CMicrobial diversityDNA extraction & sequencingUniversity of Waikato
500 mL bottleRawn/aPhysical propertiesMultiparameter field meterOn site
500 mL bottleFiltered (0.22 µm)n/aFe²⁺UV SpectroscopyOn site
330 mL rubber-sealedRaw4 °CH₂SUV SpectroscopyNZGAL
330 mL rubber-sealedRaw4 °CHCO₃⁻Automated titrationNZGAL
330 mL rubber-sealedRaw4 °CCl⁻Automated titrationNZGAL
50 mL tubeFiltered4 °CSO₄²⁻Ion chromatographyNZGAL
50 mL tubeFiltered−20 °CNH₄⁺, PO₄³⁻, NO₂⁻, NO₃⁻Flow injection analysisUniversity of Waikato
50 mL tubeFiltered4 °CBack-up samplen/an/a
15 mL tubeFiltered & acidified4 °CElements (ICP-MS)Inductively coupled plasma mass spectrometryUniversity of Waikato
15 mL tubeFiltered & alkalified4 °CAs, SbICP atomic emission spectroscopyNZGAL
50 mL syringeFilteredRoom temp.H₂, CH₄, COGas chromatographyExtremophile Research Group, GNS

Elements measured via ICP-MS

30 elements quantified by inductively coupled plasma mass spectrometry at the University of Waikato.

BoronLithiumSodiumMagnesiumAluminiumSiliconSulphurPotassiumCalciumVanadiumChromiumManganeseIronCobaltNickelCopperZincArsenicSeleniumBromineRubidiumMolybdenumSilverCadmiumCaesiumBariumMercuryThalliumLeadUranium

03: Sequencing

Microbial diversity assessment

Soil and freshwater can contain up to one billion and ten million cells respectively, making exhaustive identification impractical. The project instead extracts total microbial DNA, sequences a universal gene, and uses bioinformatics pipelines to assign taxonomy at scale using Ion Torrent Next-Generation Sequencing.

Phase 1

DNA Extraction

Summary: Chemical and physical techniques extract, purify, and concentrate DNA from microorganisms. Before cell disruption, sterile skim milk blocks possible binding sites on solids/sediments captured during filtration. Captured cells in the Sterivex filter are disrupted using cetyl trimethylammonium bromide (CTAB). Proteins, lipids, and ribonucleic acids are chemically removed using phenol-chloroform-isoamyl alcohol and chloroform-isoamyl alcohol. DNA is then concentrated via magnetic bead-adsorption.

1.

Cell disruption

Sterivex filters thaw on ice. 1.0 mL of 0.8% w/v skim milk powder solution is added and incubated for 30 min at 65 °C. 0.4 mL CTAB, 0.2 mL phosphate-buffered saline, and 0.1 mL sodium dodecyl sulfate are added and vortexed for 1 minute. The filter is placed in a 50 mL sterile capped tube and incubated at 150 rpm / 65 °C for 30 minutes. After briefly cooling on ice, the aqueous component is pushed through the syringe and collected in 2 mL tubes.

2.

DNA purification

Equal volume (1:1) of chloroform:isoamyl alcohol (24:1) is added to each tube, vortexed for 10 s, and centrifuged at 12,400 rpm for 10 minutes. The aqueous phase (typically top) is transferred to a new tube. 300 µL chloroform:isoamyl (24:1) is added to the aqueous phase, mixed on a Hula mixer for 20 minutes, then centrifuged for 10 minutes at 12,400 rpm. The aqueous phase is collected, avoiding the bottom layer.

3.

DNA concentration

Equal volume of 100% ethanol and 40 µL MoBio beads are added to the aqueous phase. Samples are agitated every 30 s–5 min, then placed on a magnetic stand for 2 minutes. Supernatant is removed and beads are washed with 500 µL of 100% ethanol, cleaning tube walls. All ethanol is removed and beads are resuspended in 30 µL Tris (pH 8.0) by gentle pipetting and brief vortexing.

Phase 2

DNA Amplification & NGS Preparation

Summary: Microbial diversity is determined by sequencing the 16S ribosomal RNA gene (~1,540 nucleotides), a short gene possessed by all bacteria and archaea. Total DNA is extracted and 16S rRNA genes are amplified by PCR modified for Ion Torrent Next-Generation Sequencing. On average, ~70,000 sequences (~270 bp each) are read per sample.

1.

PCR amplification

PCR amplifies the microbial community 16S rRNA genes using universal bacterial/archaeal primers from the Earth Microbiome Project:

F5155′-GTGCCAGCMGCCGCGGTAA-3′
R8065′-GGACTACVSGGGTATCTAAT-3′

The forward primer contains an A-adaptor, a 10–12 nucleotide sample-specific barcode, a 'GAT' barcode adaptor, and the F515 sequencing primer. The reverse primer includes a P1 sequence and R806 sequencing primer.

2.

PCR reaction composition

0.5 µL eachForward & reverse primers (10 µM)
1.0 µLBovine serum albumin (10 mg/mL)
3.0 µLdNTPs (8 mM)
3.0 µL10× PCR buffer
3.0 µLMgCl₂ (50 mM)
2.0 µLSoil community DNA (0.5 ng/µL template)

25 µL total reaction volume. Master mix is pre-treated with ethidium monoazide bromide (EMA) to remove trace DNA contaminants before primer, Taq, and template addition.

Thermocycling

Initial denaturation 94 °C / 3 min → 30 cycles: 94 °C / 45 s · 50 °C / 1 min · 72 °C / 1.5 min → Final extension 72 °C / 10 min

3.

NGS preparation & sequencing

Amplicons from triplicated PCR products are purified using SPRIselect (Beckman Coulter) to remove small nucleotide fragments. Quality and concentration are verified and adjusted to 26 pM via HS Qubit 2.0 (Life Technologies) and 9100 BioAnalyser (Agilent Technologies). Amplicons from all samples are pooled for emulsion PCR. Sequencing is conducted at the Waikato DNA Sequencing Facility using Ion Torrent PGM (Life Technologies) with Ion 318v2 chip and 400 bp chemistry.

Ion Torrent PGM sequencer.
Fig. 6. Ion Torrent PGM sequencer.
Loading the Ion 318v2 chip.
Fig. 7. Loading the Ion 318v2 chip.
Ion Torrent sequencing chips.
Fig. 8. Ion Torrent sequencing chips.
Phase 3

Post-sequencing Bioinformatics

Summary: ~70,000 raw reads per sample are processed by customised bioinformatic pipelines. Erroneous reads are removed, tagging sequences are stripped, and reads are screened for known errors and trimmed to 250 bp. Similar sequences are clustered into Operational Taxonomic Units (OTUs) at 97% similarity. OTUs are taxonomically ranked by comparison with microbial databases. For this explorer, OTUs with > 5 reads per sample are presented graphically.

Pipeline details

Raw reads in FASTQ format are processed using a custom pipeline based on Mothur and USEARCH.

Reads outside 275–345 bp are removed
Reads with any error in barcode or forward PCR primer are removed
Barcode and forward primer sequences are stripped from remaining reads
Q-score filtering removes reads with ≥ 3 expected errors
Reads are globally trimmed to 250 bp
After de-replication, unique reads are sorted by abundance
USEARCH generates OTUs via centroid-based clustering (97% similarity) and removes chimeras
OTUs are taxonomically assigned using RDP Classifier with RDP Release 11.2

04: Health & Safety

Geothermal field safety

Important safety notice

Sampling springs in geothermal environments presents serious danger and should not be attempted by the public. Ground surrounding hot springs can be hot, unstable, and undercut. Springs often emit high concentrations of toxic and/or asphyxiating gases. All project participants have extensive experience moving through and sampling geothermal ecosystems. The project attempts to minimise environmental impact from sampling.

Team equipped for geothermal field work.
Fig. 9. Team equipped for geothermal field work.
Personal protective equipment for sampling.
Fig. 10. Personal protective equipment for sampling.

Personal protective equipment & emergency kit

Comprehensive burns kit with full-body gel blanket
1st aid kit
Collapsible bucket & water bladder
Modified throwbag
Gas monitor
GPS unit
Satellite phone
Personal locator beacon
Laminated site info sheets (contacts, hazards, recommended PPE per site)

Clothing & footwear

Footwear

Leather boots topped with neoprene puttees (gaiters) are preferred over gumboots; they provide better ankle support and protection in unstable terrain.

Trousers

Locally made Cactus SuperTrousers protect against vegetation and slow water/mud infiltration if accidents occur.

General guidance for visitors

When visiting geothermal environments: do not leave formed tracks, do not approach springs or steaming soils, and always follow posted safety instructions.

References

Literature cited

Edgar, R.C. (2013). UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nature Methods 10, 996–998.

Rueckert, A. and Morgan, H.W. (2007). Removal of contaminating DNA from polymerase chain reaction using ethidium monoazide. J Microbiol Meth, 68(3), 596–600.

Schloss, P.D., Westcott, S.L., Ryabin, T., Hall, J.R., Hartmann, M., Hollister, E.B., Lesniewski, R.A., Oakley, B.B., Parks, D.H., Robinson, C.J., Sahl, J.W., Stres, B., Thallinger, G.G., Van Horn, D.J. and Weber, C.F. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol, 75(23), 7537–7541.

Wang Q., Garrity, G.M., Tiedje, J.M. and Cole, J.R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol, 73(16), 5261–5267.

Whitman, W.B., Coleman, D.C. and Wiebe, W.J. (1998). Prokaryotes, the unseen majority. Proc Natl Acad Sci USA, 95(12), 6578–6583.

Wilson, A.D. (1960). The micro-determination of ferrous iron in silicate minerals by a volumetric and a colorimetric method. Analyst, 85, 823–827.