Regional Atmospheric Modeling Research Articles

Assessing Raindrop Breakup Parameterizations Using Disdrometer Observations

Abstract An intercomparison of raindrop mean diameter frequency distribution (RDFD) is performed for numerical simulations of precipitating cloud systems using an array of models and microphysics schemes. This includes results from the Regional Atmospheric Modeling System (RAMS) double-moment microphysics, the Hebrew University Cloud Model bin microphysics (HUCM) interfaced to the RAMS parent model, and the Weather Research and Forecasting (WRF) Model with the Thompson, Morrison, double-moment 6-class (WDM6), and National Severe Storms Laboratory (NSSL) double-moment schemes. Simulations are examined with respect to the raindrop size distribution (DSD) volume-number mean diameter (Dm) and intercept parameter (Nw). When compared to a suite of disdrometer observations, the RDFD resulting from each microphysics scheme exhibits varying degrees of mean drop size constraints and peaks in the frequency distribution of Dm. A more detailed investigation of the peaked RDFD from the RAMS simulations suggests that the parameterization of raindrop collisional breakup can impose strong limitations on the evolution of simulated drop growth. As such, a summary and comparison of the drop breakup parameterizations among the aforementioned microphysics schemes is presented. While some drop breakup parameterizations are adjusted toward the observations by modifying the threshold diameter for the onset of breakup, this study explores the use of a modified maximum breakup efficiency. This method permits the parameterization to retain its threshold breakup diameter, while limiting the strength of drop breakup and permitting a broader range of drop sizes. As a result, the simulated mean drop sizes are in better agreement with observations.

Investigating the Role of Amazonian Mesoscale Wind Patterns and Strength on the Spatial Distribution of Martian Bedrock Exposures.

The Martian highlands contain Noachian-aged areally-extensive (>225km2) bedrock exposures that have been mapped using thermal and visible imaging datasets. Given their age, crater density and impact gardening should have led to the formation of decameter scale layers of regolith that would overlie and bury these outcrops if composed of competent materials like basaltic lavas. However, many of these regions lack thick regolith layers and show clear exposures of bedrock materials with elevated thermal inertia values compared to the global average. Hypothesized reasons for the lack of regolith include: (a) relatively weaker material properties than lavas, where friable materials are comminuted and deflated during wind erosion, (b) long-term protection from regolith development through burial and later exhumation through one or more surface processes, and (c) spatially concentrated aeolian erosion and wind energetics on well-lithified basaltic substrates. To test the third hypothesis, we used the Mars Regional Atmospheric Modeling System to calculate wind erosive strength at 10 regions throughout the Martian highlands and compared it to their thermophysical properties by using thermal infrared data derived from the Thermal Emission Spectrometer to understand the effect that Amazonian mesoscale wind patterns may have on the exposure of bedrock. We also investigated the effect of planet obliquity, Ls of perihelion, and atmospheric mean pressure on wind erosion potential. We found no evidence for increased aeolian activity over bedrock-containing regions relative to surrounding terrains, including at the mafic floor unit at Jezero crater (Máaz formation), supporting the first or second hypotheses for these regions.

Evaluation of future wind power potential and their projected changes in the Black Sea and possible stable locations for wind farms

This study aimed to evaluate the possible impacts of climate change on the wind power potential and assign possible stable locations in the Black Sea until the end of the 21st century. The wind fields simulated by a regional climate model (Rossby Centre regional atmospheric climate model, version RCA4) were analyzed considering two future periods (2021–2060, 2061–2100) under RCP4.5 and RCP8.5 scenarios. The temporal variability of the wind power density (WPD) was investigated using both spatial and local analyses for the future period (2021–2100). Three temporal variability indices (namely the coefficient of variation (CV), monthly (MV), and seasonal variability (SV)) were in agreement that the mean WPD in the eastern basin is significantly more variable than in other parts of the basin. Considering the impact of climate change, the future projections for both climate scenarios indicate changes, and the eastern basin will experience more changes under the RCP8.5 scenario. For the most predictable future developments of the wind power potential, 15 reference points along the Black Sea were analyzed using the intra- and inter-annual variability of the mean WPD. In the future, strong, durable, and stable wind resources in the western basin will ensure reliable, permanent, and sustainable WPD.

Effects of Non‐Photosynthetic Vegetation on Dust Emissions

