scholarly journals Lessons learned after two years of operational high-resolution (1,5 km) WRF simulations in Catalonia and a plan to increase their skill

2021 ◽  
Author(s):  
Jordi Mercader Carbó ◽  
Manel Bravo Blanco ◽  
Jordi Moré Pratdesaba ◽  
Abdelmalik Sairouni Afif
Keyword(s):  
Wind Speed ◽  
Daily Precipitation ◽  
Mesoscale Model ◽  
Initial Conditions ◽  
Skill Score ◽  
Lessons Learned ◽  
Grid Spacing ◽  
Slight Improvement ◽  
Model Topography ◽  
Score Table

<p>The WRF-ARW has been the flagship mesoscale model in the Meteorological Service of Catalonia (SMC) since 2012. Several operational runs are performed daily (initialised at 00 and 12 UTC), using both the GFS and the IFS model for initial and boundary conditions, to account for uncertainties in the synoptic evolution. To provide more accurate forecasts to end-users, a convection-allowing simulation with a grid spacing of 1,5 km was added to the operational chain, starting in the summer of 2019. </p><p>However, the verification results show that the improvement over its mother domain (a 3 km simulation with parameterised convection) is irregular because it does not happen for all the variables. For instance, the 2 m temperature forecasts are more reliable for the highest resolution domain but the wind speed at 10 m has a comparable skill. Regarding the precipitation, there is a very slight improvement only for high daily precipitation rates (50 or 80 mm) during some seasons; nevertheless, the results are worse in forecasting the occurrence of precipitation (that is, when considering low daily precipitation quantities). The comparison of the verification results among different model configurations (with various resolutions and initial conditions) can be easily performed by using a skill score table. This table and its design will also be presented in this session. </p><p>Certainly, these results help to conceive strategies to enhance the skill of the 1,5 km simulations for some of the variables that arise as more inaccurate. For instance, it is evaluated to what extent using alternative static fields (changing the model topography or the land category) improves the forecasts of temperature, humidity or wind near the surface. Furthermore, the sensitivity of precipitation forecasts to several physics schemes is tested, seeking an enhancement of their skill. </p>

scholarly journals Dynamical Downscaling of Wind Speed in Complex Terrain Prone To Bora-Type Flows

2011 ◽  
Vol 50 (8) ◽  
pp. 1676-1691 ◽  
Author(s):  
Kristian Horvath ◽  
Alica Bajić ◽  
Stjepan Ivatek-Šahdan
Keyword(s):  
Wind Speed ◽  
Complex Terrain ◽  
Dynamical Downscaling ◽  
Mesoscale Model ◽  
Spectral Power ◽  
Grid Spacing ◽  
Model Data ◽  
Mean Wind Speed ◽  
The Mean ◽  
Horizontal Grid

AbstractThe results of numerically modeled wind speed climate, a primary component of wind energy resource assessment in the complex terrain of Croatia, are given. For that purpose, dynamical downscaling of 10 yr (1992–2001) of the 40-yr ECMWF Re-Analysis (ERA-40) was performed to 8-km horizontal grid spacing with the use of a spectral, prognostic full-physics model Aire Limitée Adaptation Dynamique Développement International (ALADIN; the “ALHR” version). Then modeled data with a 60-min frequency were refined to 2-km horizontal grid spacing with a simplified and cost-effective model version, the so-called dynamical adaptation (DADA). The statistical verification of ERA-40-, ALHR-, and DADA-modeled wind speed on the basis of data from measurement stations representing different regions of Croatia suggests that downscaling was successful and that model accuracy generally improves as horizontal resolution is increased. The areas of the highest mean wind speeds correspond well to locations of frequent and strong bora flow as well as to the prominent mountain peaks. The best results are achieved with DADA and contain bias of 1% of the mean wind speed for eastern Croatia while reaching 10% for complex coastal terrain, mainly because of underestimation of the strongest winds. Root-mean-square errors for DADA are significantly smaller for flat terrain than for complex terrain, with relative values close to 12% of the mean wind speed regardless of the station location. Spectral analyses suggest that the shape of the kinetic energy spectra generally relaxes from k−3 at the upper troposphere to the shape of orographic spectra near the surface and shows no seasonal variability. Apart from the buildup of energy on smaller scales of motions, it is shown that mesoscale simulations contain a considerable amount of energy related to near-surface and mostly divergent meso-β-scale (20–200 km) motions. Spectral decomposition of measured and modeled data in temporal space indicates a reasonable performance of all model datasets in simulating the primary maximum of spectral power related to synoptic and larger-than-diurnal mesoscale motions, with somewhat increased accuracy of mesoscale model data. The primary improvement of dynamical adaptation was achieved for cross-mountain winds, whereas mixed results were found for along-mountain wind directions. Secondary diurnal and tertiary semidiurnal maxima are significantly better simulated with the mesoscale model for coastal stations but are somewhat more erroneous for the continental station. The mesoscale model data underestimate the spectral power of motions with less-than-semidiurnal periods.

