천리안 해양위성의 대기보정 및 대리교정 연구

본 학위논문은 세계최초의 정지궤도 해색 위성인 천리안 해양 위성 (GOCI : Geostationary Ocean Color Imager)에 표준으로 사용되는 대기보정 이론에 대하여 기술하고 있다. 타 극궤도 해색위성들이 1~2일 주기로 한 장소를 방문하며 전 지구를 관측하는 것과 달리 천리안 해양위성은 한반도를 포함한 동북아해역을 0.5 km 공간해상도로 낮 시간 동안 1시간의 시간간격으로 관측하고 있으며 (하루 8회 관측) 가시광~근적외파장대 (412, 443, 490, 555, 660, 680, 745, 865 nm) 영역에서 관측한다. 대기상층 위성궤도에서 일반적인 맑은 해역을 대상으로 관측된 가시광~근적외파장대 신호 중 90%이상은 대기신호이며, 해수신호의 크기는 10% 미만을 차지한다. 대기신호의 크기가 해수신호의 크기보다 10배 이상 크기 때문에 1%의 대기신호 추정 오차는 10%이상의 해수 광 스펙트럼 추정오류를 일으킨다. 이런 이유로 위성을 통한 해색원격탐사 임무는 높은 대기보정 정밀도를 요구하고 있으며 대기보정의 개발이 해색원격탐사 알고리즘 개발 중 가장 핵심이 된다. 천리안 해양위성 표준 대기보정은 NASA가 해색원격탐사 임무를 위해 개발한 SeaWiFS 표준 대기보정에 이론적인 기반을 두고 있다. SeaWiFS 방법은 우선 두개의 근적외 파장대 관측결과와 복사전달시뮬레이션 결과(조견표)를 서로 비교하여 대기 중 에어로졸 입자의 종류 및 농도 최적값을 추정해 내며 이 추정결과를 바탕으로 모든 가시광 파장의 에어로졸 반사도 스펙트럼을 다시 조견표를 이용하여 계산한다. 천리안 해양위성의 대기보정도 유사하게 두 근적외파장대 에어로졸 반사도 상관관계를 이용하여 에어로졸 종류 및 농도를 계산하는데, 이 연구를 통하여 SeaWiFS 및 다른 유사 대기보정 방법들과 비교하여 정확도 뿐 아니라 계산 효율 또한 개선하였다. 추가적으로 SeaWiFS에 적용된 수증기 흡광 보정 모델을 천리안 해양위성의 분광특성에 맞게 수정하여 적용하였으며, 탁도가 높은 해역에서 대기보정 오차를 줄이는 방법도 천리안 해양위성 관측영역의 해수 광 특성 및 반사도 정보들을 이용하여 개발하였다. 초기버전의 천리안 해양위성 표준 대기보정의 검보정 결과 탁도가 높은 연안해역에서는 10% 내외의 만족할 만한 오차수준을 보여주었으나, 탁도가 낮은 해역에서는 50% 이상의 오차를 발생되었다. 이는 대리교정 수행의 부재가 주된 요인이며, 본 연구에서는 이를 보완하기 위해 SeaWiFS 표준 대리교정 프로세스에 기반을 두고 천리안 해양위성에 맞게 대리교정을 수행하였다. 이 대리교정 방법에서는 특정 해역의 에어로졸 광특성이 항상 해양성 에어로졸이라 가정하고 이를 바탕으로 근적외 파장대 위성 관측 조도를 시뮬레이션 하여 두 근적외 파장대를 먼저 상대교정 한다. 이후, 상대교정된 두 근적외 파장대를 이용하면 맑은 해역에서 복사전달시뮬레이션을 통하여 가시광 파장대 대기조도를 모의 할 수 있게 되고, 여기에 맑은 해역의 현장 광 측정 자료가 추가되면 가시광파장대 위성관측조도의 시뮬레이션이 가능하다. 이 가시광파장대 모의 결과와 실제 위성관측조도와 비교하면 가시광파장대 대리교정을 완료할 수 있다. 본 대리교정 결과 대리교정 상수가 최대 3.2% 바뀌었으며 (490 nm 밴드) 새 대리교정 상수 적용 시 맑은 해역 대기보정 정확도가 최대 50% 이상 상승하였다. 본 연구에서는 천리안해양위성 대기보정의 성능을 평가하기 위해서 대기보정 결과 원격반사도 (remote-sensing reflectance: Rrs)를 한국해양과학기술원 해양위성연구센터에서 2010년 이후로 한반도 주변 해역 현장조사를 통해 수집한 원격반사도 자료들과 비교검정 하였으며, 검정결과 76, 84, 88, 90, 81, 82%의 정확도를 보여주었다. 추가로 현장자료가 아닌 시뮬레이션 자료를 통해 천리안 해양위성 알고리즘 뿐 아니라 다른 해색원격탐사 임무를 위해 개발된 주요 대기보정 알고리즘들 구현하여 함께 비교검증 하였고, 본 비교검증에서도 천리안 해양위성 표준 대기보정이 다른 대기보정 방법들과 비교하여 가장 낮은 오차율을 보여주었으며, 특히 다중산란 효과가 큰 작은 입자크기의 에어로졸 모델에서 더 좋은 성능을 보여주었다. 본 연구결과는 이론적으로 SeaWiFS 등 비슷한 밴드 특성을 가진 타 해색위성의 대기보정방법으로도 적용이 가능하며, 천리안 해양위성 자료처리시스템 (GOCI data processing system: GDPS) 1.5버전에의 적용될 예정이다.Chapter 1. Introduction 1 1.1 Ocean color remote sensing 1 1.2 Geostationary Ocean Color Imager (GOCI) 2 1.3 Atmospheric correction and vicarious calibration 5 Chapter 2. Initial atmospheric correction for the GOCI data 10 2.1 Introduction 10 2.2 Method 11 2.2.1 Correction for gaseous absorption and whitecap radiance 13 2.2.2 Solar irradiance normalization 15 2.2.3 Correction for molecular (Rayleigh) scattering 17 2.2.4 Cloud mask 18 2.2.5 Correction for multiple scattering by aerosols 19 2.2.6 Correction for atmospheric transmittance 22 2.2.7 Correction for near-infrared water reflectance over turbid waters 22 2.3 Conclusion 24 Chapter 3. Algorithm updates and vicarious calibration for the GOCI atmospheric correction 25 3.1 Backgrounds 25 3.2 Updates to the initial GOCI atmospheric correction algorithm 26 3.2.1 Correction for gaseous absorption and whitecap radiance 26 3.2.2 Sun-glint correction 28 3.2.3 Considering gravity effect for Rayleigh scattering 29 3.2.4 Correction for multiple scattering by aerosols - SRAMS 30 3.2.5 Correction for bidirectional effects for water reflectance 35 3.2.6 Correction for near-infrared water reflectance over turbid waters 39 3.2.7 Atmospheric transmittance with considering anisotropic angular distribution of water reflectance 40 3.3 Vicarious calibration of GOCI near-infrared bands 41 3.3.1 Method 44 3.3.2 Inter-calibration of GOCI near-infrared bands 45 3.3.3 Vicarious calibration of GOCI visible bands 49 Chapter 4. Validation results 51 4.1 Data 51 4.1.1 Synthetic data derived by simulations 51 4.1.2 In situ radiometric data measured from shipboard 52 4.1.3 AERONET-OC radiometric data 56 4.2 Validation of SRAMS scheme with simulation data 58 4.3 Assessment of the atmospheric correction improvements with in situ radiometric data 59 Chapter 5. Discussions 61 5.1 Impacts of water vapor correction on ocean color products 61 5.2 Stability for high solar and satellite zenith angle for diurnal observation 62 5.3 Cloud masking on fast-moving clouds and quality analysis 63 5.4 Evaluation of the GOCI aerosol correction scheme compared with other approaches 64 5.4.1 Aerosol correction approach for OCTS 64 5.4.2 Aerosol correction approach for MERIS 67 5.4.3 Evaluation results 69 5.5 Pitfalls in estimation of aerosol reflectance using 2-NIR bands 71 5.6 Issues in the vicarious calibration of GOCI VIS and NIR bands 72 5.7 Uncertainties from bidirectional effect 75 Chapter 6. Conclusion 76 Appendix. Glossary of symbols 82 Acknowledgements 86 References 88Docto

