MODIS ET (MOD16)¶

MODIS/Terra·Aqua Net Evapotranspiration MOD16A2(GF)·MYD16A2(GF)·MOD16A3GF·MYD16A3GF Version 6.1 (MODIS 위성 기반 전지구 500 m 실제·잠재 증발산)

지원 축소·대체 예정

한눈에 보기¶

기간 참고: 연말 결측 보완판(GF)이 장기 기록입니다. MOD16A2GF는 2000년(CMR 표기 2000-01-01, Terra 자료는 2000년 2월 하순부터), MYD16A2GF는 2002년(Aqua 자료는 2002-07-04부터)부터 2025년 12월까지 있습니다(2025년분은 2026-02-05 생산). 보완 전 forward 판(MOD16A2·MYD16A2)은 V6.1에서 2021-01-01부터만 남아 있으며 2026-10-13 확인 시 2026-09-22~29 구간까지 있습니다. 연자료는 GF판(MOD16A3GF·MYD16A3GF)만 CMR에 등록되어 있고 2025년분까지 있습니다.

개요¶

MOD16은 MODIS의 엽면적지수·광합성유효복사 흡수율(LAI·FPAR), 알베도, 토지피복과 GMAO 일별 기상 재분석을 Penman-Monteith 식 기반 알고리즘(Mu et al., 2011)에 넣어 전지구 식생 육지의 실제 증발산(ET)·잠열(LE)과 잠재 증발산(PET)·잠재 잠열(PLE)을 500 m, 8일 간격으로 계산한 4단계(L4) 산출물입니다. Terra판(MOD16)과 Aqua판(MYD16)이 있으며, 연말에 1년치 LAI·FPAR의 품질 나쁜 값을 시간 보간으로 메워 다시 만든 결측 보완판(GF)이 연구용 표준입니다. 물수지·가뭄·농업 용수 연구에 널리 쓰이지만, 도시·사막·수면·영구 빙설에는 값이 없습니다.

후속 자료와 무엇이 달라졌나¶

MOD16 C6 → MOD16 C6.1¶

전환 시점: 2021-03

과학 알고리즘은 그대로이고, 입력 MODIS L1B 보정 개선과 함께 결측 처리 방식이 바뀌었습니다. C6은 오염된 LAI·FPAR 화소를 채움값으로 남겨 빈 곳이 많았지만, C6.1은 기후값 LAI·FPAR와 결측 보완 알베도를 예비 입력으로 써서 빈 곳을 채웁니다. C6 자료는 2023-07-31 배포가 끝났습니다.

출처: MOD16 C6.1 User's Guide · Final Notice: MODIS Version 6.0 Land Data Products Decommissioned July 31, 2023

MODIS 증발산 MOD16/MYD16 C6.1 → VIIRS 증발산 VNP16·VJ116 V2.1 (공개 전) 예정¶

NASA는 MODIS 종료에 대비해 같은 MOD16 알고리즘을 VIIRS(Suomi NPP·NOAA-20·NOAA-21)에 적용한 VNP16·VJ116·VJ216 V2를 2025~2026년에 공개했지만, 자료 문제가 발견되어 2026-08-17 모두 철회하고 '과학 품질이 아니므로 연구에 쓰지 말라'고 공지했습니다. 문제를 고친 뒤 V2.1로 재처리할 예정이며 일정은 정해지지 않았습니다.

출처: VIIRS Actual and Potential Evapotranspiration Data Products Unavailable (LP DAAC, 2026-08-17) · VIIRS Version 2 Actual and Potential Evapotranspiration Data Products Released (LP DAAC, 2025-07-10) · VIIRS Land — NASA Evapotranspiration and Potential Evapotranspiration

활용 가능한 변수¶

LP DAAC 데이터 카탈로그와 C6.1 사용자 가이드 기준 변수(SDS)입니다. HDF-EOS2(HDF4) 파일의 실제 값 = 저장값 × scale_factor이며, 8일 자료는 저장값 32761~32767이 채움값(토지피복별 비식생 표시)이므로 반드시 먼저 걸러야 합니다. MYD16(Aqua) 파일도 변수 이름이 같습니다.

순간값 평균 누적 최대 최소 고정값 — 순간값: 표시 시각의 값 · 평균/누적/최대/최소: 표시된 구간 동안의 값 · 고정값: 시간에 따라 변하지 않음

