scholarly journals Experimental assimilation of synthetic bogus tropical cyclone pressure observations into a high-resolution rapid-update NWP model

Author(s):  
Susan Rennie ◽  
Jim Fraser

The effect of synthetic ‘bogus’ tropical cyclone (TC) central pressure observations on TC Owen was tested in a convective-scale numerical weather prediction (NWP) system with hourly 4D-Var assimilation. TC Owen traversed the Gulf of Carpentaria over 10–14 December 2018, entering from the east and briefly making landfall on the western edge before reversing course and retracing its path east to cross the northern tip of Queensland. The Australian Bureau of Meteorology runs a high-resolution NWP model centred over Darwin, which covers much of the Gulf of Carpentaria. The next-generation developmental version of this model includes data assimilation. Therefore, when TC Owen presented the opportunity to investigate the simulation of a TC within the domain, the developmental system was run as a case study. The modelled cyclone initially failed to intensify. The case study was then repeated including assimilation of bogus central pressure observations. This new run showed a large improvement in the intensity throughout the simulation; however, the TC track was not substantially improved. This demonstration of the potential impact of using synthetic observations may guide whether the development of a bogus observation source with sufficiently low latency for use in an hourly-cycling system should be prioritised.

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.


2013 ◽  
Vol 6 (4) ◽  
pp. 7425-7472
Author(s):  
U. Schumann ◽  
R. Hempel ◽  
H. Flentje ◽  
M. Garhammer ◽  
K. Graf ◽  
...  

Abstract. Photogrammetric methods and analysis results for contrails observed with wide-angle cameras are described. Four cameras of two different types (view angle < 90° or whole-sky imager) at the ground at various positions are used to track contrails and to derive their altitude, width, and horizontal speed. Camera models for both types are described to derive the observation angles for given image coordinates and their inverse. The models are calibrated with sightings of the Sun, the Moon and a few bright stars. The methods are applied and tested in a case study. Four persistent contrails crossing each other together with a short-lived one are observed with the cameras. Vertical and horizontal positions of the contrails are determined from the camera images to an accuracy of better than 200 m and horizontal speed to 0.2 m s−1. With this information, the aircraft causing the contrails are identified by comparison to traffic waypoint data. The observations are compared with synthetic camera pictures of contrails simulated with the contrail prediction model CoCiP, a Lagrangian model using air traffic movement data and numerical weather prediction (NWP) data as input. The results provide tests for the NWP and contrail models. The cameras show spreading and thickening contrails suggesting ice-supersaturation in the ambient air. The ice-supersaturated layer is found thicker and more humid in this case than predicted by the NWP model used. The simulated and observed contrail positions agree up to differences caused by uncertain wind data. The contrail widths, which depend on wake vortex spreading, ambient shear and turbulence, were partly wider than simulated.


2015 ◽  
Vol 143 (10) ◽  
pp. 4012-4037 ◽  
Author(s):  
Colin M. Zarzycki ◽  
Christiane Jablonowski

Abstract Tropical cyclone (TC) forecasts at 14-km horizontal resolution (0.125°) are completed using variable-resolution (V-R) grids within the Community Atmosphere Model (CAM). Forecasts are integrated twice daily from 1 August to 31 October for both 2012 and 2013, with a high-resolution nest centered over the North Atlantic and eastern Pacific Ocean basins. Using the CAM version 5 (CAM5) physical parameterization package, regional refinement is shown to significantly increase TC track forecast skill relative to unrefined grids (55 km, 0.5°). For typical TC forecast integration periods (approximately 1 week), V-R forecasts are able to nearly identically reproduce the flow field of a globally uniform high-resolution forecast. Simulated intensity is generally too strong for forecasts beyond 72 h. This intensity bias is robust regardless of whether the forecast is forced with observed or climatological sea surface temperatures and is not significantly mitigated in a suite of sensitivity simulations aimed at investigating the impact of model time step and CAM’s deep convection parameterization. Replacing components of the default physics with Cloud Layers Unified by Binormals (CLUBB) produces a statistically significant improvement in forecast intensity at longer lead times, although significant structural differences in forecasted TCs exist. CAM forecasts the recurvature of Hurricane Sandy into the northeastern United States 60 h earlier than the Global Forecast System (GFS) model using identical initial conditions, demonstrating the sensitivity of TC forecasts to model configuration. Computational costs associated with V-R simulations are dramatically decreased relative to globally uniform high-resolution simulations, demonstrating that variable-resolution techniques are a promising tool for future numerical weather prediction applications.


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.


