Sensitivity of high-resolution precipitation to physics parameterization options in WRF over equatorial regions

Author(s):  
Martina Messmer ◽  
Santos J. González-Rojí ◽  
Christoph C. Raible ◽  
Thomas F. Stocker

<p>Precipitation patterns and climate variability in East Africa and Western South America present high heterogeneity and complexity. This complexity is a result of large-scale and regional controls, such as surrounding oceans, lakes and topography. The combined effect of these controls has implications on precipitation and temperature, and hence, on water availability, biodiversity and ecosystem services. This study focuses on the impact of different physics parameterization in high-resolution experiments performed over equatorial regions with the Weather Research and Forecasting (WRF) model, and how these options affect the representation of precipitation in those regions.</p><p>As expected, weather and climate in equatorial regions are driven by physical processes different to those important in the mid-latitudes. Hence, it is necessary to test the parameterizations available in the WRF model. Several sensitivity simulations are performed over Kenya and Peru nesting the WRF model inside the state-of-the-art ERA5 reanalysis. A cascade of increasing grid resolutions is used in these simulations, reaching the spatial resolutions of 3 and 1 km in the innermost domains, and thus, convection permitting scales. Parameterization options of the planetary boundary layer (PBL), lake model, radiation, cumulus and microphysics schemes are changed, and their sensitivity to precipitation is tested. The year 2008 is simulated for each of the sensitivity simulations. This year is chosen as a good representative of precipitation dynamics and temperature, as it is neither abnormally wet or hot, nor dry or cold over Kenya and Peru. The simulated precipitation driven by the ERA5 reanalysis is compared against station data obtained from the WMO, and over Kenya additionally against observations from the Centre for Training and Integrated Research in ASAL Development (CETRAD).</p><p>Precipitation is strongly underestimated when adopting a typical parameterization setup for the mid-latitudes. However, results indicate that precipitation amounts and also patterns are substantially improved when changing the cumulus and PBL parameterisations. This strong increase in the simulated precipitation is obtained when using the Grell-Freitas ensemble, RRTM and the Yonsei University schemes for cumulus, long-wave radiation and planetary boundary layer, respectively. During some summer months, the accumulated precipitation is improved by up to 100 mm (80 %) compared to mid-latitudes configuration in several regions of the domains (near the Andes in Peru and over the flatlands in Kenya). Additionally, because the 1- and 2-way nesting options show a similar performance with respect to precipitation, the 1-way nesting option is preferred, as it does not overwrite the solutions in the parent domains. Hence, discontinuous solutions related to switching off the cumulus parameterization can be avoided.</p>

2017 ◽  
Vol 10 (5) ◽  
pp. 2031-2055 ◽  
Author(s):  
Thomas Schwitalla ◽  
Hans-Stefan Bauer ◽  
Volker Wulfmeyer ◽  
Kirsten Warrach-Sagi

Abstract. Increasing computational resources and the demands of impact modelers, stake holders, and society envision seasonal and climate simulations with the convection-permitting resolution. So far such a resolution is only achieved with a limited-area model whose results are impacted by zonal and meridional boundaries. Here, we present the setup of a latitude-belt domain that reduces disturbances originating from the western and eastern boundaries and therefore allows for studying the impact of model resolution and physical parameterization. The Weather Research and Forecasting (WRF) model coupled to the NOAH land–surface model was operated during July and August 2013 at two different horizontal resolutions, namely 0.03 (HIRES) and 0.12° (LOWRES). Both simulations were forced by the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analysis data at the northern and southern domain boundaries, and the high-resolution Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) data at the sea surface.The simulations are compared to the operational ECMWF analysis for the representation of large-scale features. To analyze the simulated precipitation, the operational ECMWF forecast, the CPC MORPHing (CMORPH), and the ENSEMBLES gridded observation precipitation data set (E-OBS) were used as references.Analyzing pressure, geopotential height, wind, and temperature fields as well as precipitation revealed (1) a benefit from the higher resolution concerning the reduction of monthly biases, root mean square error, and an improved Pearson skill score, and (2) deficiencies in the physical parameterizations leading to notable biases in distinct regions like the polar Atlantic for the LOWRES simulation, the North Pacific, and Inner Mongolia for both resolutions.In summary, the application of a latitude belt on a convection-permitting resolution shows promising results that are beneficial for future seasonal forecasting.


