Search Results

You are looking at 1 - 10 of 27 items for

  • Author or Editor: Kristopher M. Bedka x
  • Refine by Access: All Content x
Clear All Modify Search
Kristopher M. Bedka
and
John R. Mecikalski

Abstract

This study demonstrates methods to obtain high-density, satellite-derived atmospheric motion vectors (AMV) that contain both synoptic-scale and mesoscale flow components associated with and induced by cumuliform clouds through adjustments made to the University of Wisconsin—Madison Cooperative Institute for Meteorological Satellite Studies (UW-CIMSS) AMV processing algorithm. Operational AMV processing is geared toward the identification of synoptic-scale motions in geostrophic balance, which are useful in data assimilation applications. AMVs identified in the vicinity of deep convection are often rejected by quality-control checks used in the production of operational AMV datasets. Few users of these data have considered the use of AMVs with ageostrophic flow components, which often fail checks that assure both spatial coherence between neighboring AMVs and a strong correlation to an NWP-model first-guess wind field. The UW-CIMSS algorithm identifies coherent cloud and water vapor features (i.e., targets) that can be tracked within a sequence of geostationary visible (VIS) and infrared (IR) imagery. AMVs are derived through the combined use of satellite feature tracking and an NWP-model first guess. Reducing the impact of the NWP-model first guess on the final AMV field, in addition to adjusting the target selection and vector-editing schemes, is found to result in greater than a 20-fold increase in the number of AMVs obtained from the UW-CIMSS algorithm for one convective storm case examined here. Over a three-image sequence of Geostationary Operational Environmental Satellite (GOES)-12 VIS and IR data, 3516 AMVs are obtained, most of which contain flow components that deviate considerably from geostrophy. In comparison, 152 AMVs are derived when a tighter NWP-model constraint and no targeting adjustments were imposed, similar to settings used with operational AMV production algorithms. A detailed analysis reveals that many of these 3516 vectors contain low-level (100–70 kPa) convergent and midlevel (70–40 kPa) to upper-level (40–10 kPa) divergent motion components consistent with localized mesoscale flow patterns. The applicability of AMVs for estimating cloud-top cooling rates at the 1-km pixel scale is demonstrated with excellent correspondence to rates identified by a human expert.

Full access
Christopher S. Velden
and
Kristopher M. Bedka

Abstract

This study investigates the assignment of pressure heights to satellite-derived atmospheric motion vectors (AMVs), commonly known as cloud-drift and water vapor–motion winds. Large volumes of multispectral AMV datasets are compared with collocated rawinsonde wind profiles collected by the U.S. Department of Energy Atmospheric Radiation Measurement Program at three geographically disparate sites: the southern Great Plains, the North Slope of Alaska, and the tropical western Pacific Ocean. From a careful analysis of these comparisons, the authors estimate that mean AMV observation errors are ∼5–5.5 m s−1 and that vector height assignment is the dominant factor in AMV uncertainty, contributing up to 70% of the error. These comparisons also reveal that in most cases the RMS differences between matched AMVs and rawinsonde wind values are minimized if the rawinsonde values are averaged over specified layers. In other words, on average, the AMV values better correlate to a motion over a mean tropospheric layer rather than to a traditionally assigned discrete level. The height assignment behavioral characteristics are specifically identified according to AMV height (high cloud vs low cloud), type (spectral bands; clear vs cloudy), geolocation, height assignment method, and amount of environmental vertical wind shear present. The findings have potentially important implications for data assimilation of AMVs, and these are discussed.

Full access
Kristopher M. Bedka
and
Konstantin Khlopenkov

Abstract

Deep convective updrafts often penetrate through the surrounding cirrus anvil and into the lower stratosphere. Cross-tropopause transport of ice, water vapor, and chemicals occurs within these “overshooting tops” (OTs) along with a variety of hazardous weather conditions. OTs are readily apparent in satellite imagery, and, given the importance of OTs for weather and climate, a number of automated satellite-based detection methods have been developed. Some of these methods have proven to be relatively reliable, and their products are used in diverse Earth science applications. Nevertheless, analysis of these methods and feedback from product users indicate that use of fixed infrared temperature–based detection criteria often induces biases that can limit their utility for weather and climate analysis. This paper describes a new multispectral OT detection approach that improves upon those previously developed by minimizing use of fixed criteria and incorporating pattern recognition analyses to arrive at an OT probability product. The product is developed and validated using OT and non-OT anvil regions identified by a human within MODIS imagery. The product offered high skill for discriminating between OTs and anvils and matched 69% of human OT identifications for a particular probability threshold with a false-detection rate of 18%, outperforming previously existing methods. The false-detection rate drops to 1% when OT-induced texture detected within visible imagery is used to constrain the IR-based OT probability product. The OT probability product is also shown to improve severe-storm detection over the United States by 20% relative to the best existing method.

