Weather Research And Forecasting Research Articles

Sensitivity of WRF-Simulated 2 m Temperature and Precipitation to Physics Options over the Loess Plateau

The current paper evaluates the weather research and forecasting (WRF) model sensitivity to five different combinations of cumulus, microphysics, radiation, and planetary boundary layer (PBL) schemes over Loess Plateau for the period 2015, in terms of 2 m temperature and precipitation. The WRF configuration consists of a 10 km resolution domain nested in a coarser domain driven by European Center for Medium-Range Weather Forecasts Reanalysis (ERA-Interim) data. The model simulated 2 m temperature and precipitation have been evaluated at daily and monthly scales with gridded observational dataset. The analysis shows that all experiments reproduce well the daily 2 m temperature, with overestimation particularly in the low-temperature range. Precipitation is less well simulated, with underestimation in all range, especially for intense rainfall. Comparing with ERA-Interim, WRF shows no clear benefit in simulating daily 2 m temperature while prominent improvement in simulating daily precipitation. WRF simulations capture the annual cycle of monthly 2 m temperature and precipitation with a warm bias and wet bias for most experiments in summer. Some reasonable configurations are identified. The “best” configuration depends on the criteria.

Improving the atmospheric dispersion forecasts over Washington, D.C. using UrbanNet observations: A study with HYSPLIT model

UrbanNet data have been used to enhance the predictions of the Weather Research and Forecasting (WRF) model through observational nudging to improve temperature and wind predictions. HYSPLIT was utilized to understand the impact of using the observationally nudged WRF fields in dispersion modeling. The meteorological observations collected from the National Weather Service monitoring stations located at two major airports in Washington metropolitan area were also assimilated into WRF modeling. The results showed that observational nudging successfully adjusted WRF wind fields towards the observations and significantly reduced the forecast temperature bias at nighttime. The comparison of HYSPLIT simulations with and without the enhancement of the WRF model using UrbanNet and airport data showed significant differences in the pattern and direction of the dispersion plume especially during the early morning hours. Furthermore, ingesting the data from the closer airport to the downtown area in WRF provided HYSPLIT simulations very similar to the ones using UrbanNet data. This provides strong evidence that local data are essential to adjust weather prediction models routinely used to drive dispersion models. There was also evidence of increased mixing height when using local data collected in the downtown, mainly resulting from the increased surface heating in the city.

Influence of Foehn Winds of Truong Son Mountains on the High Temperatures Observed in North-Central Vietnam during 31 May–5 June 2017

Abstract North-central Vietnam often experiences high temperatures. Foehn winds originating from the Truong Son Mountains (also known as Laos winds) are believed to contribute to abnormally high temperatures; however, no quantitative research has focused on foehn warming in Vietnam. In this study, we conducted numerical simulations using the Weather Research and Forecasting (WRF) Model to investigate the contribution of foehn warming to abnormally high temperatures in north-central Vietnam in early June 2017. Generally, May–June is the monsoon period in Vietnam. Consequently, foehn warming during this season is thought to be mainly caused by latent heating and precipitation mechanism. However, the primary factor in the cases covered in this study was foehn warming with an isentropic drawdown mechanism. Diabatic heating with turbulent diffusion and sensible heat flux from mountain slopes also plays significant roles. The warming effect of the foehn winds on the temperatures during the events was approximately 2°–3°C. It was concluded that the high temperature events from 31 May to 5 June 2017 were caused by synoptic-scale warm advection and foehn warming. Sensitivity experiments were conducted on the WRF Model, utilizing three atmospheric boundary layer turbulence schemes (YSU, ACM2, and MYNN), consistently yielding results for simulated temperature and relative humidity. The wind speed bias for the MYNN scheme was found to be lower than that of the other schemes. However, this study did not delve into the underlying reasons for these differences. The optimal performance of each scheme remains an open question. Significance Statement It was hypothesized that north-central Vietnam often experiences high temperatures owing to foehn winds (also known as Laos winds) descending from the Truong Son Mountains. This study conducted numerical experiments to validate this hypothesis and investigate the associated mechanisms of foehn winds in this region. Surprisingly, despite the monsoon season, the isentropic drawdown mechanism without precipitation effects provides the primary explanation for the high temperatures caused by foehn winds in north-central Vietnam, and not the latent heating and precipitation mechanism with precipitation effects that researchers expected. The results of this study contribute to better prediction of high temperatures in this region and improve our understanding of foehn winds in tropical monsoon climate zones.

