scholarly journals High-Resolution Simulation of Hurricane Bonnie (1998). Part II: Water Budget

2006 ◽  
Vol 63 (1) ◽  
pp. 43-64 ◽  
Author(s):  
Scott A. Braun

Abstract The fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5) is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weaker updrafts, indicative of a larger role of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where cooling from melting enhances condensation. Unlike in squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance in situ stratiform precipitation development through additional growth by vapor deposition and aggregation. This study also examines artificial source terms for cloud and precipitation mass associated with setting to zero negative mixing ratios that arise from numerical advection errors. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation and 15%–20% of the precipitation. Thus, these terms must be accounted for to balance the model budgets, and the results suggest the need for improved model numerics.

2009 ◽  
Vol 137 (1) ◽  
pp. 414-432 ◽  
Author(s):  
F. Couvreux ◽  
F. Guichard ◽  
P. H. Austin ◽  
F. Chen

Abstract Mesoscale water vapor heterogeneities in the boundary layer are studied within the context of the International H2O Project (IHOP_2002). A significant portion of the water vapor variability in the IHOP_2002 occurs at the mesoscale, with the spatial pattern and the magnitude of the variability changing from day to day. On 14 June 2002, an atypical mesoscale gradient is observed, which is the reverse of the climatological gradient over this area. The factors causing this water vapor variability are investigated using complementary platforms (e.g., aircraft, satellite, and in situ) and models. The impact of surface flux heterogeneities and atmospheric variability are evaluated separately using a 1D boundary layer model, which uses surface fluxes from the High-Resolution Land Data Assimilation System (HRLDAS) and early-morning atmospheric temperature and moisture profiles from a mesoscale model. This methodology, based on the use of robust modeling components, allows the authors to tackle the question of the nature of the observed mesoscale variability. The impact of horizontal advection is inferred from a careful analysis of available observations. By isolating the individual contributions to mesoscale water vapor variability, it is shown that the observed moisture variability cannot be explained by a single process, but rather involves a combination of different factors: the boundary layer height, which is strongly controlled by the surface buoyancy flux, the surface latent heat flux, the early-morning heterogeneity of the atmosphere, horizontal advection, and the radiative impact of clouds.


2006 ◽  
Vol 63 (1) ◽  
pp. 19-42 ◽  
Author(s):  
Scott A. Braun ◽  
Michael T. Montgomery ◽  
Zhaoxia Pu

Abstract The fifth-generation Pennsylvania State University–National Center for Atmospheric Research (PSU–NCAR) Mesoscale Model (MM5) is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine how vertical wind shear impacts the distribution of vertical motion in the eyewall on both the storm and cloud scale. As in many previous studies, it is found here that the shear produces a wavenumber-1 asymmetry in the time-averaged vertical motion and rainfall. Several mechanisms for this asymmetry are evaluated. The vertical motion asymmetry is qualitatively consistent with an assumed balance between horizontal vorticity advection by the relative flow and stretching of vorticity, with relative asymmetric inflow (convergence) at low levels and outflow (divergence) at upper levels on the downshear side of the eyewall. The simulation results also show that the upward motion portion of the eyewall asymmetry is located in the direction of vortex tilt, consistent with the vertical motion that required to maintain dynamic balance. Variations in the direction and magnitude of the tilt are consistent with the presence of a vortex Rossby wave quasi mode, which is characterized by a damped precession of the upper vortex relative to the lower vortex. While the time-averaged vertical motion is characterized by ascent in a shear-induced wavenumber-1 asymmetry, the instantaneous vertical motion is typically associated with deep updraft towers that generally form on the downtilt-right side of the eyewall and dissipate on the downtilt-left side. The updrafts towers are typically associated with eyewall mesovortices rotating cyclonically around the eyewall and result from an interaction between the shear-induced relative asymmetric flow and the cyclonic circulations of the mesovortices. The eyewall mesovortices may persist for more than one orbit around the eyewall and, in these cases, can initiate multiple episodes of upward motion.


2016 ◽  
Vol 33 (12) ◽  
pp. 2663-2678 ◽  
Author(s):  
Douglas Lowenthal ◽  
A. Gannet Hallar ◽  
Ian McCubbin ◽  
Robert David ◽  
Randolph Borys ◽  
...  

