Improvements to an intermediate complexity atmospheric model for high-resolution downscaling in very complex terrain.

Author(s):  
Dylan Reynolds ◽  
Bert Kruyt ◽  
Ethan Gutmann ◽  
Tobias Jonas ◽  
Michael Lehning ◽  
...  

<p>            Snow deposition patterns in complex terrain are heavily dependent on the underlying topography. This topography affects precipitating clouds at the kilometer-scale and causes changes to the wind field at the sub-kilometer scale, resulting in altered advection of falling hydrometeors. Snow particles are particularly sensitive to changes in the near-surface flow field due to their low density. Atmospheric models which run at the kilometer scale cannot resolve the actual heterogeneity of the underlying terrain, resulting in precipitation maps which do not capture terrain-affected precipitation patterns. Thus, snow-atmosphere interactions such as preferential deposition are often not resolved in precipitation data used as input to snow models. To bridge this spatial gap and resolve snow-atmosphere interactions at the sub-kilometer scale, we couple an intermediate complexity atmospheric model (ICAR) to the COSMO NWP model. Applying this model to sub-kilometer terrain (horizontal resolution of 50 and 250 m) required changes to ICAR’s computational grid, atmospheric dynamics, and boundary layer flow. As a result, the near-surface flow now accounts for surface roughness and topographically induced speed up. This has been achieved by using terrain descriptors calculated once at initialization which consider a point’s exposure or sheltering relative to surrounding terrain. In particular, the use of a 3-dimensional Sx parameter allows us to simulate areas of stagnation and recirculation on the lee of terrain features. Our approach maintains the accurate large-scale precipitation patterns from COSMO but resolves the dynamics induced by terrain at the sub-kilometer scale without adding additional computational burden. We find that solid precipitation patterns at the ridge scale, such as preferential deposition of snow, are better resolved in the high-resolution version of ICAR than the current ICAR or COSMO models. This updated version of ICAR presents a new tool to dynamically downscale NWP output for snow models and enables future studies of snow-atmosphere interactions at domain scales of 100’s of kilometers.</p>

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).


2018 ◽  
Vol 12 (10) ◽  
pp. 3137-3160 ◽  
Author(s):  
Franziska Gerber ◽  
Nikola Besic ◽  
Varun Sharma ◽  
Rebecca Mott ◽  
Megan Daniels ◽  
...  

Abstract. Snow distribution in complex alpine terrain and its evolution in the future climate is important in a variety of applications including hydropower, avalanche forecasting and freshwater resources. However, it is still challenging to quantitatively forecast precipitation, especially over complex terrain where the interaction between local wind and precipitation fields strongly affects snow distribution at the mountain ridge scale. Therefore, it is essential to retrieve high-resolution information about precipitation processes over complex terrain. Here, we present very-high-resolution Weather Research and Forecasting model (WRF) simulations (COSMO–WRF), which are initialized by 2.2 km resolution Consortium for Small-scale Modeling (COSMO) analysis. To assess the ability of COSMO–WRF to represent spatial snow precipitation patterns, they are validated against operational weather radar measurements. Estimated COSMO–WRF precipitation is generally higher than estimated radar precipitation, most likely due to an overestimation of orographic precipitation enhancement in the model. The high precipitation amounts also lead to a higher spatial variability in the model compared to radar estimates. Overall, an autocorrelation and scale analysis of radar and COSMO–WRF precipitation patterns at a horizontal grid spacing of 450 m show that COSMO–WRF captures the spatial variability normalized by the domain-wide variability in precipitation patterns down to the scale of a few kilometers. However, simulated precipitation patterns systematically show a lower variability on the smallest scales of a few hundred meters compared to radar estimates. A comparison of spatial variability for different model resolutions gives evidence for an improved representation of local precipitation processes at a horizontal resolution of 50 m compared to 450 m. Additionally, differences of precipitation between 2830 m above sea level and the ground indicate that near-surface processes are active in the model.


2018 ◽  
Vol 18 (11) ◽  
pp. 2991-3006 ◽  
Author(s):  
Matthew D. K. Priestley ◽  
Helen F. Dacre ◽  
Len C. Shaffrey ◽  
Kevin I. Hodges ◽  
Joaquim G. Pinto