Forecasting High Wind Events in the HRRR Model over Wyoming and Colorado. Part II: Sensitivity of Surface Wind Speeds to Model Resolution and Physics

Abstract Strong wind events can cause severe economic loss and societal impacts. Peak winds over Wyoming and Colorado occur during the wintertime months, and in the first part of this two-part series, it was shown that the High-Resolution Rapid Refresh (HRRR) model displays large negative biases with respect to strong wind events. In this part of the study, we address two questions: 1) does increasing the horizontal resolution improve the representation of strong wind events over this region, and 2) are the biases in HRRR-forecasted winds related to the selected planetary boundary layer (PBL), surface layer (SL), and/or land surface model (LSM) parameterizations? We conduct Weather Research and Forecasting (WRF) Model simulations to address these two main questions. Increasing the horizontal resolution leads to a small improvement in the simulation of strong wind speeds over the complex terrain of Wyoming and Colorado. In general, changes in the PBL, SL, and LSM parameterizations show much larger changes in simulated wind speeds compared to increasing the model resolution alone. Specifically, changing from the Mellor–Yamada–Nakanishi–Niino scheme to the Mellor–Yamada–Janjić PBL and SL schemes results in nearly no change in the r2 values, but there is a decrease in the magnitude of the strong wind speed bias (from −12.52 to −10.16 m s−1). We attribute these differences to differences in the diagnosis of surface wind speeds and mixing in the boundary layer. Further analysis is conducted to determine the value of 1-km forecasts of strong winds compared with wind speed diagnostics commonly used by the National Weather Service. Significance Statement Motivated by prior studies showing low skill in the prediction of strong winds using state-of-the-art weather forecast models, in this study, we aim to investigate two questions: 1) does increasing the horizontal resolution improve the prediction of strong wind events over the complex terrain of Wyoming and Colorado, and 2) are the biases in High-Resolution Rapid Refresh (HRRR) forecasted winds related to the selected planetary boundary layer, surface layer, and/or land surface model parameterizations? We find that increasing the horizontal resolution provides a slight improvement in the prediction of strong winds. Further, considerable improvement in the prediction of strong winds is found for varying boundary layer, surface layer, and land surface parameterizations.

Modeling of air pollution due to marine traffic in Marseille

Marine traffic is constantly increasing and is unlikely to decrease in the years to come. However, merchant and tourist ships still use fuels with a high pollutant content, particularly in the Mediterranean sea where regulations are more flexible than those applied in northern European waters. The pollutants of interest in this work, which include particulate matter (PM10, PM2.5), NO2, SO2 and ozone, are responsible for respiratory diseases such as airway irritation, cardiopulmonary and lung cancers. The aim of this study is to estimate the impact of maritime traffic on air quality in Marseille in case of high port occupancy. The same ship traffic emission scenario will be studied under a series of typical meteorological conditions selected on the basis of a statistical seasonal analysis. This approach therefore differs from studies based on annualized ship traffic scenarios, which do not capture worst-case pollution. In order to do so, we have parameterized the numerical weather forecasting model WRF (Weather Research and Forecasting) that allows the inclusion of chemistry and atmospheric transport equations. The relative difference in concentration between a situation with and without ship traffic allows us to evaluate the influence of ships on air quality. The results show an average ship contribution of between 1.5% and 5.7% for PM2.5, with some local maximum of about 25% near the port. Contribution levels for gaseous species are much higher, with values reaching 40% for SO2 and 80% for NO2. The analysis of various meteorological conditions has highlighted a significant influence of wind speed and sea–land breeze dynamic on ship contribution.

Evaluation of a High Resolution WRF Model for Southeast Brazilian Coast: The Importance of Physical Parameterization to Wind Representation

