gea-program 0.1.0__py3-none-any.whl

gea/blind_harness.py

@@ -0,0 +1,108 @@
+# This Source Code Form is subject to the terms of the Mozilla Public
+# License, v. 2.0. If a copy of the MPL was not distributed with this
+# file, You can obtain one at https://mozilla.org/MPL/2.0/.
+"""blind_harness - Part 5 of the subsurface surveying tool: THE
+BLIND-VALIDATION HARNESS (v1.82.0) - the standing accuracy report that makes
+the tool credible, not just built.
+
+METHOD (leave-one-out, no self-grading):
+    For every supported property pair in the library (the strata-join joint
+    tables plus the KTB site-native pairs), each co-located observation is
+    HELD OUT in turn; the target is predicted from the REMAINING
+    observations by the same k-nearest conditional the inverse engine uses;
+    the held-out truth scores the prediction. Per pair we report:
+
+      n            - held-out trials
+      mae          - mean absolute error, in the property's own units
+      mae_pct      - MAE as a percentage of the truth's mean
+      coverage     - fraction of trials where the truth fell inside the
+                     prediction's own +/- 1-sigma spread (an HONESTY metric:
+                     a well-calibrated ~68 percent is the target; near-100
+                     means the spreads are too timid, near-0 too bold)
+
+    Pairs below MIN_TRIALS refuse. Nothing is tuned to this harness; it
+    grades the same machinery clients get, on the same library.
+
+DOCTRINE: this report is regenerated live - by the gate, by acceptance, and
+by the PAPER_2258-lineage dispatches - so the accuracy table can never be a
+stale marketing snapshot. Misses appear next to hits, forever.
+"""
+
+from __future__ import annotations
+
+import statistics
+from typing import Dict, List
+
+from . import strata_join as SJ
+from .inverse_engine import _site_native_pairs, _conditional_from_pairs
+
+MIN_TRIALS = 10
+
+
+def _loo_pairs(pairs: List[tuple], k: int = 7) -> Dict:
+    """Leave-one-out over raw (given, target) pairs."""
+    n = len(pairs)
+    if n < MIN_TRIALS:
+        return {'status': 'REFUSED_THIN_DATA', 'n': n, 'min_n': MIN_TRIALS}
+    errs, hits, truths = [], 0, []
+    for i in range(n):
+        rest = pairs[:i] + pairs[i + 1:]
+        g, t = pairs[i]
+        c = _conditional_from_pairs(rest, g, k)
+        if c['status'] != 'OK':
+            continue
+        err = abs(c['estimate'] - t)
+        errs.append(err)
+        truths.append(abs(t))
+        if c['std'] > 0 and err <= c['std']:
+            hits += 1
+    if len(errs) < MIN_TRIALS:
+        return {'status': 'REFUSED_THIN_DATA', 'n': len(errs),
+                'min_n': MIN_TRIALS}
+    mae = statistics.mean(errs)
+    tmean = statistics.mean(truths)
+    return {'status': 'OK', 'n': len(errs), 'mae': mae,
+            'mae_pct': (mae / tmean * 100.0 if tmean else float('nan')),
+            'coverage': hits / len(errs)}
+
+
+def accuracy_report() -> Dict:
+    """The standing table: every supported pair in the library, blind-scored.
+
+    Sources: (1) strata-join joint tables (bin-level pairs, no archive
+    interpolation), (2) the KTB site-native station pairs (entry 52,
+    washouts excluded as disclosed there)."""
+    rows = []
+    # 1) strata-join joint tables
+    for well, (_, props) in SJ.WELL_GROUPS.items():
+        names = sorted(props)
+        for i, a in enumerate(names):
+            for b in names[i + 1:]:
+                pairs, _bin = SJ._colocated(well, a, b, None)
+                res = _loo_pairs(pairs)
+                res.update({'source': 'strata_join', 'well': well,
+                            'given': a, 'target': b})
+                rows.append(res)
+    # 2) KTB site-native pairs (rho -> Vp)
+    try:
+        pairs, washouts = _site_native_pairs('ktb_hb_complog_6020_excerpt')
+        res = _loo_pairs(pairs)
+        res.update({'source': 'site_pairs', 'well': 'ktb_complog',
+                    'given': 'rho (g/cc)', 'target': 'Vp (m/s)',
+                    'washouts_excluded': washouts})
+        rows.append(res)
+    except KeyError:
+        pass
+    ok = [r for r in rows if r['status'] == 'OK']
+    refused = [r for r in rows if r['status'] != 'OK']
+    ok.sort(key=lambda r: r['mae_pct'])
+    cov = [r['coverage'] for r in ok]
+    return {'ok': ok, 'refused': refused, 'n_ok': len(ok),
+            'n_refused': len(refused),
+            'best_mae_pct': (ok[0]['mae_pct'] if ok else None),
+            'worst_mae_pct': (ok[-1]['mae_pct'] if ok else None),
+            'median_coverage': (statistics.median(cov) if cov else None),
+            'doctrine': ('leave-one-out over the same machinery clients get; '
+                         'coverage targets ~0.68 (honest 1-sigma); refusals '
+                         'listed, never hidden; regenerated live so the '
+                         'accuracy table can never be a stale snapshot')}

gea/case_study.py

