{
  "boundary_analysis": {
    "cases": 9,
    "radius_points": 27,
    "ratio_definition": "Length of OTHER subduction zones inside radius divided by new trench length; not coverage or probability.",
    "mean_at_500km_percent": 120.44444444444444,
    "max_difference_percentage_points": 0.3576177285318636,
    "rounding_note": "Printed lengths treated as rounded to nearest km; compatibility is conditional, not a reconstruction of original exact geometry."
  },
  "chosen_force_model": {
    "b": 1.0,
    "zeta": 1.0,
    "bending": 0.4,
    "s0": 1.5,
    "smin": 0.2,
    "lweak": 0.4,
    "zero_push_threshold": 0.7833937107336856,
    "scope": "Reduced, dimensionless, finite-window illustration; no Earth length/time calibration."
  },
  "cooling_model": {
    "hot_K": 1073.15,
    "cold_K": 873.15,
    "grain_exponent": 3.0,
    "examples": [
      {
        "Q_J_mol": 200000,
        "thermal_factor": 169.73597861697513,
        "grain_ratio_for_equal_viscosity": 0.18061009890248592,
        "viscosity_ratio_at_grain_ratio_point1": 0.16973597861697517
      },
      {
        "Q_J_mol": 300000,
        "thermal_factor": 2211.367196646855,
        "grain_ratio_for_equal_viscosity": 0.07675612574625754,
        "viscosity_ratio_at_grain_ratio_point1": 2.2113671966468558
      }
    ]
  },
  "geological_proxies": [
    {
      "unit": "HT",
      "sample": "KMA22-02",
      "recrystallized_RMS_um": 94,
      "subgrain_arithmetic_um": 67,
      "grain_calibration_stress_MPa": 24,
      "subgrain_calibration_stress_MPa": 5
    },
    {
      "unit": "LT",
      "sample": "KMA24-13",
      "recrystallized_RMS_um": 58,
      "subgrain_arithmetic_um": 28,
      "grain_calibration_stress_MPa": 33,
      "subgrain_calibration_stress_MPa": 12
    }
  ],
  "bending_checks": [
    {
      "thickness_m": 100.0,
      "moment_per_width_N": 29166666666.666668,
      "elastic_limit_N": 29166666666.666668,
      "fully_plastic_limit_N": 375000000000.0
    },
    {
      "thickness_m": 1000.0,
      "moment_per_width_N": 28316326530612.246,
      "elastic_limit_N": 29166666666666.668,
      "fully_plastic_limit_N": 37500000000000.0
    },
    {
      "thickness_m": 30000.0,
      "moment_per_width_N": 3.3740816326530612e+16,
      "elastic_limit_N": 7.875e+17,
      "fully_plastic_limit_N": 3.375e+16
    },
    {
      "thickness_m": 60000.0,
      "moment_per_width_N": 1.3499081632653061e+17,
      "elastic_limit_N": 6.3e+18,
      "fully_plastic_limit_N": 1.35e+17
    }
  ],
  "limits": [
    "No EBSD/microprobe or geological age inversion rerun.",
    "Crameri supplementary boundary analysis covers 9 cases, not all 13 database entries; reconstruction time labels differ from preferred event ages.",
    "Kotowski values transcribed from original-publisher text retrieved by web; full downloadable paper and SI unavailable.",
    "No original 2-D thermomechanical simulations or full Li-Gurnis analytical solution reproduced.",
    "Selected parameters illustrate mechanisms and failure conditions, not universal subduction thresholds."
  ]
}