This study assesses the performance of the Weather Research and Forecasting (WRF) model using a high-resolution spatial grid (1 km) with various combinations of physical parameterization packages to simulate a severe event in August 2021 in the southeastern Brazilian coast. After determining the optimal set of physical parameterizations for representing wind patterns during this event, a year-long evaluation was conducted, covering forecast horizons of 24, 48, and 72 h. The simulation results were compared with observational wind data from four weather stations. The findings highlight variations in the efficacy of different physical parameterization sets, with certain sets encountering challenges in accurately depicting the peak of the severe event. The most favorable results were achieved using a combination of Tiedtke (cumulus), Thompson (microphysics), TKE (boundary layer), Monin-Obukhov (surface layer), Unified-NOAH (land surface), and RRTMG (shortwave and longwave radiation). Over the one-year forecasting period, the WRF model effectively represented the overall wind pattern, including forecasts up to three days in advance (72-h forecast horizon). Generally, the statistical metrics indicate robust model performance, even for the 72-h forecast horizon, with correlation coefficients consistently exceeding 0.60 at all analyzed points. While the model proficiently captured wind distribution, it tended to overestimate northeast wind speed and gust intensities. Notably, forecast accuracy decreased as stations approached the ocean, exemplified by the ATPM station.

The 2023 National Offshore Wind data set (NOW-23)

Abstract. This article introduces the 2023 National Offshore Wind data set (NOW-23), which offers the latest wind resource information for offshore regions in the United States. NOW-23 supersedes, for its offshore component, the Wind Integration National Dataset (WIND) Toolkit, which was published a decade ago and is currently a primary resource for wind resource assessments and grid integration studies in the contiguous United States. By incorporating advancements in the Weather Research and Forecasting (WRF) model, NOW-23 delivers an updated and cutting-edge product to stakeholders. In this article, we present the new data set which underwent regional tuning and performance validation against available observations and has data available from 2000 through, depending on the region, 2019–2022. We also provide a summary of the uncertainty quantification in NOW-23, along with NOW-WAKES, a 1-year post-construction data set that quantifies expected offshore wake effects in the US Mid-Atlantic lease areas. Stakeholders can access the NOW-23 data set at https://doi.org/10.25984/1821404 (Bodini et al., 2020).

Evaluation and development of surface layer scheme representation of temperature inversions over boreal forests in Arctic wintertime conditions

Abstract. In this study, the Noah land surface model used in conjunction with the Mellor–Yamada–Janjić surface layer scheme (hereafter, Noah-MYJ) and the Noah multiphysics scheme (Noah-MP) from the Weather Research and Forecasting (WRF) 4.5.1 mesoscale model are evaluated with regard to their performance in reproducing positive temperature gradients over forested areas in the Arctic winter. First, simplified versions of the WRF schemes, recoded in Python, are compared with conceptual models of the surface layer in order to gain insight into the dependence of the temperature gradient on the wind speed at the top of the surface layer. It is shown that the WRF schemes place strong limits on the turbulent collapse, leading to lower surface temperature gradient at low wind speeds than in the conceptual models. We implemented modifications to the WRF schemes to correct this effect. The original and modified versions of Noah-MYJ and Noah-MP are then evaluated compared to long-term measurements at the Ameriflux Poker Flat Research Range, a forest site in interior Alaska. Noah-MP is found to perform better than Noah-MYJ because the former is a two-layer model which explicitly takes into account the effect of the forest canopy. Indeed, a non-negligible temperature gradient is maintained below the canopy at high wind speeds, leading to overall larger gradients than in the absence of vegetation. Furthermore, the modified versions are found to perform better than the original versions of each scheme because they better reproduce strong temperature gradients at low wind speeds.

Evaluation of the fully coupled WRF and WRF-Hydro modelling system initiated with satellite-based soil moisture data

ABSTRACT This study comparatively evaluates the skill of the (un)coupled Weather Research and Forecasting (WRF) – WRF-Hydro system supported with satellite-based soil moisture initiation in long-term simulations of hydrometeorological variables under humid and semi-arid climate conditions. Prior to the coupling, the WRF-Hydro model is calibrated using streamflow data at different data scales (i.e. hourly and daily) to discuss the influences of the changing temporal resolution of the reference dataset on the optimum parameter values. According to the results, fully coupled models correct the location and magnitude of the occurrence of hourly extreme precipitation. Changing the initial soil moisture data source affects the terrestrial water cycle elements in the semi-arid region (i.e. Mediterranean (MED)) more than in the humid area (i.e. Eastern Black Sea (EBS)), especially in spring and summer. The initiation of the coupled model with satellite-based soil moisture improves the streamflow outputs after a reasonable spin-up period for the MED region.

