Browse
Clouds and the Earth’s Radiant Energy System (CERES) Cloud Radiative Swath (CRS) Edition 4 Data ProductAbstract
Satellite observations from Clouds and the Earth’s Radiant Energy System (CERES) radiometers have produced over two decades of world-class data documenting time–space variations in Earth’s top-of-atmosphere (TOA) radiation budget. In addition to energy exchanges among Earth and space, climate studies require accurate information on radiant energy exchanges at the surface and within the atmosphere. The CERES Cloud Radiative Swath (CRS) data product extends the standard Single Scanner Footprint (SSF) data product by calculating a suite of radiative fluxes from the surface to TOA at the instantaneous CERES footprint scale using the NASA Langley Fu–Liou radiative transfer model. Here, we describe the CRS flux algorithm and evaluate its performance against a network of ground-based measurements and CERES TOA observations. CRS all-sky downwelling broadband fluxes show significant improvements in surface validation statistics relative to the parameterized fluxes on the SSF product, including a ∼30%–40% (∼20%) reduction in SW↓ (LW↓) root-mean-square error (RMSΔ), improved correlation coefficients, and the lowest SW↓ bias over most surface types. RMSΔ and correlation statistics improve over five different surface types under both overcast and clear-sky conditions. The global mean computed TOA outgoing LW radiation (OLR) remains within <1% (2–3 W m−2) of CERES observations, while the global mean reflected SW radiation (RSW) is excessive by ∼3.5% (∼9 W m−2) owing to cloudy-sky computation errors. As we highlight using data from two remote field campaigns, the CRS data product provides many benefits for studies requiring advanced surface radiative fluxes.
Abstract
Satellite observations from Clouds and the Earth’s Radiant Energy System (CERES) radiometers have produced over two decades of world-class data documenting time–space variations in Earth’s top-of-atmosphere (TOA) radiation budget. In addition to energy exchanges among Earth and space, climate studies require accurate information on radiant energy exchanges at the surface and within the atmosphere. The CERES Cloud Radiative Swath (CRS) data product extends the standard Single Scanner Footprint (SSF) data product by calculating a suite of radiative fluxes from the surface to TOA at the instantaneous CERES footprint scale using the NASA Langley Fu–Liou radiative transfer model. Here, we describe the CRS flux algorithm and evaluate its performance against a network of ground-based measurements and CERES TOA observations. CRS all-sky downwelling broadband fluxes show significant improvements in surface validation statistics relative to the parameterized fluxes on the SSF product, including a ∼30%–40% (∼20%) reduction in SW↓ (LW↓) root-mean-square error (RMSΔ), improved correlation coefficients, and the lowest SW↓ bias over most surface types. RMSΔ and correlation statistics improve over five different surface types under both overcast and clear-sky conditions. The global mean computed TOA outgoing LW radiation (OLR) remains within <1% (2–3 W m−2) of CERES observations, while the global mean reflected SW radiation (RSW) is excessive by ∼3.5% (∼9 W m−2) owing to cloudy-sky computation errors. As we highlight using data from two remote field campaigns, the CRS data product provides many benefits for studies requiring advanced surface radiative fluxes.
Abstract
Mobile weather radars at high frequencies (C, X, K, and W bands) often collect data using staggered pulse repetition time (PRT) or dual pulse repetition frequency (PRF) modes to extend the effective Nyquist velocity and mitigate velocity aliasing while maintaining a useful maximum unambiguous range. These processing modes produce widely dispersed “processor” dealiasing errors in radial velocity estimates. The errors can also occur in clusters in high shear areas. Removing these errors prior to quantitative analysis requires tedious manual editing and often produces “holes” or regions of missing data in high signal-to-noise areas. Here, data from three mobile weather radars were used to show that the staggered PRT errors are related to a summation of the two Nyquist velocities associated with each of the PRTs. Using observations taken during a mature mesoscale convective system, a landfalling tropical cyclone, and a tornadic supercell storm, an algorithm to automatically identify and correct staggered PRT processor errors has been developed and tested. The algorithm creates a smooth profile of Doppler velocities using a Savitzky–Golay filter independently in radius and azimuth and then combined. Errors are easily identified by comparing the velocity at each range gate to its smoothed counterpart and corrected based on specific error characteristics. The method improves past dual PRF correction methods that were less successful at correcting “grouped” errors. Given the success of the technique across low, moderate, and high radial shear regimes, the new method should improve research radar analyses by affording the ability to retain as much data as possible rather than manually or objectively removing erroneous velocities.
Abstract
Mobile weather radars at high frequencies (C, X, K, and W bands) often collect data using staggered pulse repetition time (PRT) or dual pulse repetition frequency (PRF) modes to extend the effective Nyquist velocity and mitigate velocity aliasing while maintaining a useful maximum unambiguous range. These processing modes produce widely dispersed “processor” dealiasing errors in radial velocity estimates. The errors can also occur in clusters in high shear areas. Removing these errors prior to quantitative analysis requires tedious manual editing and often produces “holes” or regions of missing data in high signal-to-noise areas. Here, data from three mobile weather radars were used to show that the staggered PRT errors are related to a summation of the two Nyquist velocities associated with each of the PRTs. Using observations taken during a mature mesoscale convective system, a landfalling tropical cyclone, and a tornadic supercell storm, an algorithm to automatically identify and correct staggered PRT processor errors has been developed and tested. The algorithm creates a smooth profile of Doppler velocities using a Savitzky–Golay filter independently in radius and azimuth and then combined. Errors are easily identified by comparing the velocity at each range gate to its smoothed counterpart and corrected based on specific error characteristics. The method improves past dual PRF correction methods that were less successful at correcting “grouped” errors. Given the success of the technique across low, moderate, and high radial shear regimes, the new method should improve research radar analyses by affording the ability to retain as much data as possible rather than manually or objectively removing erroneous velocities.
Abstract
Daily rainfall accumulation estimates have been derived from 1-min volume data collected via self-syphon rain gauges deployed in the Tropical Atmosphere–Ocean (TAO) array of oceanographic buoys. The underlying high-resolution volume data were obtained directly from the Global Tropical Moored Buoy Array (GTMBA) Project Office of NOAA/Pacific Marine Environmental Laboratory. The derived accumulations have been incorporated into the Pacific Rainfall (PACRAIN) database as estimated daily values to augment existing sea level oceanic rainfall records gathered using traditional rain gauges. They have also been included in the PACRAIN historical, monthly gridded rainfall product. The methodology presented, which employs differencing of least squares–regressed sensor levels about 0000 UTC and rain gauge syphon events, is shown to offer improved error characteristics over the methodology used to compute previously published GTMBA rain rates. In particular, the PACRAIN method yields larger coefficients of determination and smaller standard errors than the duplicated GTMBA method when applied to synthetic rainfall data with noise magnitude and decorrelation times encompassing those observed in the real 1-min data. These results are shown to be consistent with mathematical expectations. Sources of instrument and catchment errors, as well as evaporation, are discussed in the context of their potential effects on accumulation estimates for periods of a day or longer.
Significance Statement
In this paper, we describe the derivation of daily rainfall amounts from raw rain gauge data obtained from buoy-mounted rain gauges. These new accumulation estimates expand the store of rainfall estimates from locations approximating the open-ocean conditions of the tropical Pacific Ocean. The derivation technique we describe exhibits better performance than the method used to generate previously published estimates using the same dataset.
Abstract
Daily rainfall accumulation estimates have been derived from 1-min volume data collected via self-syphon rain gauges deployed in the Tropical Atmosphere–Ocean (TAO) array of oceanographic buoys. The underlying high-resolution volume data were obtained directly from the Global Tropical Moored Buoy Array (GTMBA) Project Office of NOAA/Pacific Marine Environmental Laboratory. The derived accumulations have been incorporated into the Pacific Rainfall (PACRAIN) database as estimated daily values to augment existing sea level oceanic rainfall records gathered using traditional rain gauges. They have also been included in the PACRAIN historical, monthly gridded rainfall product. The methodology presented, which employs differencing of least squares–regressed sensor levels about 0000 UTC and rain gauge syphon events, is shown to offer improved error characteristics over the methodology used to compute previously published GTMBA rain rates. In particular, the PACRAIN method yields larger coefficients of determination and smaller standard errors than the duplicated GTMBA method when applied to synthetic rainfall data with noise magnitude and decorrelation times encompassing those observed in the real 1-min data. These results are shown to be consistent with mathematical expectations. Sources of instrument and catchment errors, as well as evaporation, are discussed in the context of their potential effects on accumulation estimates for periods of a day or longer.
Significance Statement
In this paper, we describe the derivation of daily rainfall amounts from raw rain gauge data obtained from buoy-mounted rain gauges. These new accumulation estimates expand the store of rainfall estimates from locations approximating the open-ocean conditions of the tropical Pacific Ocean. The derivation technique we describe exhibits better performance than the method used to generate previously published estimates using the same dataset.
Bias and Trend Correction of Precipitation Datasets to Force Ocean ModelsAbstract
A novel method to adjust the precipitation produced by atmospheric reanalyses using observational constraints to force ocean models is described. The method allows the preservation of the qualities of the high-resolution and high-frequency output from the reanalyses while eliminating their bias and spurious trends. The method is shown to be robust to degradation in both space and time of the observation dataset. This method is applied to the ERA-Interim precipitation dataset using the Global Precipitation Climatology Project (GPCP) v2.3 as the observational reference in order to create a debiased dataset that can be used to force ocean models. The produced debiased dataset is then compared to ERA-Interim and GPCP in a suite of forced ice–ocean numerical experiments using the GFDL OM4 model. Ocean states obtained with the new precipitation dataset are consistent with results from GPCP-forced experiments with respect to global metrics but produces the extra sea surface salinity variability at the time scales unresolved by the observation-based dataset. Discrepancies between modeled and observed freshwater fluxes are discussed as well as the strategies to mitigate them and their impacts.
Abstract
A novel method to adjust the precipitation produced by atmospheric reanalyses using observational constraints to force ocean models is described. The method allows the preservation of the qualities of the high-resolution and high-frequency output from the reanalyses while eliminating their bias and spurious trends. The method is shown to be robust to degradation in both space and time of the observation dataset. This method is applied to the ERA-Interim precipitation dataset using the Global Precipitation Climatology Project (GPCP) v2.3 as the observational reference in order to create a debiased dataset that can be used to force ocean models. The produced debiased dataset is then compared to ERA-Interim and GPCP in a suite of forced ice–ocean numerical experiments using the GFDL OM4 model. Ocean states obtained with the new precipitation dataset are consistent with results from GPCP-forced experiments with respect to global metrics but produces the extra sea surface salinity variability at the time scales unresolved by the observation-based dataset. Discrepancies between modeled and observed freshwater fluxes are discussed as well as the strategies to mitigate them and their impacts.
Abstract
Vertical velocities knowledge is essential to study fine-scale dynamics in the surface layers of the ocean and to understand their impact on biological production mechanisms. However, these vertical velocities have long been neglected, simply parameterized, or considered as not measurable, due mainly to their order of magnitude (less than mm s−1 up to cm s−1), generally much lower than the one of the horizontal velocities (cm s−1 to dm s−1), hence the challenge of their in situ measurement. In this paper, we present an upgraded method for direct in situ measurement of vertical velocities using data from different acoustic Doppler current profilers (ADCPs) associated with CTD probes, and we perform a comparative analysis of the results obtained by this method. The analyzed data were collected during the FUMSECK cruise, from three ADCPs: two Workhorse (conventional ADCPs), one lowered on a carousel and the other deployed in free-fall mode, and one Sentinel V (a new-generation ADCP with four classical beams and a fifth vertical beam), also lowered on a carousel. Our analyses provide profiles of vertical velocities on the order of mm s−1, as expected, with standard deviations of a few mm s−1. While the fifth beam of the Sentinel V exhibits a better accuracy than conventional ADCPs, the free-fall technique provides a more accurate measurement compared to the carousel technique. Finally, this innovative study opens up the possibility to perform simple and direct in situ measurements of vertical velocities, coupling the free-fall technique with a five-beam ADCP.
Abstract
Vertical velocities knowledge is essential to study fine-scale dynamics in the surface layers of the ocean and to understand their impact on biological production mechanisms. However, these vertical velocities have long been neglected, simply parameterized, or considered as not measurable, due mainly to their order of magnitude (less than mm s−1 up to cm s−1), generally much lower than the one of the horizontal velocities (cm s−1 to dm s−1), hence the challenge of their in situ measurement. In this paper, we present an upgraded method for direct in situ measurement of vertical velocities using data from different acoustic Doppler current profilers (ADCPs) associated with CTD probes, and we perform a comparative analysis of the results obtained by this method. The analyzed data were collected during the FUMSECK cruise, from three ADCPs: two Workhorse (conventional ADCPs), one lowered on a carousel and the other deployed in free-fall mode, and one Sentinel V (a new-generation ADCP with four classical beams and a fifth vertical beam), also lowered on a carousel. Our analyses provide profiles of vertical velocities on the order of mm s−1, as expected, with standard deviations of a few mm s−1. While the fifth beam of the Sentinel V exhibits a better accuracy than conventional ADCPs, the free-fall technique provides a more accurate measurement compared to the carousel technique. Finally, this innovative study opens up the possibility to perform simple and direct in situ measurements of vertical velocities, coupling the free-fall technique with a five-beam ADCP.
Abstract
The Echo Classification from COnvectivity (ECCO) algorithm identifies convective and stratiform types of radar echo in three dimensions. It is based on the calculation of reflectivity texture—a combination of the intensity and the heterogeneity of the radar echoes on each horizontal plane in a 3D Cartesian volume. Reflectivity texture is translated into convectivity, which is designed to be a quantitative measure of the convective nature of each 3D radar grid point. It ranges from 0 (100% stratiform) to 1 (100% convective). By thresholding convectivity, a more traditional qualitative categorization is obtained, which classifies radar echoes as convective, mixed, or stratiform. In contrast to previous algorithms, these echo-type classifications are provided on the full 3D grid of the reflectivity field. The vertically resolved classifications, in combination with temperature data, allow for subclassifications into shallow, mid-, deep, and elevated convective features, and low, mid-, and high stratiform regions—again in three dimensions. The algorithm was validated using datasets collected over the U.S. Great Plains during the PECAN field campaign. An analysis of lightning counts shows ∼90% of lightning occurring in regions classified as convective by ECCO. A statistical comparison of ECCO echo types with the well-established GPM radar precipitation-type categories show 84% (88%) of GPM stratiform (convective) echo being classified as stratiform (convective) or mixed by ECCO. ECCO was applied to radar grids for the continental United States, the United Arab Emirates, Australia, and Europe, illustrating its robustness and adaptability to different radar grid characteristics and climatic regions.
Abstract
The Echo Classification from COnvectivity (ECCO) algorithm identifies convective and stratiform types of radar echo in three dimensions. It is based on the calculation of reflectivity texture—a combination of the intensity and the heterogeneity of the radar echoes on each horizontal plane in a 3D Cartesian volume. Reflectivity texture is translated into convectivity, which is designed to be a quantitative measure of the convective nature of each 3D radar grid point. It ranges from 0 (100% stratiform) to 1 (100% convective). By thresholding convectivity, a more traditional qualitative categorization is obtained, which classifies radar echoes as convective, mixed, or stratiform. In contrast to previous algorithms, these echo-type classifications are provided on the full 3D grid of the reflectivity field. The vertically resolved classifications, in combination with temperature data, allow for subclassifications into shallow, mid-, deep, and elevated convective features, and low, mid-, and high stratiform regions—again in three dimensions. The algorithm was validated using datasets collected over the U.S. Great Plains during the PECAN field campaign. An analysis of lightning counts shows ∼90% of lightning occurring in regions classified as convective by ECCO. A statistical comparison of ECCO echo types with the well-established GPM radar precipitation-type categories show 84% (88%) of GPM stratiform (convective) echo being classified as stratiform (convective) or mixed by ECCO. ECCO was applied to radar grids for the continental United States, the United Arab Emirates, Australia, and Europe, illustrating its robustness and adaptability to different radar grid characteristics and climatic regions.
Abstract
Using data from the airborne HIAPER Cloud Radar (HCR), a partitioning algorithm (ECCO-V) that provides vertically resolved convectivity and convective versus stratiform radar-echo classification is developed for vertically pointing radars. The algorithm is based on the calculation of reflectivity and radial velocity texture fields that measure the horizontal homogeneity of cloud and precipitation features. The texture fields are translated into convectivity, a numerical measure of the convective or stratiform nature of each data point. The convective–stratiform classification is obtained by thresholding the convectivity field. Subcategories of low, mid-, and high stratiform, shallow, mid-, deep, and elevated convective, and mixed echoes are introduced, which are based on the melting-layer and divergence-level altitudes. As the algorithm provides vertically resolved classifications, it is capable of identifying different types of vertically layered echoes, and convective features that are embedded in stratiform cloud layers. Its robustness was tested on data from four HCR field campaigns that took place in different meteorological and climatological regimes. The algorithm was adapted for use in spaceborne and ground-based radars, proving its versatility, as it is adaptable not only to different radar types and wavelengths, but also different research applications.
Abstract
Using data from the airborne HIAPER Cloud Radar (HCR), a partitioning algorithm (ECCO-V) that provides vertically resolved convectivity and convective versus stratiform radar-echo classification is developed for vertically pointing radars. The algorithm is based on the calculation of reflectivity and radial velocity texture fields that measure the horizontal homogeneity of cloud and precipitation features. The texture fields are translated into convectivity, a numerical measure of the convective or stratiform nature of each data point. The convective–stratiform classification is obtained by thresholding the convectivity field. Subcategories of low, mid-, and high stratiform, shallow, mid-, deep, and elevated convective, and mixed echoes are introduced, which are based on the melting-layer and divergence-level altitudes. As the algorithm provides vertically resolved classifications, it is capable of identifying different types of vertically layered echoes, and convective features that are embedded in stratiform cloud layers. Its robustness was tested on data from four HCR field campaigns that took place in different meteorological and climatological regimes. The algorithm was adapted for use in spaceborne and ground-based radars, proving its versatility, as it is adaptable not only to different radar types and wavelengths, but also different research applications.
The Estimation of the Total Freeboard of Arctic Sea Ice in Winter Using Passive Microwave Satellite MeasurementsAbstract
A method for estimating the total freeboard hf of the sea ice in the Arctic basin was developed in this study. To utilize the dielectric properties of microwave measurements on the sea ice freeboard, we adopted the spectral difference between the microwave frequencies (e.g., the gradient ratio). Satellite lidar altimetry data were utilized as a reference, and two pairs of gradient ratios [GR(36.5, 18.7) and GR(10.7, 6.9)] and an optional brightness temperature for first-year ice [that is, TB H (6.9)] were converted from passive microwave sensors into hf using the multiple linear regression equation. Using this method, we estimated hf without direct altimetry measurements. The developed method was evaluated using Operation IceBridge data and the relationship between the regressed hf and Operation IceBridge hf had a correlation coefficient R of 0.761 and a nearly unbiased (approximately −0.4 cm) pattern. Because passive microwave measurements are taken in the Arctic daily, the approach presented in this study has the potential to enable daily hf estimations for the Arctic.
Significance Statement
Arctic sea ice is one of the most critical indicators when monitoring climate change. Precise and continuous observations of sea ice thickness are essential to understand Arctic sea ice. This study attempts to estimate sea ice thickness using passive microwave satellite data. Passive microwave satellite observations are advantageous because of their wide spatial coverage and long-term records. Therefore, the suggested method in this study can be used for filling in gaps in coverage between sea ice thickness estimates from L-band radiometry and radar/lidar altimetry. The total freeboard is proportional to the thickness of sea ice, which is converted into thickness using the hydrostatic equation. The estimated total freeboard during two winter periods (2018/19 and 2019/20) demonstrates a plausible geographical distribution over the Arctic and indicates good agreement with airborne measurements.
Abstract
A method for estimating the total freeboard hf of the sea ice in the Arctic basin was developed in this study. To utilize the dielectric properties of microwave measurements on the sea ice freeboard, we adopted the spectral difference between the microwave frequencies (e.g., the gradient ratio). Satellite lidar altimetry data were utilized as a reference, and two pairs of gradient ratios [GR(36.5, 18.7) and GR(10.7, 6.9)] and an optional brightness temperature for first-year ice [that is, TB H (6.9)] were converted from passive microwave sensors into hf using the multiple linear regression equation. Using this method, we estimated hf without direct altimetry measurements. The developed method was evaluated using Operation IceBridge data and the relationship between the regressed hf and Operation IceBridge hf had a correlation coefficient R of 0.761 and a nearly unbiased (approximately −0.4 cm) pattern. Because passive microwave measurements are taken in the Arctic daily, the approach presented in this study has the potential to enable daily hf estimations for the Arctic.
Significance Statement
Arctic sea ice is one of the most critical indicators when monitoring climate change. Precise and continuous observations of sea ice thickness are essential to understand Arctic sea ice. This study attempts to estimate sea ice thickness using passive microwave satellite data. Passive microwave satellite observations are advantageous because of their wide spatial coverage and long-term records. Therefore, the suggested method in this study can be used for filling in gaps in coverage between sea ice thickness estimates from L-band radiometry and radar/lidar altimetry. The total freeboard is proportional to the thickness of sea ice, which is converted into thickness using the hydrostatic equation. The estimated total freeboard during two winter periods (2018/19 and 2019/20) demonstrates a plausible geographical distribution over the Arctic and indicates good agreement with airborne measurements.
Abstract
Micropulse differential absorption lidars (MPD) for water vapor, temperature, and aerosol profiling have been developed, demonstrated, and are addressing the needs of the atmospheric science community for low-cost ground-based networkable instruments capable of long-term monitoring of the lower troposphere. The MPD instruments use a diode-laser-based (DLB) architecture that can easily be adapted for a wide range of applications. In this study, a DLB direct-detection Doppler lidar based on the current MPD architecture is modeled to better understand the efficacy of the instrument for vertical wind velocity measurements, with the long-term goal of incorporating these measurements into the current network of MPD instruments. The direct-detection Doppler lidar is based on a double-edge receiver that utilizes two Fabry–Pérot interferometers and a vertical velocity retrieval that requires the ancillary measurement of the backscatter ratio, which is the ratio of the total backscatter coefficient to the molecular backscatter coefficient. The modeling in this paper accounts for the major sources of error. It indicates that the vertical velocity can be retrieved with an error of less than 0.56 m s−1 below 4 km with a 150-m range resolution and an averaging time of 5 min.
Significance Statement
Monitoring the temperature, relative humidity, and winds in the lower atmosphere is important for improving weather forecasting, particularly for severe weather such as thunderstorms. Cost-effective micropulse differential absorption lidar (MPD) instrumentation for continuous temperature and humidity monitoring has been developed and demonstrated, and its effects on weather forecasting are currently being evaluated. The modeling study described in this paper studies the feasibility of using a similar cost-effective MPD instrument architecture for monitoring vertical wind velocity in the lower atmosphere. Modeling indicates that wind velocities can be measured with less than 0.56 m s−1 accuracy and demonstrates the feasibility of adding vertical wind velocity measurements to the MPD instruments.
Abstract
Micropulse differential absorption lidars (MPD) for water vapor, temperature, and aerosol profiling have been developed, demonstrated, and are addressing the needs of the atmospheric science community for low-cost ground-based networkable instruments capable of long-term monitoring of the lower troposphere. The MPD instruments use a diode-laser-based (DLB) architecture that can easily be adapted for a wide range of applications. In this study, a DLB direct-detection Doppler lidar based on the current MPD architecture is modeled to better understand the efficacy of the instrument for vertical wind velocity measurements, with the long-term goal of incorporating these measurements into the current network of MPD instruments. The direct-detection Doppler lidar is based on a double-edge receiver that utilizes two Fabry–Pérot interferometers and a vertical velocity retrieval that requires the ancillary measurement of the backscatter ratio, which is the ratio of the total backscatter coefficient to the molecular backscatter coefficient. The modeling in this paper accounts for the major sources of error. It indicates that the vertical velocity can be retrieved with an error of less than 0.56 m s−1 below 4 km with a 150-m range resolution and an averaging time of 5 min.
Significance Statement
Monitoring the temperature, relative humidity, and winds in the lower atmosphere is important for improving weather forecasting, particularly for severe weather such as thunderstorms. Cost-effective micropulse differential absorption lidar (MPD) instrumentation for continuous temperature and humidity monitoring has been developed and demonstrated, and its effects on weather forecasting are currently being evaluated. The modeling study described in this paper studies the feasibility of using a similar cost-effective MPD instrument architecture for monitoring vertical wind velocity in the lower atmosphere. Modeling indicates that wind velocities can be measured with less than 0.56 m s−1 accuracy and demonstrates the feasibility of adding vertical wind velocity measurements to the MPD instruments.
Abstract
Observations made by weather radars play a central role in many aspects of meteorological research and forecasting. These applications commonly require that radar data be supplied on a Cartesian grid, necessitating a coordinate transformation and interpolation from the radar’s native spherical geometry using a process known as gridding. In this study, we introduce a variational gridding method and, through a series of theoretical and real data experiments, show that it outperforms existing methods in terms of data resolution, noise filtering, spatial continuity, and more. Known problems with existing gridding methods (Cressman weighted average and nearest neighbor/linear interpolation) are also underscored, suggesting the potential for substantial improvements in many applications involving gridded radar data, including operational forecasting, hydrological retrievals, and three-dimensional wind retrievals.
Abstract
Observations made by weather radars play a central role in many aspects of meteorological research and forecasting. These applications commonly require that radar data be supplied on a Cartesian grid, necessitating a coordinate transformation and interpolation from the radar’s native spherical geometry using a process known as gridding. In this study, we introduce a variational gridding method and, through a series of theoretical and real data experiments, show that it outperforms existing methods in terms of data resolution, noise filtering, spatial continuity, and more. Known problems with existing gridding methods (Cressman weighted average and nearest neighbor/linear interpolation) are also underscored, suggesting the potential for substantial improvements in many applications involving gridded radar data, including operational forecasting, hydrological retrievals, and three-dimensional wind retrievals.
