Predicting Localized Flooding during the 2015 North Alabama Christmas Flood Event Using WRF-Hydro Simulated Streamflow
No abstract available
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
No abstract available
Spatial correlations of two of the main rain drop-size distribution (DSD) parameters - namely the median-volume diameter (Do) and the normalized intercept parameter (Nw) - as well as rainfall rate (R) are determined from polarimetric radar measurements, with added information from 2D video disdrometer (2DVD) data. Two cases have been considered, (i) a widespread, long-duration rain event in Huntsville, Alabama, and (ii) an event with localized intense rain-cells within a convection line which occurred during the MC3E campaign. For the first case, data from a C-band polarimetric radar (ARMOR) were utilized, with two 2DVDs acting as ground-truth , both being located at the same site 15 km from the radar. The radar was operated in a special near-dwelling mode over the 2DVDs. In the second case, data from an S-band polarimetric radar (NPOL) data were utilized, with at least five 2DVDs located between 20 and 30 km from the radar. In both rain event cases, comparisons of Do, log10(Nw) and R were made between radar derived estimates and 2DVD-based measurements, and were found to be in good agreement, and in both cases, the radar data were subsequently used to determine the spatial correlations For the first case, the spatial decorrelation distance was found to be smallest for R (4.5 km), and largest fo Do (8.2 km). For log10(Nw) it was 7.2 km (Fig. 1). For the second case, the corresponding decorrelation distances were somewhat smaller but had a directional dependence. In Fig. 2, we show an example of Do comparisons between NPOL based estimates and 1-minute DSD based estimates from one of the five 2DVDs.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.
A fully calibrated low profile 2D video disdrometer (2DVD) has been recording many different rainfall events in Northern Alabama (USA) since June 2007. An earlier publication reported drop shapes and axis ratio distributions determined for some of the events. For one of the cases examined, a noticeable shift in the 3.5 - 3.75 mm drop axis ratio distribution was noted. In this paper, we extend the earlier work by separating the 2DVD measurements into stratiform and convective rain. The separation is made possible by using the minute-by-minute drop size distribution (DSD) measured by the 2DVD. The 1-minute DSDs are fitted to a gamma distribution, and using a simple indexing technique which involves two of the fitted parameters, periods of convective and stratiform rain are separated for a given event. The output of the DSD indexing technique is qualitatively confirmed by comparing with simultaneous time series observations from a co-located UHF profiler which continuously records height profiles of reflectivity, Doppler mean and spectral width, all of which enable the identification of bright-band periods and, furthermore, periods of moderate and deep convection. Excellent consistency is found between the output of the DSD-based separation method and the profiler observations. Next, we utilize the output of DSD index-based separation method to flag the periods of severe convection for a given event. Drop axis ratios during the flagged periods are derived and compared with those during stratiform rain periods. Five cases have been considered. Axis ratio distributions do not show appreciable differences between stratiform and convective periods for four of the cases. The fifth case (the same case as reported earlier) shows a shift in the 3.5 - 3.75 mm drop axis ratios during a prolonged period of convection. The contoured shapes for these drops determined from the 2DVD camera data indicate the possibility of non-axisymmetric oscillations, compared with the contoured images for other events which fit well to our reference drop shapes. For all of above cases, observations from a C-band polarimetric radar - situated 15 km away are examined. The variations between the co-polar radar reflectivity and the differential reflectivity as well as the specific differential phase are compared with the 2DVD data based scattering calculations for the 5 events. The implications will be discussed.
An impact-type Joss-Waldvogel disdrometer (JWD), a two-dimensional video disdrometer (2DVD), and a laser optical OTT Particle Size and Velocity (PARSIVEL) disdrometer (PD) were used to measure the raindrop size distribution (DSD) over a 6-month period in Huntsville, Alabama. Comparisons indicate event rain totals for all three disdrometers that were in reasonable agreement with a reference rain gauge. In a relative sense, hourly composite DSDs revealed that the JWD was more sensitive to small drops (,1 mm), while the PD appeared to severely underestimate small drops less than 0.76mm in diameter. The JWD and 2DVD measured comparable number concentrations of midsize drops (1-3mm) and large drops (3-5 mm), while the PD tended to measure relatively higher drop concentrations at sizes larger than 2.44mm in diameter. This concentration disparity tended to occur when hourly rain rates and drop counts exceeded 2.5mm/h and 400/min, respectively. Based on interactions with the PD manufacturer, the partially inhomogeneous laser beam is considered the cause of the PD drop count overestimation. PD drop fall speeds followed the expected terminal fall speed relationship quite well, while the 2DVD occasionally measured slower drops for diameters larger than 2.4mm, coinciding with events where wind speeds were greater than 4m/s. The underestimation of small drops by the PD had a pronounced effect on the intercept and shape of parameters of gamma-fitted DSDs, while the overestimation of midsize and larger drops resulted in higher mean values for PD integral rain parameters
This study evaluates the effectiveness of spectrophotometric measurements from space in revealing properties of lightning flash. The multicolor optical waveform data obtained by FORMOSAT-2/Imager of Sprites and Upper Atmospheric Lightning (ISUAL) were analyzed in relation to National Lightning Detection Network (NLDN), North Alabama Lightning Mapping Array (LMA). As of July 2011, we found six lightning events which were observed by ISUAL and North Alabama LMA. In two of these events, NLDN showed clear positive cloud-to-ground (CG) discharges with peak current of +139.9 kA and +41.6 kA and, around that time, LMA showed continuous intra-cloud (IC) leader activities at 4-6 km altitudes. ISUAL also observed consistent optical waveforms of the IC and CG components and, interestingly, it was found that the blue/red spectral ratio clearly decreased by a factor of 1.5-2.5 at the time of CG discharges. Other four lightning events in which NLDN did not detect any CG discharges were also investigated, but such a feature was not found in any of these cases. These results suggest that the optical color of CG component is more reddish than that of IC component and we explain this as a result of more effective Rayleigh scattering in blue light emissions coming from lower-altitude light source. This finding suggests that spectral measurements could be a new useful technique to characterize ICs and CGs from space. In this talk, we will also present a result from lightning statistical analysis of ISUAL spectrophotometric data and ULF magnetic data.
The upgrade of the National Weather Service (NWS) network of S ]band dual-polarization radars is currently underway, and the incorporation of polarimetric information into the real ]time forecasting process will enhance the forecaster fs ability to assess thunderstorms and their near ]storm environments. Recent research has suggested that the combination of polarimetric variables differential reflectivity (ZDR) and specific differential phase (KDP) can be useful in the assessment of low level wind shear within a thunderstorm. In an environment with strong low ]level veering of the wind, ZDR values will be largest along the right inflow edge of the thunderstorm near a large gradient in horizontal reflectivity (indicative of large raindrops falling with a relative lack of smaller drops), and take the shape of an arc. Meanwhile, KDP values, which are proportional to liquid water content and indicative of a large number of smaller drops, are maximized deeper into the forward flank precipitation shield than the ZDR arc as the smaller drops are being advected further from the updraft core by the low level winds than the larger raindrops. Using findings from previous work, three severe weather events that occurred in North Alabama were examined in order to assess the utility of these signatures in determining the potential for tornadic activity. The first case is from October 26, 2010, where a large number of storms indicated tornadic potential from a standard reflectivity and velocity analysis but very few storms actually produced tornadoes. The second event is from February 28, 2011, where tornadic storms were present early on in the event, but as the day progressed, the tornado threat transitioned to a high wind threat. The third case is from April 27, 2011, where multiple rounds of tornadic storms ransacked the Tennessee Valley. This event provides a dataset including multiple modes of tornadic development, including QLCS and supercell structures. The overarching goal of examining these three events is to compare dual ]polarization features from this larger dataset to previous work and to determine if these signatures can be a useful indication of the potential for tornadic activity associated with the amount of low ]level wind shear in the near ]storm environment.
In the past year, the primary objectives were to show the usefulness of total lightning as compared to traditional cloud-to-ground (CG) networks, test the lightning jump algorithm configurations in other regions of the country, increase the number of thunderstorms within our thunderstorm database, and to pinpoint environments that could prove difficult for any lightning jump configuration. A total of 561 thunderstorms have been examined in the past year (409 non-severe, 152 severe) from four regions of the country (North Alabama, Washington D.C., High Plains of CO/KS, and Oklahoma). Results continue to indicate that the 2 lightning jump algorithm configuration holds the most promise in terms of prospective operational lightning jump algorithms, with a probability of detection (POD) at 81%, a false alarm rate (FAR) of 45%, a critical success index (CSI) of 49% and a Heidke Skill Score (HSS) of 0.66. The second best performing algorithm configuration was the Threshold 4 algorithm, which had a POD of 72%, FAR of 51%, a CSI of 41% and an HSS of 0.58. Because a more complex algorithm configuration shows the most promise in terms of prospective operational lightning jump algorithms, accurate thunderstorm cell tracking work must be undertaken to track lightning trends on an individual thunderstorm basis over time. While these numbers for the 2 configuration are impressive, the algorithm does have its weaknesses. Specifically, low-topped and tropical cyclone thunderstorm environments are present issues for the 2 lightning jump algorithm, because of the suppressed vertical depth impact on overall flash counts (i.e., a relative dearth in lightning). For example, in a sample of 120 thunderstorms from northern Alabama that contained 72 missed events by the 2 algorithm 36% of the misses were associated with these two environments (17 storms).
This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation events. Soil samples were collected across four sampling events along terrestrial-aquatic gradients at three wetland sites located within the Tanglewood Biological Station in Alabama from April 2024 to June 2025. The contents in this data package include soil in-situ nitrogen mineralization rates (measurements of net nitrification, net ammonification, and net mineralization), nutrient stocks (total carbon, total nitrogen, and organic matter), stable isotopes (carbon and nitrogen), and water volumetric measurements (water-filled pore space). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024), related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025), and related surface water sediment chemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3377325 (Molina Serpas et al., 2026). In addition to this readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) net nitrification rate, (2) net ammonification rate, (3) areal net mineralization rate, (4) percent organic matter, (5) water-filled pore space, (6) total carbon content, (7) total nitrogen content, (8) stable carbon isotope (delta carbon-13), and (9) stable nitrogen isotope (delta nitrogen-15), and (10) methods codes. All files are .csv or .pdf.
This dataset supports a broader study investigating wetland hydrologic and biogeochemical responses to inundation disturbances. Bimonthly surface water and sediment sampling events were conducted at nine wetland sites situated within the Tanglewood Biological Station in Alabama from April 2023 to February 2024. The contents included in the data package include surface water nitrogen (nitrogen oxides and ammonium) and sediment potential nitrate reduction rates (measured as potential denitrification and dissimilatory nitrate reduction to ammonium processing), nutrient stocks (total carbon, total nitrogen, and organic matter), and stable isotopes (carbon and nitrogen). Water level data related to each wetland location can be found at https://data.ess-dive.lbl.gov/view/doi:10.15485/2530253 (Kirker et al., 2024) and related water geochemistry data can be found at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/3001967 (Forbes et al., 2025). In addition to a readme, this data package also includes a file-level metadata (FLMD) file that describes each file and a data dictionary (DD) that describes all column/row headers and variable definitions. This dataset is comprised of (1) file-level metadata; (2) data dictionary; (3) field metadata and international generic sample numbers (IGSNs); (4) readme; (5) the field protocol; and (6) a subfolder with sample data. The sample data subfolder contains (1) sediment potential denitrification rate, (2) sediment potential dissimilatory nitrate reduction to ammonium (DNRA) rate, (3) sediment total carbon and nitrogen content, (4) sediment stable isotopes (delta nitrogen-15 and delta carbon-13), (5) sediment percent organic matter, (6) surface water nitrous oxides, (7) surface water ammonium, and (8) methods codes. All files are .csv or .pdf.
New particle formation (NPF) and growth govern cloud condensation nuclei (CCN) concentrations in many regions. The mechanisms governing the nucleation of molecular clusters vary substantially in different regions of the atmosphere. Additionally, the growth of these clusters from ~2 to 20 nm sizes is often governed by the availability of extremely low volatility organic vapours (ELVOCs). While the pathways to ELVOC formation from the oxidation of biogenic monoterpenes with ozone is better understood, the chemical and mechanistic pathways for ELVOC formation from oxidation of anthropogenic organics are not well understood. We integrate measurements and three-dimensional regional model simulations with the Weather Research and Forecasting Model coupled to chemistry (WRF-Chem) to understand the processes governing new particle formation and growth and secondary organic aerosol (SOA) formation during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems (HI-SCALE) field campaign at the Southern Great Plains (SGP) observatory in Oklahoma, and contrast it with a site within the Bankhead National Forest (BNF), Alabama in Southeast USA, where 5-year long measurements will begin in 2024. Simulations show that nucleation rates are at least an order of magnitude higher at SGP compared to BNF during the springtime days (April 28 and May 14, 2016), largely due to lower H2SO4 concentrations at BNF, which are needed for nucleation. In addition, the larger CS at BNF (compared to SGP) increase the loss of molecular clusters by coagulation to pre-existing particles. Among the 8 different nucleation mechanisms in WRF-Chem, we find that the amine+H2SO4 nucleation mechanism dominates at the SGP site, while the pure organic ion induced nucleation mechanism dominates over BNF. Through various WRF-Chem sensitivity simulations, we find that anthropogenic ELVOCs are critical for explaining the growth of newly formed particles and the resulting number size distribution observed near the surface at the SGP site during the daytime. In addition, we show that treating organic particles as semisolid, with strong diffusion-limited uptake of organic vapours, brings model predictions into closer agreement with the observed evolution of particle size distribution. Simulations also predict that anthropogenic SOA, formed by the oxidation of aromatic volatile organic compounds (VOCs), is the dominant organic aerosol component at SGP, while biogenic SOA dominates particle composition at the BNF site in Southeast USA on these days.
The objective of this analysis is to process Lincoln Laboratory Doppler radar data obtained during FLOWS testing at Huntsville, Alabama, in the summer of 1986, to characterize windshear events. The processing includes plotting velocity and F-factor profiles, histogram analysis to summarize statistics, and correlation analysis to demonstrate any correlation between different data fields.
A three-dimensional nine element HF-CW Doppler sounder array has been used to detect ionospheric disturbances which may be due to tornadoes. The typical events chosen in the present study occurred on November 20 and 27, 1973. Both events are apparently associated with tornadoes sighted in the Huntsville, Alabama area. The Doppler records show S-shaped waves rather than the quasi-sinusoidal waves observed in conjunction with and apparently due to thunderstorms. The wave-periods are in the range of 6 to 8 minutes instead of the 3 to 5 minute periods associated with thunderstorms. Dissipation of waves is mostly due to the evanescent effect and they cannot propagate very far from the path of the tornado center. A theory is presented which is in good agreement with the observations.
Drop size distribution (DSD) measurements using ground-based disdrometers (point measurements) have often been used to derive equations to relate radar observations to the integral rainfall parameters (Atlas et al. 1999, Bringi et al., 2003, Kozu et al., 2006, Tokay and Short, 1996, Ajayi and Owolabi, 1987, Battan, 1973). Disdrometers such as JWD, MRR and several others have a major limitation in measuring drops with equi-volume diameters (D(sub eq)) larger than 5 mm because they often rely on the velocity-diameter relationship which plateaus beyond this diameter range (Atlas et al., 1973, Gunn & Kinzer, 1949). Other disdrometers such as Parsivel also lack accuracy beyond this diameter range. The 2D video disdrometer (2DVD: Schönhuber et al., 2008) on the other hand gives drop-shape contours and velocities for each individual drop/hydrometeor falling through its sensor area; this provides a unique opportunity to study the role of very-large drops on radar measurements in particular those with polarimetric radar capability where DSDs with a significant component of very large drops may require special consideration given that the differential reflectivity and other polarimetric radar parameters including attenuation-correction methods will be sensitive to the concentrations of these large drops. A recent study on the occurrence of large drops by Gatlin et al. (2014) has compiled a large and diverse set of measurements made with the 2D video disdrometers from many locations around the globe. Some of the largest drops found in this study were 9 mm D(sub eq) and larger, and in this paper, we report on three such events, with maximum D(sub eq's) of 9.0, 9.1 and 9.7 mm, which occurred in Colorado, Northern Alabama, and Oklahoma, respectively. Detailed examination of the 2DVD data - in terms of shapes and fall velocities - has confirmed that these are fully-melted hydrometeors, although for the last case in Oklahoma, a bigger and non-fully-melted hydrometeor was also observed. All three events were also captured by polarimetric radars, namely the S-band CHILL radar operated by Colorado State University (Brunkow et al., 2000), the C-band ARMOR radar (Petersen et al., 2007) operated by University of Alabama in Huntsville, and NEXRADKVNX, operated by the US National Weather Service, respectively. For the last event, several other radar observations were also made, including two X-band radars operated by the US Dept. of Energy. Analyses of 2DVD data in conjunction with the corresponding radar observations are presented, along with some discussion on sampling issues related to the measurements of such large rain drops. The latter is addressed using maximum diameter D(sub max) measurements from 1-minute DSDs using two collocated 2DVDs for 37 events in Huntsville.
In 1986 the MIST (Microburst and Severe - Thunderstorm) field study took place over northern Alabama. It was noted that the time of day when a peak in microburst activity was observed usually took place around 2000 UTC. It was further confirmed through the examination of the various PAM (Portable Automated Mesonet) stations that this was also the time of day where temperatures would reach their peak. On the afternoon of 18 August 2002, as the temperature approached a daily high temperature of 32 C, the first radar indication of the HSV microburst-producing storm initiated at 2020 UTC.
Mesoscale data collected during both the satellite precipitation and cloud experiment and the microburst and severe thunderstorm program are analyzed in order to describe features associated with two distinct mesoscale precipitation events that occurred about 10 hours apart over the region of northern Alabama to central Tennessee in June 1986. Data sets used include mesobeta-scale rawinsonde data, surface mesonet data, RADAP data, and GOES images. The present mesoscale environment involved the merger of Hurricane Bonnie remnants with a preexisting midlatitude short-wave trough.