TYPHOON KARDING (NORU) STORM SURGE ANALYSIS USING THE COAWST MODELING SYSTEM

Abstract. The Philippines, frequently affected by typhoons, faces the hazard of storm surges. This study examined the Coupled-Ocean-Atmosphere-Wave-Sediment Transport Modeling System (COAWST) to simulate Typhoon Karding's September 2022 storm surge. COAWST integrates the Regional Ocean Modeling System (ROMS), Simulating WAves Nearshore (SWAN), and Weather Research and Forecasting (WRF) coupled using the Model Coupling Toolkit (MCT). Four setups were analyzed: i. ROMS only, ii. ROMS-SWAN, iii. ROMS-WRF, and iv. ROMS-SWAN-WRF, focusing on four variables: a. Surface air pressure, b. Wind speed, c. Free-surface elevation, and d. Significant wave height. Results show that the ROMS-WRF and ROMS-SWAN-WRF setups accurately simulated Typhoon Karding's track with minimal positional error and wind speed. However, the models overestimated the typhoon's minimum air pressure with p-biases of 7.74% (i and ii), 4.1% (iii), and 3.9% (iv), and RMSE values of 68.529 hPa (i and ii), 36.744 hPa (iii), and 36.789 hPa (iv). Additionally, water levels were underestimated, with RMSE ranging from 0.31 to 0.35 meters and p-biases from −72.56% to −154.89% at Baler, Aurora validation point. At the Real, Quezon validation point, RMSE and p-bias ranged from 0.30 to 0.34 meters and −84.80% to −166.44%, respectively. Nonetheless, the models were able to simulate the storm surge and significant wave height at Baler and Real points similar to recorded data, with setup iv performing best in storm surge simulation. In summary, COAWST may be employed for typhoon simulations, with coupling being able to increase accuracy.

Investigation on the Sensitivity of Precipitation Simulation to Model Parameterization and Analysis Nudging over Hebei Province, China

The physical parameterizations have important influence on model performance in precipitation simulation and prediction; however, previous investigations are seldom conducted at very high resolution over Hebei Province, which is often influenced by extreme events such as droughts and floods. In this paper, the influence of parameterization schemes and analysis nudging on precipitation simulation is investigated using the WRF (weather research and forecasting) model with many sensitivity experiments at the cumulus “gray-zone” resolution (5 km). The model performance of different sensitivity simulations is determined by a comparison with the local high-quality observational data. The results indicate that the WRF model generally reproduces the distribution of precipitation well, and the model tends to underestimate precipitation compared with the station observations. The sensitivity simulation with the Tiedtke cumulus parameterization scheme combined with the Thompson microphysics scheme shows the best model performance, with the highest temporal correlation coefficient (0.45) and lowest root mean square error (0.34 mm/day). At the same time, analysis nudging, which incorporates observational information into simulation, can improve the model performance in precipitation simulation. Further analysis indicates that the negative bias in precipitation may be associated with the negative bias in relative humidity, which in turn is associated with the positive bias in temperature and wind speed. This study highlights the role of parameterization schemes and analysis nudging in precipitation simulation and provides a valuable reference for further investigations on precipitation forecasting applications.

Evaluation of wind farm parameterizations in the WRF model under different atmospheric stability conditions with high-resolution wake simulations

