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

본 학위논문은 세계최초의 정지궤도 해색 위성인 천리안 해양 위성 (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

Chlorophyll-a Algorithms for Oligotrophic Oceans: A Novel Approach Based on Three-Band Reflectance Difference

A new empirical algorithm is proposed to estimate surface chlorophyll-a concentrations (Chl) in the global ocean for Chl less than or equal to 0.25 milligrams per cubic meters (approximately 77% of the global ocean area). The algorithm is based on a color index (CI), defined as the difference between remote sensing reflectance (R(sub rs), sr(sup -1) in the green and a reference formed linearly between R(sub rs) in the blue and red. For low Chl waters, in situ data showed a tighter (and therefore better) relationship between CI and Chl than between traditional band-ratios and Chl, which was further validated using global data collected concurrently by ship-borne and SeaWiFS satellite instruments. Model simulations showed that for low Chl waters, compared with the band-ratio algorithm, the CI-based algorithm (CIA) was more tolerant to changes in chlorophyll-specific backscattering coefficient, and performed similarly for different relative contributions of non-phytoplankton absorption. Simulations using existing atmospheric correction approaches further demonstrated that the CIA was much less sensitive than band-ratio algorithms to various errors induced by instrument noise and imperfect atmospheric correction (including sun glint and whitecap corrections). Image and time-series analyses of SeaWiFS and MODIS/Aqua data also showed improved performance in terms of reduced image noise, more coherent spatial and temporal patterns, and consistency between the two sensors. The reduction in noise and other errors is particularly useful to improve the detection of various ocean features such as eddies. Preliminary tests over MERIS and CZCS data indicate that the new approach should be generally applicable to all existing and future ocean color instruments

Modeling Atmosphere-Ocean Radiative Transfer: A PACE Mission Perspective

The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth's sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates

An investigation into the dynamical and statistical properties of dominant ocean surface waves using close-range remote sensing

Denne avhandlingen er basert på forskningsresultat som behandler statistiske og dynamiske egenskaper av dominante vinddrevne overflatebølger i åpent hav. Med uttrykket dominante bølger refererer vi her til de største bølgene, med størst energi, i en gitt sjøtilstand. Bølgedrevne prosesser er viktige både i klimasammenheng via atmosfære--hav interaksjon som drives i stor grad av bølgebrytning, samt for kommersiell og rekreasjonell offshorevirksomhet p.g.a. risikoen for å bli utsatt for f.eks. ekstreme enkeltbølger. Både bølgebrytning og ekstrembølgestatistikk er i skrivende stund ufullstendig representert i teoretiske og numeriske modeller. Arbeidet som presenteres i denne avhandlingen undersøker de ovennevnte temaene ved bruk av bølgeobservasjoner som er primært samlet inn på Ekofiskfeltet i den sentrale delen av Nordsjøen. Observasjonsdatasettene består av en langtidstidsserie av laser-altimetermålinger og stereoskopiske videodata fra Ekofisk, samt videomålinger av brytende bølger fra et forskningstokt i nordre Stillehavet. Forskningsresultatene er presentert i artikkelform med to publiserte verk og ett innlevert manuskript. Det blir påvist en tydelig forbindelse mellom økt bølgebrytning og dominante bølgegrupper, et resultat som tidligere har blitt påvist i laboratorie- og modelleksperiment, men sjeldent ved bruk av feltobservasjoner. Tredimensjonale stereo-rekonstruksjoner viser også at ekstreme bølgekammer, både brytende og ikke-brytende, følger nylig utviklet teori om ikke-lineær bølgegruppedynamikk. Dette funnet har konsekvenser f.eks. for estimering av geometriske og kinematiske bølgeegenskaper såsom steilhet og kamhastighet fra endimensjonale tidsseriemålinger. Som følge av en langtidsanalyse av endimensjonal bølgestatistikk blir det vist at enrettet, langkammet og bratt sjø mest sannsynlig leder til ekstreme enkeltbølger med statistiske egenskaper som avviker systematisk fra ordinære statistiske modeller. Tredimensjonal, kortsiktig tid-rom-statistikk av ekstreme bølgekammer blir også undersøkt v.h.a. stereomålingene fra Ekofisk. Her blir det vist at statistiske modeller utvidet fra endimensjonale til tredimensjonale bølgefelt i snitt er velegnet til å beskrive forekomsten av de høyeste bølgekammene, spesielt for relativt store tid-rom segment.The research presented in this thesis characterizes statistical and dynamical aspects of dominant wind-generated surface gravity waves inferred from field observations in intermediate-to-deep water. Dominant waves are the most energetic waves in a sea state, and as such, understanding their behavior is important in both engineering and geophysical contexts. Large waves impart considerable impact forces on marine structures such as oil and gas platforms and offshore wind turbines, and these forces may multiply manyfold when waves break. Wave breaking in deep water, often referred to as whitecapping, is also a key, though incompletely understood, process regulating the transfer of momentum, gas and heat across the air-sea interface, and must thus be accurately parameterized in large-scale weather and climate models. Current theory holds that the wave breaking process is closely linked kinematically and dynamically to the group structure inherent in ocean surface wave fields. Wave group dynamics is also believed to govern the characteristic shape and motion of so-called extreme or rogue waves, whose correct statistical description is central to many offshore activities. The work presented herein shows, using state-of-the-art stereoscopic imaging techniques employed at the Ekofisk platform complex in the central North Sea, that large-scale wave breaking activity in the open ocean is strongly enhanced in dominant wave groups. The topic of wave group-modulated wave breaking has received considerable attention in the past two decades from theoretical, numerical and laboratory perspectives; however, quantitative field studies of the phenomenon remain comparatively rare. The current results also support the general notion that the dominant waves in a given sea state regulate the breaking of shorter waves. The statistics of extreme wave crest elevations is investigated using a novel long-term laser altimeter data set, also located at the Ekofisk field. The validity of the extreme values is verified using a newly developed despiking methodology, and the quality controlled data set, which covers storm events over an 18-year period, is used to investigate the effects of wave steepness and directionality on crest height statistics. Narrow directional spread combined with high wave steepness is found to lead to crest height statistics that deviate the most from standard linear and second-order formulations. Finally, geometric wave shape and crest speed dynamics are analyzed for the highest wave crests encountered in three-dimensional, spatially and temporally resolved segments of the stereo-reconstructed sea surface fields. The directly measured crest steepness is found to conform to the classical breaking limit of Stokes, whereas crest steepness estimated from one-dimensional time series measurements using the linear gravity-wave dispersion relation are systematically higher. This may be at least in part explained by the observation that the directly measured crest speed just before, during and after the moment of maximum crest elevation slows down compared to the linear gravity-wave phase speed estimate. For the first time, the crest speed slowdown is shown with field measurements to apply to both breaking and non-breaking dominant wave crests.Doktorgradsavhandlin

The variability and forcing of ocean whitecaps and their impact on air-sea fluxes

The breaking of ocean surface gravity waves is an important phenomenon that affects the dynamics of the upper ocean, development of the wave field, and air-sea exchange processes. As the surface expression of this process, whitecaps provide a visible signature of wave breaking; their areal extent per unit area sea surface - known as the whitecap fraction, W - can be used to quantify the amount and scale of wave breaking. W is traditionally estimated using digital images of the ocean surface and is widely used to represent whitecaps in remote sensing applications, and in the parameterisation of a host of air-sea processes in models. These parameterisations - generally functions of wind speed alone - are based on limited amounts of data, and fail to take into account the known influence of secondary factors on whitecaps. A novel approach to estimating W using satellite observations has recently been developed, based on passive radiometric measurements of brightness temperature at microwave frequencies. The satellite-based approach enables measurement of W on a global scale, and in a variety of conditions. In this work, the basic characteristics of W estimates at two different radiometric frequencies, W10 (10 GHz) and W37 (37 GHz), is investigated. The wind speed dependence, global distribution, and seasonal dependence of the estimates are investigated. Comparison is made against estimates obtained from simple, but widely used, wind speed only parameterisations formulated from in situ data. A more direct comparison of radiometric and photographic W estimates, based on ship-satellite matchups, is also made. Both comparisons indicate that satellite-based W has a different wind speed dependence, resulting in estimates that are, on average, higher at low wind speeds and lower at higher wind speeds than parameterisations formulated from in situ, photographic measurements. On a global scale, this results in satellite-based W being more uniform latitudinally than predictions from traditional formulations. A dataset comprising estimates of W10 and W37, together with collocated and concurrent estimates for a variety of forcing parameters, is used to investigate the the influence on W10 and W37 of secondary forcings, such as the wave field and environmental factors. It is found that on a global scale wind speed describes much of the variability in both W10 and W37 though the influence of secondary factors on W can be appreciable (especially for W37). Based on the magnitude of the influence of secondary forcing factors on W10 and W37, it is concluded that much of the variability in whitecap fraction is likely due to the behavior of the thinner, decaying foam patches, variability that is not captured by the retrieval using the 10 GHz channel. Though whitecap fraction offers a pragmatic approach to inferring the magnitude of processes associated with breaking surface waves, it remains an indirect measure with inherent limitations. More fundamental questions regarding the interpretation and use of W are considered. A dynamical model that relates whitecap fraction to breaking wave statistics is used to illustrate the contribution to whitecap fraction due to whitecaps in different lifetime stages. Such a model provides a framework for better relating whitecap fraction to the dynamic, active part of the wave breaking process which is likely more closely linked to processes such as breaking-induced energy dissipation, turbulent mixing, and bubble-mediated gas exchange. Finally, the implications of use of radiometric estimates for quantifying air-sea processes - specifically production of sea spray aerosol and bubble-mediated gas exchange - is discussed. It is shown that difference between the satellite-based W estimates and those predicted using traditional parameterisations provides an explanation for the consistent geographical biases in sea spray aerosol concentration found in a number of large scale models. The benefit of these novel observations will also extend to predictions of other air-sea processes, and remote sensing applications, that require estimation of W; these benefits will be enhanced if whitecaps and their radiometric signature are more closely related to the physical processes which they are used to quantify

Modeling Atmosphere-Ocean Radiative Transfer: A PACE Mission Perspective

The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth’s sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates

Modeling atmosphere-ocean radiative transfer: A PACE mission perspective

The research frontiers of radiative transfer (RT) in coupled atmosphere-ocean systems are explored to enable new science and specifically to support the upcoming Plankton, Aerosol, Cloud ocean Ecosystem (PACE) satellite mission. Given (i) the multitude of atmospheric and oceanic constituents at any given moment that each exhibits a large variety of physical and chemical properties and (ii) the diversity of light-matter interactions (scattering, absorption, and emission), tackling all outstanding RT aspects related to interpreting and/or simulating light reflected by atmosphere-ocean systems becomes impossible. Instead, we focus on both theoretical and experimental studies of RT topics important to the science threshold and goal questions of the PACE mission and the measurement capabilities of its instruments. We differentiate between (a) forward (FWD) RT studies that focus mainly on sensitivity to influencing variables and/or simulating data sets, and (b) inverse (INV) RT studies that also involve the retrieval of atmosphere and ocean parameters. Our topics cover (1) the ocean (i.e., water body): absorption and elastic/inelastic scattering by pure water (FWD RT) and models for scattering and absorption by particulates (FWD RT and INV RT); (2) the air-water interface: variations in ocean surface refractive index (INV RT) and in whitecap reflectance (INV RT); (3) the atmosphere: polarimetric and/or hyperspectral remote sensing of aerosols (INV RT) and of gases (FWD RT); and (4) atmosphere-ocean systems: benchmark comparisons, impact of the Earth’s sphericity and adjacency effects on space-borne observations, and scattering in the ultraviolet regime (FWD RT). We provide for each topic a summary of past relevant (heritage) work, followed by a discussion (for unresolved questions) and RT updates

Perspectives and Integration in SOLAS Science

Why a chapter on Perspectives and Integration in SOLAS Science in this book? SOLAS science by its nature deals with interactions that occur: across a wide spectrum of time and space scales, involve gases and particles, between the ocean and the atmosphere, across many disciplines including chemistry, biology, optics, physics, mathematics, computing, socio-economics and consequently interactions between many different scientists and across scientific generations. This chapter provides a guide through the remarkable diversity of cross-cutting approaches and tools in the gigantic puzzle of the SOLAS realm. Here we overview the existing prime components of atmospheric and oceanic observing systems, with the acquisition of ocean–atmosphere observables either from in situ or from satellites, the rich hierarchy of models to test our knowledge of Earth System functioning, and the tremendous efforts accomplished over the last decade within the COST Action 735 and SOLAS Integration project frameworks to understand, as best we can, the current physical and biogeochemical state of the atmosphere and ocean commons. A few SOLAS integrative studies illustrate the full meaning of interactions, paving the way for even tighter connections between thematic fields. Ultimately, SOLAS research will also develop with an enhanced consideration of societal demand while preserving fundamental research coherency. The exchange of energy, gases and particles across the air-sea interface is controlled by a variety of biological, chemical and physical processes that operate across broad spatial and temporal scales. These processes influence the composition, biogeochemical and chemical properties of both the oceanic and atmospheric boundary layers and ultimately shape the Earth system response to climate and environmental change, as detailed in the previous four chapters. In this cross-cutting chapter we present some of the SOLAS achievements over the last decade in terms of integration, upscaling observational information from process-oriented studies and expeditionary research with key tools such as remote sensing and modelling. Here we do not pretend to encompass the entire legacy of SOLAS efforts but rather offer a selective view of some of the major integrative SOLAS studies that combined available pieces of the immense jigsaw puzzle. These include, for instance, COST efforts to build up global climatologies of SOLAS relevant parameters such as dimethyl sulphide, interconnection between volcanic ash and ecosystem response in the eastern subarctic North Pacific, optimal strategy to derive basin-scale CO2 uptake with good precision, or significant reduction of the uncertainties in sea-salt aerosol source functions. Predicting the future trajectory of Earth’s climate and habitability is the main task ahead. Some possible routes for the SOLAS scientific community to reach this overarching goal conclude the chapter

Hyperspectral Measurements, Parameterizations, and Atmospheric Correction of Whitecaps and Foam From Visible to Shortwave Infrared for Ocean Color Remote Sensing

Breaking waves are highly reflective features on the sea surface that change the spectral properties of the ocean surface in both magnitude and spectral shape. Here, hyperspectral reflectance measurements of whitecaps from 400 to 2,500 nm were taken in Long Island Sound, USA of natural and manufactured breaking waves to explore new methods to estimate whitecap contributions to ocean color imagery. Whitecap reflectance was on average ~40% in visible wavelengths and decreased significantly into the near infrared and shortwave infrared following published trends. The spectral shape was well-characterized by a third order polynomial function of liquid water absorption that can be incorporated into coupled ocean-atmospheric models and spectral optimization routines. Localized troughs in whitecap reflectance correspond to peaks in liquid water absorption and depths of the troughs are correlated to the amount and intensity of the breaking waves. Specifically, baseline-corrected band depths at 980 and 1,200 nm explained 77 and 90% of the whitecap-enhanced reflectance on a logarithmic scale, respectively. Including these wavebands into future ocean color sensors could potentially provide new tools to estimate whitecap contributions to reflectance more accurately than with wind speed. An effective whitecap factor was defined as the optical enhancements within a pixel due to whitecaps and foam independent of spatial scale. A simple mixed-pixel model of whitecap and background reflectance explained as much of the variability in measured reflectance as more complex models incorporating semi-transparent layers of foam. Using an example atmosphere, enhanced radiance from whitecaps was detectable at the top of the atmosphere and a multiple regression of at-sensor radiance at 880, 1,038, 1,250, and 1,615 nm explained 99% of the variability in whitecap factor. A proposed model of whitecap-free reflectance includes contributions from water-leaving radiance, glint, and diffuse reflected skylight. The epsilon ratio at 753 and 869 nm commonly used for aerosol model selection is nearly invariant with whitecap factor compared to the ratio at shortwave infrared bands. While more validation data is needed, this research suggests several promising avenues to retrieve estimates of the whitecap reflectance and to use ocean color to further elucidate the physics of wave breaking and gas exchange