Browse
Multisensor Agile Adaptive Sampling (MAAS): A Methodology to Collect Radar Observations of Convective Cell Life CycleAbstract
Multisensor Agile Adaptive Sampling (MAAS), a smart sensing framework, was adapted to increase the likelihood of observing the vertical structure (with little to no gaps), spatial variability (at subkilometer scale), and temporal evolution (at ∼2-min resolution) of convective cells. This adaptation of MAAS guided two mechanically scanning C-band radars (CSAPR2 and CHIVO) by automatically analyzing the latest NEXRAD data to identify, characterize, track, and nowcast the location of all convective cells forming in the Houston domain. MAAS used either a list of predetermined rules or real-time user input to select a convective cell to be tracked and sampled by the C-band radars. The CSAPR2 tracking radar was first tasked to collect three sector plan position indicator (PPI) scans toward the selected cell. Edge computer processing of the PPI scans was used to identify additional targets within the selected cell. In less than 2 min, both the CSAPR2 and CHIVO radars were able to collect bundles of three to six range–height indicator (RHI) scans toward different targets of interest within the selected cell. Bundles were successively collected along the path of cell advection for as long as the cell met a predetermined set of criteria. Between 1 June and 30 September 2022 over 315 000 vertical cross-section observations were collected by the C-band radars through ∼1300 unique isolated convective cells, most of which were observed for over 15 min of their life cycle. To the best of our knowledge, this dataset, collected primarily through automatic means, constitutes the largest dataset of its kind.
Abstract
Multisensor Agile Adaptive Sampling (MAAS), a smart sensing framework, was adapted to increase the likelihood of observing the vertical structure (with little to no gaps), spatial variability (at subkilometer scale), and temporal evolution (at ∼2-min resolution) of convective cells. This adaptation of MAAS guided two mechanically scanning C-band radars (CSAPR2 and CHIVO) by automatically analyzing the latest NEXRAD data to identify, characterize, track, and nowcast the location of all convective cells forming in the Houston domain. MAAS used either a list of predetermined rules or real-time user input to select a convective cell to be tracked and sampled by the C-band radars. The CSAPR2 tracking radar was first tasked to collect three sector plan position indicator (PPI) scans toward the selected cell. Edge computer processing of the PPI scans was used to identify additional targets within the selected cell. In less than 2 min, both the CSAPR2 and CHIVO radars were able to collect bundles of three to six range–height indicator (RHI) scans toward different targets of interest within the selected cell. Bundles were successively collected along the path of cell advection for as long as the cell met a predetermined set of criteria. Between 1 June and 30 September 2022 over 315 000 vertical cross-section observations were collected by the C-band radars through ∼1300 unique isolated convective cells, most of which were observed for over 15 min of their life cycle. To the best of our knowledge, this dataset, collected primarily through automatic means, constitutes the largest dataset of its kind.
Abstract
The third edition of the CM SAF Cloud, Albedo and Surface Radiation dataset from AVHRR data (CLARA-A3) contains for the first time the top-of-atmosphere products reflected solar flux (RSF) and outgoing longwave radiation (OLR), which are presented and validated using CERES, HIRS, and ERA5 reference data. The products feature an unprecedented resolution (0.25°) and time span (4 decades) and offer synergy and compatibility with other CLARA-A3 products. The RSF is relatively stable; its bias with respect to (w.r.t.) ERA5 remains mostly within ±2 W m−2. Deviations are predominantly caused by absence of either morning or afternoon satellite, mostly during the first decade. The radiative impact of the Pinatubo volcanic eruption is estimated at 3 W m−2. The OLR is stable w.r.t. ERA5 and HIRS, except during 1979–80. OLR regional uncertainty w.r.t. HIRS is quantified by the mean absolute bias (MAB) and correlates with observation density and time (satellite orbital configuration), which is optimal during 2002–16, with monthly and daily MAB of approximately 1.5 and 3.5 W m−2, respectively. Daily OLR uncertainty is higher (MAB +40%) during periods with only morning or only afternoon observations (1979–87). During the CERES era (2000–20), the OLR uncertainties w.r.t. CERES-EBAF, CERES-SYN, and HIRS are very similar. The RSF uncertainty achieves optimal results during 2002–16 with a monthly MAB w.r.t. CERES-EBAF of ∼2 W m−2 and a daily MAB w.r.t. CERES-SYN of ∼5 W m−2, and it is more sensitive to orbital configuration than is OLR. Overall, validation results are satisfactory for this first release of TOA flux products in the CLARA-A3 portfolio.
Abstract
The third edition of the CM SAF Cloud, Albedo and Surface Radiation dataset from AVHRR data (CLARA-A3) contains for the first time the top-of-atmosphere products reflected solar flux (RSF) and outgoing longwave radiation (OLR), which are presented and validated using CERES, HIRS, and ERA5 reference data. The products feature an unprecedented resolution (0.25°) and time span (4 decades) and offer synergy and compatibility with other CLARA-A3 products. The RSF is relatively stable; its bias with respect to (w.r.t.) ERA5 remains mostly within ±2 W m−2. Deviations are predominantly caused by absence of either morning or afternoon satellite, mostly during the first decade. The radiative impact of the Pinatubo volcanic eruption is estimated at 3 W m−2. The OLR is stable w.r.t. ERA5 and HIRS, except during 1979–80. OLR regional uncertainty w.r.t. HIRS is quantified by the mean absolute bias (MAB) and correlates with observation density and time (satellite orbital configuration), which is optimal during 2002–16, with monthly and daily MAB of approximately 1.5 and 3.5 W m−2, respectively. Daily OLR uncertainty is higher (MAB +40%) during periods with only morning or only afternoon observations (1979–87). During the CERES era (2000–20), the OLR uncertainties w.r.t. CERES-EBAF, CERES-SYN, and HIRS are very similar. The RSF uncertainty achieves optimal results during 2002–16 with a monthly MAB w.r.t. CERES-EBAF of ∼2 W m−2 and a daily MAB w.r.t. CERES-SYN of ∼5 W m−2, and it is more sensitive to orbital configuration than is OLR. Overall, validation results are satisfactory for this first release of TOA flux products in the CLARA-A3 portfolio.
Swish Floats: An Inexpensive Neutrally Buoyant Float to Monitor Dispersion in Coastal SeasAbstract
Neutrally buoyant floats have been widely used to measure flows in the ocean, but deploying them in large numbers can be costly and impractical. This is particularly true near coastlines due to the elevated risk of instrument grounding or vessel collisions, resulting in a lack of subsurface Lagrangian measurements in coastal regions. Here, we describe an inexpensive neutrally buoyant satellite-tracked float (named “Swallow-ish,” or “Swish” floats) that has been designed and tested as a cost-effective strategy to measure subsurface dispersion in coastal areas on time scales up to a month. These autonomous instruments are inexpensive, constructed at a material cost of CAD $300 per unit; lightweight, with a mass of 5.4 kg; isopycnal; and constructed from commercially available components, using recently available global navigation satellite system technology to provide the user with a point-to-point measure of subsurface transport. We describe the float design, ballasting techniques, and the governing equations that determine their behavior. Further, through 29 deployments in two coastal seas, we calculate an uncertainty budget and determine a ballasting error of ±1.6 g, corresponding to a local depth targeting error of 16–30 m, analyze the float resurfacing data to calculate subsurface dispersion coefficients, and examine the float depth records to quantify the local internal wave field. Finally, we evaluate surface dispersion using the postresurfacing trajectories. Our findings indicate that Swish floats offer a cost-effective alternative for Lagrangian measurements of subsurface flows in coastal regions.
Abstract
Neutrally buoyant floats have been widely used to measure flows in the ocean, but deploying them in large numbers can be costly and impractical. This is particularly true near coastlines due to the elevated risk of instrument grounding or vessel collisions, resulting in a lack of subsurface Lagrangian measurements in coastal regions. Here, we describe an inexpensive neutrally buoyant satellite-tracked float (named “Swallow-ish,” or “Swish” floats) that has been designed and tested as a cost-effective strategy to measure subsurface dispersion in coastal areas on time scales up to a month. These autonomous instruments are inexpensive, constructed at a material cost of CAD $300 per unit; lightweight, with a mass of 5.4 kg; isopycnal; and constructed from commercially available components, using recently available global navigation satellite system technology to provide the user with a point-to-point measure of subsurface transport. We describe the float design, ballasting techniques, and the governing equations that determine their behavior. Further, through 29 deployments in two coastal seas, we calculate an uncertainty budget and determine a ballasting error of ±1.6 g, corresponding to a local depth targeting error of 16–30 m, analyze the float resurfacing data to calculate subsurface dispersion coefficients, and examine the float depth records to quantify the local internal wave field. Finally, we evaluate surface dispersion using the postresurfacing trajectories. Our findings indicate that Swish floats offer a cost-effective alternative for Lagrangian measurements of subsurface flows in coastal regions.
Abstract
We present a new global oxygen atlas. This atlas uses all of the available full water column profiles of oxygen, salinity, and temperature available as part of the World Ocean Database released in 2018. Instead of optimal interpolation, we use the Data Interpolating Variational Analysis (DIVA) approach to map the available profiles onto 108 depth levels between the surface and 6800 m, covering more than 99% of ocean volume. This 1/2° × 1/2° atlas covers the period 1955–2018 in 1-yr intervals. The DIVA method has significant benefits over traditional optimal interpolation. It allows the explicit inclusion of advection and boundary constraints, thus offering improvements in the representations of oxygen, salinity, and temperature in regions of strong flow and near coastal boundaries. We demonstrate these benefits of this mapping approach with some examples from this atlas. We can explore the regional and temporal variations of oxygen in the global oceans. Preliminary analyses confirm earlier analyses that the oxygen minimum zone in the eastern Pacific Ocean has expanded and intensified. Oxygen inventory changes between 1970 and 2010 are assessed and compared against prior studies. We find that the full ocean oxygen inventory decreased by 0.84% ± 0.42%. For this period, temperature-driven solubility changes explain about 21% of the oxygen decline over the full water column; in the upper 100 m, solubility changes can explain all of the oxygen decrease; for the 100–600 m depth range, it can explain only 29%, 19% between 600 and 1000 m, and just 11% in the deep ocean.
Significance Statement
The purpose of this study is to create a new oxygen atlas of the world’s oceans using a technique that better represents the effects of ocean currents and topographic boundaries, and to investigate how oxygen in the ocean has changed over recent decades. We find the total quantity of oxygen in the world’s oceans has decreased by 0.84% since 1970, similar to previous studies. We also examine how much of this change can be explained by changes in water temperature; we find that this can explain all the changes in the upper 100 m but only 21% of the oxygen decline over the whole water column.
Abstract
We present a new global oxygen atlas. This atlas uses all of the available full water column profiles of oxygen, salinity, and temperature available as part of the World Ocean Database released in 2018. Instead of optimal interpolation, we use the Data Interpolating Variational Analysis (DIVA) approach to map the available profiles onto 108 depth levels between the surface and 6800 m, covering more than 99% of ocean volume. This 1/2° × 1/2° atlas covers the period 1955–2018 in 1-yr intervals. The DIVA method has significant benefits over traditional optimal interpolation. It allows the explicit inclusion of advection and boundary constraints, thus offering improvements in the representations of oxygen, salinity, and temperature in regions of strong flow and near coastal boundaries. We demonstrate these benefits of this mapping approach with some examples from this atlas. We can explore the regional and temporal variations of oxygen in the global oceans. Preliminary analyses confirm earlier analyses that the oxygen minimum zone in the eastern Pacific Ocean has expanded and intensified. Oxygen inventory changes between 1970 and 2010 are assessed and compared against prior studies. We find that the full ocean oxygen inventory decreased by 0.84% ± 0.42%. For this period, temperature-driven solubility changes explain about 21% of the oxygen decline over the full water column; in the upper 100 m, solubility changes can explain all of the oxygen decrease; for the 100–600 m depth range, it can explain only 29%, 19% between 600 and 1000 m, and just 11% in the deep ocean.
Significance Statement
The purpose of this study is to create a new oxygen atlas of the world’s oceans using a technique that better represents the effects of ocean currents and topographic boundaries, and to investigate how oxygen in the ocean has changed over recent decades. We find the total quantity of oxygen in the world’s oceans has decreased by 0.84% since 1970, similar to previous studies. We also examine how much of this change can be explained by changes in water temperature; we find that this can explain all the changes in the upper 100 m but only 21% of the oxygen decline over the whole water column.
Abstract
Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇ HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇ HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇ HN| below 8 km.
Significance Statement
These results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.
Abstract
Estimation of uncertainties (random error statistics) of radio occultation (RO) observations is important for their effective assimilation in numerical weather prediction (NWP) models. Average uncertainties can be estimated for large samples of RO observations and these statistics may be used for specifying the observation errors in NWP data assimilation. However, the uncertainties of individual RO observations vary, and so using average uncertainty estimates will overestimate the uncertainties of some observations and underestimate those of others, reducing their overall effectiveness in the assimilation. Several parameters associated with RO observations or their atmospheric environments have been proposed to estimate individual RO errors. These include the standard deviation of bending angle (BA) departures from either climatology in the upper stratosphere and lower mesosphere (STDV) or the sample mean between 40 and 60 km (STD4060), the local spectral width (LSW), and the magnitude of the horizontal gradient of refractivity (|∇ HN|). In this paper we show how the uncertainties of two RO datasets, COSMIC-2 and Spire BA, as well as their combination, vary with these parameters. We find that the uncertainties are highly correlated with STDV and STD4060 in the stratosphere, and with LSW and |∇ HN| in the lower troposphere. These results suggest a hybrid error model for individual BA observations that uses an average statistical model of RO errors modified by STDV or STD4060 above 30 km, and LSW or |∇ HN| below 8 km.
Significance Statement
These results contribute to the understanding of the sources of uncertainties in radio occultation observations. They could be used to improve the effectiveness of these observations in their assimilation into numerical weather prediction and reanalysis models by improving the estimation of their observational errors.
Abstract
Turbulence is what we want to avoid the most during flight. Numerical weather prediction (NWP) model–based methods for diagnosing turbulence have offered valuable guidance for pilots. NWP-based turbulence diagnostics show high accuracy in detecting turbulence in general. However, there is still room for improvements such as capturing convectively induced turbulence. In such cases, observation data can be beneficial to correctly locate convective regions and help provide corresponding turbulence information. Geostationary satellite data are commonly used for upper-level turbulence detection by utilizing its water vapor band information. The Geostationary Operational Environmental Satellite (GOES)-16 carries the Advanced Baseline Imager (ABI), which enables us to observe further down into the atmosphere with improved spatial, temporal, and spectral resolutions. Its three water vapor bands allow us to observe different vertical parts of the atmosphere, and from its infrared window bands, convective activity can be inferred. Such multispectral information from ABI can be helpful in inferring turbulence intensity at different vertical levels. This study develops U-Net based machine learning models that take ABI imagery as inputs to estimate turbulence intensity at three vertical levels: 10–18, 18–24, and above 24 kft (1 kft ≈ 300 m). Among six different U-Net-based models, U-Net3+ model with a filter size of three showed the best performance against the pilot report (PIREP). Two case studies are presented to show the strengths and weaknesses of the U-Net3+ model. The results tend to be overestimated above 24 kft, but estimates of 10–18 and 18–24 kft agree well with the PIREP, especially near convective regions.
Significance Statement
Turbulence is directly related to aviation safety as well as cost-effective aircraft operation. To avoid turbulence, turbulence diagnostics are calculated from numerical weather prediction (NWP) model outputs and are provided to pilots. The goal of this study is to develop a satellite data–driven machine learning model that estimates turbulence intensity in three different vertical layers to provide additional information along with the NWP-based turbulence diagnostics. Validation results against pilot reports show that the machine learning model performs comparable to NWP-based turbulence diagnostics. Furthermore, results with different channel selections reveal that using multiple water vapor channels can help extract additional information for estimating turbulence intensity at lower levels.
Abstract
Turbulence is what we want to avoid the most during flight. Numerical weather prediction (NWP) model–based methods for diagnosing turbulence have offered valuable guidance for pilots. NWP-based turbulence diagnostics show high accuracy in detecting turbulence in general. However, there is still room for improvements such as capturing convectively induced turbulence. In such cases, observation data can be beneficial to correctly locate convective regions and help provide corresponding turbulence information. Geostationary satellite data are commonly used for upper-level turbulence detection by utilizing its water vapor band information. The Geostationary Operational Environmental Satellite (GOES)-16 carries the Advanced Baseline Imager (ABI), which enables us to observe further down into the atmosphere with improved spatial, temporal, and spectral resolutions. Its three water vapor bands allow us to observe different vertical parts of the atmosphere, and from its infrared window bands, convective activity can be inferred. Such multispectral information from ABI can be helpful in inferring turbulence intensity at different vertical levels. This study develops U-Net based machine learning models that take ABI imagery as inputs to estimate turbulence intensity at three vertical levels: 10–18, 18–24, and above 24 kft (1 kft ≈ 300 m). Among six different U-Net-based models, U-Net3+ model with a filter size of three showed the best performance against the pilot report (PIREP). Two case studies are presented to show the strengths and weaknesses of the U-Net3+ model. The results tend to be overestimated above 24 kft, but estimates of 10–18 and 18–24 kft agree well with the PIREP, especially near convective regions.
Significance Statement
Turbulence is directly related to aviation safety as well as cost-effective aircraft operation. To avoid turbulence, turbulence diagnostics are calculated from numerical weather prediction (NWP) model outputs and are provided to pilots. The goal of this study is to develop a satellite data–driven machine learning model that estimates turbulence intensity in three different vertical layers to provide additional information along with the NWP-based turbulence diagnostics. Validation results against pilot reports show that the machine learning model performs comparable to NWP-based turbulence diagnostics. Furthermore, results with different channel selections reveal that using multiple water vapor channels can help extract additional information for estimating turbulence intensity at lower levels.
Abstract
Accurate, timely, and accessible meteorological and soil moisture measurements are essential for a number of applications including weather forecasting, agricultural decision-making, and flood and drought prediction. Such data are becoming increasingly available globally, but the large number of networks and various data reporting formats often make utilization of such data difficult. The TexMesonet is a “network of networks” developed within the state of Texas to collect, process, and make public data collected from more than 1700 monitoring stations throughout the state. This paper describes the TexMesonet, with special attention paid to monitoring sites installed and managed by the Texas Water Development Board. It also provides a case study exemplifying how these data may be used and gives recommendations for future data applications.
Abstract
Accurate, timely, and accessible meteorological and soil moisture measurements are essential for a number of applications including weather forecasting, agricultural decision-making, and flood and drought prediction. Such data are becoming increasingly available globally, but the large number of networks and various data reporting formats often make utilization of such data difficult. The TexMesonet is a “network of networks” developed within the state of Texas to collect, process, and make public data collected from more than 1700 monitoring stations throughout the state. This paper describes the TexMesonet, with special attention paid to monitoring sites installed and managed by the Texas Water Development Board. It also provides a case study exemplifying how these data may be used and gives recommendations for future data applications.
Abstract
Radiative transfer (RT) is a crucial but computationally expensive process in numerical weather/climate prediction. We develop neural networks (NN) to emulate a common RT parameterization called the Rapid Radiative Transfer Model (RRTM), with the goal of creating a faster parameterization for the Global Forecast System (GFS) v16. In previous work we emulated a highly simplified version of the shortwave RRTM only—excluding many predictor variables, driven by Rapid Refresh forecasts interpolated to a consistent height grid, using only 30 sites in the Northern Hemisphere. In this work we emulate the full shortwave and longwave RRTM—with all predictor variables, driven by GFSv16 forecasts on the native pressure–sigma grid, using data from around the globe. We experiment with NNs of widely varying complexity, including the U-net++ and U-net3+ architectures and deeply supervised training, designed to ensure realistic and accurate structure in gridded predictions. We evaluate the optimal shortwave NN and optimal longwave NN in great detail—as a function of geographic location, cloud regime, and other weather types. Both NNs produce extremely reliable heating rates and fluxes. The shortwave NN has an overall RMSE/MAE/bias of 0.14/0.08/−0.002 K day−1 for heating rate and 6.3/4.3/−0.1 W m−2 for net flux. Analogous numbers for the longwave NN are 0.22/0.12/−0.0006 K day−1 and 1.07/0.76/+0.01 W m−2. Both NNs perform well in nearly all situations, and the shortwave (longwave) NN is 7510 (90) times faster than the RRTM. Both will soon be tested online in the GFSv16.
Significance Statement
Radiative transfer is an important process for weather and climate. Accurate radiative transfer models exist, such as the RRTM, but these models are computationally slow. We develop neural networks (NNs), a type of machine learning model that is often computationally fast after training, to mimic the RRTM. We wish to accelerate the RRTM by orders of magnitude without sacrificing much accuracy. We drive both the NNs and RRTM with data from the GFSv16, an operational weather model, using locations around the globe during all seasons. We show that the NNs are highly accurate and much faster than the RRTM, which suggests that the NNs could be used to solve radiative transfer inside the GFSv16.
Abstract
Radiative transfer (RT) is a crucial but computationally expensive process in numerical weather/climate prediction. We develop neural networks (NN) to emulate a common RT parameterization called the Rapid Radiative Transfer Model (RRTM), with the goal of creating a faster parameterization for the Global Forecast System (GFS) v16. In previous work we emulated a highly simplified version of the shortwave RRTM only—excluding many predictor variables, driven by Rapid Refresh forecasts interpolated to a consistent height grid, using only 30 sites in the Northern Hemisphere. In this work we emulate the full shortwave and longwave RRTM—with all predictor variables, driven by GFSv16 forecasts on the native pressure–sigma grid, using data from around the globe. We experiment with NNs of widely varying complexity, including the U-net++ and U-net3+ architectures and deeply supervised training, designed to ensure realistic and accurate structure in gridded predictions. We evaluate the optimal shortwave NN and optimal longwave NN in great detail—as a function of geographic location, cloud regime, and other weather types. Both NNs produce extremely reliable heating rates and fluxes. The shortwave NN has an overall RMSE/MAE/bias of 0.14/0.08/−0.002 K day−1 for heating rate and 6.3/4.3/−0.1 W m−2 for net flux. Analogous numbers for the longwave NN are 0.22/0.12/−0.0006 K day−1 and 1.07/0.76/+0.01 W m−2. Both NNs perform well in nearly all situations, and the shortwave (longwave) NN is 7510 (90) times faster than the RRTM. Both will soon be tested online in the GFSv16.
Significance Statement
Radiative transfer is an important process for weather and climate. Accurate radiative transfer models exist, such as the RRTM, but these models are computationally slow. We develop neural networks (NNs), a type of machine learning model that is often computationally fast after training, to mimic the RRTM. We wish to accelerate the RRTM by orders of magnitude without sacrificing much accuracy. We drive both the NNs and RRTM with data from the GFSv16, an operational weather model, using locations around the globe during all seasons. We show that the NNs are highly accurate and much faster than the RRTM, which suggests that the NNs could be used to solve radiative transfer inside the GFSv16.
Abstract
The interaction of airflow with complex terrain has the potential to significantly amplify extreme precipitation events and modify the structure and intensity of precipitating cloud systems. However, understanding and forecasting such events is challenging, in part due to the scarcity of direct in situ measurements. Doppler radar can provide the capability to monitor extreme rainfall events over land, but our understanding of airflow modulated by orographic interactions remains limited. The SAMURAI software is a three-dimensional variational data assimilation (3DVAR) technique that uses the finite element approach to retrieve kinematic and thermodynamic fields. The analysis has high fidelity to observations when retrieving flows over a flat surface, but the capability of imposing topography as a boundary constraint is not previously implemented. Here, we implement the immersed boundary method (IBM) as pseudo-observations at their native coordinates in SAMURAI to represent the topographic forcing and surface impermeability. In this technique, neither data interpolation onto a Cartesian grid nor explicit physical constraint integration during the cost function minimization is needed. Furthermore, the physical constraints are treated as pseudo-observations, offering the flexibility to adjust the strength of the boundary condition. A series of observing simulation sensitivity experiments (OSSEs) using a full-physics model and radar emulator simulating rainfall from Typhoon Chanthu (2021) over Taiwan are conducted to evaluate the retrieval accuracy and parameter settings. The OSSE results show that the strength of the IBM constraints can impact the overall wind retrievals. Analysis from real radar observations further demonstrates that the improved retrieval technique can advance scientific analyses for the underlying dynamics of orographic precipitation using radar observations.
Abstract
The interaction of airflow with complex terrain has the potential to significantly amplify extreme precipitation events and modify the structure and intensity of precipitating cloud systems. However, understanding and forecasting such events is challenging, in part due to the scarcity of direct in situ measurements. Doppler radar can provide the capability to monitor extreme rainfall events over land, but our understanding of airflow modulated by orographic interactions remains limited. The SAMURAI software is a three-dimensional variational data assimilation (3DVAR) technique that uses the finite element approach to retrieve kinematic and thermodynamic fields. The analysis has high fidelity to observations when retrieving flows over a flat surface, but the capability of imposing topography as a boundary constraint is not previously implemented. Here, we implement the immersed boundary method (IBM) as pseudo-observations at their native coordinates in SAMURAI to represent the topographic forcing and surface impermeability. In this technique, neither data interpolation onto a Cartesian grid nor explicit physical constraint integration during the cost function minimization is needed. Furthermore, the physical constraints are treated as pseudo-observations, offering the flexibility to adjust the strength of the boundary condition. A series of observing simulation sensitivity experiments (OSSEs) using a full-physics model and radar emulator simulating rainfall from Typhoon Chanthu (2021) over Taiwan are conducted to evaluate the retrieval accuracy and parameter settings. The OSSE results show that the strength of the IBM constraints can impact the overall wind retrievals. Analysis from real radar observations further demonstrates that the improved retrieval technique can advance scientific analyses for the underlying dynamics of orographic precipitation using radar observations.
Abstract
The core Argo array has operated with the design goal of uniform spatial distribution of 3° in latitude and longitude. Recent studies have acknowledged that spatial and temporal scales of variability in some parts of the ocean are not resolved by 3° sampling and have recommended increased core Argo density in the equatorial region, boundary currents, and marginal seas with an integrated vision of other Argo variants. Biogeochemical (BGC) Argo floats currently observe the ocean from a collection of pilot arrays, but recently funded proposals will transition these pilot arrays to a global array. The current BGC Argo implementation plan recommends uniform spatial distribution of BGC Argo floats. For the first time, we estimate the effectiveness of the existing BGC Argo array to resolve the anomaly from the mean using a subset of modeled, full-depth BGC fields. We also study the effectiveness of uniformly distributed BGC Argo arrays with varying float densities at observing the ocean. Then, using previous Argo trajectories, we estimate the Argo array’s future distribution and quantify how well it observes the ocean. Finally, using a novel technique for sequentially identifying the best deployment locations, we suggest the optimal array distribution for BGC Argo floats to minimize objective mapping uncertainty in a subset of BGC fields and to best constrain BGC temporal variability.
Abstract
The core Argo array has operated with the design goal of uniform spatial distribution of 3° in latitude and longitude. Recent studies have acknowledged that spatial and temporal scales of variability in some parts of the ocean are not resolved by 3° sampling and have recommended increased core Argo density in the equatorial region, boundary currents, and marginal seas with an integrated vision of other Argo variants. Biogeochemical (BGC) Argo floats currently observe the ocean from a collection of pilot arrays, but recently funded proposals will transition these pilot arrays to a global array. The current BGC Argo implementation plan recommends uniform spatial distribution of BGC Argo floats. For the first time, we estimate the effectiveness of the existing BGC Argo array to resolve the anomaly from the mean using a subset of modeled, full-depth BGC fields. We also study the effectiveness of uniformly distributed BGC Argo arrays with varying float densities at observing the ocean. Then, using previous Argo trajectories, we estimate the Argo array’s future distribution and quantify how well it observes the ocean. Finally, using a novel technique for sequentially identifying the best deployment locations, we suggest the optimal array distribution for BGC Argo floats to minimize objective mapping uncertainty in a subset of BGC fields and to best constrain BGC temporal variability.