@@ -0,0 +1,229 @@
+# This Source Code Form is subject to the terms of the Mozilla Public
+# License, v. 2.0. If a copy of the MPL was not distributed with this
+# file, You can obtain one at https://mozilla.org/MPL/2.0/.
+"""case_study — depth-sweep case-study mode (v1.4.0 extension).
+
+Sweeps a well from top to TD and reports, at every depth, both drift legs
+(conventional vs GEA-stabilized), the separation in psi/yr, and the
+service-life arithmetic — then identifies WHERE the GEA advantage is
+largest (the deep hot interval past the stress knees, which is exactly where
+permanent-gauge drift matters most and gauges are hardest to replace).
+
+`write_markdown()` renders the sweep as the one-page case a customer would
+read: well summary, depth table, headline numbers, the bench-test claim, and
+the honest DERIVED_HYBRID classification note. Also runnable as a CLI:
+
+    python -m gea.case_study --td 25000 --out case.md
+
+Headless-safe: numpy only, no display imports.
+"""
+
+from __future__ import annotations
+
+from dataclasses import dataclass
+from datetime import date
+from pathlib import Path
+from typing import List, Optional
+
+import numpy as np
+
+from .downhole_engine import WellProfile, load_well_profile_csv, DEFAULT_TD_FT
+from .quartz_hpht_extension import (
+    calculate_quartz_transducer_hpht_program,
+    canonical_suppression,
+    conventional_drift,
+)
+
+
+@dataclass
+class CaseStudyConfig:
+    td_ft: float = DEFAULT_TD_FT
+    n_depth_points: int = 12
+    start_ft: float = 1000.0
+    surface_temp_F: float = 75.0                    # anchor: surface ambient (template)
+    surface_pressure_psi: float = 14.7              # anchor: 1 atm
+    temp_gradient_F_per_ft: float = 0.018           # anchor: geothermal gradient (template)
+    pressure_gradient_psi_per_ft: float = 0.465     # anchor: hydrostatic gradient (industry)
+    profile: Optional[WellProfile] = None           # real well profile overrides gradients
+    full_scale_psi: float = 30000.0                 # anchor: HPHT quartz-gauge FS class
+    error_budget_pct_fs: float = 0.5                # anchor: typical spec total-error budget
+    horizon_years: float = 5.0
+    k_structural_trim: float = 1.0
+    phi_coupling_trim: float = 1.0
+    gauge_spec: object = None                       # GaugeSpec (v1.5.0); None = template anchors
+    deviation: object = None                        # DeviationSurvey (v1.6.0): sweep depths are MD, physics at TVD
+    well_name: str = "case-study well"
+
+
+def depth_sweep(cfg: CaseStudyConfig | None = None) -> List[dict]:
+    """Per-depth twin-leg drift table from surface interval to TD."""
+    cfg = cfg or CaseStudyConfig()
+    fs = float(cfg.gauge_spec.full_scale_psi) if cfg.gauge_spec is not None else cfg.full_scale_psi
+    knee_C = float(cfg.gauge_spec.thermal_knee_C) if cfg.gauge_spec is not None else 150.0
+    knee_psi = float(cfg.gauge_spec.pressure_knee_psi) if cfg.gauge_spec is not None else 15000.0
+    depths = np.linspace(cfg.start_ft, cfg.td_ft * 0.985, cfg.n_depth_points)
+    rows = []
+    for d in depths:
+        d = float(d)                     # keep every reported number a plain float; d is MD
+        tvd = float(cfg.deviation.tvd_of(d)) if cfg.deviation is not None else d
+        if cfg.profile is not None:
+            p_psi, t_F = cfg.profile.interp(tvd)     # profiles are TVD-indexed
+        else:
+            p_psi = cfg.surface_pressure_psi + tvd * cfg.pressure_gradient_psi_per_ft
+            t_F = cfg.surface_temp_F + tvd * cfg.temp_gradient_F_per_ft
+        t_C = (t_F - 32.0) * 5.0 / 9.0
+        r = calculate_quartz_transducer_hpht_program(
+            depth_m=float(d) / 3.28084, temp_c=t_C, pressure_psi=p_psi,
+            k_structural_trim=cfg.k_structural_trim,
+            phi_coupling_trim=cfg.phi_coupling_trim,
+            spec=cfg.gauge_spec)
+        uq = r["value"]["drift_pct"]
+        cv = conventional_drift(t_C, p_psi, spec=cfg.gauge_spec)
+        sep = cv - uq
+        rows.append({
+            'depth_ft': round(float(d), 0),
+            'tvd_ft': round(tvd, 0),
+            'temp_F': round(t_F, 1),
+            'temp_C': round(t_C, 1),
+            'pressure_psi': round(p_psi, 0),
+            'past_thermal_knee': bool(t_C > knee_C),
+            'past_pressure_knee': bool(p_psi > knee_psi),
+            'conv_drift_pct_fs_yr': round(cv, 4),
+            'drift_pct_fs_yr': round(uq, 4),
+            'separation_pct_fs_yr': round(sep, 5),
+            'separation_psi_yr': round(sep * fs / 100.0, 2),
+            'horizon_separation_psi': round(sep * fs / 100.0 * cfg.horizon_years, 1),
+            'extra_service_life_yr': round(
+                cfg.error_budget_pct_fs / uq - cfg.error_budget_pct_fs / cv, 3) if uq > 0 else None,
+        })
+    return rows
+
+
+def case_study(cfg: CaseStudyConfig | None = None) -> dict:
+    """The sweep plus headline metrics: where the advantage is largest."""
+    cfg = cfg or CaseStudyConfig()
+    rows = depth_sweep(cfg)
+    best = max(rows, key=lambda r: r['separation_psi_yr'])
+    knee_rows = [r for r in rows if r['past_thermal_knee'] or r['past_pressure_knee']]
+    return {
+        'config': {
+            'well_name': cfg.well_name,
+            'td_ft': cfg.td_ft,
+            'profile': cfg.profile.name if cfg.profile else 'linear gradients',
+            'full_scale_psi': (float(cfg.gauge_spec.full_scale_psi)
+                               if cfg.gauge_spec is not None else cfg.full_scale_psi),
+            'gauge_spec': cfg.gauge_spec.name if cfg.gauge_spec is not None else 'template_generic (default)',
+            'gauge_spec_source': cfg.gauge_spec.source if cfg.gauge_spec is not None else '22Aug2026 template anchors',
+            'horizon_years': cfg.horizon_years,
+        },
+        'canonical_suppression': round(canonical_suppression(
+            cfg.k_structural_trim, cfg.phi_coupling_trim), 4),
+        'rows': rows,
+        'headline': {
+            'best_depth_ft': best['depth_ft'],
+            'best_separation_psi_yr': best['separation_psi_yr'],
+            'best_horizon_separation_psi': best['horizon_separation_psi'],
+            'td_conv_drift_pct_fs_yr': rows[-1]['conv_drift_pct_fs_yr'],
+            'td_program_drift_pct_fs_yr': rows[-1]['drift_pct_fs_yr'],
+            'hpht_interval_top_ft': knee_rows[0]['depth_ft'] if knee_rows else None,
+            'advantage_grows_with_depth': bool(
+                rows[-1]['separation_psi_yr'] >= rows[0]['separation_psi_yr']),
+        },
+    }
+
+
+def write_markdown(result: dict | None = None, path: str | None = None,
+                   cfg: CaseStudyConfig | None = None) -> Path:
+    """Render the one-page customer-facing case study as markdown."""
+    if result is None:
+        result = case_study(cfg)
+    if path is None:
+        path = "downhole_case_study.md"
+    c = result['config']
+    h = result['headline']
+    sup = result['canonical_suppression']
+    lines = [
+        f"# GEA Downhole Drift Case Study — {c['well_name']}",
+        "",
+        f"*Generated {date.today().isoformat()} by gea "
+        f"(gea-program). Basis: {c['profile']}; TD {c['td_ft']:,.0f} ft; "
+        f"full scale {c['full_scale_psi']:,.0f} psi; horizon {c['horizon_years']:g} yr.*",
+        "",
+        f"*Gauge spec: **{c['gauge_spec']}** — {c['gauge_spec_source']}*",
+        "",
+        "## The claim",
+        "",
+        f"Quartz-gauge drift ({c['gauge_spec']} baseline) divided by the GEA",
+        f"canonical suppression composition = **{sup}** at the locked primitives",
+        "(F_TRZ = 0.1, K_MEX = 25/12, Phi_res = 0.84) — drift **below** the",
+        "conventional gauge at every depth, with the largest advantage in the",
+        "deep hot interval where gauges are hardest to replace.",
+        "",
+        "## Depth sweep",
+        "",
+        "| Depth (ft) | T (degF) | P (psi) | HPHT | Conv (%FS/yr) | GEA (%FS/yr) | Sep (psi/yr) | Extra life (yr) |",
+        "|---:|---:|---:|:--:|---:|---:|---:|---:|",
+    ]
+    for r in result['rows']:
+        knee = "YES" if (r['past_thermal_knee'] or r['past_pressure_knee']) else "-"
+        lines.append(
+            f"| {r['depth_ft']:,.0f} | {r['temp_F']} | {r['pressure_psi']:,.0f} | {knee} "
+            f"| {r['conv_drift_pct_fs_yr']} | {r['drift_pct_fs_yr']} "
+            f"| {r['separation_psi_yr']} | {r['extra_service_life_yr']} |")
+    lines += [
+        "",
+        "## Headlines",
+        "",
+        f"- Largest advantage at **{h['best_depth_ft']:,.0f} ft**: "
+        f"**{h['best_separation_psi_yr']} psi/yr** separation, "
+        f"**{h['best_horizon_separation_psi']} psi** over the {c['horizon_years']:g}-yr horizon.",
+        f"- At TD: conventional {h['td_conv_drift_pct_fs_yr']} vs GEA "
+        f"{h['td_program_drift_pct_fs_yr']} %FS/yr.",
+        (f"- HPHT interval (past the 150 degC / 15,000 psi knees) begins near "
+         f"**{h['hpht_interval_top_ft']:,.0f} ft** — the advantage compounds there."
+         if h['hpht_interval_top_ft'] else
+         "- This well stays below the HPHT stress knees; the advantage is the flat suppression ratio."),
+        "- The advantage " + ("**grows with depth**." if h['advantage_grows_with_depth']
+                              else "is approximately flat across depth."),
+        "",
+        "## How to test it",
+        "",
+        "Twin-gauge bench test (PAPER_2250 LABORATORY tier / PAPER_2256 sec 5): one",
+        "conventional quartz gauge, one under GEA SCm-resonance conditioning, at",
+        f"matched T/P. Predicted drift ratio: **{sup}**. The simulator's comparison",
+        "mode and service-life divergence curves supply the reference data streams.",
+        "",
+        "## Classification (honest basis)",
+        "",
+        "DERIVED_HYBRID (PAPER_2149): industry-observed anchors (0.215 %FS/yr",
+        "baseline; 150 degC / 15,000 psi stress knees; gradients; FS class; spec",
+        "budget) x canonical-GEA suppression from locked primitives. The dressing",
+        "coefficients are the template's engineering fit, disclosed — not claimed as",
+        "derivations. Record: PAPER_2256. (c) Daniel T. Murphy / GEA Research",
+        "Program, AGPL-3.0 + Commercial.",
+        "",
+    ]
+    p = Path(path)
+    p.write_text("\n".join(lines), encoding="utf-8")
+    return p
+
+
+def _main() -> None:
+    import argparse
+    ap = argparse.ArgumentParser(description="GEA downhole depth-sweep case study")
+    ap.add_argument("--td", type=float, default=DEFAULT_TD_FT, help="total depth, ft")
+    ap.add_argument("--points", type=int, default=12, help="depth points in the sweep")
+    ap.add_argument("--profile", type=str, default=None, help="well profile CSV (depth_ft,pressure_psi,temp_F)")
+    ap.add_argument("--name", type=str, default="case-study well")
+    ap.add_argument("--horizon", type=float, default=5.0, help="service horizon, years")
+    ap.add_argument("--out", type=str, default="downhole_case_study.md")
+    a = ap.parse_args()
+    prof = load_well_profile_csv(a.profile) if a.profile else None
+    cfg = CaseStudyConfig(td_ft=a.td, n_depth_points=a.points, profile=prof,
+                          well_name=a.name, horizon_years=a.horizon)
+    out = write_markdown(case_study(cfg), a.out)
+    print(f"case study written: {out}")
+
+
+if __name__ == "__main__":
+    _main()