Abstract. Extratropical cyclones are the most damaging natural hazard to affect western Europe. Serial clustering occurs when many intense cyclones affect one specific geographic region in a short period of time which can potentially lead to very large seasonal losses. Previous studies have shown that intense cyclones may be more likely to cluster than less intense cyclones. We revisit this topic using a high-resolution climate model with the aim to determine how important clustering is for windstorm-related losses. The role of windstorm clustering is investigated using a quantifiable metric (storm severity index, SSI) that is based on near-surface meteorological variables (10 m wind speed) and is a good proxy for losses. The SSI is used to convert a wind footprint into losses for individual windstorms or seasons. 918 years of a present-day ensemble of coupled climate model simulations from the High-Resolution Global Environment Model (HiGEM) are compared to ERA-Interim reanalysis. HiGEM is able to successfully reproduce the wintertime North Atlantic/European circulation, and represent the large-scale circulation associated with the serial clustering of European windstorms. We use two measures to identify any changes in the contribution of clustering to the seasonal windstorm loss as a function of return period. Above a return period of 3 years, the accumulated seasonal loss from HiGEM is up to 20 % larger than the accumulated seasonal loss from a set of random resamples of the HiGEM data. Seasonal losses are increased by 10 %–20 % relative to randomized seasonal losses at a return period of 200 years. The contribution of the single largest event in a season to the accumulated seasonal loss does not change with return period, generally ranging between 25 % and 50 %. Given the realistic dynamical representation of cyclone clustering in HiGEM, and comparable statistics to ERA-Interim, we conclude that our estimation of clustering and its dependence on the return period will be useful for informing the development of risk models for European windstorms, particularly for longer return periods.


2016 ◽  
Author(s):  
R. J. Haarsma ◽  
M. Roberts ◽  
P. L. Vidale ◽  
C. A. Senior ◽  
A. Bellucci ◽  
...  

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest the possibility for significant changes in both large-scale aspects of circulation, as well as improvements in small-scale processes and extremes. However, such high resolution global simulations at climate time scales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centers and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other MIPs. Increases in High Performance Computing (HPC) resources, as well as the revised experimental design for CMIP6, now enables a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950-2050, with the possibility to extend to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulation. HighResMIP thereby focuses on one of the CMIP6 broad questions: “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.


Author(s):  
Antonio Parodi ◽  
Martina Lagasio ◽  
Agostino N. Meroni ◽  
Flavio Pignone ◽  
Francesco Silvestro ◽  
...  

AbstractBetween the 4th and the 6th of November 1994, Piedmont and the western part of Liguria (two regions in north-western Italy) were hit by heavy rainfalls that caused the flooding of the Po, the Tanaro rivers and several of their tributaries, causing 70 victims and the displacement of over 2000 people. At the time of the event, no early warning system was in place and the concept of hydro-meteorological forecasting chain was in its infancy, since it was still limited to a reduced number of research applications, strongly constrained by coarse-resolution modelling capabilities both on the meteorological and the hydrological sides. In this study, the skills of the high-resolution CIMA Research Foundation operational hydro-meteorological forecasting chain are tested in the Piedmont 1994 event. The chain includes a cloud-resolving numerical weather prediction (NWP) model, a stochastic rainfall downscaling model, and a continuous distributed hydrological model. This hydro-meteorological chain is tested in a set of operational configurations, meaning that forecast products are used to initialise and force the atmospheric model at the boundaries. The set consists of four experiments with different options of the microphysical scheme, which is known to be a critical parameterisation in this kind of phenomena. Results show that all the configurations produce an adequate and timely forecast (about 2 days ahead) with realistic rainfall fields and, consequently, very good peak flow discharge curves. The added value of the high resolution of the NWP model emerges, in particular, when looking at the location of the convective part of the event, which hit the Liguria region.


2018 ◽  
Author(s):  
Toni Mitovski ◽  
Jason N. S. Cole ◽  
Norman A. McFarlane ◽  
Knut von Salzen ◽  
Guang J. Zhang

Abstract. Changes in the large-scale environment during convective precipitation events in the Tropical Western Pacific simulated by version 4.3 of the Canadian Atmospheric Model (CanAM4.3) is compared against those simulated by version 5.0 of the super parameterized Community Atmosphere Model (spCAM5). This is done by compositing sub-hourly output of convective rainfall, convective available potential energy (CAPE), CAPE generation due to large-scale forcing in the free troposphere (dCAPELSFT), and near surface vertical velocity (ω) over the time period May–July 1997. Compared to spCAM5, CanAM4.3 tends to produce more frequent light convective precipitation ( 2 mm h−1). In spCAM5 5 % of convective precipitation events lasted less than 1.5 h and 75 % lasted between 1.5 and 3.0 h while in CanAM4.3 80 % of the events lasted less than 1.5 h. Convective precipitation in spCAM5 is found to be a function of dCAPELSFT and the large-scale near surface ω with variations in ω slightly leading variations in convective precipitation. Convective precipitation in CanAM4.3 does not have the same dependency and instead is found to be a function of CAPE.


2020 ◽  
Author(s):  
Xavier Lapillonne ◽  
William Sawyer ◽  
Philippe Marti ◽  
Valentin Clement ◽  
Remo Dietlicher ◽  
...  

