An improved FSM method is used in geophysical and environmental corrections to enhance the final product of the along track Jason-2 SLA data and extend it near the Red Sea borders. In this study the ionospheric correction range, wet tropospheric correction range, sea state bias correction range and dry tropospheric correction range are enhanced and improved using FSM01, which helped to retrieve three more tracks (106, 170 and 234), earlier neglected by the distribution centers, and extend the tracks towards the coast. The FSM01 SLA is compared with Jason-2 SLA and AVISO SLA for the available 5 tracks, in which the FSM01 SLA show a good agreement and higher correlation with the Jason-2 SLA compared with that of AVISO, in addition to that it fills the gaps in the times series of all tracks. The new retrieved tracks also compared with those retrieved by AVISO, both data show similar variability, with FSM01 SLA show no gaps in the time series. The FSM01 SLA also extended towards the coast and show high correlation with the coastal tide measurements.

Three Dimensional Variational data assimilation or analysis (3DVAR) is one of most classical methods for providing the initial values for numerical models. In this method, the dimensions of the background error covariance and the observational error covariance matrices are large. Therefore, it is difficult to get the inverse of the covariance matrices and to reduce the orders of these matrices without information loss. With the use of the Sylvester Equation, on the basis of a new linear regression, a new cost function for 3DVAR was given. For the first-guess m×n field, there is an approximate 1−(m²+n²)/(mn×mn) reduction with m>1 & n>1 by using the cost function. The results of the numerical experiments show that the effect of this algorithm is no worse than that of the old cost function for 3DVAR.