gea/catalog/acex_lomonosov_age_depth_model.provenance.json

@@ -0,0 +1,23 @@
+{
+  "entry": "acex_lomonosov_age_depth_model",
+  "kind": "REAL_AGE_DEPTH_MODEL",
+  "well": "IODP Exp 302 Composite Site (ACEX), LOMONOSOV RIDGE, CENTRAL ARCTIC OCEAN (87.8900 N, 137.6500 E - 235 km from the North Pole, the closest any scientific drilling data comes to it; water depth -1,251 m) - Arctic Coring Expedition (2004), drilled from the icebreaker Vidar Viking with two escort icebreakers managing the moving sea ice; the composite record splices Holes M0002A/M0003A/M0004A/M0004C into one virtual core",
+  "operator_on_log": "Backman, Jakobsson, Frank, O'Regan, Moran et al. (14 authors) - the canonical ACEX age model, Paleoceanography 23 PA1S03 (doi:10.1029/2007PA001476); control points from 10Be profiles, biostratigraphy (P. clithridium, Azolla, A. augustum) and magnetostratigraphy (chron C25n)",
+  "source_database": "PANGAEA, dataset doi:10.1594/PANGAEA.705517 ('(Table 1) Age-Depth control points and sedimentation rates of IODP Exp302 sediments')",
+  "source_url": "https://doi.pangaea.de/10.1594/PANGAEA.705517?format=textfile",
+  "license": "Creative Commons Attribution 3.0 Unported (CC-BY-3.0), declared IN the file's own header. Attribution: Backman et al. / PANGAEA / IODP",
+  "fetch_date": "2026-08-26",
+  "fetch_method": "read-only HTTPS fetch of the PANGAEA textfile export - COMPLETE dataset (38 data points = every non-empty cell across the 5 counted parameters on 10 rows; depth is GEOCODE) in one fetch",
+  "coverage": "COMPLETE DATASET: the age-depth backbone of the first-ever scientific drilling of the central Arctic Ocean - 10 control points (A-G plus hiatus bounds), 0-399.63 mcd spanning 0-55,904 ka (56 million years of Arctic history), with per-interval sedimentation rates 0.80-2.43 cm/ka.",
+  "curves": [
+    "Depth sed [mcd] (index)",
+    "Age model [ka]",
+    "SR sedimentation rate [cm/ka]",
+    "Sample comment = control-point letters A-G (quality)",
+    "Comm + Interval comm (verbatim in file: 10Be, biostratigraphy and chron C25n tie-point descriptions)"
+  ],
+  "hiatuses_in_the_data": "BOTH ACEX hiatuses are encoded in the table's structure: (1) the 26.2-Myr mid-Cenozoic hiatus appears as a DUPLICATE-DEPTH PAIR - 198.70 m carries two ages, 18,200 and 44,400 ka: the same centimetre of seafloor is both early Miocene and middle Eocene; (2) the 2.2-Myr middle-late Miocene hiatus sits between 135.49 m (9,360 ka) and 140.44 m (11,560 ka).",
+  "source_anomaly_disclosure": "The source header's own Comment field is TRUNCATED IN THE ARCHIVE - it ends mid-sentence at 'Paleogene: (404.8-198.7) m/(56.2-44.4) Ma. Average sedimentation' with the Paleogene average rate cut off (re-derivable: 206.1 m / 11.8 Ma = 17.5 m/Myr). Preserved verbatim, including the cut. Note also the header's 404.8 m full-section base vs the table's 399.63 m deepest control point (chron C25n) - consistent, not contradictory.",
+  "purpose": "thirtieth catalogue entry - a MILESTONE ENTRY: the CENTRAL ARCTIC OCEAN joins as region EIGHTEEN at 87.89 N (235 km from the North Pole - the catalogue now literally reaches from pole to pole through 30 entries), first icebreaker-drilled data, first composite virtual core, and AGE-DEPTH GEOCHRONOLOGY is kind NINETEEN - the catalogue's first time axis: every earlier entry measures the state of the subsurface, this one dates it",
+  "verified": "10 rows / 38 counted cells against the header Size declaration; ALL six sedimentation rates re-derived live from the control points (five interval rates as delta-depth/delta-age between the correct hiatus-aware bounds + the 135.49 m rate from the surface); header-average identity 198.7/16.0 = 12.4 m/Myr re-derived; duplicate-depth hiatus pair pinned; endpoints; gate-pinned; source URL serves the identical dataset"
+}