AbstractMineral dust is among the top contributors to global aerosol loads. Ability of non‐photosynthetic vegetation (NPV) to suppress dust emission has been widely acknowledged but a realistic representation of NPV has not been tested with regional‐to‐global scale models. In this study, we implemented a satellite‐based total vegetation data set, which included NPV, into a regional atmospheric chemistry model and conducted simulations for the year 2016 over the conterminous United States. To test the response of dust simulations to the NPV coverage, we conducted a control simulation incorporating only the photosynthetic vegetation (PV). Simulated dust emissions decrease by 10%–70% over most of the southwestern US from spring to autumn due to NPV. Reductions in dust concentrations are the largest in spring, which attenuate the overpredictions of fine soil concentrations, but accentuate the underpredictions in summer. Overall, the mean errors and correlations of annual simulations are slightly improved with NPV. NPV modulates dust emissions mainly by sheltering the surface and increasing the threshold velocity through drag partitioning. Moreover, we investigated the effect of vegetation height and addressed its uncertainties through a series of sensitivity tests. We observed that a 50% variation in predefined vegetation heights results in small changes in soil concentrations over majority of southwestern US, but causes up to 30% changes at local hotspots. This study highlights the significance of including NPV into the dust model and points out the importance of validation of total vegetation datasets as well as more realistic representation of vegetation heights and seasonality.

Impacts of high-resolution atmospheric forcing and air-sea coupling on coastal ocean circulation off the Pearl River Estuary

We investigated the impacts of high-resolution atmospheric forcing and ocean-atmosphere coupling on the estuary-shelf ocean simulation off the Pearl River Estuary. We conducted process and dynamics analyses of the wind-driven coastal ocean circulation under atmospheric flux forcing from (1) global reanalysis data, (2) a high-resolution regional atmospheric model, and (3) an air-sea coupled model during an upwelling-favorable wind. The results revealed that the high-resolution atmospheric model significantly improved the representations of the near-surface wind field and air temperature. The air-sea coupled model outperformed the uncoupled model in simulating coastal currents, water temperature, and salinity. The high-resolution uncoupled model strengthened the surface wind stress and along-isobath pressure gradient force (PGF), resulting in an intensified cross-isobath transport. The improved wind forcing from the air-sea coupled model modulated the spatial variation of the net stress curl and vorticity advection and enhanced the along-isobath PGF for a stronger cross-isobath transport. The lower sea surface temperature (SST) in the air-sea coupled model reduced the air temperature and wind stress. Adjusting the SST in the air-sea coupled model improved the momentum fluxes and the associated ocean transport dynamics.

Wind-wave forecasting in enclosed basins using statistically downscaled global wind forcing

Accurate wind-wave forecasting in enclosed and semi-enclosed basins is a challenging task, demanding primarily for high-resolution wind forcing at regional scale. This is generally obtained with dynamical downscaling from a low-to-mid resolution atmospheric model. In this context, a new wave forecasting system for the marginal Adriatic Sea is herein presented aimed at proposing an alternative strategy for accurate wind-wave forecasting in (semi-) enclosed basins that does not require an ad-hoc regional atmospheric model. The system is based on the state-of-the-art WAVEWATCH III® spectral wave model forced by the global IFS-ECMWF forecast. At first, wind speed is quantile-corrected to account for the systematic underestimation over the Adriatic Sea. Then, the significant wave height in the target region and for regimes associated with marine storms is calibrated following standard procedure. Wind and wave observations from different sources are used for calibration and validation of the wave forecasts, which achieve satisfactory scores. We also compare results with those of other forecasting systems in the area, highlighting the importance of the wind forcing accuracy and the wave model calibration. Doing so, we discuss the challenges that characterise (semi-) enclosed environments in order to propose effective solutions for them and future developments.

A Daily 1‐km Resolution Greenland Rainfall Climatology (1958–2020) From Statistical Downscaling of a Regional Atmospheric Climate Model