Deterministic-stochastic empirical mode decomposition (EMD) is used to obtain low-frequency (non-diffusive, i.e., background velocity) and high-frequency (diffusive, i.e., eddies) components from a Lagrangian drifter‘s trajectory. Eddy characteristics are determined from the time series of eddy trajectories from individual Lagrangian drifter such as the eddy radial scale, eddy velocity scale, eddy Rossby number, and eddy-background kinetic energy ratio. A long-term dataset of the SOund Fixing And Ranging float time series obtained near the California coast by the Naval Postgraduate School from 1992 to 2004 at depth between 150 and 600 m (http://www.oc.nps.edu/npsRAFOS/) is used as an example to demonstrate the capability of the deterministic-stochastic EMD.

The aim of this study is to model hydrodynamic processes of the Istanbul Strait with its stratified flow characteristic and calibrate the most important parameters using local and global search algorithms. For that two open boundary conditions are defined, which are in the North and South part of the Strait. Observed bathymetric, hydrographic, meteorological and water level data are used to set up the Delft3D-FLOW model. First, the sensitivities of model parameters on the numerical model outputs are assessed using PEST toolbox. Then, the model is calibrated based on the objective functions focusing on the flowrates of upper and lower layers. The salinity and temperature profiles of the Strait are only used for model validation. The results show that the calibrated model outputs of Istanbul Strait are reliable and consistent with the in-situ measurements. The sensitivity analysis reveals that the Spatial Low-Pass Filter Coefficient, Horizontal Eddy Viscosity, Prandtl-Schmidt Number, Slope in log-log Spectrum and Manning Roughness Coefficient are most sensitive parameters affecting flowrate performance of the model. The agreement between observed salinity profiles and simulated model outputs is promising whereas the match between observed and simulated temperature profiles is weak showing that the model can be improved particularly for simulating the mixing layer.

This study utilizes repeated geoacoustic mapping to quantify the morphodynamic response of the nearshore to storm-induced changes. The aim of this study was to quantitatively map the nearshore zone of Assateague Island National Seashore (ASIS) to determine what changes in bottom sediments, benthic fauna and fish habitat are attributable to storm events including hurricane Sandy and the passage of hurricane Joaquin. Specifically, (1) the entire domain of the National Parks Service offshore area was mapped with side-scan sonar and multibeam bathymetry at a resolution comparable to that of the existing pre-storm survey, (2) a subset of the benthic stations were resampled that represented all sediment strata previously identified, and (3) newly obtained data were compared to that from the pre-storm survey to determined changes that could be attributed to specific storms such as Sandy and Joaquin. Capturing event specific dynamics requires rapid response surveys in close temporal association of the before and after period. The time-lapse between the pre-storm surveys for Sandy and our study meant that only a time and storm integrated signature for that storm could be obtained whereas with hurricane Joaquin we could identify impacts to the habitat type and geomorphology more directly related to that particular storm. This storm impacts study provides for the National Park Service direct documentation of storm-related changes in sediments and marine habitats on multiple scales: from large scale, side-scan sonar maps and interpretation of acoustic bottom types, to characterize as fully as possible habitats from 1 to 10 m up to many kilometer scales, as well as from point benthic samples within each sediment stratum and these results can help guide management of the island resources.

The effect of river fronts on oil slick transport has been demonstrated using high resolution forcing models and a fully fledged oil drift model, OpenOil. The model system is used to simulate the 2010 DeepWater Horizon oil spill. Metocean forcing data are taken from the GoM-HYCOM 1/50° ocean model with realistic river input and ECMWF global forecast products of wind and wave parameters with 1/8° resolution. The simulations are initialized from satellite observations of the surface oil patch. OpenOil includes most of the relevant processes, such as emulsification, evaporation, wave entrainment, stranding and droplet formation. The model takes account of the actual oil type and properties, using the ADIOS oil weathering database of NOAA. The effect of using a newly developed parameterization for oil droplet size distribution is studied and compared to a traditional algorithm. Although the algorithms provide different distributions for a single wave breaking event, it is found that the net difference after long simulations is negligible, indicating that the outcome is robust regarding the choice of parameterization. That indicates that the wave entrainment, vertical mixing and re-surfacing mechanisms that are part of OpenOil are more important for determining the final droplet size spectrum than the spectrum prescribed for individual wave breaking events. In both cases, the size of the droplets controls how much oil is present at the surface and hence are subject to wind and Stokes drift. The effect of removing river outflow in the ocean model is investigated in order to showcase effects of river induced fronts on oil spreading. A consistent effect on the amount and location of stranded oil is found, and considerable impact of river induced fronts is seen on the location of the surface oil patch. During a case with large river outflow (May 20-27, 2010), the total amount of stranded oil is reduced by about 50% in the simulation with no river input. The results compares well with satellite observations of the surface oil patch.

Data from environmental satellites are underutilized in many branches of operational oceanography. Several challenges exist for potential users of satellite data that may impede or preclude employing satellite products in their work. Users outside of the satellite community often encounter difficulty in discovering the types of satellite measurements that are available, and determining which satellite products are best for operational activities. In addition, the large choice of satellite data providers, each with their own data access protocols and formats, can make data access confusing and inefficient. NOAA’s CoastWatch Program is designed to be a Value Added Provider, whose mission is to help make satellite data easier for users to access. With this objective the West Coast Node of CoastWatch developed the NOAA Ocean Satellite Course, which introduces scientists and resource managers to ocean satellite products, and provides them tools to facilitate data access when using ArcGIS or R software. These tools leverage upon the data services provided by ERDDAP, a data distribution system designed to make data access easier for both humans and machines. The course has been offered annually since 2006, and over 350 participants have taken it. Results of post-course surveys are analyzed to measure course effectiveness. Lessons learned from conducting these courses are summarized.

The aim of work is to derive an explicit expression for a function of vertical mixing induced by wind-waves. To this end, in the Navier-Stokes equations, a current is decomposed into four constituents: the mean flow, the wave-orbital motion, the wave-induced turbulent and the background turbulent currents. This decomposition allows separating the wave-induced Reynolds stress, Rw, from the background one, Rb. To make a statistical closure for Rw, the Prandtl approach for the background turbulent fluctuations is used that results in an implicit expression for the wave-induced vertical mixing function, Bv. Expression for Bv is specified based on the author’s results for the eddy viscosity found earlier in the frame of the three-layer concept for a wavy air–sea interface, used for modelling wind-drift currents [1]. Finally, the explicit parameterization for Bv(a, u*, z) is found as a linear function in both the wave amplitude at depth z, a(z), and the friction velocity in the air, u*. The linear dependence of function Bv(a) on the wave amplitude provides the enhanced vertical mixing induced by wind–waves in comparison with function Bv(a) having the cubic dependence found in [2], as far as the wind-wave amplitude a(z) decays exponentially with depth.

Traditional ways of validating satellite-derived sea surface temperature (SST) and sea surface salinity (SSS) products, using comparisons with buoy measurements, do not allow for evaluating the impact of mesoscale to submesoscale variability. Here we present the validation of remotely-sensed SST and SSS data against the unmanned surface vehicle (USV) – Saildrone – measurements from the Spring 2018 Baja deployment. More specifically, biases and root mean square differences (RMSD) were calculated between USV-derived SST and SSS values, and six satellite-derived SST (MUR, OSTIA, CMC, K10, REMSS, and DMI) and three SSS (JPLSMAP, RSS40, RSS70) products. Biases between the USV SST and OSTIA/CMC/DMI were approximately zero while MUR showed a bias of 0.2C. OSTIA showed the smallest RMSD of 0.36C while DMI had the largest RMSD of 0.5C. An RMSD of 0.4C between Saildrone SST and the satellite-derived products could be explained by the daily variability in USV SST which currently cannot be resolved by remote sensing measurements. For SSS, values from the JPLSMAP product showed saltier biases of 0.2 PSU, while RSS40 and RSS70 showed fresh biases of 0.3 PSU. An RMSD of 0.4 PSU could not be explained solely by the daily variability of the USV-derived SSS. Coherences were significant at the longer wavelengths, with a local maximum at 100 km that is most likely associated with the mesoscale turbulence in the California Current System.

This research aims to estimate the micro-motion (m-m) of ships. The problem of motion and m-m detection of targets is usually solved using synthetic aperture radar (SAR) along-track interferometry (ATI) which is observed employing two radars spatially distanced by a baseline extended in the azimuth direction. This paper is proposing a new approach where the m-m estimation of ships, occupying thousands of pixels, is measured processing the information given by sub-pixel tracking generated during the coregistration process of several re-synthesized time-domain and overlapped sub-apertures. The SAR products are generated splitting the raw data, according to a small-temporal baseline strategy, observed by one single wide-band staring spotlight (ST) SAR image. The predominant vibrational modes of different ships are estimated and results are promising to extend this application in performing surveillance also of land-based industries activities. Experiments are performed processing one ST SAR image observed by the COSMO-SkyMed satellite system.

The main objective of this paper was description of seasonal and interannual trends in secondary production and mortality rates of the three most important Copepoda taxa in the Gulf of Gdańsk (southern Baltic Sea). Samples were collected monthly from 6 stations located in the western part of the Gulf of Gdańsk during three research periods: 1998-2000, 2006-2007 and 2010-2012. Production was computed basing on copepod biomass and mortality rates estimated according to vertical life table approach. Redundancy analysis was used to investigate relationship between secondary production and environmental conditions. Considering the entire research period there was significant interannual and seasonal variability of secondary production, mortality rate as well as abundance and biomass anomalies. Conducted analysis revealed correlation between increasing temperature and production of Acartia spp. and T.longicornis developmental stages, while older copepodites of P.acuspes showed almost negative correlation with temperature. The mortality rate estimations obtained for Acartia spp. Were highest in summer, while for T.longicornis peak was usually noted in spring-summer period. Lowest mortality rate estimations were noted in autumn and winter for almost all stages of investigated taxa.

We apply deep convolutional neural networks (CNNs) to estimate wave breaking type from close-range monochrome infrared imagery of the surf zone. Image features are extracted using six popular CNN architectures developed for generic image feature extraction. Logistic regression on these features is then used to classify breaker type. The six CNN-based models are compared without and with augmentation, a process that creates larger training datasets using random image transformations. The simplest model performs optimally, achieving average classification accuracies of 89% and 93%, without and with image augmentation respectively. Without augmentation, average classification accuracies vary substantially with CNN model. With augmentation, sensitivity to model choice is minimized. A class activation analysis reveals the relative importance of image features to a given classification. During its passage, the front face and crest of a spilling breaker are more important than the back face. Whereas for a plunging breaker, the crest and back face of the wave are most important. This suggests that CNN-based models utilize the distinctive `streak' temperature patterns observed on the back face of plunging breakers for classification.

In this paper, sea surface wind speed (SSWS) retrieval from Gaofen-3 (GF-3) quad-polarization stripmap (QPS) data in vertical-vertical (VV), horizontal-horizontal (HH) and vertical-horizontal (VH) polarizations is investigated in detail based on 3,170 scenes acquired from October 2016 to May 2018. The radiometric calibration factor of the VV polarization data is examined first. This calibration factor generally meets the requirement of SSWS retrieval accuracy with an absolute bias of less than 0.5 m/s but shows highly dispersed characteristics. These results lead to SSWS retrievals with a small bias of 0.18 m/s but a rather high root mean square error (RMSE) of 2.36 m/s compared with the ERA-Interim reanalysis model data. Two refitted polarization ratio (PR) models for the QPS HH polarization data are presented. Based on a combination of the incidence angle- and azimuth angle-dependent PR model and CMOD5.N, the SSWS derived from the QPS HH data shows a bias of 0.07 m/s and an RMSE of 2.26 m/s relative to the ERA-Interim reanalysis model wind speed. A linear function relating SSWS and the normalized radar cross section (NRCS) of QPS VH data is derived. The SSWS data retrieved from the QPS VH data show good agreement with the WindSat SSWS data, with a bias of 0.1 m/s and an RMSE of 2.02 m/s. We also apply the linear function to the GF-3 Wide ScanSAR data acquired for the typhoon SOULIK, which surprisingly yields a very good agreement with the model results. A comparison of SSWS retrievals among three different polarization datasets is also presented. The current study and our previous work demonstrate that the general accuracy of the SSWS retrieval based on GF-3 QPS data has an absolute bias of less than 0.3 m/s and an RMSE of 2.0 ±0.2 m/s relative to various datasets. Further improvement will depend on dedicated radiometric calibration efforts.

This study was undertaken to derive and analyze the Advanced Microwave Scanning Radiometer - EOS (AMSR-E) sea surface temperature (SST) footprint associated with the Remote Sensing Systems (RSS) Level-2 (L2) product. The footprint, in this case, is characterized by the weight attributed to each 4  4 km square contributing to the SST value of a given AMSR-E pixel. High-resolution L2 SST fields obtained from the MODerate-resolution Imaging Spectroradiometer (MODIS), carried on the same spacecraft as AMSR-E, are used as the sub-resolution “ground truth“ from which the AMSR-E footprint is determined. Mathematically, the approach is equivalent to a linear inversion problem, and its solution is pursued by means of a constrained least square approximation based on the bootstrap sampling procedure. The method yielded an elliptic-like Gaussian kernel with an aspect ratio  1.58, very close to the AMSR-E 6.93GHz channel aspect ratio,  1.7. (The 6.93GHz channel is the primary spectral frequency used to determine SST.) The semi-major axis of the estimated footprint is found to be alignedwith the instantaneous field-of-view of the sensor as expected fromthe geometric characteristics of AMSR-E. Footprintswere also analyzed year-by-year and as a function of latitude and found to be stable – no dependence on latitude or on time. Precise knowledge of the footprint is central for any satellite-derived product characterization and, in particular, for efforts to deconvolve the heavily oversampled AMSR-E SST fields and for studies devoted to product validation and comparison. A preliminarly analysis suggests that use of the derived footprint will reduce the variance between AMSR-E and MODIS fields compared to the results obtained.

Short-term (hours) biological and biogeochemical processes cannot be fully captured by the current suite of polar-orbiting satellite ocean color sensors, as their temporal resolution is limited to potentially one clear image per day. Geostationary sensors, such as the Geostationary Ocean Color Imager (GOCI) from the Republic of Korea, allow the study of these short-term processes because their geostationary orbits permit the collection of multiple images throughout each day. To assess the capability to detect changes in water properties caused by these processes, however, requires an understaning of the uncertainties introduced by the instrument and/or geophysical retrieval algorithms. This work presents a study of the variability during the day over a water region of low-productivity with the assumption that only small changes in the water properties occur during the day over the area of study. The complete GOCI mission data were processed using the SeaDAS/l2gen package. Filtering criteria were applied to assure the quality of the data. Relative differences with respect to the midday value were calculated for each hourly observation of the day. Also, the influence of the solar zenith angle in the retrieval of remote sensing reflectances and derived products was analyzed. We determined that the uncertainties in water-leaving “remote-sensing” reflectance ($R_\text{rs}$) for the 412, 443, 490, 555, 660 and 680 nm bands on GOCI are 8.05$\times10^{-4}$, 5.49$\times10^{-4}$, 4.48$\times10^{-4}$, 2.51$\times10^{-4}$, 8.83$\times10^{-5}$, and 1.36$\times10^{-4}$ sr$^{-1}$, respectively, and 1.09$\times10^{-2}$ mg m$^{-3}$ for the chlorophyll-a concentration (Chl-{\it a}), 2.09$\times10^{-3}$ m$^{-1}$ for the absorption coefficient of chromophoric dissolved organic matter at 412 nm ($a_{\text{g}}(412)$), and 3.7 mg m$^{-3}$ for particulate organic carbon (POC). We consider these to be the floor values for detectable changes in the water properties due to biological, physical or chemical processes.

The present work focuses on fresh water signatures at the sediment-water interface (1 cm) using foraminiferal species in both austral winter and summer in eleven longitudinal transects on the Western South Atlantic continental margin between 27° and 37° S, at water depths of 11.7 to 250 m. Here we show that depth, salinity, temperature, oxygen, grain size (mud and sand percentage), suspended matter, organic matter, SiO₄, NO₂, and NO₃ in this order of importance are responsible for the distribution of foraminiferal species and thecamoebians. The presence of these microfossils indicate freshwater influx in four sectors over the continental shelf: Itajaí-Açu River, Laguna estuarine system, Patos Lagoon and RdlP (Rio de la Plata) will be explored further in detail. Our findings on freshwater signature on the continental shelf sediments through benthic species indicator are comparable to other continental systems worldwide, and a paleo record study would be useful for three South American countries (Brazil, Argentina and Uruguay). A freshwater signature in the continental shelf indicates deposition sites probably linked to anthropogenic impact since most of the pollutants and contaminants are dumped into water bodies that eventually reach and accumulate in the ocean. Therefore, the freshwater-related species on the continental shelf reflects exactly where the depositional sediment sites are, and where anthropogenic impacts accumulate. Foraminiferal microhabitat occupation within these zones is discussed in detail. And we conclude that together with the fauna, the abiotic parameters play an important role in determining the occurrence and degree of marine eutrophication induced by the input of polluted river waters, also showing possible anoxic environments on the shelf.

Using a three-dimensional, hydrostatic, primitive-equation ocean model, this study investigates the dynamics of lateral circulation in a partially stratified tidal inlet, the Barataria Pass in the Gulf of Mexico, over a 25.6-hour diurnal tidal cycle. Model performance is assessed against observational data. During flood tide, the lateral circulation exhibits the characteristics similar to that induced by differential advection, i.e., lateral flow consists of two counter-rotating cells and is convergent at the surface. The analysis of momentum balance indicates that, in addition to the pressure gradient and vertical stress divergence, nonlinear advection and horizontal stress divergence are also important contributors. During ebb phase, the lateral circulation is mostly eastward for the whole water column and persisting for almost the whole period. The surface divergence suggested by the differential advection mechanism lasts for a very short period, if it ever exists. The main momentum balance across most of the transect during ebb is between the along-channel advection of cross-channel momentum and pressure gradient. The sectional averaged lateral velocity magnitude during ebb is comparable to that during flood, which is different from the idealized numerical experiment result.

Pencil-beam Doppler scatterometers are a promising remote sensing tool for measuring ocean vector winds and currents from space. While several point designs exist in the literature, these designs have been constrained by the hardware they inherited, and the design is sub-optimal. Here, I present guidelines to optimize the design of these instruments starting from the basic sensitivity equations. Unlike conventional scatterometers or pencil-beam imagers, appropriate sampling of the Doppler spectrum and optimizing the radial velocity error lead naturally to a design that incorporates a pulse-to-pulse separation and pulse length that vary with scan angle. Including this variation can improve radial velocity performance significantly and the optimal selection of system timing and bandwidth is derived. Following this, optimization of the performance based on frequency, incidence angle, antenna length, and spatial sampling strategy are considered. It is shown that antenna length influences the performance most strongly, while the errors depend only on the square root of the transmit transmit power. Finally, a set of example designs and associated performance are presented.

Current numerical methods for simulating biophysical processes in aquatic environments are typically constructed in a grid-based Eulerian framework or as an individual-based model in a particle-based Lagrangian framework. Often, the biogeochemical processes and physical (hydrodynamic) processes occur at different time and space scales, and changes in biological processes do not affect the hydrodynamic conditions. Therefore, it is possible to develop an alternative strategy to grid-based approaches for linking hydrodynamic and biogeochemical models that can significantly improve computational efficiency for this type of linked biophysical model. In this work, we utilize a new technique which links hydrodynamic effects and biological processes through a property-carrying particle model (PCPM) in a Lagrangian/Eulerian framework. The model is tested in idealized cases and its utility is demonstrated in a practical application to Sandusky Bay. Results show the integration of Lagrangian and Eulerian approaches allows for a natural coupling of mass transport (represented by particle movements and random walk) and biological processes in water columns which is described by a nutrient-phytoplankton-zooplankton-detritus (NPZD) biological model. This method is far more efficient than traditional tracer based Eulerian biophysical models for 3-D simulation, particularly for a large domain and/or ensemble simulations.

Measurements of optical properties have been used for decades to study particle distributions in the ocean. They have been found useful to constrain suspended mass concentration as well as particle-related properties such as size, composition, packing (particle porosity or density) and settling velocity. Optical properties, however, provide measurements that are biased, as certain particles (based on size, composition, shape or packing) contribute to a specific property more than others. Here we study this issue both theoretically as well as by contrasting different optical properties collected simultaneously in a bottom boundary layer, to highlight the utility of such measurements as well as the biases we are likely to encounter using different optical properties to study suspended particles. In particular, we investigate the possibility to infer settling velocity from vertical profiles of optical measurements, finding that the effects of aggregation dynamics can seldom be ignored.

As human activities of the countries in the East Asia have been remarkably expanding over recent decades, various problems in relation to ships, such as oil spill and many other coastal marine pollution, are continuously occurring in the coastal region. In order to conserve marine resources and prepare for possible ship accidents in advance, the need for efficient ship management is increasing over time. Multi-satellite, multi-sensor, multi-wavelength or multi-frequency observations make it possible to monitor a variety of vessels in the coastal region. This study presents the results of ship detection methodology applied to multi-spectral satellite images in the seas around Korean Peninsula based on optical, hyperspectral, and microwave remote sensing. To detect ships from hyperspectral images with a few hundreds of spectral channels, spectral matching algorithms are used to investigate similarity between the spectra and in-situ measurements. In the case of SAR (Synthetic Aperture Radar) images, the Constant False Alarm Rate (CFAR) algorithm is used to discriminate the vessels from backscattering coefficients of Sentinel-1 SAR and ALOS-2 PALSAR2 images. The present ship detection methods can be extensively utilized for optical, hyperspectral, and SAR images for comprehensive coastal management purposes toward perpetual sustainability in the future.

Studying the interaction between the upper ocean and the typhoons is crucial to improve our understanding of heat and momentum exchange between the ocean and the atmosphere. In recent years, the upper ocean responses to typhoons have received considerable attention. The sea surface cooling (SSC) process has been repeatedly discussed. In the present work, case studies were examined on five strong and super typhoons that occurred in 2016—LionRock, 1610; Meranti, 1614; Malakas, 1616; Megi, 1617; and Chaba 1618—to search for more evidence and new features of typhoon’s impact on the sea surface environment. The typhoon monitoring data from the Central Meteorological Observatory, the sea surface temperature (SST) data from satellite microwave and infrared remote sensing, and the sea surface height anomaly (SSHA) data from satellite altimeters were used to analyze in detail: the SSC features caused by typhoons, the relationship between the SSC and the typhoon travelling speed, and the variations in cold and warm eddies during typhoon passage. Results showed that: (1) SSC generally occurred during typhoon passage and the degree of SSC was always determined by the strength and the travelling speed of the typhoon, as well as the initial SST. (2) One day before or on the day of typhoon passage, the SSHA slightly increased due to low surface pressure. After the typhoon passed, the SSHA obviously decreased along with the SSC. The pre-existing positive SSHAs, which always represent warm eddies, decreased or disappeared during typhoon passage, whereas negative SSHAs or cold eddies were enhanced. (3) New cold eddies were generated, especially at the turning points of the typhoon path. The presence of warm eddies is suggested to have a strengthening effect on the typhoons.

Data from NASA’s Soil Moisture Active Passive Mission (SMAP) and from the California Cooperative Oceanic Fisheries Investigations (CalCOFI) were used to examine the freshening that occurred during 2015-2016 in the Southern California Current System. Overall the freshening was found to be related to the 2014-2016 Northeast Pacific Warm Anomaly. The primary goal was to determine the feasibility of using SMAP data to observe the surface salinity signal associated with the warming. As a first step direct comparisons were done with salinity from the CalCOFI data at one-meter depth. During 2015 SMAP was saltier than CalCOFI by 0.5 PSU, but biases were reduced to < 0.1 PSU during 2016. South of 33°N, and within 100 km of the coast, SMAP was fresher in 2015 by almost 0.2 PSU. CalCOFI showed freshening of 0.1 PSU. North of 33°N SMAP and CalCOFI saw significant freshening in 2016, SMAP by 0.4 PSU and CalCOFI by 0.2 PSU. Differences between SMAP and CalCOFI are consistent with the increased stratification in 2015 and changes in the mixed layer depth.

Wave energy converter (WEC) arrays deployed in coastal regions may create physical disturbances potentially resulting in environmental stresses. Presently, limited information is available on the nature of these physical disturbance or the resultant effects. A quantitative Spatial Environmental Assessment Tool (SEAT) for evaluating potential effects of wave energy converter (WEC) arrays on nearshore hydrodynamics and sediment transport is presented for the central Oregon coast (USA) through coupled numerical model simulations of an array of WECs. Derived climatological wave conditions were used as inputs to the model to allow for the calculation of risk metrics associated with various hydrodynamic and sediment transport variables such as maximum shear stress, bottom velocity, and change in bed elevation. The risk maps provided simple, quantitative, and spatially-resolved means of evaluating physical changes in the vicinity of a hypothetical WEC array in response to varying wave conditions. Near-field risk of sediment mobility was determined to be moderate in the lee of the densely spaced array, where the potential for increased sediment deposition could result in benthic habitat alteration. Modifications to the nearshore sediment deposition and erosion patterns were observed near headlands and topographic features, which could have implications for littoral sediment transport. The results illustrate the benefits of a risk evaluation tool for facilitating coastal resource management at early market marine renewable energy sites.

There is increasing interdisciplinary interest in phytoplankton community dynamics as the growing environmental problems of water quality (particularly eutrophication) and climate change demand attention. This has led to a pressing need for improved biophysical and causal understanding of Phytoplankton Functional Type (PFT) optical signals, in order that satellite radiometry may be used to detect ecologically relevant phytoplankton assemblage changes. This understanding can best be achieved with biophysically and biogeochemically consistent phytoplankton Inherent Optical Property (IOP) models, as it is only via modelling that phytoplankton assemblage characteristics can be examined systematically in relation to the bulk optical waterleaving signal. The Equivalent Algal Populations (EAP) model is used here to investigate the source and magnitude of size- and pigment- driven PFT signals in the water-leaving reflectance, as well as the potential to detect these using satellite radiometry. This model places emphasis on explicit biophysical modelling of the phytoplankton population as a holistic determinant of IOPs, and a distinctive attribute is its comprehensive handling of the spectral and angular character of phytoplankton scattering. Selected case studies and sensitivity analyses reveal that phytoplankton spectral scattering is the primary driver of the PFT-related signal. Key findings are that the backscattering-driven signal in the 520 to 600 nm region is the critical PFT identifier at marginal biomass, and that while PFT information does appear at blue and red wavelengths, it is compromised by biomass/gelbstoff ambiguity in the blue and low signal in the red, due primarily to absorption by water. These findings are hoped to provide considerable insight into the next generation of PFT algorithms.

Ship speeds extracted from AIS data vary with ice conditions. We extrapolated this variation with SAR data to a chart of expected icegoing speed. The study is for the Gulf of Bothnia in March 2013 and for ships with ice class 1A Super that are able to navigate without icbreaker assistance. The speed was normalized to 0-10 for each ship. As the matching between AIS and SAR was complicated by ice drift during the time gap, from hours to two days, we calculated a set of local SAR statistics over several scales. We used random tree regression to estimate the speed. The accuracy was quantified by mean squared error (MSE), and the fraction of estimates close to the actual speeds. These depended strongly on the route and the day. MSE varied from 0.4 to 2.7 units² for daily routes. 65 % of the estimates deviated less than one unit and 82 % less than 1.5 units from the AIS speeds. The estimated daily mean speeds were close to the observations. Largest speed decreases were provided by the estimator in a dampened form or not at all. This improved when ice chart thickness was included as one predictor.

A method for solving the turbulence equations embedded in the sigma ocean general ocean circulation model is proposed. Like the general circulation model, the turbulence equations are solved using the splitting method by physical processes. The turbulence equations are split into two main stages describing transport-diffusion and generation-dissipation processes. Parameterization of turbulence in the framework of equations allows, at the generation-dissipation stage, to use both numerical and analytical solutions and to ensure high efficiency of the algorithm. The results of large-scale ocean dynamics simulation taking into account the parameterization of vertical turbulent exchange are considered. Numerical experiments were carried out using k-omega turbulence model embedded to the Institute of Numerical Mathematics Ocean general circulation Model (INMOM). Both the circulation and turbulence models are solved using the splitting method with respect to physical processes. The coupled model is used to simulate the hydrophysical fields of the North Atlantic and Arctic Oceans for 1948--2009. The model has a horizontal resolution of 0.25 degree and 40 sigma-levels along the vertical. The sensitivity of the solution to the changes in mixing parameterization is studied. Experiments demonstrate that taking into account the climatic annual mean buoyancy frequency improves the reproduction of large-scale ocean characteristics. There is a positive effect of Prandtl number variations for reproducing the upper mixed layer depth. The experiments also demonstrate the computational effectiveness of the proposed approach in solving the turbulence equations.

Climate change such as sea level rise, change in temperature, precipitation, and storminess are expected to impact significantly coastal lagoons. The nature and magnitude of these impacts are uncertain. The objective of the research is to determine the climate change impacts on mixing and circulation at Songkhla lagoon, Thailand. Songkhla lagoon is the largest lagoonal water resource in Thailand and Southeast Asia. The lagoon is a combined freshwater and estuarine complex of high productivity which represents an extraordinary combination of environmental resources believed to be unique in the region. This work is part of a Climate Change impact assessment framework. It is the validation phase (step 5) of the framework applying a case study. Delft 3D was used to simulate CC scenarios in the climate downscaling models, part of the previous framework steps. These results were compared to the current conditions to determine the main changes in mixing and circulation in the coastal lagoon. Three indicators were applied to quantify the impacts: flushing time, salinity intrusion and stratification. The results suggest an increase in water velocities at the inlet in future scenarios and a decrease of flushing time. Salinity and stratification showed more complex changes in futures scenarios.

Parametric cyclonic wind fields are widely used worldwide for insurance risk underwriting, coastal planning, or storm surge forecasts. They support high-stakes financial, development, and emergency decisions. Yet, there is still no consensus on the best parametric approach, or relevant guidance to choose among the great variety of published models. The aim of this paper is first and foremost to demonstrate that recent progresses on estimating extreme surface wind speeds from satellite remote sensing now makes it possible to select the best option with greater objectivity. In particular, we show that the Cyclone Global Navigation Satellite System (CYGNSS) mission of NASA is able to capture a substantial part of the tropical cyclones structure, and allows identifying systematic biases in a number of parametric models. Our results also suggest that none of the traditional empirical approaches can be considered as the best option in all cases. Rather, the choice of a parametric model depends on several criteria such as cyclone intensity and/or availability of wind radii information. The benefit of using satellite remote sensing data to better select a parametric model for a specific case study is tested here by simulating hurricane Maria (2017). The significant wave heights computed by a wave-current hydrodynamic coupled model are found to be in good accordance with the predictions given by the remote sensing data in terms of bias. The results and approach presented in this study should shed new light on how to handle parametric cyclonic wind models, and help the scientific community to conduct better wind, waves and surge analysis for tropical cyclones.

Observed sea surface Ka-band normalized radar backscatter cross section (NRCS) and Doppler velocity (DV) exhibit energy at low frequencies (LF) below the surface wave range. It is shown that non-linearity in NRCS-wave slope Modulation Transfer Function (MTF) and inherent NRCS averaging within the footprint account for the NRCS and DV LF variance with the exception of VV NRCS for which almost half of the LF variance is attributable to wind fluctuations. Although the distribution of radar DV is quasi-Gaussian suggesting virtually little impact of non-linearity, the LF DV variations arise due to footprint averaging of correlated local DV and non-linear NRCS. Numerical simulations demonstrate that MTF non-linearity weakly affects traditional linear MTF estimate (less than 10% for |MTF|< 20). Thus the linear MTF is a good approximation to evaluate the DV averaged over large footprints typical of satellite observations.

Many reported sea level records are often not a single measurement. They are then a composition of different records from several tide gauges. Sometimes, they are from the same tide gauge, but the tide gauge stability has been affected by earthquakes. This is the case of Guam as discussed in the present manuscript. The claimed sea level acceleration of Guam is only the result of two earthquakes that have compromised the stability of the tide gauge.

The use of low-cost, Remotely Operated Vehicle (ROV) and underwater photogrammetry techniques for 3D reconstruction of shallow hydrothermal vent sites around Paleochori Bay, Milos Island, Greece. Characterising venting fields through interactive bathymetry models produced from still images taken from camera onboard ROV flown over areas of interest in double raster pattern. First time the shallow venting fields on Milos have been actively surveyed using ROV. Areas were successfully surveyed and the bathymetry was reconstructed using SfM photogrammetry with a ~10 cm scale resolution. A diverse range of benthic habitats were surveyed and the resulting topographic models will act as a baseline, providing further characterisation of the vent systems and any evolving seafloor morphology associated with mineral deposition.

In Shark Bay, a large hypersaline bay in Western Australia, longitudinal density gradients force gravitational circulation that is important for Bay-ocean exchange. Observations of vertical stratification and velocity are presented, confirming a steady near-bed dense water outflow from Shark Bay’s northern Geographe Channel that persisted through all stages of the tide. Outflow velocities were 2–3 times stronger than the outflows recorded previously in Naturaliste Channel (in the west) and were more resistant to breakdown by tidal mixing. Estimates of turbulent kinetic energy production derived from ADCP data using the variance method showed a more complex structure in Geographe Channel, due to shear associated with the stratified conditions, with peak levels of turbulence occurring during reversal of tidal flows. For both channels the main source of turbulence was tidal flow along the seabed, with the bottom current speed cubed, |U_b³|, providing a reasonable proxy for tidal mixing and prediction of dense water outflows from Shark Bay the majority of the time. Orientation and deeper water of Geographe Channel along the main axis of the longitudinal density gradient provided an explanation for the predominant outflow from the Bay’s northern entrance. These density-driven currents could potentially influence recruitment of commercially fished scallops and prawns through dispersal and flushing of larvae.

Ocean surface currents and winds are tightly coupled essential climate variables, and, given their short time scales, observing them at the same time and resolution is of great interest. DopplerScatt is an airborne Ka-band scatterometer that has been developed under NASA's Instrument Incubator Program (IIP) to provide a proof of concept of the feasability of measuring these variables using pencil-beam scanning Doppler scatterometry. In the first half of this paper, we present the Doppler scatterometer measurement and processing principles, paying particular attention to deriving a complete measurement error budget. Although Doppler radars have been used for the estimation of surface currents, pencil-beam Doppler Scatterometry offers challenges and opportunities that require separate treatment. The calibration of the Doppler measurement to remove platform and instrument biases has been a traditional challenge for Doppler systems, and we introduce several new techniques to mitigate these errors when conical scanning is used. The use of Ka-band for airborne Doppler scatterometry measurements is also new, and, in the second half of the paper, we examine the phenomenology of the mapping from radar cross section and radial velocity measurements to winds and surface currents. To this end, we present new Ka-band Geophysical Model Functions (GMF's) for winds and surface currents obtained from multiple airborne campaigns. We find that the wind Ka-band GMF exhibits similar dependence to wind speed as that for Ku-band scatterometers, such as QuikSCAT, albeit with much greater upwind-crosswind modulation. The surface current GMF at Ka-band is significantly different from that at C-band, and, above 4.5 m/s has a weak dependence on wind speed, although still dependent on wind direction. We examine the effects of Bragg-wave modulation by long waves through a Modululation Transfer Function (MTF), and show that the observed surface current dependence on winds is consistent with past Ka-band MTF observations. Finally, we provide a preliminary validation of our geophysical retrievals, which will be expanded in subsequent publications. Our results indicate that Ka-band Doppler scatterometry could be a feasible method for wide-swath simultaneous measurements of winds and currents from space.

The Total Algae Peak Integration Retrieval TAPIR relates the chlorophyll-a absorption coefficient at 440 nm (a440) to the reflectance peak near 683 nm induced by chlorophyll-a properties. The two-step retrieval provides both the hyperspectral quantification of the phytoplankton fluorescence and scattering and the estimation of a440 from reflectance signals. Integrating the peak, the Total Algae Peak (TAP) accounts for the variance in the peak's magnitude, shape, and central peak wavelength. TAPIR is a solely optical approach estimating a440 and supports the application of retrieval-independent individual regional bio-optical models afterwards to retrieve the chlorophyll-a concentration. Simulations reveal the major sensitivity on the considered model chlorophyll-a absorption spectrum and its single scattering albedo. Additional water and atmosphere constituents have a low impact. An uncertainty assessment reveals uncertainties of less than 30% for TAPIR a440 greater than 0.8 m-1 and less than 38% for lower a440. In optically complex waters, first validation efforts promise the applicability of TAPIR for high chlorophyll-a concentration estimations in the presence of additional water constituents. The technique is generic and considers external conditions (sun zenith angle, number of measurement bands, surface or satellite measurements, and radiometric quantity). TAPIR applies to all kind of waters including optically complex waters, arctic to tropical regions, and inland, coastal, and open ocean waters. Among other hyperspectral satellite sensors, the Environmental Mapping and Analysis Program (EnMAP) provides sufficient sampling bands for the application of TAPIR.

Wave steepness is presented as an extension and a valuable add-on to the conventional set of sea state parameters retrieved from satellite altimetry data. Following physical model based on recent advances of weak turbulence theory wave steepness is estimated from directly measured spatial gradient of wave height. In this way the method works with altimetry trajectories rather than with point-wise data. Moreover, in contrast to widely used parametric models this approach provides us with instantaneous values of wave steepness and period. Relevance of single-track estimates of wave steepness (period) is shown for wave climate studies and confirmed by a simple probabilistic model. The approach is verified via comparison against buoy and satellite data including crossover points for standard 1 second data of Ku-band altimeters. High quality of the physical model and robustness of the parametric ones are examined in terms of global wave statistics. Prospects and relevance of both approaches in the ocean wave climate studies are discussed.

Three cases of Day-Night Band (DNB) observations of the Visible Infrared Imaging Radiometer Suite (VIIRS) onboard the Suomi National Polar-orbiting Partnership (SNPP) are explored for applications to assess the ocean environment and monitor ocean dynamics. An approach to use the DNB radiance ratio was developed in order to better continuously assess the ocean diurnal and short-term environmental changes with VIIRS DNB observations. In the La Plata River Estuary, the sediment fronts showed 20–25 km diurnal inshore-offshore movements on March 13, 2017. In the waters off the Argentina coast in the South Atlantic, VIIRS DNB measurements provided continuous observations and monitoring of the algae bloom development and migration between 24–26 March 2016. This algae bloom generally kept the same spatial patterns, but moved nearly 20 km eastward in the three-day period. In the Yangtze River Estuary and Hangzhou Bay region in China’s east coast, VIIRS DNB observations also revealed the complicated coastal dynamic changes between 12–14 April 2017. Even though there are still some challenges and limitations for monitoring the ocean environment with VIIRS DNB observations, this study shows that satellite DNB observations can provide additional data sources for ocean observations, especially observations during the nighttime.

The primary measure of the quality of sea surface temperature (SST) fields obtained from satellite-borne infrared sensors has been the bias and variance of matchups with co-located in-situ values. Because such matchups tend to be widely separated, these bias and variance estimates are not necessarily a good measure of small scale (several pixels) gradients in these fields because one of the primary contributors to the uncertainty in satellite retrievals is atmospheric contamination, which tends to have large spatial scales compared with the pixel separation of infrared sensors. Hence, there is not a good measure to use in selecting SST fields appropriate for the study of submesoscale processes and, in particular, of processes associated with near-surface fronts, both of which have recently seen a rapid increase in interest. In this study, two methods are examined to address this problem, one based on spectra of the SST data and the other on their variograms. To evaluate the methods, instrument noise was estimated in Level-2 VIIRS and AVHRR SST fields of the Sargasso Sea. The two methods provided very nearly identical results for AVHRR: along-scan values of approximately 0.18 K for both day and night and along-track values of 0.21 K also for day and night. By contrast, the instrument noise estimated for VIIRS varied by method, scan geometry and day-night. Specifically, daytime, along-scan (along-track), spectral estimates were found to be approximately 0.05 K (0.08 K) and the corresponding nighttime values of 0.02 K (0.03 K). Daytime estimates based on the variogram were found to be 0.08 K (0.10 K) with the corresponding nighttime values of 0.04 K (0.06 K). Taken together: AVHRR instrument noise is significantly larger than VIIRS instrument noise, along-track noise is larger than along-scan noise and daytime levels are higher than nighttime levels. Given the similarity of results and the less stringent preprocessing requirements, the variogram is the preferred method although there is a suggestion that this approach overestimates the noise for high quality data in dynamically quiet regions. Finally, simulations of the impact of noise on the determination of SST gradients show that on average the gradient magnitude for typical ocean gradients will be accurately estimated with VIIRS but substantially overestimated with AVHRR.

A radial basis function network (RBFN) method is proposed to reconstruct daily Sea surface temperatures (SSTs) with limited SST samples. For the purpose of evaluating the SSTs using this method, non-biased SST samples in the Pacific Ocean (10°N–30°N, 115°E–135°E) are selected when the tropical storm Hagibis arrived in June 2014, and these SST samples are obtained from the OISST products according to the distribution of AVHRR L2p SST and in-situ SST data. Furthermore, an improved nearest neighbor cluster (INNC) algorithm is designed to search the optimal hidden knots for RBFNs from both the SST samples and the background fields. Then the reconstructed SSTs from the RBFN method are compared with the results from the optimum interpolation (OI) method. The statistical results show that the RBFN method has a better performance of reconstructing SST than the OI method in the study, and the average RMSE is 0.48°C for the RBFN method, which is quite smaller than the value of 0.69°C for the OI method. Additionally, the RBFN methods with different basis functions and clustering algorithms are tested, and we discover that the INNC algorithm with multi-quadric function is quite suitable for the RBFN method to reconstruct SSTs when the SST samples are sparsely distributed.

Ziyuan-3 (ZY-3) satellite is the first civilian stereo mapping satellite in China and was designed to achieve the 1: 50000 scale mapping for land and ocean. Rigorous sensor model (RSM) is required to build the relationship between the three-dimensional (3D) object space and two-dimensional (2D) image space of ZY-3 satellite imagery. However, each satellite sensor has its own imaging system with different physical sensor models, which increase the difficulty of geometric integration of multi-source images with different sensor models. Therefore, it is critical to generate generic model especially rational polynomial coefficients (RPCs) of optical imagery. Recently, relatively a few researches have been conducted on RPCs generation to ZY-3 satellite. This paper proposes an approach to evaluate the performance of RPCs generation from RSM of ZY-3 imagery. Three scenarios experiments with different terrain features (such as ocean, city and grassland) are designed and conducted to comprehensively evaluate the replacement accuracies of this approach and analyze the RPCs fitting error. All the experimental results demonstrate that the proposed method achieved encouraging accuracy of better than 1.946E-04 pixel in both x-axis direction and y-axis direction, it indicates that the RPCs is suitable for ZY-3 imagery and can be used as a replacement for the RSM of ZY-3 imagery.

This paper introduces a new destriping algorithm for remote sensing data. The method is based on combined Haar Stationary Wavelet transform and Fourier filtering. State-of-the-Art methods based on the discrete wavelet transform (DWT) may not always be effective and may cause different artifacts. Our contribution is three-fold: i) we propose to use the Undecimated Wavelet transform (UWT) to avoid as much as possible shortcomings of the classical DWT; ii) we combine a spectral filtering and UWT using the simplest possible wavelet, the Haar basis, for a computational efficiency; iii) we handle 2D fields with missing data, as commonly observed in ocean remote sensing data due to atmospheric conditions (e.g., cloud contamination). The performances of the proposed filter are tested and validated on the suppression of horizontal strip artifacts in cloudy L2 Sea Surface Temperature (SST) and ocean color snapshots.

The relations between the kinetic energy spectrum and the second order longitudinal structure function in two dimensions are derived, and several examples are considered. The forward conversion (from spectrum to structure function) is illustrated first with idealized power law spectra, representing turbulent inertial ranges. The forward conversion is also applied to the zonal kinetic energy spectrum of Nastrom and Gage (1985) and the result agrees well with the longitudinal structure function of Lindborg (1999). The inverse conversion (from structure function to spectrum) is tested with data from 2D turbulence simulations. When applied to the theoretical structure function (derived from the forward conversion of the spectrum), the result closely resembles the original spectrum, except at the largest wavenumbers. However the inverse conversion is much less successful when applied to the structure function obtained from pairs of particles in the flow. This is because the inverse conversion favors large pair separations, which are typically noisy with particle data. Fitting the structure function to a polynomial improves the result, but not sufficiently to distinguish the correct inertial range dependencies. Furthermore the inversion of non-local spectra is largely unsuccessful. Thus it appears that focusing on structure functions with Lagrangian data is preferable to estimating spectra.