Search Results
You are looking at 31 - 40 of 61 items for
- Author or Editor: Steven E. Koch x
- Refine by Access: All Content x
Abstract
An objective analysis scheme based on the Barnes technique and designed for use on an interactive computer is described. In order to meet the specific needs of the research meteorologist, the interactive Barnes scheme allows real-time assessments both of the quality of the resulting analyses and of the impact of satellite-derived data upon various meteorological data sets. Display of a number of statistical and mapped analysis quality control indicators aid the impact assessments. Simple means for taking account of the spatially clustered nature typical of satellite data are included in the internal computations of the relative weights of data at grid point locations.
An analyst is allowed the capability of modifying values of certain input parameters to the interactive Barnes scheme within internally set limits. These constraints were objectively determined and tested in a number of different situations prior to implementation. The following constraints are employed: 1) calculation of the weights as a function of a data spacing representative of the data distribution; 2) automatic elimination of detail at wavelengths smaller than twice the representative data spacing; 3) placement of bounds upon the grid spacing by the data spacing; and 4) setting of a fixed limit on the number of passes through the data to achieve rapid and sufficient convergence of the analyzed values to the observed ones. A mathematical analysis of the convergence properties of the Barnes technique is presented to support the validity of the latter constraint.
Despite these constraints, the interactive Barnes scheme remains versatile because it accepts limited inputs to the data and grid display areas, to the data and grid spacings, and to the rate of convergence of the analysis to the observations. Input parameter values are entered through a series of questions displayed on a computer video terminal and by manipulation of display function devices. The analyst immediately sees a plot of the data, the contoured grid values, superimposed in various colors if desired, and the effects of choice of analysis options. Examples of both meteorological and satellite data analyses are presented to demonstrate the objectivity, versatility and practicality of the interactive Barnes scheme.
Abstract
An objective analysis scheme based on the Barnes technique and designed for use on an interactive computer is described. In order to meet the specific needs of the research meteorologist, the interactive Barnes scheme allows real-time assessments both of the quality of the resulting analyses and of the impact of satellite-derived data upon various meteorological data sets. Display of a number of statistical and mapped analysis quality control indicators aid the impact assessments. Simple means for taking account of the spatially clustered nature typical of satellite data are included in the internal computations of the relative weights of data at grid point locations.
An analyst is allowed the capability of modifying values of certain input parameters to the interactive Barnes scheme within internally set limits. These constraints were objectively determined and tested in a number of different situations prior to implementation. The following constraints are employed: 1) calculation of the weights as a function of a data spacing representative of the data distribution; 2) automatic elimination of detail at wavelengths smaller than twice the representative data spacing; 3) placement of bounds upon the grid spacing by the data spacing; and 4) setting of a fixed limit on the number of passes through the data to achieve rapid and sufficient convergence of the analyzed values to the observed ones. A mathematical analysis of the convergence properties of the Barnes technique is presented to support the validity of the latter constraint.
Despite these constraints, the interactive Barnes scheme remains versatile because it accepts limited inputs to the data and grid display areas, to the data and grid spacings, and to the rate of convergence of the analysis to the observations. Input parameter values are entered through a series of questions displayed on a computer video terminal and by manipulation of display function devices. The analyst immediately sees a plot of the data, the contoured grid values, superimposed in various colors if desired, and the effects of choice of analysis options. Examples of both meteorological and satellite data analyses are presented to demonstrate the objectivity, versatility and practicality of the interactive Barnes scheme.
Abstract
The impact of satellite-derived cloud motion vectors (CMVs) on analysts of winds measured by rawinsondes during the 1979 SESAME Experiment is studied in two case studies (10 April and 9 May 1979). Cloud motion vectors are both arbitrarily assigned and vertically interpolated to typical “low” levels of 825 mb and &sigma = 0.9 before being combined with the rawinsonde-measured winds at these levels. Magnitudes of vector differences between the combined winds and rawinsonde-measured winds are computed to show the effect of single-level assignment of CMVs on the rawinsonde-measured wind fields. The impact of the existence of horizontal and vertical gradients of wind and moisture on these results is also examined. In addition, divergence and relative vorticity fields are derived to determine whether the addition of CMVs increase the amount of useful information to the kinematic computations.
The results show that the standard method of arbitrarily assigning wind vectors to a “low-level” coordinate surface yields systematic differences between the rawinsonde-measured winds and combined wind fields. Arbitrary assignment of cloud motions to the 0.9 sigma surface produces smaller magnitudes of vector differences than assignment to the 825-mb pressure surface. Additionally, systematic differences in the wind occur near moisture discontinuities and in regions of horizontal and vertical wind shears. If forced to make arbitrary assignments, the use of the terrain-following sigma surface yields more consistent results than arbitrary assignment to a pressure surface in the lower troposphere. The differences between the combined and rawinsonde-measured wind fields are reduced by vertical interpolation to either a pressure or sigma surface. However, the accuracy of these interpolated fields depends to a large extent on the methods used to determine cloud-base levels and the vertical wind shear.
Abstract
The impact of satellite-derived cloud motion vectors (CMVs) on analysts of winds measured by rawinsondes during the 1979 SESAME Experiment is studied in two case studies (10 April and 9 May 1979). Cloud motion vectors are both arbitrarily assigned and vertically interpolated to typical “low” levels of 825 mb and &sigma = 0.9 before being combined with the rawinsonde-measured winds at these levels. Magnitudes of vector differences between the combined winds and rawinsonde-measured winds are computed to show the effect of single-level assignment of CMVs on the rawinsonde-measured wind fields. The impact of the existence of horizontal and vertical gradients of wind and moisture on these results is also examined. In addition, divergence and relative vorticity fields are derived to determine whether the addition of CMVs increase the amount of useful information to the kinematic computations.
The results show that the standard method of arbitrarily assigning wind vectors to a “low-level” coordinate surface yields systematic differences between the rawinsonde-measured winds and combined wind fields. Arbitrary assignment of cloud motions to the 0.9 sigma surface produces smaller magnitudes of vector differences than assignment to the 825-mb pressure surface. Additionally, systematic differences in the wind occur near moisture discontinuities and in regions of horizontal and vertical wind shears. If forced to make arbitrary assignments, the use of the terrain-following sigma surface yields more consistent results than arbitrary assignment to a pressure surface in the lower troposphere. The differences between the combined and rawinsonde-measured wind fields are reduced by vertical interpolation to either a pressure or sigma surface. However, the accuracy of these interpolated fields depends to a large extent on the methods used to determine cloud-base levels and the vertical wind shear.
Abstract
The effects of sensible heating and momentum mixing on the low-level structure and dynamics of a two-dimensional cold front are studied with a hydrostatic primitive equation model. Effects of inhomogeneous heating arising from a contrast in low-level cloud cover across the front are emphasized. The relative importance of grid resolution and the choice of method for parameterizing planetary boundary layer (PBL) processes in the model are also examined. Frontal updraft dynamics are studied in terms of the following inquiries: (a) the relative importance of turbulent momentum transport, differential sensible heating, and the reduction in static stability in the heated region ahead of the front; (b) the nature of the interaction between the adiabatic, semigeostrophic frontal circulation and the thermally forced circulation; and (c) possible roles played by dry symmetric instability and density current dynamics. The terms in the frontogenesis and divergence budget equations are computed to determine the relative roles played by the various physical and dynamical processes in generating the frontal secondary circulation system.
A strong, narrow updraft jet forms in the presence of uniform sensible beating across the front. Although the greatest impact on frontogenesis occurs as a response to the reduction in static stability resulting from uniform sensible heating, additional forcing results from the nonlinear interaction between the adiabatic frontal circulation and the thermally forced circulation arising from a cross-front gradient in heating (due to the introduction of an overcast low cloud deck behind the front). The relative importance of inhomogeneous heating, however, increases with the grid resolution and the use of a multilevel treatment in place of bulk mixed-layer PBL models.
Numerical experiments reveal that symmetric instability does not create the updraft jet, despite the development of negative potential vorticity ahead of the surface cold front. Highly unbalanced dynamics and a density current-like “feeder flow” behind the cold front are strongly indicated in the presence of sensible heating effects. Budget analyses show that the frontogenetical effect of sensible heating is only indirectly important through its strengthening of the confluence (convergence) field. The nonlinear and unbalanced ageostrophic vorticity terms in the divergence budget equation exert the strongest controls on the development of the updraft jet when sensible heating is nonuniform.
These results suggest that differential cloud cover across cold fronts may promote the development of frontal squall lines. Nonhydrostatic models that include explicit prognostic equations for microphysics and use improved parameterization of boundary layer fluxes in the presence of clouds are needed to more fully address this possibility.
Abstract
The effects of sensible heating and momentum mixing on the low-level structure and dynamics of a two-dimensional cold front are studied with a hydrostatic primitive equation model. Effects of inhomogeneous heating arising from a contrast in low-level cloud cover across the front are emphasized. The relative importance of grid resolution and the choice of method for parameterizing planetary boundary layer (PBL) processes in the model are also examined. Frontal updraft dynamics are studied in terms of the following inquiries: (a) the relative importance of turbulent momentum transport, differential sensible heating, and the reduction in static stability in the heated region ahead of the front; (b) the nature of the interaction between the adiabatic, semigeostrophic frontal circulation and the thermally forced circulation; and (c) possible roles played by dry symmetric instability and density current dynamics. The terms in the frontogenesis and divergence budget equations are computed to determine the relative roles played by the various physical and dynamical processes in generating the frontal secondary circulation system.
A strong, narrow updraft jet forms in the presence of uniform sensible beating across the front. Although the greatest impact on frontogenesis occurs as a response to the reduction in static stability resulting from uniform sensible heating, additional forcing results from the nonlinear interaction between the adiabatic frontal circulation and the thermally forced circulation arising from a cross-front gradient in heating (due to the introduction of an overcast low cloud deck behind the front). The relative importance of inhomogeneous heating, however, increases with the grid resolution and the use of a multilevel treatment in place of bulk mixed-layer PBL models.
Numerical experiments reveal that symmetric instability does not create the updraft jet, despite the development of negative potential vorticity ahead of the surface cold front. Highly unbalanced dynamics and a density current-like “feeder flow” behind the cold front are strongly indicated in the presence of sensible heating effects. Budget analyses show that the frontogenetical effect of sensible heating is only indirectly important through its strengthening of the confluence (convergence) field. The nonlinear and unbalanced ageostrophic vorticity terms in the divergence budget equation exert the strongest controls on the development of the updraft jet when sensible heating is nonuniform.
These results suggest that differential cloud cover across cold fronts may promote the development of frontal squall lines. Nonhydrostatic models that include explicit prognostic equations for microphysics and use improved parameterization of boundary layer fluxes in the presence of clouds are needed to more fully address this possibility.
Abstract
Nocturnal mesoscale convective systems (MCSs) frequently develop over the Great Plains in the presence of a nocturnal low-level jet (LLJ), which contributes to convective maintenance by providing a source of instability, convergence, and low-level vertical wind shear. Although these nocturnal MCSs often dissipate during the morning, many persist into the following afternoon despite the cessation of the LLJ with the onset of solar heating. The environmental factors enabling the postsunrise persistence of nocturnal convection are currently not well understood. A thorough investigation into the processes supporting the longevity and daytime persistence of an MCS was conducted using routine observations, RAP analyses, and a WRF-ARW simulation. Elevated nocturnal convection developed in response to enhanced frontogenesis, which quickly grew upscale into a severe quasi-linear convective system (QLCS). The western portion of this QLCS reorganized into a bow echo with a pronounced cold pool and ultimately an organized leading-line, trailing-stratiform MCS as it moved into an increasingly unstable environment. Differential advection resulting from the interaction of the nocturnal LLJ with the topography of west Texas established considerable heterogeneity in moisture, CAPE, and CIN, which influenced the structure and evolution of the MCS. An inland-advected moisture plume significantly increased near-surface CAPE during the nighttime over central Texas, while the environment over southeastern Texas abruptly destabilized following the commencement of surface heating and downward moisture transport. The unique topography of the southern plains and the close proximity to the Gulf of Mexico provided an environment conducive to the postsunrise persistence of the organized MCS.
Abstract
Nocturnal mesoscale convective systems (MCSs) frequently develop over the Great Plains in the presence of a nocturnal low-level jet (LLJ), which contributes to convective maintenance by providing a source of instability, convergence, and low-level vertical wind shear. Although these nocturnal MCSs often dissipate during the morning, many persist into the following afternoon despite the cessation of the LLJ with the onset of solar heating. The environmental factors enabling the postsunrise persistence of nocturnal convection are currently not well understood. A thorough investigation into the processes supporting the longevity and daytime persistence of an MCS was conducted using routine observations, RAP analyses, and a WRF-ARW simulation. Elevated nocturnal convection developed in response to enhanced frontogenesis, which quickly grew upscale into a severe quasi-linear convective system (QLCS). The western portion of this QLCS reorganized into a bow echo with a pronounced cold pool and ultimately an organized leading-line, trailing-stratiform MCS as it moved into an increasingly unstable environment. Differential advection resulting from the interaction of the nocturnal LLJ with the topography of west Texas established considerable heterogeneity in moisture, CAPE, and CIN, which influenced the structure and evolution of the MCS. An inland-advected moisture plume significantly increased near-surface CAPE during the nighttime over central Texas, while the environment over southeastern Texas abruptly destabilized following the commencement of surface heating and downward moisture transport. The unique topography of the southern plains and the close proximity to the Gulf of Mexico provided an environment conducive to the postsunrise persistence of the organized MCS.
Abstract
To assist in optimizing a mixed-physics ensemble for warm season mesoscale convective system rainfall forecasting, the impact of various physical schemes as well as their interactions on rainfall when different initializations were used has been investigated. For this purpose, high-resolution Weather Research and Forecasting (WRF) model simulations of eight International H2O Project events were performed. For each case, three different treatments of convection, three different microphysical schemes, and two different planetary boundary layer (PBL) schemes were used. All cases were initialized with both Local Analyses and Prediction System (LAPS) “hot” start analyses and 40-km Eta Model analyses. To evaluate the impacts of the variation of two different physical schemes and their interaction on the simulated rainfall under the two different initial conditions, the factor separation method was used. The sensitivity to the use of various physical schemes and their interactions was found to be dependent on the initialization dataset. Runs initialized with Eta analyses appeared to be influenced by the use of the Betts–Miller–Janjić scheme in that model’s assimilation system, which tended to reduce the WRF’s sensitivity to changes in the microphysical scheme compared with that present when LAPS analyses were used for initialization. In addition, differences in initialized thermodynamics resulted in changes in sensitivity to PBL and convective schemes. With both initialization datasets, the greatest sensitivity to the simulated rain rate was due to changes in the convective scheme. However, for rain volume, substantial sensitivity was present due to changes in both the physical parameterizations and the initial datasets.
Abstract
To assist in optimizing a mixed-physics ensemble for warm season mesoscale convective system rainfall forecasting, the impact of various physical schemes as well as their interactions on rainfall when different initializations were used has been investigated. For this purpose, high-resolution Weather Research and Forecasting (WRF) model simulations of eight International H2O Project events were performed. For each case, three different treatments of convection, three different microphysical schemes, and two different planetary boundary layer (PBL) schemes were used. All cases were initialized with both Local Analyses and Prediction System (LAPS) “hot” start analyses and 40-km Eta Model analyses. To evaluate the impacts of the variation of two different physical schemes and their interaction on the simulated rainfall under the two different initial conditions, the factor separation method was used. The sensitivity to the use of various physical schemes and their interactions was found to be dependent on the initialization dataset. Runs initialized with Eta analyses appeared to be influenced by the use of the Betts–Miller–Janjić scheme in that model’s assimilation system, which tended to reduce the WRF’s sensitivity to changes in the microphysical scheme compared with that present when LAPS analyses were used for initialization. In addition, differences in initialized thermodynamics resulted in changes in sensitivity to PBL and convective schemes. With both initialization datasets, the greatest sensitivity to the simulated rain rate was due to changes in the convective scheme. However, for rain volume, substantial sensitivity was present due to changes in both the physical parameterizations and the initial datasets.
An assessment of the value of data from the NOAA Profiler Network (NPN) on weather forecasting is presented. A series of experiments was conducted using the Rapid Update Cycle (RUC) model/assimilation system in which various data sources were denied in order to assess the relative importance of the profiler data for short-range wind forecasts. Average verification statistics from a 13-day cold-season test period indicate that the profiler data have a positive impact on short-range (3–12 h) forecasts over the RUC domain containing the lower 48 United States, which are strongest at the 3-h projection over a central U.S. subdomain that includes most of the profiler sites, as well as downwind of the profiler observations over the eastern United States. Overall, profiler data reduce wind forecast errors at all levels from 850 to 150 hPa, especially below 300 hPa where there are relatively few automated aircraft observations. At night when fewer commercial aircraft are flying, profiler data also contribute strongly to more accurate 3-h forecasts, including near-tropopause maximum wind levels. For the test period, the profiler data contributed up to 20%–30% (at 700 hPa) of the overall reduction of 3-h wind forecast error by all data sources combined. Inclusion of wind profiler data also reduced 3-h errors for height, relative humidity, and temperature by 5%-15%, averaged over different vertical levels. Time series and statistics from large-error events demonstrate that the impact of profiler data may be much larger in peak error situations.
Three data assimilation case studies from cold and warm seasons are presented that illustrate the value of the profiler observations for improving weather forecasts. The first case study indicates that inclusion of profiler data in the RUC model runs for the 3 May 1999 Oklahoma tornado outbreak improved model guidance of convective available potential energy (CAPE), 300-hPa wind, and precipitation in southwestern Oklahoma at the onset of the event. In the second case study, inclusion of profiler data led to better RUC precipitation forecasts associated with a severe snow and ice storm that occurred over the central plains of the United States in February 2001. A third case study describes the effect of profiler data for a tornado event in Oklahoma on 8 May 2003. Summaries of National Weather Service (NWS) forecaster use of profiler data in daily operations, although subjective, support the results from these case studies and the statistical forecast model impact study in the broad sense that profiler data contribute significantly to improved short-range forecasts over the central United States where these observations currently exist.
An assessment of the value of data from the NOAA Profiler Network (NPN) on weather forecasting is presented. A series of experiments was conducted using the Rapid Update Cycle (RUC) model/assimilation system in which various data sources were denied in order to assess the relative importance of the profiler data for short-range wind forecasts. Average verification statistics from a 13-day cold-season test period indicate that the profiler data have a positive impact on short-range (3–12 h) forecasts over the RUC domain containing the lower 48 United States, which are strongest at the 3-h projection over a central U.S. subdomain that includes most of the profiler sites, as well as downwind of the profiler observations over the eastern United States. Overall, profiler data reduce wind forecast errors at all levels from 850 to 150 hPa, especially below 300 hPa where there are relatively few automated aircraft observations. At night when fewer commercial aircraft are flying, profiler data also contribute strongly to more accurate 3-h forecasts, including near-tropopause maximum wind levels. For the test period, the profiler data contributed up to 20%–30% (at 700 hPa) of the overall reduction of 3-h wind forecast error by all data sources combined. Inclusion of wind profiler data also reduced 3-h errors for height, relative humidity, and temperature by 5%-15%, averaged over different vertical levels. Time series and statistics from large-error events demonstrate that the impact of profiler data may be much larger in peak error situations.
Three data assimilation case studies from cold and warm seasons are presented that illustrate the value of the profiler observations for improving weather forecasts. The first case study indicates that inclusion of profiler data in the RUC model runs for the 3 May 1999 Oklahoma tornado outbreak improved model guidance of convective available potential energy (CAPE), 300-hPa wind, and precipitation in southwestern Oklahoma at the onset of the event. In the second case study, inclusion of profiler data led to better RUC precipitation forecasts associated with a severe snow and ice storm that occurred over the central plains of the United States in February 2001. A third case study describes the effect of profiler data for a tornado event in Oklahoma on 8 May 2003. Summaries of National Weather Service (NWS) forecaster use of profiler data in daily operations, although subjective, support the results from these case studies and the statistical forecast model impact study in the broad sense that profiler data contribute significantly to improved short-range forecasts over the central United States where these observations currently exist.
Abstract
Mesoscale model simulations are performed in order to provide insight into the complex role of jet streak adjustments in establishing an environment favorable to the generation of gravity waves on 11–12 July 1981. This wave event was observed in unprecedented detail downstream of the Rocky Mountains in Montana during the Cooperative Convective Precipitation Experiment. The high-resolution model simulations employ a variety of terrain treatments in the absence of the complicating effects of precipitation physics in order to examine the complex interactions between orography and adiabatic geostrophic adjustment processes.
Results indicate that prior to gravity wave formation, a four-stage geostrophic adjustment process modified the structure of the mid- to upper-tropospheric jet streak by creating secondary mesoscale jet streaks (jetlets) to the southeast of the polar jet streak in proximity to the gravity wave generation region (WGR). During stage I, a strong rightward-directed ageostrophic flow in the right exit region of the polar jet streak (J1) developed over west-central Montana. This thermally indirect transverse secondary circulation resulted from inertial-advective adjustments wherein momentum was transported downstream and to the right of J1 as air parcels decelerated through the exit region.
During stage II, a highly unbalanced jetlet (J2) formed just northwest of the WGR in response to the inertial-advective forcing accompanying the ageostrophic circulation associated with J1. The mass field adjusted to this ageostrophic wind field. An adiabatic cooling and warming dipole resulting from this thermally indirect secondary circulation was the cause for frontogenesis and a rightward shift in the midtropospheric pressure gradients. Since this secondary circulation associated with J2 occurred above a dramatic vertical variation in the thermal wind, the vertical transport of potentially colder air from below was larger ahead of and to the right of J1, thus shifting the new jetlet (J2) well away from J1 as the mass field adjusted to the new wind field.
Stage III was established when the new mass field, which developed in association with J2 during stage II, set up a dynamically unbalanced circulation oriented primarily across the stream, and directly over the WGR. This new leftward-directed ageostrophic cross-stream flow (A) formed between jetlet J2 and the original exit region of the polar jet streak J1.
Finally, a midlevel mesoscale jetlet (J3) is simulated to have developed in stage IV over the WGR in response to the integrated mass flux divergence associated with both the stage II and III adjustment processes. This lower-level return branch circulation to jetlet J2 was further enhanced by velocity divergence accompanying the localized cross-stream ageostrophic wind maximum (A), which develops during stage III. The entire multistage geostrophic adjustment process required about 12 h to complete over a region encompassing approximately 400 km × 400 km.
Abstract
Mesoscale model simulations are performed in order to provide insight into the complex role of jet streak adjustments in establishing an environment favorable to the generation of gravity waves on 11–12 July 1981. This wave event was observed in unprecedented detail downstream of the Rocky Mountains in Montana during the Cooperative Convective Precipitation Experiment. The high-resolution model simulations employ a variety of terrain treatments in the absence of the complicating effects of precipitation physics in order to examine the complex interactions between orography and adiabatic geostrophic adjustment processes.
Results indicate that prior to gravity wave formation, a four-stage geostrophic adjustment process modified the structure of the mid- to upper-tropospheric jet streak by creating secondary mesoscale jet streaks (jetlets) to the southeast of the polar jet streak in proximity to the gravity wave generation region (WGR). During stage I, a strong rightward-directed ageostrophic flow in the right exit region of the polar jet streak (J1) developed over west-central Montana. This thermally indirect transverse secondary circulation resulted from inertial-advective adjustments wherein momentum was transported downstream and to the right of J1 as air parcels decelerated through the exit region.
During stage II, a highly unbalanced jetlet (J2) formed just northwest of the WGR in response to the inertial-advective forcing accompanying the ageostrophic circulation associated with J1. The mass field adjusted to this ageostrophic wind field. An adiabatic cooling and warming dipole resulting from this thermally indirect secondary circulation was the cause for frontogenesis and a rightward shift in the midtropospheric pressure gradients. Since this secondary circulation associated with J2 occurred above a dramatic vertical variation in the thermal wind, the vertical transport of potentially colder air from below was larger ahead of and to the right of J1, thus shifting the new jetlet (J2) well away from J1 as the mass field adjusted to the new wind field.
Stage III was established when the new mass field, which developed in association with J2 during stage II, set up a dynamically unbalanced circulation oriented primarily across the stream, and directly over the WGR. This new leftward-directed ageostrophic cross-stream flow (A) formed between jetlet J2 and the original exit region of the polar jet streak J1.
Finally, a midlevel mesoscale jetlet (J3) is simulated to have developed in stage IV over the WGR in response to the integrated mass flux divergence associated with both the stage II and III adjustment processes. This lower-level return branch circulation to jetlet J2 was further enhanced by velocity divergence accompanying the localized cross-stream ageostrophic wind maximum (A), which develops during stage III. The entire multistage geostrophic adjustment process required about 12 h to complete over a region encompassing approximately 400 km × 400 km.
Abstract
The most significant precipitation events in California occur during the winter and are often related to synoptic-scale storms from the Pacific Ocean. Because of the terrain characteristics and the fact that the urban and infrastructural expansion is concentrated in lower elevation areas of the California Central Valley, a high risk of flooding is usually associated with these events. In the present study, the area of interest was the American River basin (ARB). The main focus of the present study was to investigate methods for Quantitative Precipitation Forecast (QPF) improvement by estimating the impact that various microphysical schemes, planetary boundary layer (PBL) schemes, and initialization methods have on cold season precipitation, primarily orographically induced. For this purpose, 3-km grid spacing Weather Research and Forecasting (WRF) model simulations of four Hydrometeorological Test bed (HMT) events were used. For each event, four different microphysical schemes and two different PBL schemes were used. All runs were initialized with both a diabatic Local Analysis and Prediction System (LAPS) “hot” start and 40-km eta analyses.
To quantify the impact of physical schemes, their interactions, and initial conditions upon simulated rain volume, the factor separation methodology was used. The results showed that simulated rain volume was particularly affected by changes in microphysical schemes for both initializations. When the initialization was changed from the LAPS to the eta analysis, the change in the PBL scheme and corresponding synergistic terms (which corresponded to the interactions between different microphysical and PBL schemes) resulted in a statistically significant impact on rain volume. In addition, by combining model runs based on the knowledge about their impact on simulated rain volume obtained through the factor separation methodology, the bias in simulated rain volume was reduced.
Abstract
The most significant precipitation events in California occur during the winter and are often related to synoptic-scale storms from the Pacific Ocean. Because of the terrain characteristics and the fact that the urban and infrastructural expansion is concentrated in lower elevation areas of the California Central Valley, a high risk of flooding is usually associated with these events. In the present study, the area of interest was the American River basin (ARB). The main focus of the present study was to investigate methods for Quantitative Precipitation Forecast (QPF) improvement by estimating the impact that various microphysical schemes, planetary boundary layer (PBL) schemes, and initialization methods have on cold season precipitation, primarily orographically induced. For this purpose, 3-km grid spacing Weather Research and Forecasting (WRF) model simulations of four Hydrometeorological Test bed (HMT) events were used. For each event, four different microphysical schemes and two different PBL schemes were used. All runs were initialized with both a diabatic Local Analysis and Prediction System (LAPS) “hot” start and 40-km eta analyses.
To quantify the impact of physical schemes, their interactions, and initial conditions upon simulated rain volume, the factor separation methodology was used. The results showed that simulated rain volume was particularly affected by changes in microphysical schemes for both initializations. When the initialization was changed from the LAPS to the eta analysis, the change in the PBL scheme and corresponding synergistic terms (which corresponded to the interactions between different microphysical and PBL schemes) resulted in a statistically significant impact on rain volume. In addition, by combining model runs based on the knowledge about their impact on simulated rain volume obtained through the factor separation methodology, the bias in simulated rain volume was reduced.