2010 ◽  
Vol 25 (4) ◽  
pp. 1293-1306 ◽  
Author(s):  
Charles R. Sampson ◽  
Paul A. Wittmann ◽  
Hendrik L. Tolman

Abstract A new algorithm to generate wave heights consistent with tropical cyclone official forecasts from the Joint Typhoon Warning Center (JTWC) has been developed. The process involves generating synthetic observations from the forecast track and the 34-, 50-, and 64-kt wind radii. The JTWC estimate of the radius of maximum winds is used in the algorithm to generate observations for the forecast intensity (wind), and the JTWC-estimated radius of the outermost closed isobar is used to assign observations at the outermost extent of the tropical cyclone circulation. These observations are then interpolated to a high-resolution latitude–longitude grid covering the entire extent of the circulation. Finally, numerical weather prediction (NWP) model fields are obtained for each forecast time, the NWP model forecast tropical cyclone is removed from these fields, and the new JTWC vortex is inserted without blending zones between the vortex and the background. These modified fields are then used as input into a wave model to generate waves consistent with the JTWC forecasts. The algorithm is applied to Typhoon Yagi (2006), in anticipation of which U.S. Navy ships were moved from Tokyo Bay to an area off the southeastern coast of Kyushu. The decision to move (sortie) the ships was based on NWP model-driven long-range wave forecasts that indicated high seas impacting the coast in the vicinity of Tokyo Bay. The sortie decision was made approximately 84 h in advance of the high seas in order to give ships time to steam the approximately 500 n mi to safety. Results from the new algorithm indicate that the high seas would not affect the coast near Tokyo Bay within 84 h. This specific forecast verifies, but altimeter observations show that it does not outperform, the NWP model-driven wave analysis and forecasts for this particular case. Overall, the performance of the new algorithm is dependent on the JTWC tropical cyclone forecast performance, which has generally outperformed those of the NWP model over the last several years.


Author(s):  
Michał Z. Ziemiański ◽  
Damian K. Wójcik ◽  
Bogdan Rosa ◽  
Zbigniew P. Piotrowski

AbstractThis paper presents the semi-implicit compressible EULAG as a new dynamical core for convective-scale numerical weather prediction. The core is implemented within the infrastructure of the operational model of the Consortium for Small Scale Modeling (COSMO), forming the NWP COSMO-EULAG model (CE). This regional high-resolution implementation of the dynamical core complements its global implementation in the Finite-Volume Module of ECMWF’s Integrated Forecasting System. The paper documents the first operational-like application of the dynamical core for realistic weather forecasts. After discussing the formulation of the core and its coupling with the host model, the paper considers several high-resolution prognostic experiments over complex Alpine orography. Standard verification experiments examine the sensitivity of the CE forecast to the choice of the advection routine and assess the forecast skills against those of the default COSMO Runge-Kutta dynamical core at 2.2 km grid size showing a general improvement. The skills are also compared using satellite observations for a weak-flow convective Alpine weather case-study, showing favorable results. Additional validation of the new CE framework for partly convection-resolving forecasts using 1.1 km, 0.55 km, 0.22 km, and 0.1 km grids, designed to challenge its numerics and test the dynamics-physics coupling, demonstrates its high robustness in simulating multi-phase flows over complex mountain terrain, with slopes reaching 85 degrees, and the flow’s realistic representation.


2021 ◽  
Author(s):  
Alberto Caldas-Alvarez ◽  
Samiro Khodayar ◽  
Peter Knippertz

