Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dataset use case study”

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.

190 records · Page 11

Analysis of Ultra High Resolution Sea Surface Temperature Level 4 Datasets

Sea surface temperature (SST) studies are often focused on improving accuracy, or understanding and quantifying uncertainties in the measurement, as SST is a leading indicator of climate change and represents the longest time series of any ocean variable observed from space. Over the past several decades SST has been studied with the use of satellite data. This allows a larger area to be studied with much more frequent measurements being taken than direct measurements collected aboard ship or buoys. The Group for High Resolution Sea Surface Temperature (GHRSST) is an international project that distributes satellite derived sea surface temperatures (SST) data from multiple platforms and sensors. The goal of the project is to distribute these SSTs for operational uses such as ocean model assimilation and decision support applications, as well as support fundamental SST research and climate studies. Examples of near real time applications include hurricane and fisheries studies and numerical weather forecasting. The JPL group has produced a new 1 km daily global Level 4 SST product, the Multiscale Ultrahigh Resolution (MUR), that blends SST data from 3 distinct NASA radiometers: the Moderate Resolution Imaging Spectroradiometer (MODIS), the Advanced Very High Resolution Radiometer (AVHRR), and the Advanced Microwave Scanning Radiometer ? Earth Observing System(AMSRE). This new product requires further validation and accuracy assessment, especially in coastal regions.We examined the accuracy of the new MUR SST product by comparing the high resolution version and a lower resolution version that has been smoothed to 19 km (but still gridded to 1 km). Both versions were compared to the same data set of in situ buoy temperature measurements with a focus on study regions of the oceans surrounding North and Central America as well as two smaller regions around the Gulf Stream and California coast. Ocean fronts exhibit high temperature gradients (Roden, 1976), and thus satellite data of SST can be used in the detection of these fronts. In this case, accuracy is less of a concern because the primary focus is on the spatial derivative of SST. We calculated the gradients for both versions of the MUR data set and did statistical comparisons focusing on the same regions.

sea surface temperature SST)↗

Updates of MERRA-2 Data and Services at NASA GES DISC

Over40 years of NASA climate reanalysis datasets from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) are available at the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). In addition to being used in traditional weather and climate research, MERRA-2 is also widely used in application studies of, e.g., wind and solar energy, air quality and health, food and drought, and heat waves. Two new MERRA-2 datasets were recently added at the GES DISC: (1) climate statistics derived fromMERRA-2 daily data to assist in the analysis of extreme temperature and precipitation events and of large-scale meteorological patterns from 1980 to the present and (2) gridded satellite and conventional observations processed in the MERRA-2 system, along with key statistics derived from the data assimilation, to help better understand how the quality of observations directly affect there analysis data. The GES DISC focuses its efforts on continually improving existing data services and to develop new data tools to satisfy various user communities. The newly added features include the following: Time series service: This is a new service for MERRA-2 data, which enables the easy and fast access of long-term hourly or daily time series at a location for popular parameters. The data is saved in a single le in Ascii format with a user- friendly structure. New analytic functions in the subsetter interface: Options for downloading daily minimum and maximum values have been added into the subsetter interface, in addition to the existing daily mean option, for all MERRA-2 and MERRA sub-daily products. Data format conversion to GeoTIFF has been implemented. New variables in Giovanni: Most monthly variables have been integrated into Giovanni, GESDISC’s online visualization and analysis tool. Due to the large data volume, hourly variables were selected based on user requests. More online information: New MERRA-2 documentation has been added: Data How-To, Data in Action, and FAQ. This presentation overviews two new MERRA-2 datasets and illustrates the new features of data services through a number of case studies. MERRA-2 data and services can be found at: https://disc.gsfc.gov/datasets?

Data management↗

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

Introduction The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗

Recent Developments to the Porous Microstructure Analysis (PuMA) Software