<p>The ICON modelling framework is a unified numerical weather and climate model used for applications ranging from operational numerical weather prediction to low and high resolution climate projection. In view of further pushing the frontier of possible applications and to make use of the latest evolution in hardware technologies, parts of the model were recently adapted to run on heterogeneous GPU system. This initial GPU port focus on components required for high-resolution climate application, and allow considering multi-years simulations at 2.8 km on the Piz Daint heterogeneous supercomputer. These simulations are planned as part of the QUIBICC project “The Quasi-Biennial Oscillation (QBO) in a changing climate”, which propose to investigate effects of climate change on the dynamics of the QBO.</p><p>Because of the low compute intensity of atmospheric model the cost of data transfer between CPU and GPU at every step of the time integration would be prohibitive if only some components would be ported to the accelerator. We therefore present a full port strategy where all components required for the simulations are running on the GPU. For the dynamics, most of the physical parameterizations and infrastructure code the OpenACC compiler directives are used. For the soil parameterization, a Fortran based domain specific language (DSL) the CLAW-DSL has been considered. We discuss the challenges associated to port a large community code, about 1 million lines of code, as well as to run simulations on large-scale system at 2.8 km horizontal resolution in terms of run time and I/O constraints. We show performance comparison of the full model on CPU and GPU, achieving a speed up factor of approximately 5x, as well as scaling results on up to 2000 GPU nodes. Finally we discuss challenges and planned development regarding performance portability and high level DSL which will be used with the ICON model in the near future.</p>


2015 ◽  
Vol 8 (8) ◽  
pp. 2645-2653 ◽  
Author(s):  
C. G. Nunalee ◽  
Á. Horváth ◽  
S. Basu

Abstract. Recent decades have witnessed a drastic increase in the fidelity of numerical weather prediction (NWP) modeling. Currently, both research-grade and operational NWP models regularly perform simulations with horizontal grid spacings as fine as 1 km. This migration towards higher resolution potentially improves NWP model solutions by increasing the resolvability of mesoscale processes and reducing dependency on empirical physics parameterizations. However, at the same time, the accuracy of high-resolution simulations, particularly in the atmospheric boundary layer (ABL), is also sensitive to orographic forcing which can have significant variability on the same spatial scale as, or smaller than, NWP model grids. Despite this sensitivity, many high-resolution atmospheric simulations do not consider uncertainty with respect to selection of static terrain height data set. In this paper, we use the Weather Research and Forecasting (WRF) model to simulate realistic cases of lower tropospheric flow over and downstream of mountainous islands using the default global 30 s United States Geographic Survey terrain height data set (GTOPO30), the Shuttle Radar Topography Mission (SRTM), and the Global Multi-resolution Terrain Elevation Data set (GMTED2010) terrain height data sets. While the differences between the SRTM-based and GMTED2010-based simulations are extremely small, the GTOPO30-based simulations differ significantly. Our results demonstrate cases where the differences between the source terrain data sets are significant enough to produce entirely different orographic wake mechanics, such as vortex shedding vs. no vortex shedding. These results are also compared to MODIS visible satellite imagery and ASCAT near-surface wind retrievals. Collectively, these results highlight the importance of utilizing accurate static orographic boundary conditions when running high-resolution mesoscale models.


2016 ◽  
Vol 9 (11) ◽  
pp. 4185-4208 ◽  
Author(s):  
Reindert J. Haarsma ◽  
Malcolm J. Roberts ◽  
Pier Luigi Vidale ◽  
Catherine A. Senior ◽  
Alessio Bellucci ◽  
...  

Abstract. Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.


2008 ◽  
Vol 2 (1) ◽  
pp. 81-87 ◽  
Author(s):  
O. Jorba ◽  
C. Marrero ◽  
E. Cuevas ◽  
J. M. Baldasano

Abstract. On 28–29 November 2005 an extratropical storm affected the Canary Islands causing significant damage related to high average wind speeds and intense gusts over some islands of the archipelago. Delta was the twenty-sixth tropical or subtropical storm of the 2005 Atlantic hurricane season. It represents an unusual meteorological phenomenon for that region, and its impacts were underestimated by the different operational meteorological forecasts during the previous days of the arrival of the low near Canary Islands. The aim of this study is to reproduce the local effects of the flow that were observed over the Canary Islands during the travel of the Delta storm near the region using high-resolution mesoscale meteorological simulations. The Advanced Research Weather Research & Forecasting Model (WRF-ARW) is applied at 9, 3 and 1 km horizontal resolution using ECMWF forecasts as initial and boundary conditions. The high-resolution simulation will outline the main features that contributed to the high wind speeds observed in the archipelago. Variations in vertical static stability, vertical windshear and the intense synoptic winds of the southwestern part of Delta with a warm core at 850 hPa were the main characteristics that contributed to the development and amplification of intense gravity waves while the large-scale flow interacted with the complex topography of the islands.


Sign in / Sign up

Export Citation Format

Share Document