AbstractThe Isotopic Fractionation in Snow (IFRACS) study was conducted at Storm Peak Laboratory (SPL) in northwestern Colorado during the winter of 2014 to elucidate snow growth processes in mixed-phase clouds. The isotopic composition (δ18O and δD) of water vapor, cloud water, and snow in mixed-phase orographic clouds were measured simultaneously for the first time. The depletion of heavy isotopes [18O and deuterium (D)] was greatest for vapor, followed by snow, then cloud. The vapor, cloud, and snow compositions were highly correlated, suggesting similar cloud processes throughout the experiment. The isotopic composition of the water vapor was directly related to its concentration. Isotopic fractionation during condensation of vapor to cloud drops was accurately reproduced assuming equilibrium fractionation. This was not the case for snow, which grows by riming and vapor deposition. This implies stratification of vapor with altitude. The relationship between temperature at SPL and δ18O was used to show that the snow gained most of its mass within 922 m above SPL. Relatively invariant deuterium excess (d) in vapor, cloud water, and snow from day to day suggests a constant vapor source and Rayleigh fractionation during transport. The diurnal variation of vapor d reflected the differences between surface and free-tropospheric air during the afternoon and early morning hours, respectively. These observations will be used to validate simulations of snow growth using an isotope-enabled mesoscale model with explicit microphysics.


2010 ◽  
Vol 10 (6) ◽  
pp. 1129-1149 ◽  
Author(s):  
M. Milelli ◽  
M. Turco ◽  
E. Oberto

Abstract. The forecast in areas of very complex topography, as for instance the Alpine region, is still a challenge even for the new generation of numerical weather prediction models which aim at reaching the km-scale. The problem is enhanced by a general lack of standard observations, which is even more evident over the southern side of the Alps. For this reason, it would be useful to increase the performance of the mathematical models by locally assimilating non-conventional data. Since in ARPA Piemonte there is the availability of a great number of non-GTS stations, it has been decided to assimilate the 2 m temperature, coming from this dataset, in the very-high resolution version of the COSMO model, which has a horizontal resolution of about 3 km, more similar to the average resolution of the thermometers. Four different weather situations have been considered, ranging from spring to winter, from cloudy to clear sky. The aim of the work is to investigate the effects of the assimilation of non-GTS data in order to create an operational very high-resolution analysis, but also to test the option of running in the future a very short-range forecast starting from these analyses (RUC or Rapid Update Cycle). The results, in terms of Root Mean Square Error, Mean Error and diurnal cycle of some surface variables such as 2 m temperature, 2 m relative humidity and 10 m wind intensity show a positive impact during the assimilation cycle which tends to dissipate a few hours after the end of it. Moreover, the 2 m temperature assimilation has a slightly positive or neutral impact on the vertical profiles of temperature, eventhough some calibration is needed for the precipitation field which is too much perturbed during the assimilation cycle, while it is unaffected in the forecast period. So the stability of the planetary boundary layer, on the one hand, has not been particularly improved by the new-data assimilation, but, on the other hand, it has not been destroyed. It has to be pointed out that a correct description of the planetary boundary layer, even only the lowest part of it, could be helpful to the forecasters and, in general, to the users, in order to deal with meteorological hazards such as snow (in particular snow/rain limit definition), or fog (description of temperature inversions).


2006 ◽  
Vol 63 (12) ◽  
pp. 3421-3436 ◽  
Author(s):  
Marat Khairoutdinov ◽  
David Randall

Results are presented from a high-resolution three-dimensional simulation of shallow-to-deep convection transition based on idealization of observations made during the Large-Scale Biosphere–Atmosphere (LBA) experiment in Amazonia, Brazil, during the Tropical Rainfall Measuring Mission (TRMM)-LBA mission on 23 February. The doubly periodic grid has 1536 × 1536 × 256 grid cells with horizontal grid spacing of 100 m, thus covering an area of 154 × 154 km2. The vertical resolution varies from 50 m in the boundary layer to 100 m in the free troposphere and gradually coarsens to 250 m near the domain top at 25.4 km. The length of the simulation is 6 h, starting from an early morning sounding corresponding to 0730 local time. Convection is forced by prescribed surface latent and sensible heat fluxes and prescribed horizontally uniform radiative heating Despite a considerable amount of convective available potential energy (CAPE) in the range of 1600–2400 J kg−1, and despite virtually no convective inhibition (CIN) in the mean sounding throughout the simulation, the cumulus convection starts as shallow, gradually developing into congestus, and becomes deep only toward the end of simulation. Analysis shows that the reason is that the shallow clouds generated by the boundary layer turbulence are too small to penetrate deep into the troposphere, as they are quickly diluted by mixing with the environment. Precipitation and the associated cold pools are needed to generate thermals big enough to support the growth of deep clouds. This positive feedback involving precipitation is supported by a sensitivity experiment in which the cold pools are effectively eliminated by artificially switching off the evaporation of precipitation; in the experiment, the convection remains shallow throughout the entire simulation, with a few congestus but no deep clouds. The probability distribution function (PDF) of cloud size during the shallow, congestus, and deep phases is analyzed using a new method. During each of the three phases, the shallow clouds dominate the mode of the PDFs at about 1-km diameter. During the deep phase, the PDFs show cloud bases as wide as 4 km. Analysis of the joint PDFs of cloud size and in-cloud variables demonstrates that, as expected, the bigger clouds are far less diluted above their bases than their smaller counterparts. Also, thermodynamic properties at cloud bases are found to be nearly identical for all cloud sizes, with the moist static energy exceeding the mean value by as much as 4 kJ kg−1. The width of the moist static energy distribution in the boundary layer is mostly due to variability of water vapor; therefore, clouds appear to grow from the air with the highest water vapor content available. No undiluted cloudy parcels are found near the level of neutral buoyancy. It appears that a simple entraining-plume model explains the entrainment rates rather well. The least diluted plumes in the simulation correspond to an entrainment parameter of about 0.1 km−1.