scholarly journals Use of Multiple Verification Methods to Evaluate Forecasts of Convection from Hot- and Cold-Start Convection-Allowing Models

2013 ◽  
Vol 28 (1) ◽  
pp. 119-138 ◽  
Author(s):  
Derek R. Stratman ◽  
Michael C. Coniglio ◽  
Steven E. Koch ◽  
Ming Xue
Keyword(s):  
Mesoscale Model ◽  
Initial Conditions ◽  
Spatial Scales ◽  
Cold Start ◽  
Skill Score ◽  
Test Bed ◽  
Hot Start ◽  
Verification Methods ◽  
Skill Scores ◽  
Convective Cells

Abstract This study uses both traditional and newer verification methods to evaluate two 4-km grid-spacing Weather Research and Forecasting Model (WRF) forecasts: a “cold start” forecast that uses the 12-km North American Mesoscale Model (NAM) analysis and forecast cycle to derive the initial and boundary conditions (C0) and a “hot start” forecast that adds radar data into the initial conditions using a three-dimensional variational data assimilation (3DVAR)/cloud analysis technique (CN). These forecasts were evaluated as part of 2009 and 2010 NOAA Hazardous Weather Test Bed (HWT) Spring Forecasting Experiments. The Spring Forecasting Experiment participants noted that the skill of CN’s explicit forecasts of convection estimated by some traditional objective metrics often seemed large compared to the subjectively determined skill. The Gilbert skill score (GSS) reveals CN scores higher than C0 at lower thresholds likely due to CN having higher-frequency biases than C0, but the difference is negligible at higher thresholds, where CN’s and C0’s frequency biases are similar. This suggests that if traditional skill scores are used to quantify convective forecasts, then higher (>35 dBZ) reflectivity thresholds should be used to be consistent with expert’s subjective assessments of the lack of forecast skill for individual convective cells. The spatial verification methods show that both CN and C0 generally have little to no skill at scales <8–12Δx starting at forecast hour 1, but CN has more skill at larger spatial scales (40–320 km) than C0 for the majority of the forecasting period. This indicates that the hot start provides little to no benefit for forecasts of convective cells, but that it has some benefit for larger mesoscale precipitation systems.

Two mechanisms for accelerated diffusion of COVID-19 outbreaks in regions with high intensity of population and polluting industrialization: the air pollution-to-human and human-to-human transmission dynamics (Preprint)

2020 ◽  
Author(s):  
Mario Coccia
Keyword(s):  
Air Pollution ◽  
Wind Speed ◽  
Critical Role ◽  
Meteorological Condition ◽  
Meteorological Conditions ◽  
Lessons Learned ◽  
Transmission Dynamics ◽  
Critical Threshold ◽  
Viral Infectivity ◽  
Low Wind Speed