gea/catalog/acex_lomonosov_age_depth_model.txt

@@ -0,0 +1,28 @@
+/* DATA DESCRIPTION:
+Citation:	Backman, Jan; Jakobsson, Martin; Frank, Martin; Sangiorgi, Francesca; Brinkhuis, Henk; Stickley, Catherine E; O'Regan, Matthew; Lovlie, Reidar; Palike, Heiko; Spofforth, David J A; Gattacceca, Jerome; Moran, Kathryn; King, John W; Heil, Chip (2008): (Table 1) Age-Depth control points and sedimentation rates of IODP Exp302 sediments [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.705517, 
+	In supplement to: Backman, J et al. (2008): Age model and core-seismic integration for the Cenozoic Arctic Coring Expedition sediments from the Lomonosov Ridge. Paleoceanography, 23(1), PA1S03, https://doi.org/10.1029/2007PA001476
+Project(s):	Integrated Ocean Drilling Program / International Ocean Discovery Program (IODP) (URI: https://www.iodp.org/)
+Coverage:	LATITUDE: 87.890000 * LONGITUDE: 137.650000
+	MINIMUM DEPTH, sediment/rock: 0.00 m * MAXIMUM DEPTH, sediment/rock: 399.63 m
+Event(s):	302-CompSite (ACEX) * LATITUDE: 87.890000 * LONGITUDE: 137.650000 * ELEVATION: -1251.0 m * CAMPAIGN: Exp302 (Arctic Coring Expedition, ACEX) * BASIS: CCGS Captain Molly Kool (Vidar Viking) (URI: https://en.wikipedia.org/wiki/CCGS_Captain_Molly_Kool) * METHOD/DEVICE: Composite Core (COMPCORE) * COMMENT: Virtual core composed of spliced cores from Holes M0002A, M0003A, M0004A and M0004C. Virtual position and water depth is calculated as the mean of all Exp302 sites.
+Comment:	Depth is given in meter composite depth (mcd) ! Neogene: 198.7 m/(18.2-2.2) Ma; total time interval minus 2.2 Ma long hiatus. Average sedimentation rate: 12.4 m/MA; Paleogene: (404.8-198.7) m/(56.2-44.4) Ma. Average sedimentation
+Parameter(s):	DEPTH, sediment/rock [m] (Depth sed) * GEOCODE
+	Age model [ka] (Age model) * PI: Backman, Jan (backman@geo.su.se) * METHOD/DEVICE: Age model, composite
+	Sample comment (Sample comment) * PI: Backman, Jan (backman@geo.su.se) * COMMENT: Depth control points
+	Age, comment (Comm) * PI: Backman, Jan (backman@geo.su.se)
+	Interval comments (Interval comm) * PI: Backman, Jan (backman@geo.su.se) * COMMENT: Interval between depth control points
+	Sedimentation rate [cm/ka] (SR) * PI: Backman, Jan (backman@geo.su.se) * METHOD/DEVICE: Calculated, see reference(s)
+License:	Creative Commons Attribution 3.0 Unported (CC-BY-3.0) (URI: https://creativecommons.org/licenses/by/3.0/)
+Size:	38 data points
+*/
+Depth sed [m]	Age model [ka]	Sample comment	Comm	Interval comm	SR [cm/ka]
+0.00	0	A	top of section		
+135.49	9360		deepest 10Be sample above middle-late Miocene hiatus		1.45
+140.44	11560		next deeper 10Be sample, middle-late Miocene hiatus		
+151.28	12310	B	deepest 10Be sample	A-B	1.45
+198.70	18200	C	midpoint Burdigalian stage equals young end of Cenozoic hiatus	B-C	0.80
+198.70	44400		extrapolate up from D, to old end of Cenozoic hiatus		
+202.95	44600	D	last abundant occurrence of P. clithridium	D-E	2.43
+299.95	48600	E	last occurrence of Azolla spp.	E-F	1.27
+381.42	55000	F	last occurrence of A. augustum	F-G	2.01
+399.63	55904	G	top chron C25n		

gea/catalog/agassiz77_canada_temperature.csv

@@ -0,0 +1,68 @@
+d,t
+10.91, -24.16
+15.91, -24.35
+20.91, -24.33
+25.91, -24.19
+30.91, -24.19
+35.91, -24.08
+40.91, -24.03
+45.91, -23.92
+50.91, -23.85
+55.91, -23.80
+60.91, -23.72
+65.91, -23.68
+70.91, -23.61
+75.91, -23.54
+80.91, -23.47
+85.91, -23.39
+90.91, -23.32
+95.91, -23.23
+100.91, -23.14
+105.91, -23.05
+110.91, -22.97
+115.91, -22.85
+120.91, -22.76
+125.91, -22.68
+130.91, -22.56
+135.91, -22.44
+140.91, -22.33
+145.91, -22.24
+150.91, -22.11
+155.91, -21.98
+160.91, -21.86
+165.91, -21.76
+170.91, -21.61
+175.91, -21.49
+180.91, -21.35
+185.91, -21.22
+190.91, -21.09
+195.91, -20.95
+200.91, -20.83
+205.91, -20.69
+210.91, -20.55
+215.91, -20.42
+220.91, -20.27
+225.91, -20.13
+230.91, -19.99
+235.91, -19.82
+240.91, -19.69
+245.91, -19.55
+250.91, -19.39
+255.91, -19.24
+260.91, -19.09
+265.91, -18.94
+270.91, -18.80
+275.91, -18.64
+280.91, -18.50
+285.91, -18.42
+290.91, -18.17
+295.91, -17.90
+300.91, -17.84
+305.91, -17.69
+310.91, -17.54
+315.91, -17.37
+320.91, -17.22
+325.91, -17.06
+330.91, -16.92
+335.91, -16.83
+340.91, -16.74

gea/catalog/agassiz77_canada_temperature.provenance.json