AbstractA daily, gridded 1‐km rainfall climatology (1958–2020) for Greenland is presented, including the Greenland ice sheet (GrIS), the peripheral glaciers and ice caps (GIC), and ice‐free tundra. It is obtained by statistically downscaling the 5.5 km output of the regional atmospheric climate model version 2 to a resolution of 1 km, using the elevation dependence of snowfall fraction. Based on this new product, the average total annual rainfall in Greenland during 1958–2020 is estimated to be 111.4 ± 11.2 Gt/year, of which 28.6 ± 6.1 Gt/year falls on the GrIS, 11.4 ± 1.4 Gt/year on the GIC, and 71.4 ± 9.0 Gt/year on the tundra. The downscaled 1 km product better resolves especially the relatively small GIC, more than doubling (+124%) their rainfall compared to the 5.5 km product. The relatively warm southern regions of Greenland have the highest rainfall amounts, with annual values locally exceeding 1,000 mm. We confirm a significant positive trend in Greenland rainfall (>40 mm/decade), notably in the northwest and mainly due to an increase in rainfall fraction (>3.5%/decade) during July and August. For the whole of Greenland, during 1991–2020 the seasonal rainfall peak has shifted from July to August, with significant increases in September, which receives more rain than June. An analysis of rainfall fraction and near‐surface temperature shows that local warming rates are a good predictor of recent rainfall changes.

Fully Dynamic High–Resolution Model for Dispersion of Icelandic Airborne Mineral Dust

Icelandic topsoil sediments, as confirmed by numerous scientific studies, represent the largest and the most important European source of mineral dust. Strong winds, connected with the intensive cyclonic circulation in the North Atlantic, induce intense emissions of mineral dust from local sources all year and carry away these fine aerosol particles for thousands of kilometers. Various impacts of airborne mineral dust particles on local air quality, human health, transportation, climate and marine ecosystems motivated us to design a fully dynamic coupled atmosphere–dust numerical modelling system in order to simulate, predict and quantify the Icelandic mineral dust process including: local measurements and source specification over Iceland. In this study, we used the Dust Regional Atmospheric Model (DREAM) with improved Icelandic high resolution dust source specification and implemented spatially variable particle size distribution, variable snow cover and soil wetness. Three case studies of intense short- and long-range transport were selected to evaluate the model performance. Results demonstrated the model’s capability to forecast major transport features, such as timing, and horizontal and vertical distribution of the processes. This modelling system can be used as an operational forecasting system, but also as a reliable tool for assessing climate and environmental Icelandic dust impacts.

Urban ventilation assessment with improved vertical wind profile in high-density cities – Comparisons between LiDAR and conventional methods

The vertical wind speed profile is crucial to urban ventilation assessment and urban planning/design. This study uses Light Detection and Ranging (LiDAR) observation as benchmark to evaluate accuracy of wind profiles estimated by conventional methods. The conventional methods include Boundary Layer Wind Tunnel (BLWT), Regional Atmospheric Modeling System (RAMS), Power Law (PL), and Weather Research and Forecasting (WRF). The evaluation involves two typical urban sites with different densities under summer weak-wind conditions. Large Eddie Simulations (LES) are conducted to further investigate the sensitivity of urban ventilation assessment results to the deviations of wind profiles. The results indicate significant deviations in LES caused by conventional methods. The largest deviations of wind velocity ratio are found in mesoscale meteorological models (RAMS and WRF (>65%)). Deviations caused by physical and empirical models are smaller but still significant (BLWT (>25%) and PL (>40%)). Consequently, large deviations (>100%) of wind-relevant criterion for outdoor thermal comfort are observed. Finally, to balance accuracy and data availability, we recommend power law method as the optimal method to provide inflow boundary condition for numerical simulations when LiDAR observation is not available. We provide new and valuable understandings to improve urban ventilation assessment in high-density cities.

Precipitation extremes over the tropical Americas under RCP4.5 and RCP8.5 climate change scenarios: Results from dynamical downscaling simulations

AbstractThe Regional Atmospheric Modelling System (RAMS) driven by data from the CMIP5 Earth System Model HadGEM2‐ES was used to simulate daily precipitation over the tropical Americas for both current and future climate, including distinct scenarios (representative concentration pathways 4.5 and 8.5) and different time horizons (short‐term, mid‐term and long‐term changes). The major objective was to evaluate possible future changes in extreme events, with emphasis on the intensity of precipitation events and the duration of wet and dry spells. According to RAMS, longer dry spells are expected over most regions of the tropical Americas in the future, with indications for Northeast Brazil, Caribbean, Northern Amazon, and shorter wet spells over Central America and Amazon. With the exception of the Caribbean, there is a general tendency towards the increased frequency of intense precipitation in the tropical Americas.