8일 자료 (MOD16A2GF·MYD16A2GF 연말 보완판 / MOD16A2·MYD16A2 forward판)¶

연자료 (MOD16A3GF·MYD16A3GF, 연말 보완판만)¶

전체 변수 목록은 공식 문서에서 확인하세요.

받는 방법¶

• NASA Earthdata Search (LP DAAC, Earthdata Cloud) — HTTPS 다운로드 / earthaccess / Earthdata Search
  • short name은 MOD16A2GF·MYD16A2GF(8일 보완판), MOD16A3GF·MYD16A3GF(연 보완판), MOD16A2·MYD16A2(8일 forward판)이고 version은 061입니다. 파일은 data.lpdaac.earthdatacloud.nasa.gov 아래 HTTPS로 받으며 무료 Earthdata 로그인이 필요합니다. 예전 LP DAAC Data Pool 배포는 2025-06-30에 끝났습니다.
• AppEEARS (영역·지점 부분 추출) — 웹 / API — 지점 시계열·영역을 GeoTIFF·NetCDF4로 추출
  • 필요한 영역·기간·변수만 골라 받을 수 있고 위경도 투영으로 바꿔 받을 수 있어 한반도 같은 작은 영역 분석에 편합니다.
• AWS us-west-2 직접 접근 (S3) — 클라우드 저장소(S3, us-west-2 안에서만)
  • s3://lp-prod-protected/MOD16A2GF.061/ 처럼 저장되어 있으며 임시 자격증명이 필요합니다.

예시 코드¶

earthaccess로 2023년 한반도 주변 타일의 Terra 8일 보완판을 받아 pyhdf로 읽는 예입니다(HDF4 형식이라 xarray 기본 엔진으로는 열리지 않을 수 있습니다). 채움값(32761 이상)을 먼저 걸러야 합니다.

# pip install earthaccess pyhdf numpy
import earthaccess
import numpy as np
from pyhdf.SD import SD, SDC

earthaccess.login()  # 무료 NASA Earthdata 계정
results = earthaccess.search_data(
    short_name="MOD16A2GF", version="061",
    temporal=("2023-07-01", "2023-07-31"),
    bounding_box=(124, 33, 132, 39),   # (서, 남, 동, 북)
)
files = earthaccess.download(results, "mod16")

hdf = SD(files[0], SDC.READ)
raw = hdf.select("ET_500m").get().astype("float64")
raw[raw > 32700] = np.nan          # 32761~32767: 물·사막·도시·빙설 등 채움값
et = raw * 0.1                      # mm / 8일(연말 구간은 5~6일)
qc = hdf.select("ET_QC_500m").get()
print(files[0], np.nanmean(et))

실전 팁¶

• 연구에는 연말 결측 보완판(MOD16A2GF·MYD16A2GF·MxD16A3GF)을 쓰세요. forward판(MOD16A2)은 구름 낀 LAI·FPAR 대신 최근 5년 기후값을 넣어 만든 준실시간 자료이고 V6.1에서는 2021년부터만 남아 있습니다.
• ET_500m 저장값 32761~32767(연자료 65529~65535)은 결측이 아니라 비식생 토지피복 표시입니다(32766 물, 32765 사막·희박 식생, 32764 영구 빙설, 32763 습지, 32762 도시, 32761 미분류). scale_factor를 곱하기 전에 걸러야 터무니없는 값이 섞이지 않습니다.
• 단위가 다릅니다. ET·PET는 8일 합계(mm), LE·PLE는 하루 평균(J m-2 day-1, scale 10000)입니다. 연말 마지막 구간은 5~6일 합계입니다.
• MODIS 토지피복(IGBP 17종, 3년 평활)에 따라 정한 생물군계별 고정 매개변수를 쓰므로, 토지피복 분류 오류가 그대로 증발산 오차가 되고 같은 생물군계 안의 지역 차이는 반영되지 않습니다. 한 해 내내 유효한 알베도가 없는 화소는 0.4를 넣어 계산합니다.
• 기상 입력은 GMAO 재분석(0.5°×0.6° 또는 1°×1.25°)을 500 m로 보간한 것이라, 산악·해안처럼 기상이 급변하는 곳에서는 500 m 해상도만큼의 정보가 없습니다.
• 2000년(Terra)·2002년(Aqua) 첫해 GF 자료는 연초 관측이 없어 시간 보간이 불완전하므로 덜 믿을 만합니다.
• MODIS Sinusoidal 투영 타일(한반도는 대부분 h27v05·h28v05)입니다. 위경도 격자로 바꾸려면 GDAL·rioxarray로 재투영하거나 AppEEARS를 쓰세요.
• Terra·Aqua는 2026년 말 이후 운용 종료가 예정되어 있고, 후속 VIIRS 증발산은 2026-08-17 철회되었습니다. 장기 시계열을 이어 가야 한다면 MODIS 종료 시점과 VIIRS V2.1 재공개 공지를 함께 확인하세요.