A semi-analytical algorithm for deriving the particle size distribution slope of turbid inland water based on OLCI data: A case study in Lake Hongze

The particle size distribution (PSD) slope (ξ) can indicate the predominant particle size, material composition, and inherent optical properties (IOPs) of inland waters. However, few semi-analytical methods have been proposed for deriving ξ from the surface remote sensing reflectance due to the variable optical state of inland waters. A semi-analytical algorithm was developed for inland waters having a wide range of turbidity and ξ in this study. Application of the proposed model to Ocean and Land Color Instrument (OLCI) imagery of the water body resulted in several important observations: (1) the proposed algorithm (754 nm and 779 nm combination) was capable of retrieving ξ with R2 being 0.72 (p < 0.01, n = 60), and MAPE and RMSE being 4.37% and 0.22 (n = 30) respectively; (2) the ξ in HZL was lower in summer than other seasons during the period considered, this variation was driven by the phenological cycle of algae and the runoff caused by rainfall; (3) the band optimization proposed in this study is important for calculating the particle backscattering slope (η) and deriving ξ because it is feasible for both algae dominant and sediment governed turbid inland lakes. These observations help improve our understanding of the relationship between IOPs and ξ, which are affected by different bio-optic processes and algal phenology in the lake environment

Developments in Earth observation for the assessment and monitoring of inland, transitional, coastal and shelf-sea waters