BACKGROUND Coronavirus disease 2019 (COVID-19) is viral infection that generates a severe acute respiratory syndrome with serious pneumonia that may result in progressive respiratory failure and death. OBJECTIVE This study has two goals. The first is to explain the main factors determining the diffusion of COVID-19 that is generating a high level of deaths. The second is to suggest a strategy to cope with future epidemic threats with of accelerated viral infectivity in society. METHODS Correlation and regression analyses on on data of N=55 Italian province capitals, and data of infected individuals at as of April 2020. RESULTS The main results are: o The accelerate and vast diffusion of COVID-19 in North Italy has a high association with air pollution. o Hinterland cities have average days of exceeding the limits set for PM10 (particulate matter 10 micrometers or less in diameter) equal to 80 days, and an average number of infected more than 2,000 individuals as of April 1st, 2020, coastal cities have days of exceeding the limits set for PM10 equal to 60 days and have about 700 infected in average. o Cities that average number of 125 days exceeding the limits set for PM10, last year, they have an average number of infected individual higher than 3,200 units, whereas cities having less than 100 days (average number of 48 days) exceeding the limits set for PM10, they have an average number of about 900 infected individuals. o The results reveal that accelerated transmission dynamics of COVID-19 in specific environments is due to two mechanisms given by: air pollution-to-human transmission and human-to-human transmission; in particular, the mechanisms of air pollution-to-human transmission play a critical role rather than human-to-human transmission. o The finding here suggests that to minimize future epidemic similar to COVID-19, the max number of days per year in which cities can exceed the limits set for PM10 or for ozone, considering their meteorological condition, is less than 50 days. After this critical threshold, the analytical output here suggests that environmental inconsistencies because of the combination between air pollution and meteorological conditions (with high moisture%, low wind speed and fog) trigger a take-off of viral infectivity (accelerated epidemic diffusion) with damages for health of population, economy and society. CONCLUSIONS Considering the complex interaction between air pollution, meteorological conditions and biological characteristics of viral infectivity, lessons learned for COVID-19 have to be applied for a proactive socioeconomic strategy to cope with future epidemics, especially an environmental policy based on reduction of air pollution mainly in hinterland zones of countries, having low wind speed, high percentage of moisture and fog that create an environment that can damage immune system of people and foster a fast transmission of viral infectivity similar to the COVID-19. CLINICALTRIAL not applicable

scholarly journals Using Fractal Downscaling of Satellite Precipitation Products for Hydrometeorological Applications

2010 ◽  
Vol 27 (3) ◽  
pp. 409-427 ◽  
Author(s):  
Kun Tao ◽  
Ana P. Barros
Keyword(s):  
Spatial Resolution ◽  
Stage Iv ◽  
Tropical Rainfall Measuring Mission ◽  
Skill Score ◽  
Probability Of Detection ◽  
Grid Spacing ◽  
Satellite Precipitation ◽  
Target Field ◽  
Skill Scores ◽  
Central United States

Abstract The objective of spatial downscaling strategies is to increase the information content of coarse datasets at smaller scales. In the case of quantitative precipitation estimation (QPE) for hydrological applications, the goal is to close the scale gap between the spatial resolution of coarse datasets (e.g., gridded satellite precipitation products at resolution L × L) and the high resolution (l × l; L ≫ l) necessary to capture the spatial features that determine spatial variability of water flows and water stores in the landscape. In essence, the downscaling process consists of weaving subgrid-scale heterogeneity over a desired range of wavelengths in the original field. The defining question is, which properties, statistical and otherwise, of the target field (the known observable at the desired spatial resolution) should be matched, with the caveat that downscaling methods be as a general as possible and therefore ideally without case-specific constraints and/or calibration requirements? Here, the attention is focused on two simple fractal downscaling methods using iterated functions systems (IFS) and fractal Brownian surfaces (FBS) that meet this requirement. The two methods were applied to disaggregate spatially 27 summertime convective storms in the central United States during 2007 at three consecutive times (1800, 2100, and 0000 UTC, thus 81 fields overall) from the Tropical Rainfall Measuring Mission (TRMM) version 6 (V6) 3B42 precipitation product (∼25-km grid spacing) to the same resolution as the NCEP stage IV products (∼4-km grid spacing). Results from bilinear interpolation are used as the control. A fundamental distinction between IFS and FBS is that the latter implies a distribution of downscaled fields and thus an ensemble solution, whereas the former provides a single solution. The downscaling effectiveness is assessed using fractal measures (the spectral exponent β, fractal dimension D, Hurst coefficient H, and roughness amplitude R) and traditional operational scores statistics scores [false alarm rate (FR), probability of detection (PD), threat score (TS), and Heidke skill score (HSS)], as well as bias and the root-mean-square error (RMSE). The results show that both IFS and FBS fractal interpolation perform well with regard to operational skill scores, and they meet the additional requirement of generating structurally consistent fields. Furthermore, confidence intervals can be directly generated from the FBS ensemble. The results were used to diagnose errors relevant for hydrometeorological applications, in particular a spatial displacement with characteristic length of at least 50 km (2500 km2) in the location of peak rainfall intensities for the cases studied.