공지 및 변경 이력¶

논문에 넣을 레퍼런스¶

데이터셋 인용¶

MOD16A2GF V6.1 (Terra 8일 보완판) — DOI 10.5067/MODIS/MOD16A2GF.061

@misc{mod16a2gf_v061,
  author    = {Running, Steve and Mu, Qiaozhen and Zhao, Maosheng and Moreno, Alvaro},
  title     = {{MODIS/Terra} Net Evapotranspiration Gap-Filled 8-Day {L4} Global 500m {SIN} Grid {V061}},
  year      = {2021},
  publisher = {NASA Land Processes Distributed Active Archive Center},
  doi       = {10.5067/MODIS/MOD16A2GF.061}
}

MOD16A3GF V6.1 (Terra 연 보완판) — DOI 10.5067/MODIS/MOD16A3GF.061

@misc{mod16a3gf_v061,
  author    = {Running, Steve and Mu, Qiaozhen and Zhao, Maosheng and Moreno, Alvaro},
  title     = {{MODIS/Terra} Net Evapotranspiration Gap-Filled Yearly {L4} Global 500m {SIN} Grid {V061}},
  year      = {2021},
  publisher = {NASA Land Processes Distributed Active Archive Center},
  doi       = {10.5067/MODIS/MOD16A3GF.061}
}

MOD16A2 V6.1 (Terra 8일 forward판) — DOI 10.5067/MODIS/MOD16A2.061

@misc{mod16a2_v061,
  author    = {Running, Steve and Mu, Qiaozhen and Zhao, Maosheng},
  title     = {{MODIS/Terra} Net Evapotranspiration 8-Day {L4} Global 500m {SIN} Grid {V061}},
  year      = {2021},
  publisher = {NASA Land Processes Distributed Active Archive Center},
  doi       = {10.5067/MODIS/MOD16A2.061}
}

논문 인용¶

MOD16 알고리즘 (Mu et al., 2011) — DOI 10.1016/j.rse.2011.02.019

@article{mu2011mod16,
  author  = {Mu, Qiaozhen and Zhao, Maosheng and Running, Steven W.},
  title   = {Improvements to a {MODIS} global terrestrial evapotranspiration algorithm},
  journal = {Remote Sensing of Environment},
  volume  = {115},
  number  = {8},
  pages   = {1781--1800},
  year    = {2011},
  doi     = {10.1016/j.rse.2011.02.019}
}

사사(Acknowledgements) 작성¶

쓴 제품(Terra·Aqua, 8일·연, 보완판·forward판)에 맞는 LP DAAC 데이터셋 DOI를 인용하세요. Aqua판 DOI는 10.5067/MODIS/MYD16A2GF.061, 10.5067/MODIS/MYD16A3GF.061, 10.5067/MODIS/MYD16A2.061입니다. 알고리즘은 Mu et al. (2011)을 함께 인용합니다.

Information on the datasets used in this study was compiled with the help of the Korean Climate Data Navigator (KCDN) (Sim, 2026; https://doi.org/10.5281/zenodo.23300923).

참고 링크¶

• MOD16A2GF V061 데이터 카탈로그 (LP DAAC)
• MOD16 C6.1 User's Guide
• MOD16 ATBD
• MOD16 검증 상태 (MODIS Land)
• MODIS/VIIRS 지면 품질평가(LDOPE) 알려진 문제
• MODIS Land News (Terra·Aqua 궤도 표류 안내)
• VIIRS 증발산 제품 소개 (VIIRS Land)