2020 ◽  
Vol 12 (4) ◽  
pp. 3097-3112
Author(s):  
Emily Collier ◽  
Thomas Mölg

Abstract. Climate impact assessments require information about climate change at regional and ideally also local scales. In dendroecological studies, this information has traditionally been obtained using statistical methods, which preclude the linkage of local climate changes to large-scale drivers in a process-based way. As part of recent efforts to investigate the impact of climate change on forest ecosystems in Bavaria, Germany, we developed a high-resolution atmospheric modelling dataset, BAYWRF, for this region over the thirty-year period of September 1987 to August 2018. The atmospheric model employed in this study, the Weather Research and Forecasting (WRF) model, was configured with two nested domains of 7.5 and 1.5 km grid spacing centred over Bavaria and forced at the outer lateral boundaries by ERA5 reanalysis data. Using an extensive network of observational data, we evaluate (i) the impact of using grid analysis nudging for a single-year simulation of the period of September 2017 to August 2018 and (ii) the full BAYWRF dataset generated using nudging. The evaluation shows that the model represents variability in near-surface meteorological conditions generally well, although there are both seasonal and spatial biases in the dataset that interested users should take into account. BAYWRF provides a unique and valuable tool for investigating climate change in Bavaria with high interdisciplinary relevance. Data from the finest-resolution WRF domain are available for download at daily temporal resolution from a public repository at the Open Science Framework (Collier, 2020; https://doi.org/10.17605/OSF.IO/AQ58B).


2016 ◽  
Author(s):  
Thomas Schwitalla ◽  
Hans-Stefan Bauer ◽  
Volker Wulfmeyer ◽  
Kirsten Warrach-Sagi

Abstract. The impact of a convection permitting (CP) northern hemisphere latitude-belt simulation with the Weather Research and Forecasting (WRF) model was investigated during the July and August 2013. For this application, the WRF model together with the NOAH land-surface model (LSM) was applied at two different horizontal resolutions, 0.03° (HIRES) and 0.12° (LOWRES). The set-up as a latitude-belt domain avoids disturbances that originate from the western and eastern boundaries and therefore allows to study the impact of model resolution and physical parameterizations on the results. Both simulations were forced by ECMWF operational analysis data at the northern and southern domain boundaries and the high-resolution Operational Sea Surface Temperature and Sea Ice Analysis (OSTIA) data at the sea surface. The simulations are compared to the operational ECMWF analysis for the representation of large scale features. To compare the simulated precipitation, the operational ECMWF forecast, the CPC MORPHing (CMORPH), and the ENSEMBLES gridded observation precipitation data set (E-OBS) were used. Compared to the operational high-resolution ECMWF analysis, both simulations are able to capture the large scale circulation pattern though the strength of the Pacific high is considerably overestimated in the LOWRES simulation. Major differences between ECMWF and WRF occur during July 2013 when the lower resolution simulation shows a significant negative bias over the North Atlantic which is not observed in the CP simulation. The analysis indicates deficiencies in the applied combinations of cloud microphysics and convection parametrization on the coarser grid scale in subpolar regions. The overall representation of the 500 hPa geopotential height surface is also improved by the CP simulation compared to the LOWRES simulation apart across Newfoundland where the geopotential height is higher than in the LOWRES simulation due to a northward shift of the location of the Atlantic high pressure system. Both simulations show higher wind speeds in the boundary layer by about 1.5 m s−1 compared to the the ECMWF analysis. Due to the higher surface evaporation, this results in a moist bias of 0.5 g kg−1 at 925 hPa in the planetary boundary layer compared to the ECMWF analysis. Major differences between ECMWF and WRF occur in the simulation of the 2-m temperatures over the Asian desert and steppe regions. They are significantly higher in WRF by about 5 K both during day- and night-time presumably as a result of different soil hydraulic parameters used in the NOAH land surface model for steppe regions. The precipitation of the HIRES simulation shows a better spatial agreement with CMORPH especially over mountainous terrain. The overall bias reduces from 80 mm at the coarser resolution to 50 mm in the HIRES simulation and the root mean square error is reduced by about 35 % when compared to the CMORPH precipitation analysis. The precipitation distribution agrees much better with the CMORPH data than the LOWRES simulation which tends to overestimate precipitation, mainly caused by the convection parametrization. Especially over Europe the CP resolution reduces the precipitation bias by about 30 % to 20 mm as a result of a better terrain representation and due to the avoidance of the convection parameterization.


2021 ◽  
Vol 22 (5) ◽  
pp. 1169-1186
Author(s):  
Ioannis Sofokleous ◽  
Adriana Bruggeman ◽  
Silas Michaelides ◽  
Panos Hadjinicolaou ◽  
George Zittis ◽  
...  

ABSTRACTA stepwise evaluation method and a comprehensive scoring approach are proposed and implemented to select a model setup and physics parameterizations of the Weather Research and Forecasting (WRF) Model for high-resolution precipitation simulations. The ERA5 reanalysis data were dynamically downscaled to 1-km resolution for the topographically complex domain of the eastern Mediterranean island of Cyprus. The performance of the simulations was examined for three domain configurations, two model initialization approaches and 18 combinations of atmospheric physics parameterizations. Two continuous and two categorical scores were used for the evaluation. A new extreme event score, which combines hits and frequency bias, was introduced as a complementary evaluator of extremes. A composite scaled score was used to identify the overall best performing parameterizations. The least errors in mean daily and monthly precipitation amounts and daily extremes were found for the domain configuration with the largest extent and three nested domains. A 5-day initialization frequency did not improve precipitation, relative to 30-day continuous simulations. The parameterization type with the largest impact on precipitation was microphysics. The cumulus parameterization was also found to have an impact on the 1-km nested domain, despite that it was only activated in the coarser “parent” domains. Comparison of simulations with 12-, 4-, and 1-km resolution revealed the better skill of the model at 1 km. The impact of the various model configurations in the small-sized domain was different from the impact in larger model domains; this could be further explored for other atmospheric variables.


2018 ◽  
Vol 57 (11) ◽  
pp. 2679-2696 ◽  
Author(s):  
Jennifer D. Hegarty ◽  
Jasper Lewis ◽  
Erica L. McGrath-Spangler ◽  
John Henderson ◽  
Amy Jo Scarino ◽  
...  

AbstractThe daytime planetary boundary layer (PBL) was examined for the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) Baltimore (Maryland)–Washington, D.C., campaign of July 2011 using PBL height (PBLH) retrievals from aerosol backscatter measurements from ground-based micropulse lidar (MPL), the NASA Langley Research Center airborne High Spectral Resolution Lidar-1 (HSRL-1), and the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) on the Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite. High-resolution Weather Research and Forecasting (WRF) Model simulations with horizontal grid spacing of 1 km and different combinations of PBL schemes, urban parameterization, and sea surface temperature inputs were evaluated against PBLHs derived from lidars, ozonesondes, and radiosondes. MPL and WRF PBLHs depicted a growing PBL in the morning that reached a peak height by midafternoon. WRF PBLHs calculated from gridded output profiles generally showed more rapid growth and higher peak heights than did the MPLs, and all WRF–lidar differences were dependent on model configuration, PBLH calculation method, and synoptic conditions. At inland locations, WRF simulated an earlier descent of the PBL top in the afternoon relative to the MPL retrievals and radiosonde PBLHs. At Edgewood, Maryland, the influence of the Chesapeake Bay breeze on the PBLH was captured by both the ozonesonde and WRF data but generally not by the MPL PBLH retrievals because of generally weaker gradients in the aerosol backscatter profile and limited normalized relative backscatter data near the top height of the marine layer.