scholarly journals Impact of physical parameterizations and initial conditions on simulated atmospheric transport and CO<sub>2</sub> mole fractions in the US Midwest

2018 ◽  
Vol 18 (20) ◽  
pp. 14813-14835 ◽  
Author(s):  
Liza I. Díaz-Isaac ◽  
Thomas Lauvaux ◽  
Kenneth J. Davis
Keyword(s):  
Boundary Layer ◽  
Wind Speed ◽  
Boundary Conditions ◽  
Planetary Boundary Layer ◽  
Initial Conditions ◽  
Atmospheric Transport ◽  
Atmospheric Co2 ◽  
The Us ◽  
Physics Parameterizations ◽  
Physical Parameterizations

Abstract. Atmospheric transport model errors are one of the main contributors to the uncertainty affecting CO2 inverse flux estimates. In this study, we determine the leading causes of transport errors over the US upper Midwest with a large set of simulations generated with the Weather Research and Forecasting (WRF) mesoscale model. The various WRF simulations are performed using different meteorological driver datasets and physical parameterizations including planetary boundary layer (PBL) schemes, land surface models (LSMs), cumulus parameterizations and microphysics parameterizations. All the different model configurations were coupled to CO2 fluxes and lateral boundary conditions from the CarbonTracker inversion system to simulate atmospheric CO2 mole fractions. PBL height, wind speed, wind direction, and atmospheric CO2 mole fractions are compared to observations during a month in the summer of 2008, and statistical analyses were performed to evaluate the impact of both physics parameterizations and meteorological datasets on these variables. All of the physical parameterizations and the meteorological initial and boundary conditions contribute 3 to 4 ppm to the model-to-model variability in daytime PBL CO2 except for the microphysics parameterization which has a smaller contribution. PBL height varies across ensemble members by 300 to 400 m, and this variability is controlled by the same physics parameterizations. Daily PBL CO2 mole fraction errors are correlated with errors in the PBL height. We show that specific model configurations systematically overestimate or underestimate the PBL height averaged across the region with biases closely correlated with the choice of LSM, PBL scheme, and cumulus parameterization (CP). Domain average PBL wind speed is overestimated in nearly every model configuration. Both planetary boundary layer height (PBLH) and PBL wind speed biases show coherent spatial variations across the Midwest, with PBLH overestimated averaged across configurations by 300–400 m in the west, and PBL winds overestimated by about 1 m s−1 on average in the east. We find model configurations with lower biases averaged across the domain, but no single configuration is optimal across the entire region and for all meteorological variables. We conclude that model ensembles that include multiple physics parameterizations and meteorological initial conditions are likely to be necessary to encompass the atmospheric conditions most important to the transport of CO2 in the PBL, but that construction of such an ensemble will be challenging due to ensemble biases that vary across the region.

scholarly journals Uncertainty Propagation in Coupled Atmosphere–Wave–Ocean Prediction System: A Study of Hurricane Earl (2010)

2019 ◽  
Vol 147 (1) ◽  
pp. 221-245 ◽  
Author(s):  
Guotu Li ◽  
Milan Curcic ◽  
Mohamed Iskandarani ◽  
Shuyi S. Chen ◽  
Omar M. Knio
Keyword(s):  
Wind Speed ◽  
Initial Conditions ◽  
Uncertainty Propagation ◽  
Storm Track ◽  
Coupled Model ◽  
Surrogate Models ◽  
Maximum Wind Speed ◽  
Track Forecast ◽  
Initial Strength ◽  
Maximum Wind