2006 ◽  
Vol 63 (1) ◽  
pp. 127-150 ◽  
Author(s):  
Greg M. McFarquhar ◽  
Henian Zhang ◽  
Gerald Heymsfield ◽  
Jeffrey B. Halverson ◽  
Robbie Hood ◽  
...  

Abstract Fine-resolution simulations of Hurricane Erin are conducted using the fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5) to investigate roles of thermodynamic, boundary layer, and microphysical processes on Erin’s structure and evolution. Choice of boundary layer scheme has the biggest impact on simulations, with the minimum surface pressure (Pmin) averaged over the last 18 h (when Erin is relatively mature) varying by over 20 hPa. Over the same period, coefficients used to describe graupel fall speeds (Vg) affect Pmin by up to 7 hPa, almost equivalent to the maximum 9-hPa difference between microphysical parameterization schemes; faster Vg and schemes with more hydrometeor categories generally give lower Pmin. Compared to radar reflectivity factor (Z) observed by the NOAA P-3 lower fuselage radar and the NASA ER-2 Doppler radar (EDOP) in Erin, all simulations overpredict the normalized frequency of occurrence of Z larger than 40 dBZ and underpredict that between 20 and 40 dBZ near the surface; simulations overpredict Z larger than 25 to 30 dBZ and underpredict that between 15 and 25 or 30 dBZ near the melting layer, the upper limit depending on altitude. Brightness temperatures (Tb) computed from modeled fields at 37.1- and 85.5-GHz channels that respond to scattering by graupel-size ice show enhanced scattering, mainly due to graupel, compared to observations. Simulated graupel mixing ratios are about 10 times larger than values observed in other hurricanes. For the control run at 6.5 km averaged over the last 18 simulated hours, Doppler velocities computed from modeled fields (Vdop) greater than 5 m s−1 make up 12% of Erin’s simulated area for the base simulation but less than 2% of the observed area. In the eyewall, 5% of model updrafts above 9 km are stronger than 10 m s−1, whereas statistics from other hurricanes show that 5% of updrafts are stronger than only 5 m s−1. Variations in distributions of Z, vertical motion, and graupel mixing ratios between schemes are not sufficient to explain systematic offsets between observations and models. A new iterative condensation scheme, used with the Reisner mixed-phase microphysics scheme, limits unphysical increases of equivalent potential temperature associated with many condensation schemes and reduces the frequency of Z larger than 50 dBZ, but has minimal effect on Z below 50 dBZ, which represent 95% of the modeled hurricane rain area. However, the new scheme changes the Erin simulations in that 95% of the updrafts are weaker than 5 m s−1 and Pmin is up to 12 hPa higher over the last 18 simulated hours.


2018 ◽  
Author(s):  
Chang-Gong Shan ◽  
Wei Wang ◽  
Cheng Liu ◽  
You-Wen Sun ◽  
Yuan Tian ◽  
...  

Abstract. Observations of stable isotopes of water vapor provide important information for water cycle. The volume mixing ratios (VMR) of H2O (XH2O) and HDO (XHDO) have been retrieved based on a high-resolution ground-based Fourier transform infrared spectroscopy (FTIR) at Hefei site, and the isotopic composition δD was calculated. Time series of XH2O were compared with the Greenhouse gases Observing Satellite (GOSAT) data, showing a good agreement. The daily averaged δD ranges from −17.02 ‰ to −282.3 ‰ between September 2015 and September 2016. Also, the relationships of meteorological parameters with stable isotopologue were analyzed. δD values showed an obvious positive correlation with temperature and ln(XH2O) and a weak correlation with relative humidity. Further, 51.35 % of airmass at Hefei site comes from the southeast of China, and the main potential sources of δD are in the east of China over the observation period based on the back trajectories model. Furthermore, the δD values of evapotranspiration were calculated based on Keeling plot. Observations of the stable isotopes of water vapor by high-resolution ground-based FTIR provide information on study of the variation of the atmospheric water vapor at Hefei site.