Full access
John R. Mecikalski
and
Kristopher M. Bedka

Abstract

This study identifies the precursor signals of convective initiation within sequences of 1-km-resolution visible (VIS) and 4–8-km infrared (IR) imagery from the Geostationary Operational Environmental Satellite (GOES) instrument. Convective initiation (CI) is defined for this study as the first detection of Weather Surveillance Radar-1988 Doppler (WSR-88D) reflectivities ≥35 dBZ produced by convective clouds. Results indicate that CI may be forecasted ∼30–45 min in advance through the monitoring of key IR fields for convective clouds. This is made possible by the coincident use of three components of GOES data: 1) a cumulus cloud “mask” at 1-km resolution using VIS and IR data, 2) satellite-derived atmospheric motion vectors (AMVs) for tracking individual cumulus clouds, and 3) IR brightness temperature (TB ) and multispectral band-differencing time trends. In effect, these techniques isolate only the cumulus convection in satellite imagery, track moving cumulus convection, and evaluate various IR cloud properties in time. Convective initiation is predicted by accumulating information within a satellite pixel that is attributed to the first occurrence of a ≥35 dBZ radar echo. Through the incorporation of satellite tracking of moving cumulus clouds, this work represents a significant advance in the use of routinely available GOES data for monitoring aspects of cumulus clouds important for nowcasting CI (0–1-h forecasts). Once cumulus cloud tracking is established, eight predictor fields based on Lagrangian trends in IR data are used to characterize cloud conditions consistent with CI. Cumulus cloud pixels for which ≥7 of the 8 CI indicators are satisfied are labeled as having high CI potential, assuming an extrapolation of past trends into the future. Comparison to future WSR-88D imagery then measures the method's predictive skill. Convective initiation predictability is demonstrated using several convective events—one during IHOP_2002—that occur over a variety of synoptic and mesoscale forcing regimes.

Full access
Sarah M. Griffin
,
Kristopher M. Bedka
, and
Christopher S. Velden

Abstract

Assigning accurate heights to convective cloud tops that penetrate into the upper troposphere–lower stratosphere (UTLS) region using infrared (IR) satellite imagery has been an unresolved issue for the satellite research community. The height assignment for the tops of optically thick clouds is typically accomplished by matching the observed IR brightness temperature (BT) with a collocated rawinsonde or numerical weather prediction (NWP) profile. However, “overshooting tops” (OTs) are typically colder (in BT) than any vertical level in the associated profile, leaving the height of these tops undetermined using this standard approach. A new method is described here for calculating the heights of convectively driven OTs using the characteristic temperature lapse rate of the cloud top as it ascends into the UTLS region. Using 108 MODIS-identified OT events that are directly observed by the CloudSat Cloud Profiling Radar (CPR), the MODIS-derived brightness temperature difference (BTD) between the OT and anvil regions can be defined. This BTD is combined with the CPR- and NWP-derived height difference between these two regions to determine the mean lapse rate, −7.34 K km−1, for the 108 events. The anvil height is typically well known, and an automated OT detection algorithm is used to derive BTD, so the lapse rate allows a height to be calculated for any detected OT. An empirical fit between MODIS and geostationary imager IR BT for OTs and anvil regions was performed to enable application of this method to coarser-spatial-resolution geostationary data. Validation indicates that ~75% (65%) of MODIS (geostationary) OT heights are within ±500 m of the coincident CPR-estimated heights.

Full access
Kristopher M. Bedka
,
Richard Dworak
,
Jason Brunner
, and
Wayne Feltz

Abstract