Evaluation of regional climate models and future wave characteristics in an enclosed sea: A case study of the Black Sea

The main objective of this study is to evaluate a possible climate change impact on wave climate in an enclosed sea (Black Sea) until the end of the 21st century. To achieve the goal, the wave data are generated from a third-generation wave model forced with the widely used reanalysis data (Climate Forecast System Reanalysis), and the hindcast wind fields simulated by a Regional Climate Model, namely Rossby Centre regional atmospheric model (version RCA4). The pattern of change in the wave climate is discussed by comparing the historical (1970–2005) with the future for periods 2021–2060 and 2061–2100 under two climate scenarios. The projected changes in the mean wave characteristics are characterized by an increase in the near future for both RCP scenarios that is more prominent under RCP8.5 scenario in the eastern basin. The middle future under RCP8.5 scenario suggests more significant changes in annual mean values than those obtained under RCP4.5 scenario. In the context of climate change, the possible effects are more prominent on extreme wave heights. For example, the increases in the wave heights for the 100-year return period reach the value of 5% and 25% in the western and eastern Black Sea, respectively.

High-Resolution Ammonia Emissions from Nitrogen Fertilizer Application in China during 2005–2020

The accurate estimation of ammonia emission is essential for quantifying secondary inorganic aerosol formation and reactive nitrogen deposition. During the last decades, both fertilizer type and the total amount of nitrogen fertilizer in China have changed, while the resulting changes in ammonia emissions and their spatio-temporal variations are unclear. In this study, we compile a long-term (2005–2020) high-resolution ammonia emission inventory for synthetic fertilizer in China with bottom-up method. We parameterized emissions factors (EFs) considering the impacts of soil properties, method of fertilizer application, fertilizer type, crop type, ambient temperature and wind speed. Meanwhile, the monthly nitrogen fertilizer application is calculated by detailed information on crop-specific fertilizer application practices. For the spatial distribution, the ammonia emissions from fertilizer mostly concentrate in eastern and southwestern China, coincident with the high density of agriculture activity and population in these regions. For the seasonal variation, the ammonia emissions from fertilizer application peak in spring and summer because of dense fertilizer application and high ambient temperature. For the long-term trend, we estimate that the emissions from synthetic fertilizer increased from 5.38 Tg in 2005 to 5.53 Tg in 2008 and remained nearly unchanged during 2008–2012, then decreased to 3.96 Tg in 2020. Urea, ammonium bicarbonate (ABC) and nitrogenous compound fertilizer are major fertilizer types used in China. Despite the increased use of nitrogen fertilizer, ammonia emissions remained stable throughout 2008–2012 with the declined use of ABC. This stable period also reflects ammonia emission increases in western China, offsetting the decreases in eastern China. Furthermore, our emission inventory provides a monthly estimation at a spatial resolution of 0.1 degrees, which can be applied to global and regional atmospheric chemistry model simulations.

Assessing the contribution of dynamical downscaling to austral autumn Northeast Brazil seasonal precipitation prediction performance

This study aims to assess the contribution to prediction performance of dynamically downscaling seasonal precipitation for the Brazilian Northeast region (NEB) produced in February with a global climate model for the austral autumn (March-April-May). This season includes the peak precipitation period of a large portion of this region. To achieve this aim, 32 years of hindcasts for the period 1979–2010 were produced with a regional model known as the Brazilian developments on the Regional Atmospheric Modeling System (BRAMS) run at 30 km spatial resolution, which was unidirectional nested to a global model run at approximately 200 km spatial resolution. The used global model was the Atmospheric General Circulation Models of the Brazilian Center for Weather Forecast and Climate Studies (CPTEC/AGCM), which provided initial and lateral boundary conditions for BRAMS. Hindcast ensembles of 10 members were produced with both the global and the regional models, simulating an operational prediction environment using persisted sea surface temperature (SST) anomalies observed in the previous January of each of the 32 hindcast years as oceanic boundary conditions. A comparative deterministic and probabilistic performance assessment of the two models was conducted. The deterministic assessment included the investigation of the mean error (bias) of the produced seasonal precipitation hindcasts, the degree of linear association between the observed and predicted precipitation anomalies measuring phase error through the correlation of the ensemble means with the corresponding observations, and amplitude error investigated through the ratio between the hindcast ensemble mean and observed variances. The probabilistic assessment was performed investigation the reliability, resolution, and discrimination characteristics of tercile probability predictions for the above normal and below normal categories, and also for probabilistic predictions for the occurrence of the symmetric event positive (or negative) precipitation anomaly. The two investigated modeling approaches showed similar performance, with an indication for dynamically downscaled predictions to have reduced amplitude error and slightly improved reliability and resolution compared to the predictions produced with the global model.