2010 ◽  
Vol 3 (2) ◽  
pp. 1263-1301 ◽  
Author(s):  
L. Yurganov ◽  
W. McMillan ◽  
C. Wilson ◽  
M. Fischer ◽  
S. Biraud

Abstract. CO mixing ratios weighted over the bottom 2-km thick atmospheric layer between 2002 and 2009 were retrieved from downwelling infrared (IR) radiance spectra of the clear sky measured by a zenith-viewing Atmospheric Emitted Radiance Interferometer (AERI) deployed at the Southern Great Plains (SGP) observatory of the Atmospheric Radiation Measurements (ARM) Program near Lamont, Oklahoma. A version of the algorithm proposed by He at al. (2001) was significantly improved and validated. Essentially, the new algorithm retrieves a CO mixing ratio that is determined by the convolution of the a priori profile (assumed to be constant with altitude), the true profile, and the averaging kernel which maximizes near the surface. Approximately 70% of the CO signal comes from the boundary layer and the remaining 30% come from the lower part of the free troposphere. Archived temperature and water vapor profiles retrieved from the same AERI spectra through automated ARM processing were used as input data for the CO retrievals. We found the archived water vapor profiles required additional constraint using SGP Microwave Radiometer retrievals of total precipitable water vapor. Additionally, a correction for scattered solar light was developed. The retrieved CO was validated using simultaneous independently measured CO profiles. An aircraft supplied in situ CO measurements at altitudes up to 4572 m above sea level once or twice a week between March 2006 and December 2008. The aircraft measurements were supplemented with ground-based CO measurements at the SGP and retrievals from the Atmospheric IR Sounder (AIRS) above 5 km to create full tropospheric CO profiles. Comparison of the convolved profiles to the AERI CO retrievals found a squared correlation coefficient of 0.57, a standard deviation of ±11.7 ppbv, a bias of 16 ppbv, and a slope of 0.92. Averaged seasonal and diurnal cycles measured by AERI are compared with those measured continuously in situ at the SGP in the boundary layer. Monthly mean CO values measured by AERI between 2002 and 2009 are compared with those measured by AIRS over North America, the Northern Hemisphere mid-latitudes, and over the tropics.


2018 ◽  
Vol 33 (6) ◽  
pp. 1547-1566 ◽  
Author(s):  
Clark Evans ◽  
Steven J. Weiss ◽  
Israel L. Jirak ◽  
Andrew R. Dean ◽  
David S. Nevius

Abstract This study evaluates forecast vertical thermodynamic profiles and derived thermodynamic parameters from two regional/convection-allowing model pairs, the North American Mesoscale Forecast System and the North American Mesoscale Nest model pair and the Rapid Refresh and High Resolution Rapid Refresh model pair, in warm-season, thunderstorm-supporting environments. Differences in bias and mean absolute error between the regional and convection-allowing models in each of the two pairs, while often statistically significant, are practically small for the variables, parameters, and vertical levels considered, such that the smaller-scale variability resolved by convection-allowing models does not degrade their forecast skill. Model biases shared by the regional and convection-allowing models in each pair are documented, particularly the substantial cool and moist biases in the planetary boundary layer arising from the Mellor–Yamada–Janjić planetary boundary layer parameterization used by the North American Mesoscale model and the Nest version as well as the middle-tropospheric moist bias shared by the Rapid Refresh and High Resolution Rapid Refresh models. Bias and mean absolute errors typically have larger magnitudes in the evening, when buoyancy is a significant contributor to turbulent vertical mixing, than in the morning. Vertical thermodynamic profile biases extend over a deep vertical layer in the western United States given strong sensible heating of the underlying surface. The results suggest that convection-allowing models can fulfill the use cases typically and historically met by regional models in operations at forecast entities such as the Storm Prediction Center, a fruitful finding given the proposed elimination of regional models with the Next-Generation Global Prediction System initiative.


2017 ◽  
Author(s):  
Maria Filioglou ◽  
Anna Nikandrova ◽  
Sami Niemelä ◽  
Holger Baars ◽  
Tero Mielonen ◽  
...  

Abstract. We present tropospheric water vapor profiles measured with a Raman lidar during three field campaigns held in Finland. Co-located radio soundings are available throughout the period for the calibration of the lidar signals. We investigate the possibility of calibrating the lidar water vapor profiles in the absence of co-existing on site soundings using water vapor profiles from the combined Advanced Infrared Sounder (AIRS) and the Advanced Microwave Radiometer (AMSU) satellite product; the Aire Limitee Adaptation dynamique Development International and High Resolution Limited Area Model (ALADIN HIRLAM) numerical weather prediction (NWP) system, and the nearest radio sounding station located 100 km away from the lidar site (only for the permanent location of the lidar). The uncertainties of the calibration factor derived from the soundings, the satellite and the model data are


Sign in / Sign up

Export Citation Format

Share Document