Abstract. Wind farm parameterizations (WFPs) are used in mesoscale models for predicting wind farm power production and its impact on wind resources while considering the variability of the regional wind climate. However, the performance of WFPs is influenced by various factors including atmospheric stability. In this study, we compared two widely used WFPs in the Weather Research and Forecasting (WRF) model to large-eddy simulations (LES) of turbine wakes performed with the same model. The Fitch WFP and the explicit wake parameterization were evaluated for their ability to represent wind speed and turbulent kinetic energy (TKE) in a two-turbine wind farm layout under neutral, unstable, and stable atmospheric stability conditions. To ensure a fair comparison, the inflow conditions were kept as close as possible between the LES and mesoscale simulations for each type of stability condition, and the LES results were spatially aggregated to align with the mesoscale grid spacing. Our findings indicate that the performance of WFPs varies depending on the specific variable (wind speed or TKE) and the area of interest downwind of the turbine when compared to the LES reference. The WFPs can accurately depict the vertical profiles of the wind speed deficit for either the grid cell containing the wind turbines or the grid cells in the far wake, but not both simultaneously. The WFPs with an explicit source of TKE overestimate TKE values at the first grid cell containing the wind turbine; however, for downwind grid cells, agreement improves. On the other hand, WFPs without a TKE source underestimate TKE in all downwind grid cells. These agreement patterns between the WFPs and the LES reference are consistent under the three atmospheric stability conditions. However, the WFPs resemble less the wind speed and TKE from the LES reference under stable conditions than that under neutral or unstable conditions.

Challenges of constructing and selecting the “perfect” boundary conditions for the large-eddy simulation model PALM

Abstract. We present the process of and difficulties in acquiring the proper boundary conditions (BCs) for the state-of-the-art large-eddy simulation (LES)-based PALM model system. We use the mesoscale Weather Research and Forecasting (WRF) model as a source of inputs for the PALM preprocessor and investigate the influence of the mesoscale model on the performance of the PALM model. A total of 16 different WRF configurations were used as a proxy for a multi-model ensemble. We developed a technique for selecting suitable sets of BCs, performed PALM model simulations driven by these BCs, and investigated the consequences of selecting a sub-optimal WRF configuration. The procedure was tested for four episodes in different seasons of the year 2019, during which WRF and PALM outputs were evaluated against the atmospheric radiosounding observations. We show that the PALM model outputs are heavily dependent on the imposed BCs and have different responses at different times of the day and in different seasons. We demonstrate that the main driver of errors is the mesoscale model and that the PALM model is capable of attenuating but not fully correcting them. The PALM model attenuates the impact of errors in BCs in wind speed, while for the air temperature, PALM shows variable behavior with respect to driving conditions. This study stresses the importance of high-quality driving BCs and the complexity of the process of their construction and selection.

Modelling wind farm effects in HARMONIE–AROME (cycle 43.2.2) – Part 1: Implementation and evaluation

Abstract. With increasing number and proximity of wind farms, it becomes crucial to consider wind farm effects (WFEs) in the numerical weather prediction (NWP) models used to forecast power production. Furthermore, these WFEs are also expected to affect other weather-related parameters at least locally. Thus, we implement the explicit wake parameterization (EWP) in the NWP model HARMONIE–AROME (hereafter HARMONIE) along-side the existing wind farm parameterization (WFP) by Fitch et al. (2012) (FITCH). We evaluate and compare the two WFPs against research flight measurements as well as against similar simulations performed with the Weather Research and Forecasting (WRF) model using case studies. The case studies include a case for WFEs above a wind farm as well as two cases for WFEs at hub height in the wake of farms. The results show that EWP and FITCH have been correctly implemented in HARMONIE. For the simulated cases, EWP underestimates the WFEs on wind speed and strongly underestimates the effect on turbulent kinetic energy (TKE). FITCH agrees better with the observations, and WFEs on TKE are particularly well captured by HARMONIE–FITCH. After this successful evaluation, simulations with all wind turbines in Europe will be performed with HARMONIE and presented in the second part of this paper series.

Impacts of projected urban growth on simulated near-surface temperature in Mexico City Metropolitan Area: Implications for urban vulnerability