The Porous Microstructure Analysis (PuMA) software is an open source framework for image-based simulation, primarily used to determine effective properties based on material microstructure. PuMA was originally developed for the study of NASA thermal protection materials; however, many of the solvers in PuMA have applicability to a broad range of materials science applications. PuMA version 3.2 computes material surface area, pore diameters, effective thermal conductivity, continuum and rarefied tortuosity, and permeability. For anisotropic materials, PuMA can estimate material orientation and compute anisotropic thermal conductivity and elasticity. In this talk, a brief overview of the PuMA software and underlying methods will be presented, as well as some recent and ongoing developments, including the use of immersed boundary methods for image-based simulation and the development of a new weave segmentation tool, called TomoSAM. Cut-Cell method for heat and mass transfer For simulations on complex microstructures, traditional unstructured meshing techniques often prove to be difficult and time-intensive. Voxel-based solvers, which represent the surface as a staircase structure, are relatively simple to implement but can lose accuracy when feature resolution is poor. In this work, we present a novel 3D cut-cell method for solving the variable coefficient Poisson equation on complex microstructures, suitable for the determination of effective thermal conductivity or tortuosity of a material. The method uses a Marching Cubes/Marching Squares surface reconstruction to create cut-cells and determine geometric quantities. A flux-correction method is extended to 3D, with least squares gradient reconstruction, to solve for the boundary fluxes in the cut-cells. Verification cases show the solver achieves globally 2nd order accuracy on complex microstructures. TomoSAM TomoSAM, a module of the PuMA software, has been developed as a plugin for 3D Slicer, a software platform used for 3D image processing and visualization. It utilizes the Segment Anything Model (SAM), a deep learning model capable of identifying objects and generating image masks based on minimal user input. This feature enables efficient segmentation of complex 3D datasets, particularly of woven materials, from tomography or similar imaging methods, reducing the need for manual segmentation.

Tomography↗

The Impact of Microphysics and Planetary Boundary Layer Physics on Model Simulation of U.S. Deep South Summer Convection

