2022 개정 교육과정 연계 온디바이스 해양생지화학·탄소순환 STEM 시리즈 #63

해양 생물학적 탄소 펌프(BCP) POC 플럭스 서보 엔진

마틴 곡선(Martin Curve) 멱법칙 침강 감쇄 모델 및 0-Alloc C# 바이오-ARGO 광학 역산 제어 루프

ZeliDesk 해양탄소순환·생지화학팀 바이오-ARGO플로트·심해탄소격리 분과

2026년 9월 12일

1. 마틴(Martin) 멱함수 침강 모델 및 광학 후방산란 변환식

유광층 하부 경계면($z_0 = 100\text{ m}$)을 통과한 POC 침강 플럭스 $F(z)$는 수심 $z$에 따라 마틴 멱함수 곡선을 따릅니다:

여기서 $b$는 무차원 감쇄 지수(Attenuation Exponent, $0.6 \sim 1.4$)로, 미생물 분해 활성도와 수온에 강하게 의존합니다. 해수 입자의 $700\text{ nm}$ 광학 후방산란 계수 $b_{bp}(700)$로부터 현장 POC 농도 $C_{POC}$는 다음과 같이 선형 환산됩니다:

침강 속도 $w_s$와 농도의 곱으로 정의되는 순 플럭스 $F_{net} = C_{POC} \cdot w_s$로부터 수심별 탄소 손실률(재무기화 플럭스)은 공간 미분 $\frac{dF}{dz}$로 직결됩니다:

2. 0-Alloc C# 스택 기반 바이오-ARGO 탄소 플럭스 서보

플로트 하강/상승 센서 스트림으로부터 마틴 $b$ 지수와 심해 격리율을 실시간 산출하는 ref struct 루프입니다.

using System;
using System.Runtime.CompilerServices;

namespace ZeliDesk.Stem.MarineBiogeochemistry
{
    public readonly ref struct BioArgoTelemetry
    {
        public readonly double CurrentDepthM;
        public readonly double Backscattering700MInv;
        public readonly double ChlorophyllFluorescenceMgM3;
        public readonly double DissolvedOxygenUmolKg;
        public readonly double WaterTempC;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public BioArgoTelemetry(double depthM, double bbp700, double chl, double o2, double tempC)
        {
            CurrentDepthM = depthM;
            Backscattering700MInv = bbp700;
            ChlorophyllFluorescenceMgM3 = chl;
            DissolvedOxygenUmolKg = o2;
            WaterTempC = tempC;
        }
    }

    public readonly ref struct PocFluxAnalysisResult
    {
        public readonly double EstimatedPocConcentrationMgPerM3;
        public readonly double MartinAttenuationExponentB;
        public readonly double DeepCarbonSequestrationFluxMgM2Day;
        public readonly bool IsTwilightZoneActive;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PocFluxAnalysisResult(double poc, double bExp, double seqFlux, bool active)
        {
            EstimatedPocConcentrationMgPerM3 = poc;
            MartinAttenuationExponentB = bExp;
            DeepCarbonSequestrationFluxMgM2Day = seqFlux;
            IsTwilightZoneActive = active;
        }
    }

    public ref struct BioArgoPocFluxServo
    {
        private const double AlphaOpt = 37500.0; // bbp -> POC 변환 계수
        private readonly double _referenceFlux100m;
        private double _prevDepthM;
        private double _prevPoc;

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public BioArgoPocFluxServo(double refFlux100m)
        {
            _referenceFlux100m = refFlux100m;
            _prevDepthM = 100.0;
            _prevPoc = 0.0;
        }

        [MethodImpl(MethodImplOptions.AggressiveInlining)]
        public PocFluxAnalysisResult Step(in BioArgoTelemetry telemetry)
        {
            // 1. 광학 후방산란 기반 POC 농도 산출 (mg C/m^3)
            double pocConc = Math.Max(0.1, AlphaOpt * telemetry.Backscattering700MInv);

            // 2. 수온 의존적 마틴 감쇄 지수 b 모델링 (수온이 높을수록 미생물 분해 활발하여 b 상승)
            double bExp = 0.65 + 0.035 * Math.Max(0.0, telemetry.WaterTempC);

            // 3. 현재 수심에서의 마틴 침강 플럭스 계산: F(z) = F100 * (z / 100)^(-b)
            double normalizedDepth = Math.Max(1.0, telemetry.CurrentDepthM / 100.0);
            double depthFlux = _referenceFlux100m * Math.Pow(normalizedDepth, -bExp);

            // 4. 심해저(1000m 기준) 영구 탄소 격리 플럭스 추정
            double deepFlux1000m = _referenceFlux100m * Math.Pow(10.0, -bExp);

            // 5. 중층 박명대(Twilight Zone: 100m ~ 1000m) 활성 판별
            bool isTwilight = telemetry.CurrentDepthM >= 100.0 && telemetry.CurrentDepthM <= 1000.0;

            _prevDepthM = telemetry.CurrentDepthM;
            _prevPoc = pocConc;

            return new PocFluxAnalysisResult(pocConc, bExp, deepFlux1000m, isTwilight);
        }
    }
}

3. 교육적 가치 및 중고교 실험 시뮬레이션 적용