The Earth's surface waters are a fundamental resource and encompass a broad range of ecosystems that are core to global biogeochemical cycling and food and energy production. Despite this, the Earth's surface waters are impacted by multiple natural and anthropogenic pressures and drivers of environmental change. The complex interaction between physical, chemical and biological processes in surface waters poses significant challenges for in situ monitoring and assessment and often limits our ability to adequately capture the dynamics of aquatic systems and our understanding of their status, functioning and response to pressures. Here we explore the opportunities that Earth observation (EO) has to offer to basin-scale monitoring of water quality over the surface water continuum comprising inland, transition and coastal water bodies, with a particular focus on the Danube and Black Sea region. This review summarises the technological advances in EO and the opportunities that the next generation satellites offer for water quality monitoring. We provide an overview of algorithms for the retrieval of water quality parameters and demonstrate how such models have been used for the assessment and monitoring of inland, transitional, coastal and shelf-sea systems. Further, we argue that very few studies have investigated the connectivity between these systems especially in large river-sea systems such as the Danube-Black Sea. Subsequently, we describe current capability in operational processing of archive and near real-time satellite data. We conclude that while the operational use of satellites for the assessment and monitoring of surface waters is still developing for inland and coastal waters and more work is required on the development and validation of remote sensing algorithms for these optically complex waters, the potential that these data streams offer for developing an improved, potentially paradigm-shifting understanding of physical and biogeochemical processes across large scale river-sea continuum including the Danube-Black Sea is considerable

Assessment of Normalized Water-Leaving Radiance Derived From Goci Using AERONET-OC Data

The geostationary ocean color imager (GOCI), as the world’s first operational geostationary ocean color sensor, is aiming at monitoring short-term and small-scale changes of waters over the northwestern Pacific Ocean. Before assessing its capability of detecting subdiurnal changes of seawater properties, a fundamental understanding of the uncertainties of normalized water-leaving radiance (nLw) products introduced by atmospheric correction algorithms is necessarily required. This paper presents the uncertainties by accessing GOCI-derived nLw products generated by two commonly used operational atmospheric algorithms, the Korea Ocean Satellite Center (KOSC) standard atmospheric algorithm adopted in GOCI Data Processing System (GDPS) and the NASA standard atmospheric algorithm implemented in Sea-Viewing Wide Field-of-View Sensor Data Analysis System (SeaDAS/l2gen package), with Aerosol Robotic Network Ocean Color (AERONET-OC) provided nLw data. The nLw data acquired from the GOCI sensor based on two algorithms and four AERONET-OC sites of Ariake, Ieodo, Socheongcho, and Gageocho from October 2011 to March 2019 were obtained, matched, and analyzed. The GDPS-generated nLw data are slightly better than that with SeaDAS at visible bands; however, the mean percentage relative errors for both algorithms at blue bands are over 30%. The nLw data derived by GDPS is of better quality both in clear and turbid water, although underestimation is observed at near-infrared (NIR) band (865 nm) in turbid water. The nLw data derived by SeaDAS are underestimated in both clear and turbid water, and the underestimation worsens toward short visible bands. Moreover, both algorithms perform better at noon (02 and 03 Universal Time Coordinated (UTC)), and worse in the early morning and late afternoon. It is speculated that the uncertainties in nLw measurements arose from aerosol models, NIR water-leaving radiance correction method, and bidirectional reflectance distribution function (BRDF) correction method in corresponding atmospheric correction procedure

Exploring Himawari-8 geostationary observations for the advanced coastal monitoring of the Great Barrier Reef

Larissa developed an algorithm to enable water-quality assessment within the Great Barrier Reef (GBR) using weather satellite observations collected every 10 minutes. This unprecedented temporal resolution records the dynamic nature of water quality fluctuations for the entire GBR, with applications for improved monitoring and management