Atmosphere ◽  
2016 ◽  
Vol 26 (4) ◽  
pp. 673-686 ◽  
Author(s):  
Misun Kang ◽  
Yun-Kyu Lim ◽  
Changbum Cho ◽  
Kyu Rang Kim ◽  
Jun Sang Park ◽  
...  

2005 ◽  
Vol 20 (6) ◽  
pp. 1048-1060 ◽  
Author(s):  
Isidora Jankov ◽  
William A. Gallus ◽  
Moti Segal ◽  
Brent Shaw ◽  
Steven E. Koch

Abstract In recent years, a mixed-physics ensemble approach has been investigated as a method to better predict mesoscale convective system (MCS) rainfall. For both mixed-physics ensemble design and interpretation, knowledge of the general impact of various physical schemes and their interactions on warm season MCS rainfall forecasts would be useful. Adopting the newly emerging Weather Research and Forecasting (WRF) model for this purpose would further emphasize such benefits. To pursue this goal, a matrix of 18 WRF model configurations, created using different physical scheme combinations, was run with 12-km grid spacing for eight International H2O Project (IHOP) MCS cases. For each case, three different treatments of convection, three different microphysical schemes, and two different planetary boundary layer schemes were used. Sensitivity to physics changes was determined using the correspondence ratio and the squared correlation coefficient. The factor separation method was also used to quantify in detail the impacts of the variation of two different physical schemes and their interaction on the simulated rainfall. Skill score measures averaged over all eight cases for all 18 configurations indicated that no one configuration was obviously best at all times and thresholds. The greatest variability in forecasts was found to come from changes in the choice of convective scheme, although notable impacts also occurred from changes in the microphysics and planetary boundary layer (PBL) schemes. Specifically, changes in convective treatment notably impacted the forecast of system average rain rate, while forecasts of total domain rain volume were influenced by choices of microphysics and convective treatment. The impact of interactions (synergy) of different physical schemes, although occasionally of comparable magnitude to the impacts from changing one scheme alone (compared to a control run), varied greatly among cases and over time, and was typically not statistically significant.