This study focuses on understanding the evolution of Hurricane Earl (2010) with respect to random perturbations in the storm’s initial strength, size, and asymmetry in wind distribution. We rely on the Unified Wave Interface-Coupled Model (UWIN-CM), a fully coupled atmosphere–wave–ocean system to generate a storm realization ensemble, and use polynomial chaos (PC) expansions to build surrogate models for time evolution of both the maximum wind speed and minimum sea level pressure in Earl. The resulting PC surrogate models provide statistical insights on probability distributions of model responses throughout the simulation time span. Statistical analysis of rapid intensification (RI) suggests that initial perturbations having intensified and counterclockwise-rotated winds are more likely to undergo RI. In addition, for the range of initial conditions considered RI seems mostly sensitive to azimuthally averaged maximum wind speed and asymmetry orientation, rather than storm size and asymmetry magnitude; this is consistent with global sensitivity analysis of PC surrogate models. Finally, we combine initial condition perturbations with a stochastic kinetic energy backscatter scheme (SKEBS) forcing in the UWIN-CM simulations and conclude that the storm tracks are substantially influenced by the SKEBS forcing perturbations, whereas the perturbations in initial conditions alone had only limited impact on the storm-track forecast.

scholarly journals Impact of Parameterization of Physical Processes on Simulation of Track and Intensity of Tropical Cyclone Nargis (2008) with WRF-NMM Model

2012 ◽  
Vol 2012 ◽  
pp. 1-18 ◽  
Author(s):  
Sujata Pattanayak ◽  
U. C. Mohanty ◽  
Krishna K. Osuri
Keyword(s):  
Land Surface ◽  
Mesoscale Model ◽  
Initial Conditions ◽  
Surface Wind ◽  
India Meteorological Department ◽  
Heat Fluxes ◽  
Cumulus Convection ◽  
Cyclonic Storm ◽  
Parameterization Schemes ◽  
Cyclone Nargis

The present study is carried out to investigate the performance of different cumulus convection, planetary boundary layer, land surface processes, and microphysics parameterization schemes in the simulation of a very severe cyclonic storm (VSCS) Nargis (2008), developed in the central Bay of Bengal on 27 April 2008. For this purpose, the nonhydrostatic mesoscale model (NMM) dynamic core of weather research and forecasting (WRF) system is used. Model-simulated track positions and intensity in terms of minimum central mean sea level pressure (MSLP), maximum surface wind (10 m), and precipitation are verified with observations as provided by the India Meteorological Department (IMD) and Tropical Rainfall Measurement Mission (TRMM). The estimated optimum combination is reinvestigated with six different initial conditions of the same case to have better conclusion on the performance of WRF-NMM. A few more diagnostic fields like vertical velocity, vorticity, and heat fluxes are also evaluated. The results indicate that cumulus convection play an important role in the movement of the cyclone, and PBL has a crucial role in the intensification of the storm. The combination of Simplified Arakawa Schubert (SAS) convection, Yonsei University (YSU) PBL, NMM land surface, and Ferrier microphysics parameterization schemes in WRF-NMM give better track and intensity forecast with minimum vector displacement error.

scholarly journals Ensemble forecasting with a stochastic convective parametrization based on equilibrium statistics

2011 ◽  
Vol 11 (11) ◽  
pp. 30457-30485 ◽  
Author(s):  
P. Groenemeijer ◽  
G. C. Craig
Keyword(s):  
Large Scale ◽  
Mesoscale Model ◽  
Initial Conditions ◽  
Deep Convection ◽  
Horizontal Resolution ◽  
Ensemble Forecasting ◽  
Total Variance ◽  
Ensemble Prediction ◽  
Ensemble Prediction System ◽  
The Impact

Abstract. The stochastic Plant-Craig scheme for deep convection was implemented in the COSMO mesoscale model and used for ensemble forecasting. Ensembles consisting of 100 48 h forecasts at 7 km horizontal resolution were generated for a 2000 × 2000 km domain covering central Europe. Forecasts were made for seven case studies and characterized by different large-scale meteorological environments. Each 100 member ensemble consisted of 10 groups of 10 members, with each group driven by boundary and initial conditions from a selected member from the global ECMWF Ensemble Prediction System. The precipitation variability within and among these groups of members was computed, and it was found that the relative contribution to the ensemble variance introduced by the stochastic convection scheme was substantial, amounting to as much as 76% of the total variance in the ensemble in one of the studied cases. The impact of the scheme was not confined to the grid scale, and typically contributed 25–50% of the total variance even after the precipitation fields had been smoothed to a resolution of 35 km. The variability of precipitation introduced by the scheme was approximately proportional to the total amount of convection that occurred, while the variability due to large-scale conditions changed from case to case, being highest in cases exhibiting strong mid-tropospheric flow and pronounced meso- to synoptic scale vorticity extrema. The stochastic scheme was thus found to be an important source of variability in precipitation cases of weak large-scale flow lacking strong vorticity extrema, but high convective activity.

