Browse
The Mount Washington Observatory Regional Mesonet: A Technical Overview of a Mountain-Based MesonetAbstract
The Mount Washington Observatory Regional Mesonet (MWRM) is a network of 18 remote meteorological monitoring stations (as of 2022), including the Auto Road Vertical Profile (ARVP), located across the White Mountains of northern New Hampshire. Each station measures temperature and relative humidity, with additional variables at many locations. All stations need to withstand the frequent combination of intense cold, high precipitation amounts, icing, and hurricane-force winds in a mountain environment. Due to these challenges, the MWRM employs rugged instrumentation, an innovative radio-communications relay approach, and carefully selected sites that balance ideal measuring environments with station survivability. Data collected from the MWRM are used operationally by forecasters (including Mount Washington Observatory and National Weather Service staff) to validate model guidance, by alpine and climate scientists, by recreationalists accessing conditions in the backcountry, by groups operating on the mountain (Cog Railway, toll Auto Road), and by search and rescue organizations. This paper provides a detailed description of the network, with emphasis on how the challenging climate and terrain of this mountain region impacts sensor selection, site maintenance, and overall operation.
Significance Statement
The mountain environment is a heterogeneous landscape, and interactions between the atmosphere and terrain can cause a wide variety of conditions across time and space. Our network of remote stations at different elevations across the White Mountains allows data users to understand how the weather varies spatially across the mountain range where conditions on higher peaks can be drastically, and dangerously, different. Sharing information about the MWRM can help other groups establish networks in similar challenging environments, and broaden our understanding of weather and climate in mountainous regions.
Abstract
The Mount Washington Observatory Regional Mesonet (MWRM) is a network of 18 remote meteorological monitoring stations (as of 2022), including the Auto Road Vertical Profile (ARVP), located across the White Mountains of northern New Hampshire. Each station measures temperature and relative humidity, with additional variables at many locations. All stations need to withstand the frequent combination of intense cold, high precipitation amounts, icing, and hurricane-force winds in a mountain environment. Due to these challenges, the MWRM employs rugged instrumentation, an innovative radio-communications relay approach, and carefully selected sites that balance ideal measuring environments with station survivability. Data collected from the MWRM are used operationally by forecasters (including Mount Washington Observatory and National Weather Service staff) to validate model guidance, by alpine and climate scientists, by recreationalists accessing conditions in the backcountry, by groups operating on the mountain (Cog Railway, toll Auto Road), and by search and rescue organizations. This paper provides a detailed description of the network, with emphasis on how the challenging climate and terrain of this mountain region impacts sensor selection, site maintenance, and overall operation.
Significance Statement
The mountain environment is a heterogeneous landscape, and interactions between the atmosphere and terrain can cause a wide variety of conditions across time and space. Our network of remote stations at different elevations across the White Mountains allows data users to understand how the weather varies spatially across the mountain range where conditions on higher peaks can be drastically, and dangerously, different. Sharing information about the MWRM can help other groups establish networks in similar challenging environments, and broaden our understanding of weather and climate in mountainous regions.
Understanding Differences in Sea Surface Temperature IntercomparisonsAbstract
Our study shows that the intercomparison among sea surface temperature (SST) products is influenced by the choice of SST reference, and the interpolation of SST products. The influence of reference SST depends on whether the reference SSTs are averaged to a grid or in pointwise in situ locations, including buoy or Argo observations, and filtered by first-guess or climatology quality control (QC) algorithms. The influence of the interpolation depends on whether SST products are in their original grids or preprocessed into common coarse grids. The impacts of these factors are demonstrated in our assessments of eight widely used SST products (DOISST, MUR25, MGDSST, GAMSSA, OSTIA, GPB, CCI, CMC) relative to buoy observations: (i) when the reference SSTs are averaged onto 0.25° × 0.25° grid boxes, the magnitude of biases is lower in DOISST and MGDSST (<0.03°C), and magnitude of root-mean-square differences (RMSDs) is lower in DOISST (0.38°C) and OSTIA (0.43°C); (ii) when the same reference SSTs are evaluated at pointwise in situ locations, the standard deviations (SDs) are smaller in DOISST (0.38°C) and OSTIA (0.39°C) on 0.25° × 0.25° grids; but the SDs become smaller in OSTIA (0.34°C) and CMC (0.37°C) on products’ original grids, showing the advantage of those high-resolution analyses for resolving finer-scale SSTs; (iii) when a loose QC algorithm is applied to the reference buoy observations, SDs increase; and vice versa; however, the relative performance of products remains the same; and (iv) when the drifting-buoy or Argo observations are used as the reference, the magnitude of RMSDs and SDs become smaller, potentially due to changes in observing intervals. These results suggest that high-resolution SST analyses may take advantage in intercomparisons.
Significance Statement
Intercomparisons of gridded SST products be affected by how the products are compared with in situ observations: whether the products are in coarse (0.25°) or original (0.05°–0.10°) grids, whether the in situ SSTs are in their reported locations or gridded and how they are quality controlled, and whether the biases of satellite SSTs are corrected by localized matchups or large-scale patterns. By taking all these factors into account, our analyses indicate that the NOAA DOISST is among the best SST products for the long period (1981–present) and relatively coarse (0.25°) resolution that it was designed for.
Abstract
Our study shows that the intercomparison among sea surface temperature (SST) products is influenced by the choice of SST reference, and the interpolation of SST products. The influence of reference SST depends on whether the reference SSTs are averaged to a grid or in pointwise in situ locations, including buoy or Argo observations, and filtered by first-guess or climatology quality control (QC) algorithms. The influence of the interpolation depends on whether SST products are in their original grids or preprocessed into common coarse grids. The impacts of these factors are demonstrated in our assessments of eight widely used SST products (DOISST, MUR25, MGDSST, GAMSSA, OSTIA, GPB, CCI, CMC) relative to buoy observations: (i) when the reference SSTs are averaged onto 0.25° × 0.25° grid boxes, the magnitude of biases is lower in DOISST and MGDSST (<0.03°C), and magnitude of root-mean-square differences (RMSDs) is lower in DOISST (0.38°C) and OSTIA (0.43°C); (ii) when the same reference SSTs are evaluated at pointwise in situ locations, the standard deviations (SDs) are smaller in DOISST (0.38°C) and OSTIA (0.39°C) on 0.25° × 0.25° grids; but the SDs become smaller in OSTIA (0.34°C) and CMC (0.37°C) on products’ original grids, showing the advantage of those high-resolution analyses for resolving finer-scale SSTs; (iii) when a loose QC algorithm is applied to the reference buoy observations, SDs increase; and vice versa; however, the relative performance of products remains the same; and (iv) when the drifting-buoy or Argo observations are used as the reference, the magnitude of RMSDs and SDs become smaller, potentially due to changes in observing intervals. These results suggest that high-resolution SST analyses may take advantage in intercomparisons.
Significance Statement
Intercomparisons of gridded SST products be affected by how the products are compared with in situ observations: whether the products are in coarse (0.25°) or original (0.05°–0.10°) grids, whether the in situ SSTs are in their reported locations or gridded and how they are quality controlled, and whether the biases of satellite SSTs are corrected by localized matchups or large-scale patterns. By taking all these factors into account, our analyses indicate that the NOAA DOISST is among the best SST products for the long period (1981–present) and relatively coarse (0.25°) resolution that it was designed for.
Abstract
The Ocean Color Component of the Aerosol Robotic Network (AERONET-OC) aims at supporting the assessment of satellite ocean color radiometric products with in situ reference data derived from automated above-water measurements. This study, applying metrology principles and taking advantage of recent technology and science advances, revisits the uncertainty estimates formerly provided for AERONET-OC normalized water-leaving radiances L WN. The new uncertainty values are quantified for a number of AERONET-OC sites located in marine regions representative of chlorophyll-a-dominated waters (i.e., Case 1) and a variety of optically complex waters. Results show uncertainties typically increasing with the optical complexity of water and wind speed. Relative and absolute uncertainty values are provided for the various sites together with contributions from each source of uncertainty affecting measurements. In view of supporting AERONET-OC data users, the study also suggests practical solutions to quantify uncertainties for L WN from its spectral values. Additionally, results from an evaluation of the temporal variability characterizing L WN at various AERONET-OC sites are presented to address the impact of temporal mismatches between in situ and satellite data in matchup analysis.
Abstract
The Ocean Color Component of the Aerosol Robotic Network (AERONET-OC) aims at supporting the assessment of satellite ocean color radiometric products with in situ reference data derived from automated above-water measurements. This study, applying metrology principles and taking advantage of recent technology and science advances, revisits the uncertainty estimates formerly provided for AERONET-OC normalized water-leaving radiances L WN. The new uncertainty values are quantified for a number of AERONET-OC sites located in marine regions representative of chlorophyll-a-dominated waters (i.e., Case 1) and a variety of optically complex waters. Results show uncertainties typically increasing with the optical complexity of water and wind speed. Relative and absolute uncertainty values are provided for the various sites together with contributions from each source of uncertainty affecting measurements. In view of supporting AERONET-OC data users, the study also suggests practical solutions to quantify uncertainties for L WN from its spectral values. Additionally, results from an evaluation of the temporal variability characterizing L WN at various AERONET-OC sites are presented to address the impact of temporal mismatches between in situ and satellite data in matchup analysis.
Abstract
Marine air temperatures recorded on ships during the daytime are known to be biased warm on average due to energy storage by the superstructure of the vessels. This makes unadjusted daytime observations unsuitable for many applications including for the monitoring of long-term temperature change over the oceans. In this paper a physics-based approach is used to estimate this heating bias in ship observations from ICOADS. Under this approach, empirically determined coefficients represent the energy transfer terms of a heat budget model that quantifies the heating bias and is applied as a function of cloud cover and the relative wind speed over individual ships. The coefficients for each ship are derived from the anomalous diurnal heating relative to nighttime air temperature. Model coefficients, cloud cover, and relative wind speed are then used to estimate the heating bias ship by ship and generate nighttime-equivalent time series. A variety of methodological approaches were tested. Application of this method enables the inclusion of some daytime observations in climate records based on marine air temperatures, allowing an earlier start date and giving an increase in spatial coverage compared to existing records that exclude daytime observations.
Significance Statement
Currently, the longest available record of air temperature over the oceans starts in 1880. We present an approach that enables observations of air temperatures over the oceans to be used in the creation of long-term climate records that are presently excluded. We do this by estimating the biases inherent in daytime temperature reports from ships, and adjust for these biases by implementing a numerical heat-budget model. The adjustment can be applied to the variety of ship types present in observational archives. The resulting adjusted temperatures can be used to create a more spatially complete record over the oceans, that extends further back in time, potentially into the late eighteenth century.
Abstract
Marine air temperatures recorded on ships during the daytime are known to be biased warm on average due to energy storage by the superstructure of the vessels. This makes unadjusted daytime observations unsuitable for many applications including for the monitoring of long-term temperature change over the oceans. In this paper a physics-based approach is used to estimate this heating bias in ship observations from ICOADS. Under this approach, empirically determined coefficients represent the energy transfer terms of a heat budget model that quantifies the heating bias and is applied as a function of cloud cover and the relative wind speed over individual ships. The coefficients for each ship are derived from the anomalous diurnal heating relative to nighttime air temperature. Model coefficients, cloud cover, and relative wind speed are then used to estimate the heating bias ship by ship and generate nighttime-equivalent time series. A variety of methodological approaches were tested. Application of this method enables the inclusion of some daytime observations in climate records based on marine air temperatures, allowing an earlier start date and giving an increase in spatial coverage compared to existing records that exclude daytime observations.
Significance Statement
Currently, the longest available record of air temperature over the oceans starts in 1880. We present an approach that enables observations of air temperatures over the oceans to be used in the creation of long-term climate records that are presently excluded. We do this by estimating the biases inherent in daytime temperature reports from ships, and adjust for these biases by implementing a numerical heat-budget model. The adjustment can be applied to the variety of ship types present in observational archives. The resulting adjusted temperatures can be used to create a more spatially complete record over the oceans, that extends further back in time, potentially into the late eighteenth century.
Abstract
The Polar Radiant Energy in the Far Infrared Experiment (PREFIRE) will fill a gap in our understanding of polar processes and the polar climate by offering widespread, spectrally resolved measurements through the far-infrared (FIR) with two identical CubeSat spacecraft. While the polar regions are typically difficult for skillful cloud identification due to cold surface temperatures, the reflection by bright surfaces, and frequent temperature inversions, the inclusion of the FIR may offer increased spectral sensitivity, allowing for the detection of even thin ice clouds. This study assesses the potential skill, as well as limitations, of a neural network (NN)-based cloud mask using simulated spectra mimicking what the PREFIRE mission will capture. Analysis focuses on the polar regions. Clouds are found to be detected approximately 90% of time using the derived neural network. The NN’s assigned confidence for whether a scene is “clear” or “cloudy” proves to be a skillful way in which quality flags can be attached to predictions. Clouds with higher cloud-top heights are typically more easily detected. Low-altitude clouds over polar surfaces, which are the most difficult for the NN to detect, are still detected over 80% of the time. The FIR portion of the spectrum is found to increase the detection of clear scenes and increase mid- to high-altitude cloud detection. Cloud detection skill improves through the use of the overlapping fields of view produced by the PREFIRE instrument’s sampling strategy. Overlapping fields of view increase accuracy relative to the baseline NN while simultaneously predicting on a sub-FOV scale.
Significance Statement
Clouds play an important role in defining the Arctic and Antarctic climates. The purpose of this study is to explore the potential of never-before systematically measured radiative properties of the atmosphere to aid in the detection of polar clouds, which are traditionally difficult to detect. Satellite measurements of emitted radiation at wavelengths longer than 15 μm, combined with complex machine learning methods, may allow us to better understand the occurrence of various cloud types at both poles. The occurrence of these clouds can determine whether the surface warms or cools, influencing surface temperatures and the rate at which ice melts or refreezes. Understanding the frequencies of these various clouds is increasingly important within the context of our rapidly changing climate.
Abstract
The Polar Radiant Energy in the Far Infrared Experiment (PREFIRE) will fill a gap in our understanding of polar processes and the polar climate by offering widespread, spectrally resolved measurements through the far-infrared (FIR) with two identical CubeSat spacecraft. While the polar regions are typically difficult for skillful cloud identification due to cold surface temperatures, the reflection by bright surfaces, and frequent temperature inversions, the inclusion of the FIR may offer increased spectral sensitivity, allowing for the detection of even thin ice clouds. This study assesses the potential skill, as well as limitations, of a neural network (NN)-based cloud mask using simulated spectra mimicking what the PREFIRE mission will capture. Analysis focuses on the polar regions. Clouds are found to be detected approximately 90% of time using the derived neural network. The NN’s assigned confidence for whether a scene is “clear” or “cloudy” proves to be a skillful way in which quality flags can be attached to predictions. Clouds with higher cloud-top heights are typically more easily detected. Low-altitude clouds over polar surfaces, which are the most difficult for the NN to detect, are still detected over 80% of the time. The FIR portion of the spectrum is found to increase the detection of clear scenes and increase mid- to high-altitude cloud detection. Cloud detection skill improves through the use of the overlapping fields of view produced by the PREFIRE instrument’s sampling strategy. Overlapping fields of view increase accuracy relative to the baseline NN while simultaneously predicting on a sub-FOV scale.
Significance Statement
Clouds play an important role in defining the Arctic and Antarctic climates. The purpose of this study is to explore the potential of never-before systematically measured radiative properties of the atmosphere to aid in the detection of polar clouds, which are traditionally difficult to detect. Satellite measurements of emitted radiation at wavelengths longer than 15 μm, combined with complex machine learning methods, may allow us to better understand the occurrence of various cloud types at both poles. The occurrence of these clouds can determine whether the surface warms or cools, influencing surface temperatures and the rate at which ice melts or refreezes. Understanding the frequencies of these various clouds is increasingly important within the context of our rapidly changing climate.
Abstract
Zenith hydrostatic delay (ZHD) is a crucial parameter in Global Navigation Satellite System (GNSS) navigation and positioning and GNSS meteorology. Since the Saastamoinen ZHD model has a larger error in China, it is significant to improve the Saastamoinen ZHD model. This work first estimated the Saastamoinen model using the integrated ZHD as reference values obtained from radiosonde data collected at 73 stations in China from 2012 to 2016. Then, the residuals between the reference values and the Saastamoinen modeled ZHDs were calculated, and the correlations between the residuals and meteorological parameters were explored. The continuous wavelet transform method was used to recognize the annual and semiannual characteristics of the residuals. Because of the nonlinear variation characteristics of residuals, the nonlinear least squares estimation method was introduced to establish an improved ZHD model—China Revised Zenith Hydrostatic Delay (CRZHD)—adapted for China. The accuracy of the CRZHD model was assessed using radiosonde data and International GNSS Service (IGS) data in 2017; the radiosonde data results show that the CRZHD model is superior to the Saastamoinen model with a 69.6% improvement. The three IGS stations with continuous meteorological data present that the BIAS and RMSE are decreased by 2.7 and 1.5 (URUM), 5.9 and 5.3 (BJFS), and 9.6 and 8.8 mm (TCMS), respectively. The performance of the CRZHD model retrieving PWV was discussed using radiosonde data in 2017. It is shown that the CRZHD model retrieving PWV (CRZHD-PWV) outperforms the Saastamoinen model (SAAS-PWV), in which the precision is improved by 44.4%. The BIAS ranged from −1 to 1 mm and RMSE ranged from 0 to 2 mm in CRZHD-PWV account for 89.0% and 95.9%, while SAAS-PWV account for 46.6% and 58.9%.
Significance Statement
Zenith hydrostatic delay (ZHD) is one of the most important parameters in Global Navigation Satellite System (GNSS) navigation and positioning and GNSS meteorology, which can be derived from a precise ZHD model due to its stability. This research established an improved ZHD model for China to obtain accurate ZHD, which is a prerequisite for pinpoint precipitable water vapor (PWV) retrieval. And the PWV value is beneficial to analyze the change in precipitation in some regions, forecast the short-term rainfall, and monitor the climate.
Abstract
Zenith hydrostatic delay (ZHD) is a crucial parameter in Global Navigation Satellite System (GNSS) navigation and positioning and GNSS meteorology. Since the Saastamoinen ZHD model has a larger error in China, it is significant to improve the Saastamoinen ZHD model. This work first estimated the Saastamoinen model using the integrated ZHD as reference values obtained from radiosonde data collected at 73 stations in China from 2012 to 2016. Then, the residuals between the reference values and the Saastamoinen modeled ZHDs were calculated, and the correlations between the residuals and meteorological parameters were explored. The continuous wavelet transform method was used to recognize the annual and semiannual characteristics of the residuals. Because of the nonlinear variation characteristics of residuals, the nonlinear least squares estimation method was introduced to establish an improved ZHD model—China Revised Zenith Hydrostatic Delay (CRZHD)—adapted for China. The accuracy of the CRZHD model was assessed using radiosonde data and International GNSS Service (IGS) data in 2017; the radiosonde data results show that the CRZHD model is superior to the Saastamoinen model with a 69.6% improvement. The three IGS stations with continuous meteorological data present that the BIAS and RMSE are decreased by 2.7 and 1.5 (URUM), 5.9 and 5.3 (BJFS), and 9.6 and 8.8 mm (TCMS), respectively. The performance of the CRZHD model retrieving PWV was discussed using radiosonde data in 2017. It is shown that the CRZHD model retrieving PWV (CRZHD-PWV) outperforms the Saastamoinen model (SAAS-PWV), in which the precision is improved by 44.4%. The BIAS ranged from −1 to 1 mm and RMSE ranged from 0 to 2 mm in CRZHD-PWV account for 89.0% and 95.9%, while SAAS-PWV account for 46.6% and 58.9%.
Significance Statement
Zenith hydrostatic delay (ZHD) is one of the most important parameters in Global Navigation Satellite System (GNSS) navigation and positioning and GNSS meteorology, which can be derived from a precise ZHD model due to its stability. This research established an improved ZHD model for China to obtain accurate ZHD, which is a prerequisite for pinpoint precipitable water vapor (PWV) retrieval. And the PWV value is beneficial to analyze the change in precipitation in some regions, forecast the short-term rainfall, and monitor the climate.
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.
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 R2 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 R2 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
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.