2017 ◽  
Author(s):  
Ting Yang ◽  
Zifa Wang ◽  
Wei Zhang ◽  
Alex Gbaguidi ◽  
Nubuo Sugimoto ◽  
...  

Abstract. Predicting air pollution events in low atmosphere over megacities requires thorough understanding of the tropospheric dynamic and chemical processes, involving notably, continuous and accurate determination of the boundary layer height (BLH). Through intensive observations experimented over Beijing (China), and an exhaustive evaluation existing algorithms applied to the BLH determination, persistent critical limitations are noticed, in particular over polluted episodes. Basically, under weak thermal convection with high aerosol loading, none of the retrieval algorithms is able to fully capture the diurnal cycle of the BLH due to pollutant insufficient vertical mixing in the boundary layer associated with the impact of gravity waves on the tropospheric structure. Subsequently, a new approach based on gravity wave theory (the cubic root gradient method: CRGM), is developed to overcome such weakness and accurately reproduce the fluctuations of the BLH under various atmospheric pollution conditions. Comprehensive evaluation of CRGM highlights its high performance in determining BLH from Lidar. In comparison with the existing retrieval algorithms, the CRGM potentially reduces related computational uncertainties and errors from BLH determination (strong increase of correlation coefficient from 0.44 to 0.91 and significant decrease of the root mean square error from 643 m to 142 m). Such newly developed technique is undoubtedly expected to contribute to improve the accuracy of air quality modelling and forecasting systems.


Sign in / Sign up

Export Citation Format

Share Document