scholarly journals Assessing the Potential of a Long-Term Climate Forecast for Cuba Using the WRF Model

2021 ◽  
Vol 4 (1) ◽  
Author(s):  
Lourdes Álvarez-Escudero ◽  
Yandy G. Mayor ◽  
Israel Borrajero-Montejo ◽  
Arnoldo Bezanilla-Morlot
Keyword(s):  
Mesoscale Model ◽  
Initial Conditions ◽  
Wrf Model ◽  
Reanalysis Data ◽  
Weather Research And Forecasting ◽  
Tropical Storms ◽  
Precipitation Regime ◽  
Wide Debate ◽  
One Year ◽  
Cumulative Precipitation

Seasonal climatic prediction studies are a matter of wide debate all over the world. Cuba, a mainly agricultural nation, should greatly benefit from the knowledge, which is available months in advance of the precipitation regime and allows for the proper management of water resources. In this work, a series of six experiments were made with a mesoscale model WRF (Weather Research and Forecasting Model) that produced a 15-month forecast for each month of cumulative precipitation starting at two dates, and for three non-consecutive years with different meteorological characteristics: one dry year (2004), one year that started dry and turned rainy (2005), and one year where several tropical storms occurred (2008). ERA-Interim reanalysis data were used for the initial and border conditions and experiments started 1 month before the beginning of the rainy and the dry seasons, respectively. In a general sense, the experience of using WRF indicated that it was a valid resource for seasonal forecast, since the results obtained were in the same range as those reported by the literature for similar cases. Several limitations were revealed by the results: the forecasts underestimated the monthly cumulative precipitation figures, tropical storms entering through the borders sometimes followed courses different from the real courses inside the working domain, storms that developed inside the domain were not reproduced by WRF, and differences in initial conditions led to significantly different forecasts for the corresponding time steps (nonlinearity). Changing the model parameterizations and initial conditions of the ensemble forecast experiments was recommended.

scholarly journals Global simulations of the atmosphere at 1.45 km grid-spacing with the Integrated Forecasting System

2020 ◽  
Author(s):  
Peter Düben ◽  
Nils Wedi ◽  
Sami Saarinen ◽  
Christian Zeman
Keyword(s):  
Initial Conditions ◽  
Weather Forecast ◽  
Wave Model ◽  
Shallow Convection ◽  
Grid Spacing ◽  
Time Step ◽  
Weather Forecasts ◽  
Forecasting System ◽  
Diabatic Forcing ◽  
Physical Phenomena

&lt;p&gt;Global simulations with 1.45 km grid-spacing are presented that were performed with the Integrated Forecasting System (IFS) of the European Centre for Medium-Range Weather Forecasts (ECMWF). Simulations are uncoupled (without ocean, sea-ice or wave model), using 62 or 137 vertical levels and the full complexity of weather forecast simulations including recent date initial conditions, real-world topography, and state-of-the-art physical parametrizations and diabatic forcing including shallow convection, turbulent diffusion, radiation and five categories for the water substance (vapour, liquid, ice, rain, snow). Simulations are evaluated with regard to computational efficiency and model fidelity. Scaling results are presented that were performed on the fastest supercomputer in Europe - Piz Daint (Top 500, Nov 2018). Important choices for the model configuration at this unprecedented resolution for the IFS are discussed such as the use of hydrostatic and non-hydrostatic equations or the time resolution of physical phenomena which is defined by the length of the time step.&amp;#160;&lt;/p&gt;&lt;p&gt;Our simulations indicate that the IFS model &amp;#8212; based on spectral transforms with a semi-implicit, semi-Lagrangian time-stepping scheme in contrast to more local discretization techniques &amp;#8212; can provide a meaningful baseline reference for O(1) km global simulations.&lt;/p&gt;

Sign in / Sign up

Export Citation Format

Share Document