Browse
Scripps Argo Trajectory-Based Velocity Product: Global Estimates of Absolute Velocity Derived from Core, Biogeochemical, and Deep Argo Float Trajectories at Parking DepthAbstract
Global estimates of absolute velocities can be derived from Argo float trajectories during drift at parking depth. A new velocity dataset developed and maintained at Scripps Institution of Oceanography is presented based on all Core, Biogeochemical, and Deep Argo float trajectories collected between 2001 and 2020. Discrepancies between velocity estimates from the Scripps dataset and other existing products including YoMaHa and ANDRO are associated with quality control criteria, as well as selected parking depth and cycle time. In the Scripps product, over 1.3 million velocity estimates are used to reconstruct a time-mean velocity field for the 800–1200 dbar layer at 1° horizontal resolution. This dataset provides a benchmark to evaluate the veracity of the BRAN2020 reanalysis in representing the observed variability of absolute velocities and offers a compelling opportunity for improved characterization and representation in forecast and reanalysis systems.
Significance Statement
The aim of this study is to provide observation-based estimates of the large-scale, subsurface ocean circulation. We exploit the drift of autonomous profiling floats to carefully isolate the inferred circulation at the parking depth, and combine observations from over 11 000 floats, sampling between 2001 and 2020, to deliver a new dataset with unprecedented accuracy. The new estimates of subsurface currents are suitable for assessing global models, reanalyses, and forecasts, and for constraining ocean circulation in data-assimilating models.
Abstract
Global estimates of absolute velocities can be derived from Argo float trajectories during drift at parking depth. A new velocity dataset developed and maintained at Scripps Institution of Oceanography is presented based on all Core, Biogeochemical, and Deep Argo float trajectories collected between 2001 and 2020. Discrepancies between velocity estimates from the Scripps dataset and other existing products including YoMaHa and ANDRO are associated with quality control criteria, as well as selected parking depth and cycle time. In the Scripps product, over 1.3 million velocity estimates are used to reconstruct a time-mean velocity field for the 800–1200 dbar layer at 1° horizontal resolution. This dataset provides a benchmark to evaluate the veracity of the BRAN2020 reanalysis in representing the observed variability of absolute velocities and offers a compelling opportunity for improved characterization and representation in forecast and reanalysis systems.
Significance Statement
The aim of this study is to provide observation-based estimates of the large-scale, subsurface ocean circulation. We exploit the drift of autonomous profiling floats to carefully isolate the inferred circulation at the parking depth, and combine observations from over 11 000 floats, sampling between 2001 and 2020, to deliver a new dataset with unprecedented accuracy. The new estimates of subsurface currents are suitable for assessing global models, reanalyses, and forecasts, and for constraining ocean circulation in data-assimilating models.
Abstract
The existence of significant cross-polar antenna patterns, as well as the scan-dependent measurement biases, inherent to the polarimetric phased array radar (PPAR), are among the most important risk factors for using this technology in weather observations. The cross-polar patterns on receive induce cross coupling between returns from the two orthogonal fields causing biases in polarimetric variable estimates. Furthermore, the electromagnetic coupling in hardware may exacerbate the cross-coupling effects. To address this problem, a pulse-to-pulse phase coding in either the horizontal or vertical ports of the transmission elements has been proposed. However, it does not affect the scan-dependent system biases in PPAR estimates which require corrections via calibration mechanisms. Further, the cross-coupling signals are proportional to the cross-polar pattern power levels, rendering mitigation effective only at steering angles where these levels are sufficiently low (e.g., approximately less than ~-25 dB). In that regard, any approach that augments the number of such steering angles benefits the cross-coupling mitigation effectiveness. Herein, a simple approach that has a potential to achieve this via antenna tilt is presented.
Abstract
The existence of significant cross-polar antenna patterns, as well as the scan-dependent measurement biases, inherent to the polarimetric phased array radar (PPAR), are among the most important risk factors for using this technology in weather observations. The cross-polar patterns on receive induce cross coupling between returns from the two orthogonal fields causing biases in polarimetric variable estimates. Furthermore, the electromagnetic coupling in hardware may exacerbate the cross-coupling effects. To address this problem, a pulse-to-pulse phase coding in either the horizontal or vertical ports of the transmission elements has been proposed. However, it does not affect the scan-dependent system biases in PPAR estimates which require corrections via calibration mechanisms. Further, the cross-coupling signals are proportional to the cross-polar pattern power levels, rendering mitigation effective only at steering angles where these levels are sufficiently low (e.g., approximately less than ~-25 dB). In that regard, any approach that augments the number of such steering angles benefits the cross-coupling mitigation effectiveness. Herein, a simple approach that has a potential to achieve this via antenna tilt is presented.
Abstract
The Polar Radiant Energy in the Far InfraRed Experiment (PREFIRE) mission will measure the Earth’s emission at wavelengths ranging from 3-54 µm. The pre-launch clear-sky retrieval algorithm, evaluated with simulated test data, indicates that PREFIRE measurements will be valuable for retrieving atmospheric water vapor and temperature profiles. Far infrared measurements provide unique retrieval information, indicated by the high ranking of select FIR channels as primary contributors to the total degrees of freedom for signal (DFS). In utilizing all the PREFIRE channels, the average total DFS of 4 test regions ranges from 1.90 - 4.71. The information content increases with higher column water vapor and in the presence of near surface temperature inversions. Using the DFS profiles for guidance, the retrieval concentrates information into 7 distinct layers to reduce the retrieval uncertainty per layer. Sensitivity tests indicate forward model error due to surface emissivity uncertainty results in about a 9% increase in column water vapor uncertainty. The clear-sky retrieval is sensitive to the presence of undetected ice clouds, especially those with optical depths larger than 0.2. Hence, in addition to a separate PREFIRE cloud mask, optimal estimation retrieval metrics are explored as possible indicators of cloudy scenes.
Abstract
The Polar Radiant Energy in the Far InfraRed Experiment (PREFIRE) mission will measure the Earth’s emission at wavelengths ranging from 3-54 µm. The pre-launch clear-sky retrieval algorithm, evaluated with simulated test data, indicates that PREFIRE measurements will be valuable for retrieving atmospheric water vapor and temperature profiles. Far infrared measurements provide unique retrieval information, indicated by the high ranking of select FIR channels as primary contributors to the total degrees of freedom for signal (DFS). In utilizing all the PREFIRE channels, the average total DFS of 4 test regions ranges from 1.90 - 4.71. The information content increases with higher column water vapor and in the presence of near surface temperature inversions. Using the DFS profiles for guidance, the retrieval concentrates information into 7 distinct layers to reduce the retrieval uncertainty per layer. Sensitivity tests indicate forward model error due to surface emissivity uncertainty results in about a 9% increase in column water vapor uncertainty. The clear-sky retrieval is sensitive to the presence of undetected ice clouds, especially those with optical depths larger than 0.2. Hence, in addition to a separate PREFIRE cloud mask, optimal estimation retrieval metrics are explored as possible indicators of cloudy scenes.
Application of Machine Learning Techniques to Ocean Mooring Time Series DataAbstract
In situ observations are vital to improving our understanding of the variability and dynamics of the ocean. A critical component of the ocean circulation is the strong, narrow, and highly variable western boundary currents. Ocean moorings that extend from the seafloor to the surface remain the most effective and efficient method to fully observe these currents. For various reasons, mooring instruments may not provide continuous records. Here we assess the application of the Iterative Completion Self-Organizing Maps (ITCOMPSOM) machine learning technique to fill observational data gaps in a 7.5 yr time series of the East Australian Current. The method was validated by withholding parts of fully known profiles, and reconstructing them. For 20% random withholding of known velocity data, validation statistics of the u- and υ-velocity components are R 2 coefficients of 0.70 and 0.88 and root-mean-square errors of 0.038 and 0.064 m s−1, respectively. Withholding 100 days of known velocity profiles over a depth range between 60 and 700 m has mean profile residual differences between true and predicted u and υ velocity of 0.009 and 0.02 m s−1, respectively. The ITCOMPSOM also reproduces the known velocity variability. For 20% withholding of salinity and temperature data, root-mean-square errors of 0.04 and 0.38°C, respectively, are obtained. The ITCOMPSOM validation statistics are significantly better than those obtained when standard data filling methods are used. We suggest that machine learning techniques can be an appropriate method to fill missing data and enable production of observational-derived data products.
Significance Statement
Moored observational time series of ocean boundary currents monitor the full-depth variability and change of these dynamic currents and are used to understand their influence on large-scale ocean climate, regional shelf–coastal processes, extreme weather, and seasonal climate. In this study we apply a machine learning technique, Iterative Completion Self-Organizing Maps (ITCOMPSOM), to fill data gaps in a boundary current moored observational data record. The ITCOMPSOM provides an improved method to fill data gaps in the mooring record and if applied to other observational data records may improve the reconstruction of missing data. The derived gridded data product should improve the accessibility and potentially increase the use of these data.
Abstract
In situ observations are vital to improving our understanding of the variability and dynamics of the ocean. A critical component of the ocean circulation is the strong, narrow, and highly variable western boundary currents. Ocean moorings that extend from the seafloor to the surface remain the most effective and efficient method to fully observe these currents. For various reasons, mooring instruments may not provide continuous records. Here we assess the application of the Iterative Completion Self-Organizing Maps (ITCOMPSOM) machine learning technique to fill observational data gaps in a 7.5 yr time series of the East Australian Current. The method was validated by withholding parts of fully known profiles, and reconstructing them. For 20% random withholding of known velocity data, validation statistics of the u- and υ-velocity components are R 2 coefficients of 0.70 and 0.88 and root-mean-square errors of 0.038 and 0.064 m s−1, respectively. Withholding 100 days of known velocity profiles over a depth range between 60 and 700 m has mean profile residual differences between true and predicted u and υ velocity of 0.009 and 0.02 m s−1, respectively. The ITCOMPSOM also reproduces the known velocity variability. For 20% withholding of salinity and temperature data, root-mean-square errors of 0.04 and 0.38°C, respectively, are obtained. The ITCOMPSOM validation statistics are significantly better than those obtained when standard data filling methods are used. We suggest that machine learning techniques can be an appropriate method to fill missing data and enable production of observational-derived data products.
Significance Statement
Moored observational time series of ocean boundary currents monitor the full-depth variability and change of these dynamic currents and are used to understand their influence on large-scale ocean climate, regional shelf–coastal processes, extreme weather, and seasonal climate. In this study we apply a machine learning technique, Iterative Completion Self-Organizing Maps (ITCOMPSOM), to fill data gaps in a boundary current moored observational data record. The ITCOMPSOM provides an improved method to fill data gaps in the mooring record and if applied to other observational data records may improve the reconstruction of missing data. The derived gridded data product should improve the accessibility and potentially increase the use of these data.
Abstract
Accurate ocean wave measurements are needed for the safe design and operation of offshore facilities, but despite many ocean wave measurements, the accuracy of wave measurement systems remains an ongoing issue. Of paramount importance are measurements during extreme sea states. This paper examines wave measurements made with an Optech Laser, a Rosemount WaveRadar, and a Datawell Waverider buoy at North Rankin A platform (NRA), Australia; Ekofisk, North Sea; and several South China Sea locations. We evaluate the relative performance of these instruments based upon various frequency domain comparisons, including comparisons of their 1-D frequency spectra using spectrograms, spectral moments, high-frequency tail slopes, and significant wave heights derived from their wave spectra. A spectral relationship (transfer function) in terms of mean spectral ratio of the instruments is developed, which can be used for spectral calibration. On average, Laser and Waverider spectral estimates agree well at all sea states. However, at low wind speeds, the higher frequency spectral levels of the Laser are relatively high and noisy compared with the other two instruments. Radar higher frequency spectral estimates are relatively low compared to the other two instruments, particularly at lower sea states. In additionally, the higher frequency tail slopes of all three instruments vary between f ‒4 and f ‒5. However, at higher sea states, the Waverider tail slopes become steeper than f ‒5. The Radar produces the lowest significant wave heights (H m0) compared to the Laser and Waverider, but its second moment period (T m02) estimates are longer than the Laser and Waverider.
Abstract
Accurate ocean wave measurements are needed for the safe design and operation of offshore facilities, but despite many ocean wave measurements, the accuracy of wave measurement systems remains an ongoing issue. Of paramount importance are measurements during extreme sea states. This paper examines wave measurements made with an Optech Laser, a Rosemount WaveRadar, and a Datawell Waverider buoy at North Rankin A platform (NRA), Australia; Ekofisk, North Sea; and several South China Sea locations. We evaluate the relative performance of these instruments based upon various frequency domain comparisons, including comparisons of their 1-D frequency spectra using spectrograms, spectral moments, high-frequency tail slopes, and significant wave heights derived from their wave spectra. A spectral relationship (transfer function) in terms of mean spectral ratio of the instruments is developed, which can be used for spectral calibration. On average, Laser and Waverider spectral estimates agree well at all sea states. However, at low wind speeds, the higher frequency spectral levels of the Laser are relatively high and noisy compared with the other two instruments. Radar higher frequency spectral estimates are relatively low compared to the other two instruments, particularly at lower sea states. In additionally, the higher frequency tail slopes of all three instruments vary between f ‒4 and f ‒5. However, at higher sea states, the Waverider tail slopes become steeper than f ‒5. The Radar produces the lowest significant wave heights (H m0) compared to the Laser and Waverider, but its second moment period (T m02) estimates are longer than the Laser and Waverider.
Abstract
The present work details the measurement capabilities of Wave Glider autonomous surface vehicles (ASVs) for research-grade meteorology, wave, and current data. Methodologies for motion compensation are described and tested, including a correction technique to account for Doppler shifting of the wave signal. Wave Glider measurements are evaluated against observations obtained from World Meteorological Organization (WMO)-compliant moored buoy assets located off the coast of Southern California. The validation spans a range of field conditions and includes multiple deployments to assess the quality of vehicle-based observations. Results indicate that Wave Gliders can accurately measure wave spectral information, bulk wave parameters, water velocities, bulk winds, and other atmospheric variables with the application of appropriate motion compensation techniques. Measurement errors were found to be comparable to those from reference moored buoys and within WMO operational requirements. The findings of this study represent a step toward enabling the use of ASV-based data for the calibration and validation of remote observations and assimilation into forecast models.
Abstract
The present work details the measurement capabilities of Wave Glider autonomous surface vehicles (ASVs) for research-grade meteorology, wave, and current data. Methodologies for motion compensation are described and tested, including a correction technique to account for Doppler shifting of the wave signal. Wave Glider measurements are evaluated against observations obtained from World Meteorological Organization (WMO)-compliant moored buoy assets located off the coast of Southern California. The validation spans a range of field conditions and includes multiple deployments to assess the quality of vehicle-based observations. Results indicate that Wave Gliders can accurately measure wave spectral information, bulk wave parameters, water velocities, bulk winds, and other atmospheric variables with the application of appropriate motion compensation techniques. Measurement errors were found to be comparable to those from reference moored buoys and within WMO operational requirements. The findings of this study represent a step toward enabling the use of ASV-based data for the calibration and validation of remote observations and assimilation into forecast models.
Abstract
The assimilation of hyperspectral infrared sounders (HIS) observations aboard Earth-observing satellites has become vital to numerical weather prediction, yet this assimilation is predicated on the assumption of clear-sky observations. Using collocated assimilated observations from the Atmospheric Infrared Sounder (AIRS) and the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP), it is found that nearly 7.7% of HIS observations assimilated by the Naval Research Laboratory Variational Data Assimilation System–Accelerated Representer (NAVDAS-AR) are contaminated by cirrus clouds. These contaminating clouds primarily exhibit visible cloud optical depths at 532 nm (COD532nm) below 0.10 and cloud-top temperatures between 240 and 185 K as expected for cirrus clouds. These contamination statistics are consistent with simulations from the Radiative Transfer for TOVS (RTTOV) model showing a cirrus cloud with a COD532nm of 0.10 imparts brightness temperature differences below typical innovation thresholds used by NAVDAS-AR. Using a one-dimensional variational (1DVar) assimilation system coupled with RTTOV for forward and gradient radiative transfer, the analysis temperature and moisture impact of assimilating cirrus-contaminated HIS observations is estimated. Large differences of 2.5 K in temperature and 11 K in dewpoint are possible for a cloud with COD532nm of 0.10 and cloud-top temperature of 210 K. When normalized by the contamination statistics, global differences of nearly 0.11 K in temperature and 0.34 K in dewpoint are possible, with temperature and dewpoint tropospheric root-mean-squared errors (RMSDs) as large as 0.06 and 0.11 K, respectively. While in isolation these global estimates are not particularly concerning, differences are likely much larger in regions with high cirrus frequency.
Abstract
The assimilation of hyperspectral infrared sounders (HIS) observations aboard Earth-observing satellites has become vital to numerical weather prediction, yet this assimilation is predicated on the assumption of clear-sky observations. Using collocated assimilated observations from the Atmospheric Infrared Sounder (AIRS) and the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP), it is found that nearly 7.7% of HIS observations assimilated by the Naval Research Laboratory Variational Data Assimilation System–Accelerated Representer (NAVDAS-AR) are contaminated by cirrus clouds. These contaminating clouds primarily exhibit visible cloud optical depths at 532 nm (COD532nm) below 0.10 and cloud-top temperatures between 240 and 185 K as expected for cirrus clouds. These contamination statistics are consistent with simulations from the Radiative Transfer for TOVS (RTTOV) model showing a cirrus cloud with a COD532nm of 0.10 imparts brightness temperature differences below typical innovation thresholds used by NAVDAS-AR. Using a one-dimensional variational (1DVar) assimilation system coupled with RTTOV for forward and gradient radiative transfer, the analysis temperature and moisture impact of assimilating cirrus-contaminated HIS observations is estimated. Large differences of 2.5 K in temperature and 11 K in dewpoint are possible for a cloud with COD532nm of 0.10 and cloud-top temperature of 210 K. When normalized by the contamination statistics, global differences of nearly 0.11 K in temperature and 0.34 K in dewpoint are possible, with temperature and dewpoint tropospheric root-mean-squared errors (RMSDs) as large as 0.06 and 0.11 K, respectively. While in isolation these global estimates are not particularly concerning, differences are likely much larger in regions with high cirrus frequency.
Abstract
This manuscript presents several improvements to methods for despiking and measuring turbulent dissipation values with acoustic Doppler velocimeters (ADVs). This includes an improved inertial subrange fitting algorithm relevant for all experimental conditions as well as other modifications designed to address failures of existing methods in the presence of large infragravity (IG) frequency bores and other intermittent, nonlinear processes. We provide a modified despiking algorithm, wavenumber spectrum calculation algorithm, and inertial subrange fitting algorithm that together produce reliable dissipation measurements in the presence of IG frequency bores, representing turbulence over a 30 min interval. We use a semi-idealized model to show that our spectrum calculation approach works substantially better than existing wave correction equations that rely on Gaussian-based velocity distributions. We also find that our inertial subrange fitting algorithm provides more robust results than existing approaches that rely on identifying a single best fit and that this improvement is independent of environmental conditions. Finally, we perform a detailed error analysis to assist in future use of these algorithms and identify areas that need careful consideration. This error analysis uses error distribution widths to find, with 95% confidence, an average systematic uncertainty of ±15.2% and statistical uncertainty of ±7.8% for our final dissipation measurements. In addition, we find that small changes to ADV despiking approaches can lead to large uncertainties in turbulent dissipation and that further work is needed to ensure more reliable despiking algorithms.
Significance Statement
Turbulent mixing is a process where the random movement of water can lead to water with different properties irreversibly mixing. This process is important to understand in estuaries because the extent of mixing of freshwater and saltwater inside an estuary alters its overall circulation and thus affects ecosystem health and the distribution of pollution or larvae in an estuary, among other things. Existing approaches to measuring turbulent dissipation, an important parameter for evaluating turbulent mixing, make assumptions that fail in the presence of certain processes, such as long-period, breaking waves in shallow estuaries. We evaluate and improve data analysis techniques to account for such processes and accurately measure turbulent dissipation in shallow estuaries. Some of our improvements are also relevant to a broad array of coastal and oceanic conditions.
Abstract
This manuscript presents several improvements to methods for despiking and measuring turbulent dissipation values with acoustic Doppler velocimeters (ADVs). This includes an improved inertial subrange fitting algorithm relevant for all experimental conditions as well as other modifications designed to address failures of existing methods in the presence of large infragravity (IG) frequency bores and other intermittent, nonlinear processes. We provide a modified despiking algorithm, wavenumber spectrum calculation algorithm, and inertial subrange fitting algorithm that together produce reliable dissipation measurements in the presence of IG frequency bores, representing turbulence over a 30 min interval. We use a semi-idealized model to show that our spectrum calculation approach works substantially better than existing wave correction equations that rely on Gaussian-based velocity distributions. We also find that our inertial subrange fitting algorithm provides more robust results than existing approaches that rely on identifying a single best fit and that this improvement is independent of environmental conditions. Finally, we perform a detailed error analysis to assist in future use of these algorithms and identify areas that need careful consideration. This error analysis uses error distribution widths to find, with 95% confidence, an average systematic uncertainty of ±15.2% and statistical uncertainty of ±7.8% for our final dissipation measurements. In addition, we find that small changes to ADV despiking approaches can lead to large uncertainties in turbulent dissipation and that further work is needed to ensure more reliable despiking algorithms.
Significance Statement
Turbulent mixing is a process where the random movement of water can lead to water with different properties irreversibly mixing. This process is important to understand in estuaries because the extent of mixing of freshwater and saltwater inside an estuary alters its overall circulation and thus affects ecosystem health and the distribution of pollution or larvae in an estuary, among other things. Existing approaches to measuring turbulent dissipation, an important parameter for evaluating turbulent mixing, make assumptions that fail in the presence of certain processes, such as long-period, breaking waves in shallow estuaries. We evaluate and improve data analysis techniques to account for such processes and accurately measure turbulent dissipation in shallow estuaries. Some of our improvements are also relevant to a broad array of coastal and oceanic conditions.