Urbanization impacts the surface temperature fields increasing the vulnerability of urban residents to heat exposure. Identifying vulnerable urban populations to extreme heat exposure is crucial to develop mitigation and adaptation strategies towards sustainability. We used an urban growth model (SLEUTH) to simulate emerging urban areas in Mexico City Metropolitan Area (MCMA) under a hypothetical land-use policy scenario projected to 2060 in which no restrictions were posed to urban growth. SLEUTH outputs were used in the numerical model Weather Research and Forecasting (WRF) to quantify expected changes in near-surface temperature within the MCMA. We calculated and mapped heat exposure as differences in average (Tmean), maximum (Tmax) and minimum (Tmin) temperatures over the diurnal cycle between future and current land cover conditions. Population vulnerability to projected increases in heat exposure was determined using a set of socioeconomic indicators. SLEUTH simulations showed an urban area expansion of nearly 4,790 km2 by 2060. Overall, changes in Tmin were greater than changes observed for Tmax and Tmean. Tmean, Tmax and Tmin increases up to 0.6°C, 1.3°C and 2.6°C, respectively, were recorded for the MCMA with greatest temperature changes observed in the State of Mexico. Results suggested the presence of socioeconomic disparities in the projected spatial exposure of urban-induced heat in MCMA. We argue that our results could be used to inform and guide locally tailored actions aimed at reducing exposure and increasing population´s capacities to cope and adapt to future threats.

Analysis of an intense O3 pollution episode on the Atlantic coast of the Iberian Peninsula using photochemical modeling: characterization of transport pathways and accumulation processes

Abstract. A tropospheric O3 pollution episode over the Atlantic coast of the Iberian Peninsula during 2–6 August 2018 has been analyzed. The episode was characterized by a permanent wind shear throughout the entire period, making the observed ozone surface distribution especially difficult to explain. A new methodology is described to analyze upper-level atmospheric parameters, such as temperature, wind direction, wind speed, and O3 concentrations, added to the traditional use of surface parameters, using the Weather Research and Forecasting (WRF) and Comprehensive Air Quality Model with Extensions (CAMx) models as well as available surface and upper-air observations. Results indicate that the episode was characterized by a first phase of a sudden increase in O3 concentrations produced by fumigation and inter-regional transport processes within the Iberian Peninsula, followed by continental O3 transport from Europe to the Atlantic coast. An Atlantic front produced the dissipation of the episode, generating an “ozone front” accompanying the cold front passage across the region.

Optimization and Evaluation of the Weather Research and Forecasting (WRF) Model for Wind Energy Resource Assessment and Mapping in Iran

This study aims to optimize the Weather Research and Forecasting (WRF) model regarding the choice of the best planetary boundary layer (PBL) physical scheme and to evaluate the model’s performance for wind energy assessment and mapping over the Iranian territory. In this initiative, five PBL and surface layer parameterization schemes were tested, and their performance was evaluated via comparison with observational wind data. The study used two-way nesting domains with spatial resolutions of 15 km and 5 km to represent atmospheric circulation patterns affecting the study area. Additionally, a seventeen-year simulation (2004–2020) was conducted, producing wind datasets for the entire Iranian territory. The accuracy of the WRF model was assessed by comparing its results with observations from multiple sites and with the high-resolution Global Wind Atlas. Statistical parameters and wind power density were calculated from the simulated data and compared with observations to evaluate wind energy potential at specific sites. The model’s performance was sensitive to the horizontal resolution of the terrain data, with weaker simulations for wind speeds below 3 m/s and above 10 m/s. The results confirm that the WRF model provides reliable wind speed data for realistic wind energy assessment studies in Iran. The model-generated wind resource map identifies areas with high wind (wind speed > 5.6 m/s) potential that are currently without wind farms or Aeolic parks for exploitation of the wind energy potential. The Sistan Basin in eastern Iran was identified as the area with the highest wind power density, while areas west of the Zagros Mountains and in southwest Iran showed high aeolian potential during summer. A novelty of this research is the application of the WRF model in an area characterized by high topographical complexities and specific geographical features. The results provide practical solutions and valuable insights for industry stakeholders, facilitating informed decision making, reducing uncertainties, and promoting the effective utilization of wind energy resources in the region.

Outlook for Offshore Wind Energy Development in Mexico from WRF Simulations and CMIP6 Projections