Two satellite infrared-based overshooting convective cloud-top (OT) detection methods have recently been described in the literature: 1) the 11-μm infrared window channel texture (IRW texture) method, which uses IRW channel brightness temperature (BT) spatial gradients and thresholds, and 2) the water vapor minus IRW BT difference (WV-IRW BTD). While both methods show good performance in published case study examples, it is important to quantitatively validate these methods relative to overshooting top events across the globe. Unfortunately, no overshooting top database currently exists that could be used in such study. This study examines National Aeronautics and Space Administration CloudSat Cloud Profiling Radar data to develop an OT detection validation database that is used to evaluate the IRW-texture and WV-IRW BTD OT detection methods. CloudSat data were manually examined over a 1.5-yr period to identify cases in which the cloud top penetrates above the tropopause height defined by a numerical weather prediction model and the surrounding cirrus anvil cloud top, producing 111 confirmed overshooting top events. When applied to Moderate Resolution Imaging Spectroradiometer (MODIS)-based Geostationary Operational Environmental Satellite-R Series (GOES-R) Advanced Baseline Imager proxy data, the IRW-texture (WV-IRW BTD) method offered a 76% (96%) probability of OT detection (POD) and 16% (81%) false-alarm ratio. Case study examples show that WV-IRW BTD > 0 K identifies much of the deep convective cloud top, while the IRW-texture method focuses only on regions with a spatial scale near that of commonly observed OTs. The POD decreases by 20% when IRW-texture is applied to current geostationary imager data, highlighting the importance of imager spatial resolution for observing and detecting OT regions.

Full access
Jason M. Apke
,
Yoo-Jeong Noh
, and
Kristopher Bedka

Abstract

This study introduces a validation technique for quantitative comparison of algorithms that retrieve winds from passive detection of cloud- and water vapor–drift motions, also known as atmospheric motion vectors (AMVs). The technique leverages airborne wind-profiling lidar data collected in tandem with 1-min refresh-rate geostationary satellite imagery. AMVs derived with different approaches are used with accompanying numerical weather prediction model data to estimate the full profiles of lidar-sampled winds, which enables ranking of feature tracking, quality control, and height-assignment accuracy and encourages mesoscale, multilayer, multiband wind retrieval solutions. The technique is used to compare the performance of two brightness motion, or “optical flow,” retrieval algorithms used within AMVs, 1) patch matching (PM; used within operational AMVs) and 2) an advanced variational optical flow (VOF) method enabled for most atmospheric motions by new-generation imagers. The VOF AMVs produce more accurate wind retrievals than the PM method within the benchmark in all imager bands explored. It is further shown that image regions with low texture and multilayer-cloud scenes in visible and infrared bands are tracked significantly better with the VOF approach, implying VOF produces representative AMVs where PM typically breaks down. It is also demonstrated that VOF AMVs have reduced accuracy where the brightness texture does not advect with the mean wind (e.g., gravity waves), where the image temporal noise exceeds the natural variability, and when the height assignment is poor. Finally, it is found that VOF AMVs have improved performance when using fine-temporal refresh-rate imagery, such as 1- versus 10-min data.

Open access
Cameron R. Homeyer
,
Joel D. McAuliffe
, and
Kristopher M. Bedka

Abstract

Expansive cirrus clouds present above the anvils of extratropical convection have been observed in satellite and aircraft-based imagery for several decades. Despite knowledge of their occurrence, the precise mechanisms and atmospheric conditions leading to their formation and maintenance are not entirely known. Here, the formation of these cirrus “plumes” is examined using a combination of satellite imagery, four-dimensional ground-based radar observations, assimilated atmospheric states from a state-of-the-art reanalysis, and idealized numerical simulations with explicitly resolved convection. Using data from 20 recent events (2013–present), it is found that convective cores of storms with above-anvil cirrus plumes reach altitudes 1–6 km above the tropopause. Thus, it is likely that these clouds represent the injection of cloud material into the lower stratosphere. Comparison of storms with above-anvil cirrus plumes and observed tropopause-penetrating convection without plumes reveals an association with large vector differences between the motion of a storm and the environmental wind in the upper troposphere and lower stratosphere (UTLS), suggesting that gravity wave breaking and/or stretching of the tropopause-penetrating cloud are/is more prevalent in plume-producing storms. A weak relationship is found between plume occurrence and the stability of the lower stratosphere (or tropopause structure), and no relationship is found with the duration of stratospheric penetration or stratospheric humidity. Idealized model simulations of tropopause-penetrating convection with small and large magnitudes of storm-relative wind in the UTLS are found to reproduce the observationally established storm-relative wind relationship and show that frequent gravity wave breaking is the primary mechanism responsible for plume formation.