Development of Forecast System for the Upwelled Coastal Cold Waters in the Eastern Coast of Korea

In this study, a forecast system for the occurrence of cold waters in the eastern coast of Korea was developed by combining the ROMS ocean model and the WRF regional atmospheric model. The 1-day hindcast and 3.5-day forecasts from April 1 to August 31, 2021, successfully reproduced and forecasted the temporal and spatial changes of coastal cold waters. Three case studies with different surface boundary conditions in the ROMS model confirmed that the wind is a primary driver for the occurrence of cold waters in the region. The performance of the ROMS was improved by high-resolution WRF outputs, and through error correction using the water temperature data from the real-time stations, the average RMSE of 9 stations for Forecast Day-3 was 0.73℃, indicating a prossibility of high-precision forecast.

Developing a new wind dataset by blending satellite data and WRF model wind predictions

This paper presents an approach to improve wind datasets developed using the regional atmospheric model Weather Research Forecasting by combining its predictions with remotely sensed wind observations in enhanced wind speed analyses that leads to blended winds. In this study, satellite data derived from scatterometers, radiometers, and synthetic aperture radar are used. The spatial and temporal features of each wind product are thoroughly analysed. For the probabilistic evaluation of their skill, comprehensive comparisons with available buoy data are carried out. The statistical analysis shows that the combined use of satellite and numerical weather prediction model data improves the agreement with buoy measurements, demonstrating the added value of using the blended product. As an application of the method, new improved satellite wind speeds are presented in the form of a wind energy assessment along the Iberian coastal area. From inspection of the provided wind power maps, northern and central regions emerge as the most promising areas for wind harnessing offshore despite some seasonal variations. Finally, potential wind farm sites are provided, along with insights into multi-year wind speed distribution. The results show how the new dataset can be used for the selection of promising areas for wind exploitation.

Simulation of an orographic gravity wave above the Amazon rainforest and its influence on gases transport near the surface

The influence of the topography at the ATTO site (Amazon Tall Tower Observatory) in the formation of Gravity Waves (GW) was investigated in this work. The role of these waves in gas transport on different heights in the nocturnal boundary layer, was also investigated. For this, were analyzed: experimental data of wind speed, temperature, humidity and CO2 to identify the GW. Furthermore satellite images, ECMWF ERA5 (European Center for Medium-Range Weather Forecasts) reanalysis, and numerical simulations from the mesoscale JULES-CCATT-BRAMS model (Brazilian Regional Atmospheric Modeling System) version 5.3 were used to perform GW simulations. The simulations involved two configurations (Control-With Topography and No Topography). The sattelite images made it clear that the GW was not associated with convective activities. The simulations showed clear ondulatory signs in air temperature, vertical component of the wind, and gases like CO and O3. Moreover, they also reveal that the wave originated by the presence of the local topography and this orograph gravity wave had and important role in the transport of gases above the ATTO site.

Mars Exploration Using Sailplanes

We present the preliminary design of sailplanes, used for Mars exploration. The sailplanes mitigate the weight and energy storage limitations traditionally associated with powered flight by instead exploiting atmospheric wind gradients for dynamic soaring, and slope/thermal updrafts for static soaring. Equations of motion for the sailplanes were combined with wind profiles from the Mars Regional Atmospheric Modeling System (MRAMS) for two representative sites: Jezero crater, Perseverance’s landing site, and over a section of the Valles Marineris canyon. Optimal flight trajectories were obtained from the constrained optimization problem, using the lift coefficient and the roll angle as control parameters. Numerical results for complete dynamic soaring cycles demonstrated that the total sailplane energy at the end of a soaring cycle increases by 6.8–11%. The absence of a propulsion system, allowing for a compact form factor, means the sailplanes can be packaged into CubeSats and deployed as secondary payloads at a relatively low cost; providing scientific data over locations inaccessible by current landers and rovers. Various sailplane deployment methods are considered, including rapid deployment during Entry, Descent, and Landing (EDL) of a Mars Science Laboratory-class (MSL) vehicle and slow deployment using a blimp.