This study presents a comprehensive assessment of the offshore wind energy potential in Mexico across 40 years (1979–2018) of numerical simulations using the Weather Research and Forecasting (WRF) model and data from the Coupled Model Intercomparison Project Phase 6 (CMIP6). The WRF identifies three regions with moderate to good wind potential: off the north coast of Tamaulipas (Zone I), the northwest coast of Yucatan (Zone II), and the Gulf of Tehuantepec (Zone III). The analysis involves comparing 47 CMIP6 climate models with the WRF results and selecting the best performing models to obtain future projections for the short term (2040–2069) and the long term (2070–2099). Two ensemble-based strategies were implemented. The first one, which uses an intersection approach from which four CMIP6 models were considered, reveals positive percentage differences in Zone II for both future projections, especially for the long-term one. In Zones I and III, positive values are also observed near the coast, mainly for the long-term projection, but they are considerably lower compared to those in Zone II. The second ensemble strategy uses weight assignment through the Mean Absolute Percentage Error, so that a greater weight is given to the model that performed better in each particular zone, potentially providing more accurate results. The findings suggest the likelihood of increased offshore wind energy in these three zones of Mexico, for both short- and long-term future projections, with positive percentage differences of up to 10% in certain areas.

The Role of Parameterizations and Model Coupling on Simulations of Energy and Water Balances – Investigations With the Atmospheric Model WRF and the Hydrologic Model WRF‐Hydro

AbstractThe distributed hydrologic model WRF‐Hydro can operate in a fully‐coupled mode with the atmospheric Weather, Research and Forecasting (WRF) model. WRF‐Hydro enhances the modeling of terrestrial hydrologic processes in coupled WRF/WRF‐Hydro by simulating lateral surface and subsurface water flows. The objectives of this study are (a) to examine the effect of WRF‐Hydro on the surface energy and water balance in fully‐coupled WRF/WRF‐Hydro simulations and (b) to examine the impact of five WRF physics parameterizations on WRF‐Hydro streamflow. The study area is the Mediterranean island of Cyprus and 31 mountainous watersheds. The domain‐average soil moisture was 20% higher in the coupled WRF/WRF‐Hydro, relative to the standalone WRF model, during a 1‐year simulation. The higher soil moisture could explain the increase in latent heat (36%) and evapotranspiration (33%). The increase in these fluxes was less with a modification in the model transpiration parameterization to represent nocturnal transpiration and the use of remote sensing leaf area index data. The simulated precipitation of the coupled model increased up to 3%, relative to WRF. Two‐year long WRF‐Hydro simulations gave a median Nash‐Sutcliffe Efficiency for daily streamflow of the 31 watersheds of 0.5 for observed precipitation forcing and between −1.9 and 0.2 for the forcing of the five WRF parameterizations. This study showed that the enhancement of the standalone WRF model with lateral water flow processes in the coupled mode with WRF‐Hydro modifies the terrestrial energy and water balance. The improved terrestrial process representation should be considered for future hydrological cycle studies with WRF.

Implications of WRF model resolutions on resolving rainfall variability with topography over East Africa

There is an increasing need to improve the accuracy of extreme weather forecasts for life-saving applications and in support of various socioeconomic sectors in East Africa, a region with remarkable mesoscale systems due to its complex topography defined by sharp gradients in elevation, inland water bodies, and landuse conversions. This study sought to investigate the impacts of the Weather Research and Forecasting (WRF) model spatial resolution on resolving rainfall variability with topography utilizing nested domains at 12 and 2.4 km resolutions. The model was driven by the National Centers for Environmental Prediction (NCEP)-Global Data Assimilation System (GDAS) Global Forecast System (GFS) final (FNL) reanalysis to simulate the weather patterns over East Africa from 3rd April 2018 to 30th April 2018, which were evaluated against several freely available gridded weather datasets alongside rainfall data from the Kenya Meteorological Department (KMD) stations. The reference datasets and the model outputs revealed that the highlands had more rainfall events and higher maximum daily rainfall intensity compared to the surrounding lowlands, attributed to orographic lifting enhancing convection. Rainfall was inversely proportional to altitude from 500 m to 1,100 m above sea level (ASL) for both coarse and fine resolutions. The convection-permitting setup was superior in three aspects: resolving the inverse altitude-rainfall relationship for altitudes beyond 3000 m ASL, simulating heavy rainfall events over the lowlands, and resolution of the diurnal cycle of low-level wind. Although the coarse resolution setup reasonably simulated rainfall over large mountains, only the convection-permitting configuration could accurately resolve rainfall variability over contrasting topographical features. The study notes that high-resolution modeling systems and topography-sensitive bias correction techniques are critical for improving the quality of operational weather forecasts in East Africa.