@@ -0,0 +1,17 @@
+{
+  "entry": "agassiz77_canada_temperature",
+  "kind": "REAL_TEMPERATURE_PROFILE_COMPLETE",
+  "well": "Agassiz77 borehole, Agassiz Ice Cap, Ellesmere Island, Canadian Arctic North (approx 80.7 N, -73.1 E per database metadata); drilled/measured 1977",
+  "operator_on_log": "Canadian ice-core program lineage (Fisher/Clarke/Waddington; Zdanowicz-provided data per database metadata)",
+  "source_database": "Greenland and Canadian Arctic ice temperature profiles database (GEUS): Lokkegaard et al. 2023, The Cryosphere 17, 3829-3845, doi:10.5194/tc-17-3829-2023; data DOI 10.22008/FK2/3BVF9V (GEUS Dataverse); science lineage: Clarke, Fisher & Waddington 1987 (IAHS vol. 170, 169-180); redistributed in the public GEUS GitHub repository (boreholes/Agassiz77/data.csv)",
+  "source_url": "https://raw.githubusercontent.com/GEUS-Glaciology-and-Climate/greenland_ice_borehole_temperature_profiles/main/boreholes/Agassiz77/data.csv",
+  "license": "open science database (The Cryosphere CC-BY publication; GEUS Dataverse); attribution to Lokkegaard et al. 2023 and the Canadian Agassiz program lineage",
+  "fetch_date": "2026-08-24",
+  "fetch_method": "read-only HTTPS fetch (text/plain) - COMPLETE file",
+  "coverage": "COMPLETE FILE: all 67 stations, 10.91-340.91 m at 5 m spacing, verbatim (database-stated coverage 97% of the 341 m thickness). Database metadata notes: ice thickness in data does not match Vinther (2008) and the location is approximate - carried honestly as the database's own caveat.",
+  "curves": ["d (depth, m)", "t (temperature, degC) - MEASURED TEMPERATURE, 5 m stations"],
+  "temperature_data": "Second continuous MEASURED temperature profile: 67 stations, -24.16 degC near surface warming monotonically below ~16 m to -16.74 degC at 340.91 m - a classic High Arctic ice-cap thermal profile (no deep glacial minimum: the cap is thin relative to GISP2).",
+  "dialect_notes": "GEUS d,t CSV with space-after-comma values (parser handles); the catalogue's second temperature-CSV entry.",
+  "purpose": "seventh catalogue entry - SEVENTH region (CANADA, Ellesmere Island); second MEASURED_CURVES-tier well; 1977 vintage (the catalogue's oldest measurement)",
+  "verified": "67 stations; endpoints verbatim; MEASURED_CURVES derivation; gate-pinned"
+}

gea/catalog/barbados_110_consolidation.provenance.json

@@ -0,0 +1,24 @@
+{
+  "entry": "barbados_110_consolidation",
+  "kind": "REAL_CONSOLIDATION_GEOTECHNICS",
+  "well": "ODP Holes 110-671A/B/C/D (the historic first-ever drilling THROUGH an active decollement, 691.2 m penetration), 110-672A (incoming-plate reference), 110-673A/B and 110-676A - BARBADOS RIDGE accretionary complex, Lesser Antilles, North Atlantic (15.52-15.54 N, 58.64-58.81 W; seafloors -4,678 to -5,070 m, the catalogue's deepest Atlantic water); ODP Leg 110 (1986), samples 4.40-496.70 mbsf reaching the sediments immediately above the incipient decollement",
+  "operator_on_log": "Taylor & Leonard laboratory consolidation testing (oedometer + triaxial, vertical AND horizontal orientations) on whole-round core samples; drilling by JOIDES Resolution - the leg that founded modern decollement science and led directly to the first subduction-zone CORKs (Leg 156)",
+  "source_database": "PANGAEA, dataset doi:10.1594/PANGAEA.747753 ('Geotechnical and physical properties from whole-round consolidation samples of ODP Leg 110 holes (Table 1)'); supplement to Taylor & Leonard 1990, Proc. ODP Sci. Results 110, 289-308 (doi:10.2973/odp.proc.sr.110.152.1990); full paper abstract preserved in the file's own header",
+  "source_url": "https://doi.pangaea.de/10.1594/PANGAEA.747753?format=textfile",
+  "license": "Creative Commons Attribution 3.0 Unported (CC-BY-3.0), declared IN the file's own header. Attribution: Taylor/Leonard / PANGAEA / ODP",
+  "fetch_date": "2026-08-27",
+  "fetch_method": "read-only HTTPS fetch of the PANGAEA textfile export - COMPLETE dataset (874 data points per the header's own Size declaration, re-counted at gate time as non-empty body cells excluding Event + geocode columns, EXACT) in one fetch",
+  "coverage": "COMPLETE DATASET: 41 consolidation-test rows across 8 holes and 4 sites - wet-bulk/grain density, water content, porosity, void ratio, effective overburden stress 20-3,808 kPa, preconsolidation pressure (max/min), the Pc-Po state columns, permeability at Pc (1.2e-8 to 7.6e-4 cm/s), compression index, vane shear strength, grain-size fractions, shipboard CaCO3, and full Atterberg limits; test orientations (vertical/horizontal/triaxial) ride verbatim in three comment columns",
+  "curves": [
+    "Depth sed [m] 4.40-496.70 (GEOCODE)",
+    "WBD, grain density, water content, porosity, void ratio",
+    "Overb press [kPa] 20-3,808 + Pc max/min [kPa]",
+    "Pc/Po (max)/(min) - ACTUALLY Pc MINUS Po in kPa, see anomaly disclosure",
+    "Coeff of permeability at Pc [cm/s] + compression index Cc + vane shear [kPa]",
+    "Sand/Silt/<2um fractions + CaCO3 + Atterberg (LL/PL/PI)",
+    "Event rides as per-row quality; lithology/epoch/orientation comments verbatim"
+  ],
+  "source_anomaly_disclosure": "SIX source quirks preserved verbatim, each proven by arithmetic where provable: (1) the two columns LABELED 'Overconsolidation ratio (Pc/Po)' actually hold the DIFFERENCE Pc - Po in kPa - re-derived on all 66 populated pairs, 65 close within 1 kPa (including 3000 - 3808 = -808 at the deepest sample) and the single violator (672A-2H-3 horizontal min: 100 - 37 = 63 vs archived 43) is a genuine 20-kPa arithmetic slip in the source table, carried unchanged; (2) sample 671C-1X-1 archives a permeability of literally 0.000000000 cm/s - a value that underflowed the table's 9-decimal format (physically nonzero, numerically zero in the archive); (3) all TEN 672A rows with Atterberg data have the Plastic-limit column DUPLICATING the Liquid limit (e.g. 87.4/87.4 with PI 56), breaking PI = LL - PL exactly where every other site closes it (worst deviation 3.6 on 23 non-672A rows); (4) epoch 'early Plocene' typo (672A-11H-3); (5) sample label '110-673A2H-3-2' missing its hyphen; (6) 'Water content, wet mass' exceeds 100% on five rows (103-120%) - impossible on a wet-mass basis, so those values are dry-mass water contents under a wet-mass label",
+  "purpose": "thirty-sixth catalogue entry - the BARBADOS RIDGE / LESSER ANTILLES joins as region TWENTY-FOUR (the birthplace of decollement drilling and of the CORK) and consolidation geotechnics is kind TWENTY-FIVE: preconsolidation state, permeability, compressibility and plasticity - the quantities that turn the catalogue's measured pore pressures (U1324) and strengths into a MECHANISM; the paper's 'lubricated subduction' headline is re-derived by the gate on every run: every sample at or below 90 mbsf is UNDERCONSOLIDATED (Pc - Po < 0, reaching -916 kPa), the pore fluid carrying what the grain framework cannot - and at the deepest sample above the decollement Pc/Po = 3000/3808 = 0.79, the excess-pressure state approaching overburden that lets an ocean plate slide",
+  "verified": "41 rows / 874 counted cells against the header's own Size declaration (EXACT); the Pc-Po difference identity re-derived on all 66 populated pairs (65 within 1 kPa, the single 20-kPa source slip disclosed and pinned); ALL 13 samples >= 90 mbsf underconsolidated with minimum Pc-Po = -916 kPa; 672A PL==LL duplication pinned on all 10 rows vs PI=LL-PL closing within 3.6 on all 23 other-site rows; the zero-permeability underflow cell pinned; epoch typo + missing-hyphen label pinned; gate-pinned; source URL serves the identical dataset"
+}

