Applying Fog Forecasting Techniques Using AWIPS and the Internet
Robert E. Cox
National Weather Service, Wichita, Kansas
Abstract
1. Introduction
Given the level of performance in forecasting fog, one can deduce a great need exists for improvement. This paper will attempt to address this need by assessing the use of available tools to forecast radiation and advection fog utilizing AWIPS and the Internet. As an operational demonstration of these tools, four case studies will be examined. It is hoped that this study will result in improved understanding and utilization of tools by forecasters when predicting dense fog.2. Radiation Fog
2.1 Contributing Factors
S= K dT/dZ |
(1)
|
where K term is the thermal conductivity of the soil and dT/dZ is the change in temperature with respect to height. Thermal conductivity is defined as the amount of heat that the soil will conduct under a specified temperature gradient ( Rosenburg et al. 1983). Thermal conductivity is very difficult to measure in real-time, because it requires a steady-state situation. However, it does play a vital role in the heat flux which is important when forecasting radiation fog. More specifically, low thermal conductivity is desired as to limit heat flow to the surface as much as possible (UCAR 2001).
There are several factors that affect thermal conductivity. They include the following; porosity, moisture content and organic matter content of the soil (Rosenburg et al. 1983). Clay and silt loams are very common soils found in the Great Plains. These soils have the ability to hold more moisture than sandy soils, because they are more porous. Therefore, clay and silt loams have a lower thermal conductivity. At similar moisture content, the clay and silt loams also will have lower thermal conductivities than sand (Figure 1) (Rosenburg et al. 1983).
As a result, clay and silt loams will usually have lower heat fluxes towards the surface, which is ideal in radiation fog. However, if the soil is dry and compacted this increases the thermal conductivity, because it decreases the volume of the insulating pore spaces (Rosenburg et al. 1983).Kansas Agricultural reports in AWIPS usually report the mean four-inch soil temperature for the day (Figure 2). However, this alone cannot give a good estimate of what the actual soil temperature is during the times of potential fog. Forecasters need to take into account how the thermal energy and moisture is being distributed in the soil. The Meteorological Services of Canada 1978 desribes the transfer of moisture through soil layers as follows: “In the presence of a temperature gradient, water tends to be distilled from warmer regions to condense in cooler regions. Because of the heat storage and relative constant temperatures of the deeper soil layers, the moisture flow due to temperature gradients is usually downward in summer and upward in winter.” The same principle also applies to temperature. The heat flows from warmer to colder regions (Rosenburg et al. 1983).
During the summer, one tends have larger ranges of soil temperature, especially near the surface due to increased insolation. This is shown in Figure 3. This means that most of the soil heat flux is being directed downward into the soil, even through much of the night (Rosenburg et al. 1983). During the summer, moisture flux is also being directed downward into the soil which will increase the thermal conductivity, by drying out the soil. Therefore, one will likely not see the radiation fog development during the summer, unless the soil is very moist with high evapotranspiration. Evapotranspiration can play a role in radiation fog development in the Central Plains during the spring and fall during the early morning hours by adding moisture to the boundary layer, but it does not play much of a role during the winter since most of the grass and plants are dormant.
2.2 Forecast Tools
2.2.a UPS Technique
|
(2)
|
Where:
Tb = boundary layer temperature ( °C)
Tsfc= shelter temperature forecast ( °C)
u= boundary layer wind speed (knots)
Their method utilized FOUS ETA/NGM data which are not readily available for many TAF sites across the Plains. However, AWIPS has the capability to compute the MRi utilizing the model data. The NWS Weather Forecast Office (WFO) in Wichita, Kansas has the MRi in their AWIPS volume browser 2which is computed utilizing the lowest sigma layer of the boundary layer data and the extrapolated 2 m temperatures. The forecaster can then create images of the MRi to display where the best decoupling may take place. The forecaster must use the 2 m temperature data with caution and compare it with real-time observations.
BUFKIT version 4 (Figure 4) capabilities incorporate utilization of the UPS technique which includes diagnosis of the hydrolapse rate, vertical humidity flux, crossover temperature, ground temperature and MRi. BUFKIT is the most flexible of the programs available as the user may alter the expected crossover temperature to explore different outcomes created by the software. This technique works very well, especially in conditions when there is little or no pattern change during the next 12 to 18 hours. Of course, this is all based on model data, so the forecaster will have to take into account all of the raw data which includes surface observations, profilers, RAOBS (Radiosonde Observations, AMDAR (Aircraft Meteorological Data Relay), MDCRS (Meteorological Data Collection and Reporting System), and ACARS (ARINC Communications Addressing and Reporting System).
2.2.b Soil Temperature and Moisture Sources
The Oklahoma Mesonet (Figure 5) is a great real-time source to view the hourly soil temperatures, but more of these sites are needed nationwide to examine soil heat fluxes.
The RUC20 (Rapid Update Cycle 20km Resolution) model is another helpful tool in determining how much soil moisture is present (Figure 6).3There are 6 levels in the RUC land-surface model, extending down to 3 m deep, but the field shown is for the top 2 cm of soil only. Therefore, this field responds quickly to recent precipitation or surface drying and may not be indicative of deep soil moisture which may move upwards during the winter. The variable displayed is the soil volumetric moisture content, the ratio of water volume to total volume in the soil. Values of 0.25 or so are relatively high values. (NOAA/Earth System Research Laboratory RUC development group). Of course, radar and rain gage data must be examined to determine areas of recent rainfall.
2.2.c Other Tools
3.0 Advection Fog
3.1 Contributing Factors
3.2 Advection Fog Forecast Tools
In order to predict when this fog will occur, use of real-time data is paramount. This includes examining surface and boundary layer data from observations, profilers, ACARS, RAOBS, and model soundings. On a surface chart, the advection fog often develops along and just north of a warm/stationary front with a tight thermal/moisture gradient. This is also usually east of a developing surface low pressure system. (FMI 2005) Moisture transport vectors and ThetaE advection in the boundary layer in AWIPS are also good tools to ascertain the degree of moisture advection. Moisture Transport Vectors are defined as: the product of the horizontal velocity vector, V, and the mixing ratio, q. Units are gm-m/kg-s; values typically range from 50-250, depending upon the level and the season.
The NOAA Profiler Network (NPN) also has some data that is not viewable in AWIPS. One of the most interesting pieces of information which can be used for advection fog forecasting is the signal power display.4 Meteorological features such as moisture advection and cloud layers are often visible in the signal power displays. See Figure 7. Unfortunately, this data is not as useful for fog prediction when the moisture layer is very shallow.
4.0 Case Studies
4.1 December 20-21, 2005 Stratus-Build-Down Radiation Fog Event
In this event, the NAM (North American Mesoscale)/RUC models struggled with clearing out the low level moisture too quickly, which made forecasting this event quite a challenge. Afternoon crossover values ranged from the lower to mid 20s F, but the models were projecting these values to lower during the course of the night mainly across Central and South Central Kansas. This could have led one to believe that strong subsidence from the surface ridge would create dew deposition. Thus, the fog may end up being delayed or patchy. Across Southeast Kansas the models were keeping crossover values about the same which meant the fog could be a problem in this region if they decouple. Figure 8 displays there was some wet snow on the ground. Figure 9 shows the lowest ceilings over Central and South Central Kansas while Southeast Kansas had ceilings mainly in the MVFR category during the afternoon. During the late afternoon and early evening, there was some clearing across the western sections of the forecast area, but the lack of isentropic downglide or rising theta surfaces displayed on the 270k theta surface really slowed down the eastward progression of the clearing line. See Figures 10 and 11.
The NAM/RUC were both showing the clearing line arriving at Wichita by midnight according to the boundary layer RH and soundings. This indeed was not the case, as one can see in Figures 12 and 13.
Around 8 pm, the forecasters made some phone calls to the western fringes of the clearing line. At that time, dense fog was forming in those areas, and the crossover values were being met. The satellite fog imagery displayed this fog layer fairly well as seen in Figures 14 and 15. This was also where the cloud layer was the thinnest which allowed for ideal fog top cooling, downward movement of heat and moisture and good decoupling across this region (Figure 16). As the night wore on, the clearing line made very slow progress to the east. It finally made it to Salina and Hutchinson around 09z, but not before the widespread dense fog developed and persisted for several hours. The dense fog mostly formed where the ceilings started out the lowest in Central Kansas. During the evening, the ceilings were mainly below 500ft and continued to lower until the dense fog settled in. South Central Kansas started out with ceilings around 800ft and then the dense fog developed between 06 and 08z when boundary layer winds started to decouple. Southeast Kansas never did receive dense fog, possibly due to thicker cloud decks which did not allow for good fog top cooling, and lack of decoupling until daybreak. However, the ceilings did lower about a 1,000 ft during the night.
4.2 January 3, 2006 Radiation Fog Event
The second case took place on January 3, 2006 in Central Missouri. This event was very similar to the case of December 20-21st, but in this case the models did a better job with the low level moisture (Figure 17). The stratus deck came to a halt just east of Springfield (Figure 18), due to lack of isentropic downglide (Figure 19). This stratus deck then began to build to the west shortly after midnight with 700ft ceilings at Springfield (Figure 20). The formation of the stratus deck at this layer could have been a result of a dewpoint inversion formed from the dew deposition process in the lower part of boundary layer (Figure 21). These ceilings continued to lower through daybreak until dense fog formed across the back edge of stratus (Figure 22). This is where the cloud deck was the thinnest which allowed for ideal fog top cooling, downward moving heat and moisture, and decoupling. The only difference between this fog event and the fog event that took place on December 20-21st was that the ground was dry and there could have been more subsidence from the surface ridge. This could have offset the timing of the fog due to the soil temperature flux upward and dew deposition. However, the fog top cooling and good decoupling towards sunrise provided the necessary ingredients to form dense fog over central Missouri between 12 and 16z. Further west, there was more isentropic downglide (Figure 23) and possibly higher deposition rates which could have depleted the boundary layer moisture around 12z.
4.3 December 3, 2001 Advection Fog Event
The 04z surface analysis (Figure 24) on December 3rd, 2001 depicted low pressure over Southeast Colorado. From this low pressure system, a trough extended northeast from the Oklahoma Panhandle through Wichita and extending northeast. There was a tight moisture gradient along this boundary mainly across South Central Kansas with temperatures in the lower 50s and dewpoints around 50 south of the boundary. North of the boundary, temperatures were generally in the mid 30s, dewpoints in the mid 20s with good radiational cooling in place (Figure 24). The 08z and 13z surface maps (Figures 25 and 26), indicated the low level moisture advecting northward into Northern Kansas. This low level moisture was moving up a gradual slope from Northern Oklahoma to Central Kansas (Figure 27). Water Vapor imagery also indicated a cirrus shield stretching from the Central Rockies to Eastern Kansas. The forecaster on duty expected the cirrus shield to suppress the fog potential, but in retrospect it did not have much affect on it. The fog initiated in Southern Kansas at Wichita between 08z-09z and areas to the north towards 12z. The visibility dropped to one quarter of a mile within 30 minutes across the area (Figure 28). The fog was not a long-lived event with most of the fog dissipating within a couple of hours of its development. The observations also show a 3 degree C temperature climb during this timeframe. The fog that moved into these areas all had one thing in common. The dewpoints that advected in were several degrees higher than the ambient air temperature. The fog also was denser over the areas that had good radiational cooling prior to the advection of higher dewpoints. It appears like the dense fog did not develop over Southeast Kansas due to lack of radiational cooling and higher low level moisture in place.
The BUFKIT data was very useful during this event. The ETA showed the low level moisture advecting into the lowest 25 mb of the sounding after 06z, in Wichita, Hutchinson and Salina. See Figures 29 and 30. The ETA also showed Russell receiving the low level moisture, but after examining the observations the low level moisture never did arrive at that point. See Figures 31 and 32. Soil temperatures were also about 5 deg c cooler than surface temperatures which could have played a role in strengthening the inversion and limiting the heat flux upwards. See Figure 33.
4.4 January 12, 2005 Advection Fog Event
At 06z on January 12, 2005, a surface warm frontal boundary extended from the Oklahoma Panhandle to Northwest Arkansas (Figure 34). This boundary then extended northeast through St. Louis into Central Ohio. Temperatures north of this boundary were primarily in the 20s and 30s with temperatures ranging from the 40s to the 60s south of the boundary (Figure 34). The visibility was around ¼ of a mile in Central Oklahoma near the frontal boundary (Figure 35). By 15z, the frontal boundary had moved north into Wichita and extended through St. Louis to near the Michigan and Northern Indiana border (Figures 36 and 37). The surface visibility and ceilings dropped to the LIFR category with the visibility generally one quarter of a mile or less, especially along and north of the front (Figures 38 and 39). You will also notice the temperatures and dewpoints showed a dramatic rise when the dense fog started to develop and continued to rise until the inversion mixed out late in the day. There was good moisture transport in the boundary layer (Figure 40). The profilers could have also a benefited in this case. These profilers displayed an increase in signal power in the lowest 1km mainly after 09z which means there was an increase in low level moisture (Figures 41 and 42). In addition, the MRi was greater than 0.01 across the region, which means that the boundary layer was well mixed (Figure 43). The RUC soundings displayed a strong inversion over Northern Indiana and Southeast Kansas(Figures 44 and 45). One interesting thing to note in this case was that the dense fog developed over both snow covered areas and non-snow covered areas (Figure 46). In the non-snow covered areas the ground temperatures were also slightly warmer than the surface temperatures (Figure 47). Therefore, this event was primarily driven by advection.
5.0 Conclusions
________
1Issue is taken here by the author given the number of advection cases, some of which will be demonstrated, where surface temperature actually exhibit warming.
2 The Volume Browser provides access to numerical models and other gridded data sources, such as RAOB, BUFR, and Profiler data. Through the Browser interface, you can choose the data source(s), field(s), plane(s), and point(s), and generate a customized list of model graphics or images for display. Forecast Systems Laboratory, “AWIPS D2D User Manual”, http://www-md.fsl.noaa.gov/~cheatwoo/AWIPS/OB4-UM/3_1.htm (September 2004)
3 Equivalent Potential Temperature - The temperature a parcel of air would have if a) it was lifted until it became saturated, b) all water vapor was condensed out, and c) it was returned adiabatically (i.e., without transfer of heat or mass) to a pressure of 1000 millibars.
4The signal power is a measure of the amount of backscattered power received from the atmosphere. High signal power values (greater the 60dB) are typically associated with high moisture content or the presence of particles, while low power values (less than 40dB) usually indicate a dry or stable atmosphere.”
Acknowledgements:
The author would like to thank Ken Cook and Jim Caruso, both at NWS Wichita, and Jeff Manion (Central Region Headquarters Kansas City, Missouri), for their review and comments on this paper. The author would also like to thank Jim Caruso for his work on a case study that was conducted on December 3, 2001 which was incorporated into this study.
References:
Baker, R., Cramer, J., Peters, J., 2002: Radiation Fog: UPS Airlines Conceptual Models and Forecast Methods. Available on the web at:
http://www.meted.ucar.edu/dlac/website/resources/upspaper.pdf.
Finnish Meteorological Institute: Manual of Synoptic Satellite Meteorology Conceptual Models and Case Studies. Available on the web at: http://www.zamg.ac.at/docu/Manual/SatManu/main.htm.
Goodwin, L. “Analysis of Weather-Related Crashes on U.S. Highways. Miretek Systems Inc.”, 2002.
Meteorological Services of Canada 1978: Manual of Standards and Obtaining Soil Temperature Data. Available on the web at: http://www.msc-smc.ec.gc.ca/msb/manuals/soiltemp/introduction_e.html.
Pearson, D. C. 2002: “VFR Flight Not Recommended” A Study of Weather-Related Fatal Aviation Accidents. Technical Attachment SR SSD 2002-18. Available on the web at:
http://www.srh.noaa.gov/topics/attach/html/ssd02-18.htm.
Rosenburg, N.J., Blad, B.L., Verma, S.B., 1983: Microclimate: The Biological Environment. 2, 94-115.
Tardiff, R., 2001: Dew and Frost Deposition and their Role in the Regulation of Radiation Fog Onset: A Numerical Sensitivity Study. Available on the web at: http://www.rap.ucar.edu/staff/tardif/Documents/CUprojects/ATOC5225/dew_frost_fog.pdf.
University Cooperation for Atmospheric Research, 2001: Forecasting Radiation Fog.
Available on the web at: http://meted.ucar.edu/mesoprim/mpradfog/frameset.htm.
Forecast Systems Laboratory, 2004: AWIPS D2D User Manual. Available on the web at: http://www-md.fsl.noaa.gov/~cheatwoo/AWIPS/OB4-UM/3_1.htm.
Other websites used for research:
AMS Glossary:
http://amsglossary.allenpress.com/glossary/browse?s=h&p=28
Atmospheric Radiation Measurement:
National Weather Service, NOAA Profiler Network:
http://www.profiler.noaa.gov/npn/aboutProfilerData.jsp
Oklahoma Mesonet:
http://www.mesonet.ou.edu/public/current.html
BUFKIT Website:
http://www.wbuf.noaa.gov/bufkit/bufkit.html