Abstract. Heavy precipitation is one of the most devastating weather extremes in the western Mediterranean region. Our capacity to prevent negative impacts from such extreme events requires advancements in numerical weather prediction, data assimilation and new observation techniques. In this paper we investigate the impact of two state-of-the-art data sets with very high resolution, Global Positioning System-Zenith Total Delays (GPS-ZTD) with a 10 min temporal resolution and radiosondes with ~700 levels, on the representation of convective precipitation in nudging experiments. Specifically, we investigate whether the high temporal resolution, quality, and coverage of GPS-ZTDs can outweigh their lack of vertical information or if radiosonde profiles are more valuable despite their scarce coverage and low temporal resolution (24 h to 6 h). The study focuses on the Intensive Observation Period 6 (IOP6) of the Hydrological Cycle in the Mediterranean eXperiment (HyMeX; 24 September 2012). This event is selected due to its severity (100 mm/12 h), the availability of observations for nudging and validation, and the large observation impact found in preliminary sensitivity experiments. We systematically compare simulations performed with the COnsortium for Small scale MOdelling (COSMO) model assimilating GPS, high- and low vertical resolution radiosoundings in model resolutions of 7 km, 2.8 km and 500 m. The results show that the additional GPS and radiosonde observations cannot compensate errors in the model dynamics and physics. In this regard the reference COSMO runs have an atmospheric moisture wet bias prior to precipitation onset but a negative bias in rainfall, indicative of deficiencies in the numerics and physics, unable to convert the moisture excess into sufficient precipitation. Nudging GPS and high-resolution soundings corrects atmospheric humidity, but even further reduces total precipitation. This case study also demonstrates the potential impact of individual observations in highly unstable environments. We show that assimilating a low-resolution sounding from Nimes (southern France) while precipitation is taking place induces a 40 % increase in precipitation during the subsequent three hours. This precipitation increase is brought about by the moistening of the 700  hPa level (7.5 g kg−1) upstream of the main precipitating systems, reducing the entrainment of dry air above the boundary layer. The moist layer was missed by GPS observations and high-resolution soundings alike, pointing to the importance of profile information and timing. However, assimilating GPS was beneficial for simulating the temporal evolution of precipitation. Finally, regarding the scale dependency, no resolution is particularly sensitive to a specific observation type, however the 2.8 km run has overall better scores, possibly as this is the optimally tuned operational version of COSMO. In follow-up experiments the Icosahedral Nonhydrostatic Model (ICON) will be investigated for this case study to assert whether its numerical and physics updates, compared to its predecessor COSMO, are able to improve the quality of the simulations.


2020 ◽  
Vol 148 (10) ◽  
pp. 4247-4265 ◽  
Author(s):  
Domingo Muñoz-Esparza ◽  
Robert D. Sharman ◽  
Stanley B. Trier

AbstractMesoscale numerical weather prediction (NWP) models are routinely exercised at kilometer-scale horizontal grid spacings (Δx). Such fine grids will usually allow at least partial resolution of small-scale gravity waves and turbulence in the upper troposphere and lower stratosphere (UTLS). However, planetary boundary layer (PBL) parameterization schemes used with these NWP model simulations typically apply explicit subgrid-scale vertical diffusion throughout the entire vertical extent of the domain, an effect that cannot be ignored. By way of an example case of observed widespread turbulence over the U.S. Great Plains, we demonstrate that the PBL scheme’s mixing in NWP model simulations of Δx = 1 km can have significant effects on the onset and characteristics of the modeled UTLS gravity waves. Qualitatively, PBL scheme diffusion is found to affect not only background conditions responsible for UTLS wave activity, but also to control the local vertical mixing that triggers or hinders the onset and propagation of these waves. Comparisons are made to a reference large-eddy simulation with Δx = 250 m to statistically quantify these effects. A significant and systematic overestimation of resolved vertical velocities, wave-scale fluxes, and kinetic energy is uncovered in the 1-km simulations, both in clear-air and in-cloud conditions. These findings are especially relevant for upper-level gravity wave and turbulence simulations using high-resolution kilometer-scale NWP models.


2020 ◽  
Vol 13 (5) ◽  
pp. 2279-2298
Author(s):  
Guillaume Thomas ◽  
Jean-François Mahfouf ◽  
Thibaut Montmerle

Abstract. This paper presents the potential of nonlinear and linear versions of an observation operator for simulating polarimetric variables observed by weather radars. These variables, deduced from the horizontally and vertically polarized backscattered radiations, give information about the shape, the phase and the distributions of hydrometeors. Different studies in observation space are presented as a first step toward their inclusion in a variational data assimilation context, which is not treated here. Input variables are prognostic variables forecasted by the AROME-France numerical weather prediction (NWP) model at convective scale, including liquid and solid hydrometeor contents. A nonlinear observation operator, based on the T-matrix method, allows us to simulate the horizontal and the vertical reflectivities (ZHH and ZVV), the differential reflectivity ZDR, the specific differential phase KDP and the co-polar correlation coefficient ρHV. To assess the uncertainty of such simulations, perturbations have been applied to input parameters of the operator, such as dielectric constant, shape and orientation of the scatterers. Statistics of innovations, defined by the difference between simulated and observed values, are then performed. After some specific filtering procedures, shapes close to a Gaussian distribution have been found for both reflectivities and for ZDR, contrary to KDP and ρHV. A linearized version of this observation operator has been obtained by its Jacobian matrix estimated with the finite difference method. This step allows us to study the sensitivity of polarimetric variables to hydrometeor content perturbations, in the model geometry as well as in the radar one. The polarimetric variables ZHH and ZDR appear to be good candidates for hydrometeor initialization, while KDP seems to be useful only for rain contents. Due to the weak sensitivity of ρHV, its use in data assimilation is expected to be very challenging.


Sign in / Sign up

Export Citation Format

Share Document