Satellite Ocean Colour: Current Status and Future Perspective

Spectrally resolved water-leaving radiances (ocean colour) and inferred chlorophyll concentration are key to studying phytoplankton dynamics at seasonal and interannual scales, for a better understanding of the role of phytoplankton in marine biogeochemistry; the global carbon cycle; and the response of marine ecosystems to climate variability, change and feedback processes. Ocean colour data also have a critical role in operational observation systems monitoring coastal eutrophication, harmful algal blooms, and sediment plumes. The contiguous ocean-colour record reached 21 years in 2018; however, it is comprised of a number of one-off missions such that creating a consistent time-series of ocean-colour data requires merging of the individual sensors (including MERIS, Aqua-MODIS, SeaWiFS, VIIRS, and OLCI) with differing sensor characteristics, without introducing artefacts. By contrast, the next decade will see consistent observations from operational ocean colour series with sensors of similar design and with a replacement strategy. Also, by 2029 the record will start to be of sufficient duration to discriminate climate change impacts from natural variability, at least in some regions. This paper describes the current status and future prospects in the field of ocean colour focusing on large to medium resolution observations of oceans and coastal seas. It reviews the user requirements in terms of products and uncertainty characteristics and then describes features of current and future satellite ocean-colour sensors, both operational and innovative. The key role of in situ validation and calibration is highlighted as are ground segments that process the data received from the ocean-colour sensors and deliver analysis-ready products to end-users. Example applications of the ocean-colour data are presented, focusing on the climate data record and operational applications including water quality and assimilation into numerical models. Current capacity building and training activities pertinent to ocean colour are described and finally a summary of future perspectives is provided

Feasibility Study for an Aquatic Ecosystem Earth Observing System Version 1.2.

International audienceMany Earth observing sensors have been designed, built and launched with primary objectives of either terrestrial or ocean remote sensing applications. Often the data from these sensors are also used for freshwater, estuarine and coastal water quality observations, bathymetry and benthic mapping. However, such land and ocean specific sensors are not designed for these complex aquatic environments and consequently are not likely to perform as well as a dedicated sensor would. As a CEOS action, CSIRO and DLR have taken the lead on a feasibility assessment to determine the benefits and technological difficulties of designing an Earth observing satellite mission focused on the biogeochemistry of inland, estuarine, deltaic and near coastal waters as well as mapping macrophytes, macro-algae, sea grasses and coral reefs. These environments need higher spatial resolution than current and planned ocean colour sensors offer and need higher spectral resolution than current and planned land Earth observing sensors offer (with the exception of several R&D type imaging spectrometry satellite missions). The results indicate that a dedicated sensor of (non-oceanic) aquatic ecosystems could be a multispectral sensor with ~26 bands in the 380-780 nm wavelength range for retrieving the aquatic ecosystem variables as well as another 15 spectral bands between 360-380 nm and 780-1400 nm for removing atmospheric and air-water interface effects. These requirements are very close to defining an imaging spectrometer with spectral bands between 360 and 1000 nm (suitable for Si based detectors), possibly augmented by a SWIR imaging spectrometer. In that case the spectral bands would ideally have 5 nm spacing and Full Width Half Maximum (FWHM), although it may be necessary to go to 8 nm wide spectral bands (between 380 to 780nm where the fine spectral features occur -mainly due to photosynthetic or accessory pigments) to obtain enough signal to noise. The spatial resolution of such a global mapping mission would be between ~17 and ~33 m enabling imaging of the vast majority of water bodies (lakes, reservoirs, lagoons, estuaries etc.) larger than 0.2 ha and ~25% of river reaches globally (at ~17 m resolution) whilst maintaining sufficient radiometric resolution

Determination of the downwelling diffuse attenuation coefficient of lakewater with the sentinel-3A OLCI