Surface Wind Responses to Mesoscale Sea Surface Temperature over Western Boundary Current Regions Assessed by Spectral Transfer Functions

Abstract Satellite observations have revealed that mesoscale sea surface temperature (SST) perturbations can exert distinct influence on sea surface wind by modifying the overlying atmospheric boundary layer. Recently, spectral transfer functions have been shown to be useful to elucidate the wind response features. Spectral transfer functions can represent spatially lagged responses, their horizontal scale dependence, and background wind speed dependence. By adopting the transfer function analysis, the present study explores seasonality and regional differences in the wind response over the major western boundary current regions. Transfer functions estimated from satellite observations are found to be largely consistent among seasons and regions, suggesting that the underlying dominant dynamics are ubiquitous. Nevertheless, the wind response exhibits statistically significant seasonal and regional differences depending on background wind speed. When background wind is stronger (weaker) than 8.5 m s−1, the wind response is stronger (weaker) in winter than in summer. The Agulhas Retroflection region exhibits stronger wind response typically by 30% than the Gulf Stream and Kuroshio Extension regions. Although observed wind distributions are reasonably reconstructed from the transfer functions and observed SST, surface wind convergence zones along the Gulf Stream and Kuroshio Extension are underrepresented. The state-of-the-art atmospheric reanalysis and regional model represent well the structure of the transfer functions in the wavenumber space. The amplitude is, however, underestimated by typically 30%. The transfer function analysis can be adapted to many other atmospheric responses besides sea surface wind, and thus provide new insights into the climatic role of the mesoscale air–sea coupling.

Assessment on the Water Vapor Flux from Atmospheric Reanalysis Data in the South China Sea on 2019 Summer

Abstract This paper assesses the water vapor flux performance of three reanalysis datasets (ERA5, JRA55, NCEP-2) on the South China Sea. The radiosonde data were from the South China Sea Scientific Expedition organized by Sun Yat-sen University in the 2019 summer (SCSEX2019). The comparison shows that all reanalyses underestimate the temperature and specific humidity under 500 hPa. As for the wind profile, the most significant difference appeared at 1800 UTC when there was no conventional radiosonde observation around the experiment area. As for the water vapor flux, ERA5 seems to give the best zonal flux but the worst meridional one. A deeper analysis shows that the bias in the wind mainly caused the difference in water vapor flux from ERA5. As for JRA55 and NCEP-2, the humidity and wind field bias coincidentally canceled each other, inducing a much smaller bias, especially in meridional water vapor flux. Therefore, to get a more realistic water vapor flux, a correction in the wind profile was most needed for ERA5. In contrast, the simultaneous improvement on both wind and humidity fields might produce a better water vapor flux for JRA55 and NCEP-2. Significance Statement This paper mainly aims to assess three atmospheric reanalyses from the viewpoint of the water vapor flux over the South China Sea during the monsoon period. The observation data contain more than 120 radiosonde profiles. Our work has given an objective comparison among the reanalyses and observations. We also tried to explain the bias in the water vapor flux over the ocean from the reanalyses. The results of our work might help understand the monsoon precipitation given by atmospheric reanalyses or regional climate models and enlighten the development of atmospheric assimilation products.

Development of the Real‐Time 30‐s‐Update Big Data Assimilation System for Convective Rainfall Prediction With a Phased Array Weather Radar: Description and Preliminary Evaluation

AbstractWe present the first ever real‐time numerical weather prediction system with 30‐s update cycles at a 500‐m grid spacing for the prediction of convective precipitation in the subsequent 30 min using a new‐generation multi‐parameter phased array weather radar. The system comprises a regional atmospheric model known as the SCALE and the local ensemble transform Kalman filter (LETKF). To accelerate the SCALE‐LETKF system, data transfer between the two aforementioned components is performed using a memory copy instead of a file I/O. A complete real‐time workflow including domain nesting and observational data transfer is constructed. A real‐time test in July and August 2020 showed that the system is fast enough for a real‐time application of 30‐s forecast‐analysis cycles and 30‐min prediction. The development includes a new thinning method considering the spatially correlated observation errors in the dense radar data. This new thinning method is effective in two past case studies in the summer of 2019.