Full access
Kristopher M. Bedka
,
Christopher S. Velden
,
Ralph A. Petersen
,
Wayne F. Feltz
, and
John R. Mecikalski

Abstract

Geostationary satellite-derived atmospheric motion vectors (AMVs) have been used over several decades in a wide variety of meteorological applications. The ever-increasing horizontal and vertical resolution of numerical weather prediction models puts a greater demand on satellite-derived wind products to monitor flow accurately at smaller scales and higher temporal resolution. The focus of this paper is to evaluate the accuracy and potential applications of a newly developed experimental mesoscale AMV product derived from Geostationary Operational Environmental Satellite (GOES) imagery. The mesoscale AMV product is derived through a variant on processing methods used within the University of Wisconsin—Madison Cooperative Institute for Meteorological Satellite Studies (UW-CIMSS) AMV algorithm and features a significant increase in vector density throughout the troposphere and lower stratosphere over current NOAA/National Environmental Satellite, Data, and Information Service (NESDIS) processing methods for GOES-12 Imager data. The primary objectives of this paper are to 1) highlight applications of experimental GOES mesoscale AMVs toward weather diagnosis and forecasting, 2) compare the coverage and accuracy of mesoscale AMVs with the NOAA/NESDIS operational AMV product, and 3) demonstrate the utility of 6-min NOAA Wind Profiler Network observations for satellite-derived AMV validation. Although the more conservative NOAA/NESDIS AMV product exhibits closer statistical agreement to rawinsonde and wind profiler observations than do the experimental mesoscale AMVs, a comparison of these two products for selected events shows that the mesoscale product better depicts the circulation center of a midlatitude cyclone, boundary layer confluence patterns, and a narrow low-level jet that is well correlated with subsequent severe thunderstorm development. Thus, while the individual experimental mesoscale AMVs may sacrifice some absolute accuracy, they show promise in providing greater temporal and spatial flow detail that can benefit diagnosis of upper-air flow patterns in near–real time. The results also show good agreement between 6-min wind profiler and rawinsonde observations within the 700–200-hPa layer, with larger differences in the stratosphere, near the mean top of the planetary boundary layer, and just above the earth’s surface. Despite these larger differences within select layers, the stability of the difference profile with height builds confidence in the use of 6-min, ∼404-MHz NOAA Wind Profiler Network observations to evaluate and better understand satellite AMV error characteristics.

Full access
Kristopher M. Bedka
,
John T. Allen
,
Heinz Jurgen Punge
,
Michael Kunz
, and
Denis Simanovic

ABSTRACT

A 10-yr geostationary (GEO) overshooting cloud-top (OT) detection database using Multifunction Transport Satellite (MTSAT) Japanese Advanced Meteorological Imager (JAMI) observations has been developed over the Australian region. GEO satellite imagers collect spatially and temporally detailed observations of deep convection, providing insight into the development and evolution of hazardous storms, particularly where surface observations of hazardous storms and deep convection are sparse and ground-based radar or lightning sensor networks are limited. Hazardous storms often produce one or more OTs that indicate the location of strong updrafts where weather hazards are typically concentrated, which can cause substantial impacts on the ground such as hail, damaging winds, tornadoes, and lightning and to aviation such as turbulence and in-flight icing. The 10-yr OT database produced using an automated OT detection algorithm is demonstrated for analysis of storm frequency, diurnally, spatially, and seasonally relative to known features such as the Australian monsoon, expected regions of hazardous storms along the southeastern coastal regions of southern Queensland and New South Wales, and the preferential extratropical cyclone track along the Indian Ocean and southern Australian coast. A filter based on atmospheric instability, deep-layer wind shear, and freezing level was used to identify OTs that could have produced hail. The filtered OT database is used to generate a hail frequency estimate that identifies a region extending from north of Brisbane to Sydney and the Goldfields–Esperance region of eastern Western Australia as the most hail-prone regions.

Full access