The Ocean and Land Color Imager (OLCI) on the Sentinel-3A satellite, which was launched by the European Space Agency in 2016, is a new-generation water color sensor with a spatial resolution of 300 m and 21 bands in the range of 400-1020 nm. The OLCI is important to the expansion of remote sensing monitoring of inland waters using water color satellite data. In this study, we developed a dual band ratio algorithm for the downwelling diffuse attenuation coefficient at 490 nm (Kd(490)) for the waters of Lake Taihu, a large shallow lake in China, based on data measured during seven surveys conducted between 2008 and 2017 in combination with Sentinel-3A-OLCI data. The results show that: (1) Compared to the available Kd(490) estimation algorithms, the dual band ratio (681 nm/560 nm and 754 nm/560 nm) algorithm developed in this study had a higher estimation accuracy (N = 26, coefficient of determination (R2) = 0.81, root-mean-square error (RMSE) = 0.99m-1and mean absolute percentage error (MAPE) = 19.55%) and validation accuracy (N = 14, R2= 0.83, RMSE = 1.06 m-1and MAPE = 27.30%), making it more suitable for turbid inland waters; (2) A comparison of the OLCI Kd(490) product and a similar Moderate Resolution Imaging Spectroradiometer (MODIS) product reveals a high consistency between the OLCI and MODIS products in terms of the spatial distribution of Kd(490). However, the OLCI product has a smoother spatial distribution and finer textural characteristics than the MODIS product and contains notably higher-quality data; (3) The Kd(490) values for Lake Taihu exhibit notable spatial and temporal variations. Kd(490) is higher in seasons with relatively high wind speeds and in open waters that are prone to wind- and wave-induced sediment resuspension. Finally, the Sentinel-3A-OLCI has a higher spatial resolution and is equipped with a relatively wide dynamic range of spectral bands suitable for inland waters. The Sentinel-3B satellite will be launched soon and, together with the Sentinel-3A satellite, will form a two-satellite network with the ability to make observations twice every three days. This satellite network will have a wider range of application and play an important role in the monitoring of inland waters with complex optical properties

An Overview of Approaches and Challenges for Retrieving Marine Inherent Optical Properties from Ocean Color Remote Sensing

Ocean color measured from satellites provides daily global, synoptic views of spectral water-leaving reflectancesthat can be used to generate estimates of marine inherent optical properties (IOPs). These reflectances, namelythe ratio of spectral upwelled radiances to spectral downwelled irradiances, describe the light exiting a watermass that defines its color. IOPs are the spectral absorption and scattering characteristics of ocean water and itsdissolved and particulate constituents. Because of their dependence on the concentration and composition ofmarine constituents, IOPs can be used to describe the contents of the upper ocean mixed layer. This informationis critical to further our scientific understanding of biogeochemical oceanic processes, such as organic carbonproduction and export, phytoplankton dynamics, and responses to climatic disturbances. Given their im-portance, the international ocean color community has invested significant effort in improving the quality of satellite-derived IOP products, both regionally and globally. Recognizing the current influx of data products intothe community and the need to improve current algorithms in anticipation of new satellite instruments (e.g., theglobal, hyperspectral spectroradiometer of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mis-sion), we present a synopsis of the current state of the art in the retrieval of these core optical properties.Contemporary approaches for obtaining IOPs from satellite ocean color are reviewed and, for clarity, separatedbased their inversion methodology or the type of IOPs sought. Summaries of known uncertainties associated witheach approach are provided, as well as common performance metrics used to evaluate them. We discuss currentknowledge gaps and make recommendations for future investment for upcoming missions whose instrumentcharacteristics diverge sufficiently from heritage and existing sensors to warrant reassessing current approaches

An Overview of Approaches and Challenges for Retrieving Marine Inherent Optical Properties from Ocean Color Remote Sensing

Ocean color measured from satellites provides daily global, synoptic views of spectral water-leaving reflectances that can be used to generate estimates of marine inherent optical properties (IOPs). These reflectances, namely the ratio of spectral upwelled radiances to spectral downwelled irradiances, describe the light exiting a water mass that defines its color. IOPs are the spectral absorption and scattering characteristics of ocean water and its dissolved and particulate constituents. Because of their dependence on the concentration and composition of marine constituents, IOPs can be used to describe the contents of the upper ocean mixed layer. This information is critical to further our scientific understanding of biogeochemical oceanic processes, such as organic carbon production and export, phytoplankton dynamics, and responses to climatic disturbances. Given their importance, the international ocean color community has invested significant effort in improving the quality of satellite-derived IOP products, both regionally and globally. Recognizing the current influx of data products into the community and the need to improve current algorithms in anticipation of new satellite instruments (e.g., the global, hyperspectral spectroradiometer of the NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission), we present a synopsis of the current state of the art in the retrieval of these core optical properties. Contemporary approaches for obtaining IOPs from satellite ocean color are reviewed and, for clarity, separated based their inversion methodology or the type of IOPs sought. Summaries of known uncertainties associated with each approach are provided, as well as common performance metrics used to evaluate them. We discuss current knowledge gaps and make recommendations for future investment for upcoming missions whose instrument characteristics diverge sufficiently from heritage and existing sensors to warrant reassessing current approaches