Inspection of output from various configurations of high-resolution, explicit convection forecast models such as the Weather Research and Forecasting (WRF) model indicates significant sensitivity to the choices of model physics pararneterizations employed. Some of the largest apparent sensitivities are related to the specifications of the cloud microphysics and planetary boundary layer physics packages. In addition, these sensitivities appear to be especially pronounced for the weakly-sheared, multicell modes of deep convection characteristic of the Deep South of the United States during the boreal summer. Possible ocean-land sensitivities also argue for further examination of the impacts of using unique ocean-land surface initialization datasets provided by the NASA Short-term Prediction Research and Transition (SPoRn Center to select NOAAlNWS weather forecast offices. To obtain better quantitative understanding of these sensitivities and also to determine the utility of the ocean-land initialization data, we have executed matrices of regional WRF forecasts for selected convective events near Mobile, AL (MOB), and Houston, TX (HGX). The matrices consist of identically initialized WRF 24-h forecasts using any of eight microphysics choices and any of three planetary boWldary layer choices. The resulting 24 simulations performed for each event within either the MOB or HGX regions are then compared to identify the sensitivities of various convective storm metrics to the physics choices. Particular emphasis is placed on sensitivities of precipitation timing, intensity, and coverage, as well as amount and coverage oflightuing activity diagnosed from storm kinematics and graupel in the mixed phase layer. The results confirm impressions gleaned from study of the behavior of variously configured WRF runs contained in the ensembles produced each spring at the Center for the Analysis and Prediction of Storms, but with the benefit of more straightforward control of the physics package choices. The design of the experiments thus allows for more direct interpretation of the sensitivities to each possible physics combination. The results should assist forecasters in their efforts to anticipate and correct for possible biases in simulated WRF convection patterns, and help the modeling community refine their model parameterizations.

McCaul, Eugene W., Jr.↗

The Impacts of Microphysics and Planetary Boundary Layer Physics on Model Simulations of U. S. Deep South Summer Convection

Inspection of output from various configurations of high-resolution, explicit convection forecast models such as the Weather Research and Forecasting (WRF) model indicates significant sensitivity to the choices of model physics parameterizations employed. Some of the largest apparent sensitivities are related to the specifications of the cloud microphysics and planetary boundary layer physics packages. In addition, these sensitivities appear to be especially pronounced for the weakly-sheared, multicell modes of deep convection characteristic of the Deep South of the United States during the boreal summer. Possible ocean-land sensitivities also argue for further examination of the impacts of using unique ocean-land surface initialization datasets provided by the NASA Short-term Prediction Research and Transition (SPoRT Center to select NOAA/NWS weather forecast offices. To obtain better quantitative understanding of these sensitivities and also to determine the utility of the ocean-land initialization data, we have executed matrices of regional WRF forecasts for selected convective events near Mobile, AL (MOB), and Houston, TX (HGX). The matrices consist of identically initialized WRF 24-h forecasts using any of eight microphysics choices and any of three planetary boundary layer choices. The resulting 24 simulations performed for each event within either the MOB or HGX regions are then compared to identify the sensitivities of various convective storm metrics to the physics choices. Particular emphasis is placed on sensitivities of precipitation timing, intensity, and coverage, as well as amount and coverage of lightning activity diagnosed from storm kinematics and graupel in the mixed phase layer. The results confirm impressions gleaned from study of the behavior of variously configured WRF runs contained in the ensembles produced each spring at the Center for the Analysis and Prediction of Storms, but with the benefit of more straightforward control of the physics package choices. The design of the experiments thus allows for more direct interpretation of the sensitivities to each possible physics combination. The results should assist forecasters in their efforts to anticipate and correct for possible biases in simulated WRF convection patterns, and help the modeling community refine their model parameterizations.

McCaul, E. W., Jr.↗

Use of Collocated KWAJEX Satellite, Aircraft, and Ground Measurements for Understanding Ambiguities in TRMM Radiometer Rain Profile Algorithm

Coordinated ground, aircraft, and satellite observations are analyzed from the 1999 TRMM Kwajalein Atoll field experiment (KWAJEX) to better understand the relationships between cloud microphysical processes and microwave radiation intensities in the context of physical evaluation of the Level 2 TRMM radiometer rain profile algorithm and uncertainties with its assumed microphysics-radiation relationships. This talk focuses on the results of a multi-dataset analysis based on measurements from KWAJEX surface, air, and satellite platforms to test the hypothesis that uncertainties in the passive microwave radiometer algorithm (TMI 2a12 in the nomenclature of TRMM) are systematically coupled and correlated with the magnitudes of deviation of the assumed 3-dimensional microphysical properties from observed microphysical properties. Re-stated, this study focuses on identifying the weaknesses in the operational TRMM 2a12 radiometer algorithm based on observed microphysics and radiation data in terms of over-simplifications used in its theoretical microphysical underpinnings. The analysis makes use of a common transform coordinate system derived from the measuring capabilities of the aircraft radiometer used to survey the experimental study area, i.e., the 4-channel AMPR radiometer flown on the NASA DC-8 aircraft. Normalized emission and scattering indices derived from radiometer brightness temperatures at the four measuring frequencies enable a 2-dimensional coordinate system that facilities compositing of Kwajalein S-band ground radar reflectivities, ARMAR Ku-band aircraft radar reflectivities, TMI spacecraft radiometer brightness temperatures, PR Ku-band spacecraft radar reflectivities, bulk microphysical parameters derived from the aircraft-mounted cloud microphysics laser probes (including liquid/ice water contents, effective liquid/ice hydrometeor radii, and effective liquid/ice hydrometeor variances), and rainrates derived from any of the individual ground, aircraft, or satellite algorithms applied to the radar or radiometer measurements, or their combination. The results support the study's underlying hypothesis, particularly in context of ice phase processes, in that the cloud regions where the 2a12 algorithm's microphysical database most misrepresents the microphysical conditions as determined by the laser probes, are where retrieved surface rainrates are most erroneous relative to other reference rainrates as determined by ground and aircraft radar. In reaching these conclusions, TMI and PR brightness temperatures and reflectivities have been synthesized from the aircraft AMPR and ARMAR measurements with the analysis conducted in a composite framework to eliminate measurement noise associated with the case study approach and single element volumes obfuscated by heterogeneous beam filling effects. In diagnosing the performance of the 2a12 algorithm, weaknesses have been found in the cloud-radiation database used to provide microphysical guidance to the algorithm for upper cloud ice microphysics. It is also necessary to adjust a fractional convective rainfall factor within the algorithm somewhat arbitrarily to achieve satisfactory algorithm accuracy.

Smith, Eric A.↗

High-Resolution Specification of the Land and Ocean Surface for Improving Regional Mesoscale Model Predictions

The exchange of energy and moisture between the Earth's surface and the atmospheric boundary layer plays a critical role in many meteorological processes. High-resolution, accurate representations of surface properties such as sea-surface temperature (SST), soil temperature and moisture content, ground fluxes, and vegetation are necessary to better understand the Earth-atmosphere interactions and improve numerical predictions of sensible weather. The NASA Short-term Prediction Research and Transition (SPoRT) Center has been conducting separate studies to examine the impacts of high-resolution land-surface initialization data from the Goddard Space Flight Center Land Information System (LIS) on subsequent WRF forecasts, as well as the influence of initializing WRF with SST composites derived from the MODIS instrument. This current project addresses the combined impacts of using high-resolution lower boundary data over both land (LIS data) and water (MODIS SSTs) on the subsequent daily WRF forecasts over Florida during May 2004. For this experiment, the WRF model is configured to run on a nested domain with 9- km and 3-kin grid spacing, centered on the Florida peninsula and adjacent coastal waters of the Gulf of Mexico and Atlantic Ocean. A control configuration of WRF is established to take all initial condition data from the NCEP Eta model. Meanwhile, two WRF experimental runs are configured to use high-resolution initialization data from (1) LIS land-surface data only, and (2) a combination of LIS data and high-resolution MODIS SST composites. The experiment involves running 24-hour simulations of the control WRF configuration, the MS-initialized WRF, and the LIS+MODIS-initialized WRF daily for the entire month of May 2004. All atmospheric data for initial and boundary conditions for the Control, LIS, and LIS+MODIS runs come from the NCEP Eta model on a 40-km grid. Verification statistics are generated at land surface observation sites and buoys, and the impacts of the high-resolution lower boundary data on the development and evolution of mesoscale circulations such as sea and land breezes are examined, This paper will present the results of these WRF modeling experiments using LIS and MODIS lower boundary datasets over the Florida peninsula during May 2004.

Case, Jonathan L.↗

Impact of MODIS High-Resolution Sea-Surface Temperatures on WRF Forecasts at NWS Miami, FL

Over the past few years,studies at the Short-term Prediction Research and Transition (SPoRT) Center have suggested that the use of Moderate Resolution Imaging Spectroradiometer (MODIS) composite sea-surface temperature (SST) products in regional weather forecast models can have a significant positive impact on short-term numerical weather prediction in coastal regions. The recent paper by LaCasse et al. (2007, Monthly Weather Review) highlights lower atmospheric differences in regional numerical simulations over the Florida offshore waters using 2-km SST composites derived from the MODIS instrument aboard the polar-orbiting Aqua and Terra Earth Observing System satellites. To help quantify the value of this impact on NWS Weather Forecast Offices (WFOs), the SPoRT Center and the NWS WFO at Miami, FL (MIA) are collaborating on a project to investigate the impact of using the high-resolution MODIS SST fields within the Weather Research and Forecasting (WRF) prediction system. The scientific hypothesis being tested is: More accurate specification of the lower-boundary forcing within WRF will result in improved land/sea fluxes and hence, more accurate evolution of coastal mesoscale circulations and the associated sensible weather elements. The NWS MIA is currently running the WRF system in real-time to support daily forecast operations, using the National Centers for Environmental Prediction Nonhydrostatic Mesoscale Model dynamical core within the NWS Science and Training Resource Center's Environmental Modeling System (EMS) software; The EMS is a standalone modeling system capable of downloading the necessary daily datasets, and initializing, running and displaying WRF forecasts in the NWS Advanced Weather Interactive Processing System (AWIPS) with little intervention required by forecasters. Twenty-seven hour forecasts are run daily with start times of 0300,0900, 1500, and 2100 UTC on a domain with 4-km grid spacing covering the southern half of Florida and the far western portions of the Bahamas, the Florida Keys, the Straights of Florida, and adjacent waters of the Gulf of Mexico and Atlantic Ocean. Each model run is initialized using the Local Analysis and Prediction System (LAPS) analyses available in AWIPS, invoking the diabatic. "hot-start" capability. In this WRF model "hot-start", the LAPS-analyzed cloud and precipitation features are converted into model microphysics fields with enhanced vertical velocity profiles, effectively reducing the model spin-up time required to predict precipitation systems. The SSTs are initialized with the NCEP Real-Time Global (RTG) analyses at l/12 degree resolution (approx. 9 km); however, the RTG product does not exhibit fine-scale details consistent with its grid resolution. SPoRT is conducting parallel WRF EMS runs identical to the operational runs at NWS MIA in every respect except for the use of MODIS SST composites in place of the RTG product as the initial and boundary conditions over water. The MODIS SST composites for initializing the SPoRT WRF runs are generated on a 2-km grid four times daily at 0400, 0700, 1600, and 1900 UTC, based on the times of the overhead passes of the Aqua and Terra satellites. The incorporation of the MODIS SST composites into the SPoRTWRF runs is staggered such that the 0400UTC composite initializes the 0900 UTC WRF, the 0700 UTC composite initializes the 1500 UTC WRF, the 1600 UTC composite initializes the 2100 UTC WRF, and the 1900 UTC composite initializes the 0300 UTC WRF. A comparison of the SPoRT and Miami forecasts is underway in 2007, and includes quantitative verification of near-surface temperature, dewpoint, and wind forecasts at surface observation locations. In addition, particular days of interest are being analyzed to determine the impact of the MODIS SST data on the development and evolution of predicted sea/land-breeze circulations, clouds, and precipitation. This paper will present verification results comparing the NWS MIA forecasts the SPoRT experimental WRF forecasts, and highlight any substantial differences noted in the predicted mesoscale phenomena.

Case, Jonathan L.↗

Martian B Storm Genesis and Evolution: Initial Analysis of Thermal Datasets.

Introduction: Dust lifting on Mars likely occurs primarily as a result of the exchange of momentum between the atmosphere and the surface via saltation. During saltation, sand-sized particles are mobilized but do not enter into suspension. When these larger particles fall back to the surface, kinetic energy is transferred to smaller dust particles which are then lofted into suspension in the atmosphere. Depending on the altitude to which dust is lofted, it can have a significant effect on atmospheric temperatures. As a strong absorber and emitter in the infrared, dust can influence atmospheric heating and modify the global circulation and weather on Mars [1,2]. Although dust is present in Mars’ atmosphere throughout the year, the atmosphere is generally dustier during the second half of the year when Mars is near perihelion. Observations reveal that episodic global-scale dust storms and fairly regular regional-scale dust storms are superimposed on a well-defined and highly repeatable seasonal cycle of dust opacity and associated mid-level atmospheric temperature responses. Kass et al. (2016) used 50 Pa temperature observations from MRO/MCS to identify three highly repeatable time periods during which regional dust storms occur, and designated them the “A”, “B” and “C” storms. While “A” and “C” storms have been studied a fair amount to-date, “B” storms have not yet been investigated in detail. This study explores the generation and evolution of the annually recurring regional dust storm known as the “B” storm, which was identified and categorized by Kass et al. (2016) based on 25 km (50 Pa) temperature observations. The B storm is a southern-hemisphere (SH) phenomenon that originates at the cap edge just after perihelion and which reaches peak intensity during the SH summer solstice, Ls 270. It may originate from the cap edge storms that spawn near the edge of the seasonal CO2 cap during retreat, but the mechanisms for B storm genesis have yet to be determined definitively [1]. Methods: We will use observational data sets and a global climate model (GCM) to investigate “B” regional storms. The data analysis component will include the analysis of imagery from MGS/MOC and MRO/MARCI, and spectroscopic data sets of dust and temperatures from MGS/TES and MRO/MCS with the goal of fully characterizing the behavior of these storms. Both MGS and TES provide data well-suited for temperature analysis at 25 km. MCS measures atmospheric temperature, dust extinction, and water ice extinction at 5 km intervals from the surface to about 80 km. TES measured atmospheric temperatures, column dust and water ice opacities, and column water vapor abundances. Measurements made by TES extended from the surface to about 40 km [1]. At the 50 Pa (25 km) level, local dust events usually confined to shallower depths are effectively filtered out of the analysis leaving the regional dust events identifiable by their temperature signatures [1]. Our preliminary analysis makes use of the fact that the brightness temperature at 15 microns (T15 temperature) is a close approximation to observed temperature at 25 km. We first reproduce the zonal mean 50 Pa level temperature plots for MY 29-32 to establish a baseline for our procedures moving forward [1]. Expanding on Kass et al. (2016), we include recent MCS data from MY 33 and 34 as well. Preliminary Analysis: The daytime (3PM) T15 temperatures in Figure 1 indicate: in MY 29, a strong A storm at Ls 240, a B storm at high southern latitudes just after Ls 270, and a C storm at Ls 320; in MY 30, a B storm at Ls 270; in MY 31 & MY 32, a B storm just before Ls 270; in MY 33, a B storm at Ls 270; and in MY 34, a strong A storm in the northern hemisphere at Ls 210, and a B storm around Ls 270 although there is a data gap. For the B storms, each is indicative of lofted dust and resultant warming. The daytime temperature structure illustrates that the B storm occurs annually around Ls 270 and is confined to high southern latitudes. It reaches its peak intensity around SH summer solstice, Ls 270, consistently for all six MY assessed. Since direct solar heating is absent overnight, the nighttime T15 temperatures (Figure 2) are often useful for differentiating the heat signature of direct solar heating from the dynamical response to that heating. However, in the southern polar latitudes at perihelion the sun does not set and direct solar heating remains present throughout the night. Importantly for our study, dust lofted in the B storm experiences this direct heating day and night for the entirety of its lifetime. The B storm expands as far north as -60 latitude and decays in latitudinal extent more gradually than it grows. This feature is less obvious in the nighttime (3AM) T15 temperatures (Figure 2). The temperature signal is stronger at night for MY 30-33. The warm pool is larger in area relative to the background at night in these four cases. This more uniform warming masks the “tail” feature somewhat, such that it is barely noticeable during these years. Unfortunately, gaps in MCS data in MY 29 and 34 prevent confirmation of the tail feature during those years, however, the B storm temperature signature follows a very different pattern than that described for MY 30-33. MY 29 and 34 appear to show smaller centers of warming at night and larger centers of warming during the day. This is in opposition to that previously described for MY 30-33. Conclusions and Future Work: We will continue investigating the heat signatures of B storms by looking at the total column heating as recorded by TES. We will also look at lower altitudes for patterns that may describe the relationship between B storms and the cap edge storms that develop while the seasonal cap is retreating. In the future, we will use GCM simulations to determine the atmospheric and thermo-dynamic conditions associated with these storms.

Courtney Marylou Batterson↗