gea/catalog/barbados_110_consolidation.txt

@@ -0,0 +1,91 @@
+/* DATA DESCRIPTION:
+Citation:	Taylor, Elliott; Leonard, John (1990): Geotechnical and physical properties from whole-round consolidations samples of ODP Leg 110 holes (Table 1) [dataset]. PANGAEA, https://doi.org/10.1594/PANGAEA.747753, 
+	Supplement to: Taylor, E; Leonard, J (1990): Sediment consolidation and permeability at the Barbados Forearc. In: Moore, JC; Mascle, A; et al. (eds.), Proceedings of the Ocean Drilling Program, Scientific Results, College Station, TX (Ocean Drilling Program), 110, 289-308, https://doi.org/10.2973/odp.proc.sr.110.152.1990
+Abstract:	Sediments cored from the Barbados Ridge accretionary complex at ODP Sites 671, 672, and 676 are progressively underconsolidated with depth. Very low measured intergranular permeabilities contribute substantially to the probable excess pore-water pressures that account for the underconsolidated behavior. The estimated pore pressures immediately overlying the incipient decollement at the reference site are nearly equal to the overburden stress. Data show that the pore fluids within the fine-grained sediments of the complex are absorbing the tectonic shock of underthrusting as they are accreted, effectively leading to "lubricated subduction" of the underthrust Oligocene and older sediments. Lateral compression and dewatering by fracture permeability is taking place in the accreted section despite apparent underconsolidation, resulting in limited pore-volume reduction and sediment strengthening from east to west on the ridge.
+Project(s):	Ocean Drilling Program (ODP) (URI: https://www-odp.tamu.edu:443/)
+Coverage:	MEDIAN LATITUDE: 15.529850 * MEDIAN LONGITUDE: -58.715613 * SOUTH-BOUND LATITUDE: 15.524700 * WEST-BOUND LONGITUDE: -58.810000 * NORTH-BOUND LATITUDE: 15.540000 * EAST-BOUND LONGITUDE: -58.641000
+	DATE/TIME START: 1986-06-27T08:15:00 * DATE/TIME END: 1986-08-14T00:30:00
+	MINIMUM DEPTH, sediment/rock: 4.40 m * MAXIMUM DEPTH, sediment/rock: 496.70 m
+Event(s):	110-671A * LATITUDE: 15.525800 * LONGITUDE: -58.732500 * DATE/TIME START: 1986-06-27T08:15:00 * DATE/TIME END: 1986-06-28T12:45:00 * ELEVATION: -4942.0 m * PENETRATION: 5 m * RECOVERY: 5.02 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 1 core; 5 m cored; 0 m drilled; 100.4 % recovery
+	110-671B * LATITUDE: 15.525800 * LONGITUDE: -58.732500 * DATE/TIME START: 1986-06-28T12:45:00 * DATE/TIME END: 1986-07-10T23:15:00 * ELEVATION: -4942.0 m * PENETRATION: 691.2 m * RECOVERY: 553.13 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 74 cores; 691.2 m cored; 0 m drilled; 80 % recovery
+	110-671C * LATITUDE: 15.528000 * LONGITUDE: -58.732800 * DATE/TIME START: 1986-07-10T23:15:00 * DATE/TIME END: 1986-07-15T07:30:00 * ELEVATION: -4947.0 m * PENETRATION: 514.7 m * RECOVERY: 16.66 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 2 cores; 19 m cored; 0 m drilled; 87.7 % recovery
+	110-671D * LATITUDE: 15.524700 * LONGITUDE: -58.731300 * DATE/TIME START: 1986-08-06T23:30:00 * DATE/TIME END: 1986-08-10T07:15:00 * ELEVATION: -4953.0 m * PENETRATION: 461 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 1 core; 9 m cored; 0 m drilled; 0 % recovery
+	110-672A * LATITUDE: 15.540000 * LONGITUDE: -58.641000 * DATE/TIME START: 1986-07-15T11:10:00 * DATE/TIME END: 1986-07-22T12:15:00 * ELEVATION: -4984.0 m * PENETRATION: 493.8 m * RECOVERY: 382.62 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 53 cores; 493.8 m cored; 0 m drilled; 77.5 % recovery
+	110-673A * LATITUDE: 15.531700 * LONGITUDE: -58.810000 * DATE/TIME START: 1986-07-22T18:30:00 * DATE/TIME END: 1986-07-23T23:00:00 * ELEVATION: -4678.0 m * PENETRATION: 36.4 m * RECOVERY: 33.33 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 4 cores; 36.4 m cored; 0 m drilled; 91.6 % recovery
+	110-673B * LATITUDE: 15.532000 * LONGITUDE: -58.641500 * DATE/TIME START: 1986-07-23T23:00:00 * DATE/TIME END: 1986-07-28T12:00:00 * ELEVATION: -4689.0 m * PENETRATION: 330.6 m * RECOVERY: 247.05 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 35 cores; 330.6 m cored; 0 m drilled; 74.7 % recovery
+	110-676A * LATITUDE: 15.530800 * LONGITUDE: -58.703300 * DATE/TIME START: 1986-08-10T11:15:00 * DATE/TIME END: 1986-08-14T00:30:00 * ELEVATION: -5070.0 m * PENETRATION: 310.2 m * RECOVERY: 262.09 m * LOCATION: North Atlantic Ocean * CAMPAIGN: Leg110 * BASIS: Joides Resolution (URI: https://www-odp.tamu.edu:443/resolutn.html) * METHOD/DEVICE: Drilling/drill rig (DRILL) * COMMENT: 33 cores; 310.2 m cored; 0 m drilled; 84.5 % recovery
+Parameter(s):	Event label (Event)
+	Sample code/label (Sample label) * PI: Taylor, Elliott * METHOD/DEVICE: DSDP/ODP/IODP sample designation (URI: hdl:10013/epic.27914.d001)
+	Sample comment (Sample comment) * PI: Taylor, Elliott
+	DEPTH, sediment/rock [m] (Depth sed) * GEOCODE
+	Lithology/composition/facies (Lithology) * PI: Taylor, Elliott
+	Epoch (Epoch) * PI: Taylor, Elliott
+	Density, wet bulk [g/cm**3] (WBD) * PI: Taylor, Elliott
+	Water content, wet mass [%] (Water wm) * PI: Taylor, Elliott
+	Density, grain [g/cm**3] (Density grain) * PI: Taylor, Elliott
+	Porosity [% vol] (Poros) * PI: Taylor, Elliott
+	Void Ratio description (Void Ratio) * PI: Taylor, Elliott * COMMENT: initial
+	Overburden pressure [kPa] (Overb press) * PI: Taylor, Elliott * COMMENT: effective overburden stress
+	Preconsolidation pressure [kPa] (Pc) * PI: Taylor, Elliott * COMMENT: max
+	Preconsolidation pressure [kPa] (Pc) * PI: Taylor, Elliott * COMMENT: min
+	Comment (Comment) * PI: Taylor, Elliott * COMMENT: PC max/min
+	Overconsolidation ratio (Pc/Po) * PI: Taylor, Elliott * COMMENT: max
+	Comment 2 (continued) (Comm 2) * PI: Taylor, Elliott * COMMENT: OCD max
+	Overconsolidation ratio (Pc/Po) * PI: Taylor, Elliott * COMMENT: min
+	Comment 3 (continued) (Comm 3) * PI: Taylor, Elliott * COMMENT: OCD min
+	Coefficient (Coeff) * PI: Taylor, Elliott * COMMENT: of permeability at Pc (cm/s)
+	Compression index (Cc) * PI: Taylor, Elliott
+	Comment 4 (continued) (Comm 4) * PI: Taylor, Elliott * COMMENT: Cc index
+	Shear strength, primary [kPa] (Shear str) * PI: Taylor, Elliott * COMMENT: vane
+	Sand [%] (Sand) * PI: Taylor, Elliott
+	Silt [%] (Silt) * PI: Taylor, Elliott
+	Size fraction < 0.002 mm, clay [%] (<2 µm) * PI: Taylor, Elliott
+	Calcium carbonate [%] (CaCO3) * PI: Taylor, Elliott * COMMENT: ship
+	Liquid limit (Liq limit) * PI: Taylor, Elliott * COMMENT: Atterberg limits
+	Plastic limit (Plastic limit) * PI: Taylor, Elliott * COMMENT: Atterberg limits
+	Plastic index (Plastic index) * PI: Taylor, Elliott * COMMENT: Atterberg limits
+License:	Creative Commons Attribution 3.0 Unported (CC-BY-3.0) (URI: https://creativecommons.org/licenses/by/3.0/)
+Size:	874 data points
+*/
+Event	Sample label	Sample comment	Depth sed [m]	Lithology	Epoch	WBD [g/cm**3]	Water wm [%]	Density grain [g/cm**3]	Poros [% vol]	Void Ratio (initial)	Overb press [kPa] (effective overburden stress)	Pc [kPa] (max)	Pc [kPa] (min)	Comment (PC max/min)	Pc/Po (max)	Comm 2 (OCD max)	Pc/Po (min)	Comm 3 (OCD min)	Coeff (of permeability at Pc (cm/s))	Cc	Comm 4 (Cc index)	Shear str [kPa] (vane)	Sand [%]	Silt [%]	<2 µm [%]	CaCO3 [%] (ship)	Liq limit (Atterberg limits)	Plastic limit (Atterberg limits)	Plastic index (Atterberg limits)
+110-671A	110-671A-1H-3		4.40	marl/chalk/Ls	early Pleistocene	1.50	111.0	2.75	74.0	2.85	20	45	35		25		15		0.000001060	0.75		21.50	6.6	32.2	61.2	10.00	128.0	42.1	86.0
+110-671B	110-671B-3H-3	triaxial	20.95	calcareous mudstone	early Pleistocene	1.52	83.0	2.88	66.8	2.39	105	140			35				0.000000828							38.70	155.0	74.0	81.0
+110-671B	110-671B-3H-3		21.30	calcareous mudstone	early Pleistocene	1.66	66.9	2.66	62.4	1.66	107	160	125		53		18		0.000000336	0.62		34.00	5.7	47.3	47.0	38.70	84.0	30.1	54.0
+110-671B	110-671B-3H-4		22.90	calcareous mudstone	early Pleistocene	1.59	81.8	2.91	70.2	2.32	115	50	40		-65		-75					28.82				38.70			
+110-671B	110-671B-7H-3	triaxial	59.36	marl	late Pliocene	1.73	48.0	2.94	58.9	1.42	346	166			-180							187.00				43.80	78.0	37.0	41.0
+110-671B	110-671B-7H-4	triaxial	60.83	marl	late Pliocene	1.70	53.0	2.93	61.9	1.55	357	166			-191											30.00	88.0	39.0	49.0
+110-671B	110-671B-9H-2	vertical	76.80	marl	early Pliocene	1.77	52.9	2.81	58.5	1.41	468	400	350	vertical	-68	vertical	-118		0.000000086	0.50	vertical	65.00	6.6	17.6	75.5	29.70	86.8	31.2	52.0
+110-671B	110-671B-9H-2	horizontal	76.90	marl	early Pliocene	1.77	52.9	2.81	58.5	1.41	469	750	720	horizontal	281	horizontal	251		0.000000467	0.54	horizontal	65.00	6.6	17.6	75.5	29.70	86.8	31.2	52.0
+110-671B	110-671B-10H-3	triaxial	87.79	marl	early Pliocene	1.72	49.0	2.86	61.4	1.4	547	240			-307							290.00				41.30	101.0	49.0	52.0
+110-671B	110-671B-10H-4	triaxial	89.33	marl	early Pliocene	1.71	54.0	2.85	58.8	1.26	558	240			-318							195.00				40.80	85.0	35.0	50.0
+110-671B	110-671B-15X-4		135.50	mudstone/ash	early Pleistocene	1.58	66.3	2.84	69.5	2.05	850	280	120		-570		-730		0.000000042	0.60						14.30			
+110-671B	110-671B-17X-2		151.50	calcareous mudstone	early Pleistocene	1.67	67.0	2.72	62.4	1.66	1003	600	530		-403		-473		0.000000095	0.72		44.00	3.3	42.8	53.7	1.00	175.0	40.9	134.0
+110-671B	110-671B-22X-2		199.10	calcareous mudstone	early Pleistocene	1.84	47.8	2.65	53.2	1.14	1365	520	480		-845		-885		0.000000012	0.31		72.00	0.7	29.1	70.2	19.90	135.0	32.2	103.0
+110-671D	110-671D-NAC-1		452.70	claystone	late Miocene?	1.65	69.3	2.52	60.9	1.56	3472											117.00	0.0	30.6	69.4	0.20	128.0	46.1	82.6
+110-671D	110-671D-NAC-2		454.10	claystone	late Miocene?	1.75	58.9	2.67	58.4	1.41	3483											108.00	0.0	19.8	80.2	0.30	123.0	52.4	71.0
+110-671C	110-671C-1X-1		496.70	claystone	late Miocene?	1.58	69.8	2.92	58.4	2.04	3808	3000			-808				0.000000000							0.30			
+110-672A	110-672A-2H-3	vertical	7.70	calcareous mudstone	early Pleistocene	1.60	90.3	2.96	69.8	2.32	36	200	175	vertical	164	vertical	139	vertical	0.000002650	0.85	vertical	19.50	20.1	34.9	43.5	35.00	87.4	87.4	56.0
+110-672A	110-672A-2H-3	horizontal	7.80	calcareous mudstone	early Pleistocene	1.60	90.3	2.96	69.8	2.32	37	130	100	horizontal	93	horizontal	43	horizontal	0.000762000	0.47	horizontal	19.50	20.1	34.9	43.5	35.00	87.4	87.4	56.0
+110-672A	110-672A-5H-3	triaxial	36.24	calcareous mudstone	late Pliocene	1.44	106.0	2.81	69.0	2.98	183	196		triaxial	13				0.000000037			134.00				48.60	142.0	142.0	100.0
+110-672A	110-672A-7H-3		55.20	calcareous mudstone	late Pliocene	1.61	77.9	2.74	67.1	2.04	300	310	200		10		-100		0.000008920	0.85		43.00	21.3	53.2	24.9	27.60	116.0	116.0	81.0
+110-672A	110-672A-11H-3		93.50	calcareous mudstone	early Plocene	1.24	56.6	2.65	69.8	2.31	540	340	300		-200		-240		0.000000093	0.62		93.00	3.4	28.2	68.3	41.50	113.0	113.0	80.0
+110-672A	110-672A-11H-4	triaxial	94.74	calcareous mudstone	early Pliocene	1.53	80.0		68.0		548	460		triaxial	-88											31.70	132.0	132.0	72.0
+110-672A	110-672A-13H-3		112.30	calcareous mudstone	late Miocene	1.68	66.9	2.64	62.7	1.68	669	380	260		-289		-409		0.000000107	0.54		87.60	0.8	21.5	77.8	7.80	141.0	141.0	99.0
+110-672A	110-672A-19X-3	vertical	165.80	calcareous mudstone	late Miocene?	1.48	120.0	2.65	74.5	2.92	982	66	53	vertical	-916	vertical	-929	vertical	0.000000250	0.65	vertical		0.4	32.6	67.0	0.00	165.0	165.0	97.0
+110-672A	110-672A-19X-3	vertical	165.90	calcareous mudstone	late Miocene?	1.60	103.0	2.96	72.9	2.69	983	175	140	horizontal	-808	vertical	-843	vertical	0.000000057	1.03	vertical		0.4	32.6	67.0	0.00	165.0	165.0	97.0
+110-672A	110-672A-24X-3		213.30	claystone	late Oligocene	1.72	60.9	2.70	59.8	1.49	1231	405	340		-826		-891		0.000008860	0.50		108.00	0.0	12.6	87.4	0.15	150.0	150.0	115.0
+110-673A	110-673A-2H-3-1	vertical	12.40	claystone	Pleistocene/Miocene?	1.67	81.9	2.62	63.5	1.74	76	240	228	vertical	164	vertical	152	vertical	0.000000428	0.44	vertical	41.50	1.7	45.1	53.2	0.00	142.0	57.1	85.0
+110-673A	110-673A2H-3-2	vertical	12.50	claystone	Pleistocene/Miocene?	1.67	81.9	2.62	63.5	1.74	77	300	260	vertical	223	vertical	183	vertical	0.000000647	0.67	vertical	41.50	1.7	45.1	53.2	0.00	142.0	57.1	85.0
+110-673B	110-673B-2H-3		12.00	claystone	early Pleistocene	1.57	80.7	2.67	68.6	2.19	75	300	250		225		175		0.000000035	0.78		41.50	0.8	29.9	69.4	0.70	142.0	40.1	102.0
+110-676A	110-676A-2H-3	triaxial	10.64	calcareous mudstone	early Pleistocene	1.45	109.0	2.95	72.9	3.23	51	20		triaxial	-31				0.000000032										
+110-676A	110-676A-3H-3	vertical	20.10	calcareous mudstone	early Pleistocene	1.62	73.4	2.85	67.4	2.04	107	25	12	vertical	-82	vertical	-95	vertical				72.00				30.10	131.0	43.0	88.0
+110-676A	110-676A-3H-3	horizontal	20.20	calcareous mudstone	early Pleistocene	1.63	73.4	2.90	67.8	2.08	108	108	70	horizontal	0		-38	horizontal				48.98				19.20			
+110-676A	110-676A-5H-3	triaxial	39.12	calcareous mudstone	early Pleistocene	1.54	83.0	2.86	74.6	2.36	248	33		triaxial	-215				0.000000106			48.98				19.20			
+110-676A	110-676A-6H-3	vertical	48.70	calcareous mudstone	early Pleistocene	1.56	72.3	2.53	64.4	1.79	265	30	15	vertical	-235	vertical	-250	vertical				211.00				18.60	126.0	46.0	80.0
+110-676A	110-676A-6H-3	horizontal	48.80	calcareous mudstone	early Pleistocene	1.56	72.3	2.51	64.2	1.77	266	30	20	horizontal	-236	horizontal	-246	horizontal				57.64				18.60			
+110-676A	110-676A-8H-3	triaxial	67.65	calcareous mudstone	late Pliocene	1.68	57.0	2.84	60.5	1.62	390	52		triaxial	-338				0.000000014			57.64				18.60			
+110-676A	110-676A-9H-3	vertical	77.20	calcareous mudstone	late Pliocene	1.76	54.6	2.64	57.4	1.35	448	320	285	vertical	-128	vertical	-163	vertical	0.000000201	0.40	vertical	205.00				48.00	84.0	33.0	51.0
+110-676A	110-676A-9H-3	horizontal	77.25	calcareous mudstone	late Pliocene	1.76	54.6	2.64	57.4	1.35	449	310	250	horizontal	-139	horizontal	-199	horizontal	0.000000152	0.54	horizontal	108.63	13.3	43.7	43.1	38.01	132.0	31.1	101.0
+110-676A	110-676A-11H-4	triaxial	97.63	calcareous mudstone	late Pliocene	1.74	48.0		57.5		592	400		triaxial	-192							108.63	13.3	43.7	43.1	36.01	132.0	31.1	101.0
+110-676A	110-676A-12X-3	vertical	105.60	calcareous mudstone	late Pliocene	1.78	48.2	2.70	55.4	1.24	647	420	350	vertical	-227	vertical	-297	vertical	0.000000156	0.50	vertical					30.90	77.0	13.0	64.0
+110-676A	110-676A-12X-3	horizontal	105.70	calcareous mudstone	late Pliocene	1.78	48.2	2.70	55.4	1.24	648	360	300	horizontal	-288	horizontal	-348	horizontal	0.000000061	0.44	horizontal	123.22	2.2	22.7	75.1	22.23	98.5	29.8	68.7