Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mesoscale modeling”

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.

At least 415 records · Page 23

High-Resolution Simulation of Hurricane Bonnie (1998): Water Budget - Part 2

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, Scott A.↗

The Water Budget of a Simulated Hurricane

The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weaker updrafts, indicative of a larger role of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where cooling from melting enhances condensation. Unlike in squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance in situ stratiform precipitation development through additional growth by vapor deposition and aggregation. We examine artificial source terms for cloud and precipitation mass associated with setting to zero negative mixing ratios that arise from numerical advection errors. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation and 15-20% of the precipitation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Braun, S.↗

Lightning Mapping and the Nowcasting of Severe Storms

This paper describes a successful research and operational collaboration between NASA scientists and NWS forecasters to improve severe stor m warnings using real-time data from a regional VHF total lightning mapping array (LMA). Key objectives of our research using LMA data ar e: a) Identification of intensifying and weakening storms using the time rate-of-change of total flash rate; b) Improved severe storm poten tial situational awareness; c) Evaluation of the potential of total f lash rate trend to improve severe storm probability of detection (POD ) and lead time; and d) Validation of mesoscale model forecasts of th understorm initiation. The LMA data are distributed for ingest and di splay in the WFO AWIPS decision support system, and archived at each WFO for case studies, event playbacks, and assessments using the NWS Warning Event Simulator. The Huntsville WFO has upgraded severe thund erstorm warnings to verified tornado warnings and avoided a false ala rm on a severe storm through the added information on storm growth, intensification, and decay that can be deduced from the magnitude and temporal trend of total flash rates. We present detailed case studies of the observed relationships between lightning activity and tornadi c storm development as determined by radar reflectivity and velocity fields, and thunderstorms forecast by the Weather Research and Foreca st (WRF) model. From these collaborative studies, forecasters can eva luate the added value of total lightning data within the forecast and warning decision-making process (http://weather.msfc.nasa.gov/sport) .

Goodman, S.↗

The Three Gorges Dam Affects Regional Precipitation

Issues regarding building large-scale dams as a solution to power generation and flood control problems have been widely discussed by both natural and social scientists from various disciplines, as well as the policy-makers and public. Since the Chinese government officially approved the Three Gorges Dam (TGD) projects, this largest hydroelectric project in the world has drawn a lot of debates ranging from its social and economic to climatic impacts. The TGD has been partially in use since June 2003. The impact of the TGD is examined through analysis of the National Aeronautics and Space Administration (NASA) Tropical Rainfall Measuring Mission (TRMM) rainfall rate and Moderate Resolution Imaging Spectroradiometer (MODIS) land surface temperature and high-resolution simulation using the Pennsylvania State University-National Center for Atmospheric Research (PSU-NCAR) fifth-generation Mesoscale Model (MM5). The independent satellite data sets and numerical simulation clearly indicate that the land use change associated with the TGD construction has increased the precipitation in the region between Daba and Qinling mountains and reduced the precipitation in the vicinity of the TGD after the TGD water level abruptly rose from 66 to 135 m in June 2003. This study suggests that the climatic effect of the TGD is on the regional scale (approx.100 km) rather than on the local scale (approx.10 km) as projected in previous studies.

Wu, Liguang↗

An Investigation of the Characterization of Cloud Contamination in Hyperspectral Radiances

In regions lacking direct observations, the assimilation of radiances from infrared and microwave sounders is the primary method for characterizing the atmosphere in the analysis process. In recent years, technological advances have led to the launching of more advanced sounders, particularly in the thermal infrared spectrum. With the advent of these hyperspectral sounders, the amount of data available for the analysis process has and will continue to be dramatically increased. However, the utilization of infrared radiances in variational assimilation can be problematic in the presence of clouds; specifically the assessment of the presence of clouds in an instantaneous field of view (IFOV) and the contamination in the individual channels within the IFOV. Various techniques have been developed to determine if a channel is contaminated by clouds. The work presented in this paper and subsequent presentation will investigate traditional techniques and compare them to a new technique, the C02 sorting technique, which utilizes the high spectral resolution of the Atmospheric Infrared Sounder (AIRS) within the framework of the Gridpoint Statistical Interpolation (GSI) 3DVAR system. Ultimately, this work is done in preparation for the assessment of short-term forecast impacts with the regional assimilation of AIRS radiances within the analysis fields of the Weather Research and Forecast Nonhydrostatic Mesoscale Model (WRF-NMM) at the NASA Short-term Prediction Research and Transition (SPORT) Center.

McCarty, William↗

Impact of Local Sensors

Forecasters at the 45th Weather Squadron (45 WS) use observations from the Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) wind tower network and the CCAFS (XMR) daily rawinsonde observations (RAOB) to issue and verify wind advisories and warnings for operations. These observations are also used by the National Weather Service (NWS) Spaceflight Meteorology Group (SMG) in Houston, Texas and the NWS Melbourne, Florida (NWS MLB) to initialize their locally-run mesoscale models. In addition, SMG uses these observations to support shuttle landings at the Shuttle Landing Facility (SLF). Due to impending budget cuts, some or all of the KSC/CCAFS wind towers on the east-central Florida mainland and the XMR RAOBs may be eliminated. The locations of the mainland towers and XMR RAOB site are shown in Figure I. The loss of these data may impact the forecast capability of the 45 WS, SMG and NWS MLB. The AMU was tasked to conduct an objective independent modeling study to help determine how important these observations are to the accuracy of the model output used by the forecasters. To accomplish this, the Applied Meteorology Unit (AMU) performed a sensitivity study using the Weather Research and Forecasting (WRF) model initialized with and without KSC/CCAFS wind tower and XMR RAOB data.

Watson, Leela R.↗

Tool for Forecasting Cool-Season Peak Winds Across Kennedy Space Center and Cape Canaveral Air Force Station (CCAFS)

Peak wind speed is important element in 24-Hour and Weekly Planning Forecasts issued by 45th Weather Squadron (45 WS). Forecasts issued for planning operations at KSC/CCAFS. 45 WS wind advisories issued for wind gusts greater than or equal to 25 kt. 35 kt and 50 kt from surface to 300 ft. AMU developed cool-season (Oct - Apr) tool to help 45 WS forecast: daily peak wind speed, 5-minute average speed at time of peak wind, and probability peak speed greater than or equal to 25 kt, 35 kt, 50 kt. AMU tool also forecasts daily average wind speed from 30 ft to 60 ft. Phase I and II tools delivered as a Microsoft Excel graphical user interface (GUI). Phase II tool also delivered as Meteorological Interactive Data Display System (MIDDS) GUI. Phase I and II forecast methods were compared to climatology, 45 WS wind advisories and North American Mesoscale model (MesoNAM) forecasts in a verification data set.

Barrett, Joe H., III↗

30 WS North Base Wind Study

The 30 Weather Squadron (30 WS) is concerned about strong winds observed at their northern towers without advance warning. They state that terrain influences along the extreme northern fringes of Vandenberg Air Force Base (VAFB) make it difficult for forecasters to issue timely and accurate high wind warnings for northeasterly wind events. These events tend to occur during the winter or early spring when they are under the influence of the Great Basin high pressure weather regime. The Launch Weather Officers (LWOs) have seen these rapid wind increases in the current northern Towers 60, 70 and 71 in excess of their 35 kt operational warning threshold. For this task, the 30 WS requested the Applied Meteorology Unit (AMU) analyze data from days when these towers reported winds in excess of 35 kt and determine if there were any precursors in the observations that would allow the LWOs to better forecast and warn their operational customers for these wind events. The 30 WS provided wind tower data for the cool season (October - March) from the period January 2004-March 20 IO. The AMU decoded and evaluated the wind tower data for 66 days identified by the 30 WS as having high-wind events. Out of the 66 event days, only 30 had wind speed observations of > or =35 kt from at least one of the three northern towers. The AMU analyzed surface and upper air charts to determine the synoptic conditions for each event day along with tower peak wind speed and direction time series and wind rose charts for all 30 event days. The analysis revealed a trend on all event days in which the tower winds shifted to the northeast for a period of time before the first recorded > or =35 kt wind speed. The time periods for the 30 event days ranged from 20 minutes to several hours, with a median value of 110 minutes. This trend, if monitored, could give the 30 WS forecasters a precursor to assist in issuing an operational warning before a high wind event occurs. The AMU recommends developing a high-wind alert capability for VAFB using a local mesoscale model to forecast these wind events. The model should incorporate all of the VAFB local data sets and have a forecast capability of between 2 to 24 hours. Such a model would allow the meteorologists at VAFB to alert the operational customers of high wind events in a timely manner so protective action could be taken.

Wheeler, Mark↗

MesoNAM Verification Phase II

The 45th Weather Squadron Launch Weather Officers use the 12-km resolution North American Mesoscale model (MesoNAM) forecasts to support launch weather operations. In Phase I, the performance of the model at KSC/CCAFS was measured objectively by conducting a detailed statistical analysis of model output compared to observed values. The objective analysis compared the MesoNAM forecast winds, temperature, and dew point to the observed values from the sensors in the KSC/CCAFS wind tower network. In Phase II, the AMU modified the current tool by adding an additional 15 months of model output to the database and recalculating the verification statistics. The bias, standard deviation of bias, Root Mean Square Error, and Hypothesis test for bias were calculated to verify the performance of the model. The results indicated that the accuracy decreased as the forecast progressed, there was a diurnal signal in temperature with a cool bias during the late night and a warm bias during the afternoon, and there was a diurnal signal in dewpoint temperature with a low bias during the afternoon and a high bias during the late night.

Watson, Leela R.↗

African Easterly Waves in 30-day High-Resolution Global Simulations: A Case Study During the 2006 NAMMA Period

In this study, extended -range (30 -day) high-resolution simulations with the NASA global mesoscale model are conducted to simulate the initiation and propagation of six consecutive African easterly waves (AEWs) from late August to September 2006 and their association with hurricane formation. It is shown that the statistical characteristics of individual AEWs are realistically simulated with larger errors in the 5th and 6th AEWs. Remarkable simulations of a mean African easterly jet (AEJ) are also obtained. Nine additional 30 -day experiments suggest that although land surface processes might contribute to the predictability of the AEJ and AEWs, the initiation and detailed evolution of AEWs still depend on the accurate representation of dynamic and land surface initial conditions and their time -varying nonlinear interactions. Of interest is the potential to extend the lead time for predicting hurricane formation (e.g., a lead time of up to 22 days) as the 4th AEW is realistically simulated.

Shen, Bo-Wen↗

Stratosphere-to-Troposphere Transport Revealed by Ground-based Lidar and Ozonesonde at a Midlatitude Site

This paper presents ozone structures measured by a ground-based ozone lidar and ozonesonde at Huntsville, Alabama, on 27-29 April 2010 originating from a stratosphere-to-troposphere transport event associated with a cutoff cyclone and tropopause fold. In this case, the tropopause reached 6 km and the stratospheric intrusion resulted in a 2-km thick elevated ozone layer with values between 70 and 85 ppbv descending from the 306-K to 298-K isentropic surface at a rate of 5 km day1. The potential temperature was provided by a collocated microwave profiling radiometer. We examine the corresponding meteorological fields and potential vorticity (PV) structures derived from the analysis data from the North American Mesoscale model. The 2-PVU (PV unit) surface, defined as the dynamic tropopause, is able to capture the variations of the ozone tropopause estimated from the ozonesonde and lidar measurements. The estimated ozone/PV ratio, from the measured ozone and model derived PV, for the mixing layer between the troposphere and stratosphere is approximately 41 ppbv/PVU with an uncertainty of approximately 33%. Within two days, the estimated mass of ozone irreversibly transported from the stratospheric into the troposphere is between 0.07 Tg (0.9 10(exp33) molecules) and 0.11 Tg (1.3 10(exp33) molecules) with an estimated uncertainty of 59%. Tropospheric ozone exhibited enormous variability due to the complicated mixing processes. Low ozone and large variability were observed in the mid-troposphere after the stratospheric intrusion due to the westerly advection including the transition from a cyclonic system to an anticyclonic system. This study using high temporal and vertical-resolution measurements suggests that, in this case, stratospheric air quickly lost its stratospheric characteristics once it is irreversibly mixed down into the troposphere.

Kuang, Shi↗

Benefits of a 4th Ice Class in the Simulated Radar Reflectivities of Convective Systems Using a Bulk Microphysics Scheme

Numerous cloud microphysical schemes designed for cloud and mesoscale models are currently in use, ranging from simple bulk to multi-moment, multi-class to explicit bin schemes. This study details the benefits of adding a 4th ice class (hail) to an already improved 3-class ice bulk microphysics scheme developed for the Goddard Cumulus Ensemble model based on Rutledge and Hobbs (1983,1984). Besides the addition and modification of several hail processes from Lin et al. (1983), further modifications were made to the 3-ice processes, including allowing greater ice super saturation and mitigating spurious evaporationsublimation in the saturation adjustment scheme, allowing graupelhail to become snow via vapor growth and hail to become graupel via riming, and the inclusion of a rain evaporation correction and vapor diffusivity factor. The improved 3-ice snowgraupel size-mapping schemes were adjusted to be more stable at higher mixing rations and to increase the aggregation effect for snow. A snow density mapping was also added. The new scheme was applied to an intense continental squall line and a weaker, loosely-organized continental case using three different hail intercepts. Peak simulated reflectivities agree well with radar for both the intense and weaker case and were better than earlier 3-ice versions when using a moderate and large intercept for hail, respectively. Simulated reflectivity distributions versus height were also improved versus radar in both cases compared to earlier 3-ice versions. The bin-based rain evaporation correction affected the squall line case more but did not change the overall agreement in reflectivity distributions.

bulk microphysics schemes↗

Intercomparison of Martian Lower Atmosphere Simulated Using Different Planetary Boundary Layer Parameterization Schemes

We use the mesoscale modeling capability of Mars Weather Research and Forecasting (MarsWRF) model to study the sensitivity of the simulated Martian lower atmosphere to differences in the parameterization of the planetary boundary layer (PBL). Characterization of the Martian atmosphere and realistic representation of processes such as mixing of tracers like dust depend on how well the model reproduces the evolution of the PBL structure. MarsWRF is based on the NCAR WRF model and it retains some of the PBL schemes available in the earth version. Published studies have examined the performance of different PBL schemes in NCAR WRF with the help of observations. Currently such assessments are not feasible for Martian atmospheric models due to lack of observations. It is of interest though to study the sensitivity of the model to PBL parameterization. Typically, for standard Martian atmospheric simulations, we have used the Medium Range Forecast (MRF) PBL scheme, which considers a correction term to the vertical gradients to incorporate nonlocal effects. For this study, we have also used two other parameterizations, a non-local closure scheme called Yonsei University (YSU) PBL scheme and a turbulent kinetic energy closure scheme called Mellor- Yamada-Janjic (MYJ) PBL scheme. We will present intercomparisons of the near surface temperature profiles, boundary layer heights, and wind obtained from the different simulations. We plan to use available temperature observations from Mini TES instrument onboard the rovers Spirit and Opportunity in evaluating the model results.

Natarajan, Murali↗

Derivation of Aerosol Profiles for MC3E Convection Studies and Use in Simulations of the 20 May Squall Line Case

Advancing understanding of deep convection microphysics via mesoscale modeling studies of well-observed case studies requires observation-based aerosol inputs. Here, we derive hygroscopic aerosol size distribution input profiles from ground-based and airborne measurements for six convection case studies observed during the Midlatitude Continental Convective Cloud Experiment (MC3E) over Oklahoma. We demonstrate use of an input profile in simulations of the only well-observed case study that produced extensive stratiform outflow on 20 May 2011. At well-sampled elevations between −11 and −23 C over widespread stratiform rain, ice crystal number concentrations are consistently dominated by a single mode near approx. 400 microm in randomly oriented maximum dimension (Dmax). The ice mass at −23 C is primarily in a closely collocated mode, whereas a mass mode near Dmax approx. 1000 microns becomes dominant with decreasing elevation to the −11 C level, consistent with possible aggregation during sedimentation. However, simulations with and without observation-based aerosol inputs systematically overpredict mass peak Dmax by a factor of 3-5 and underpredict ice number concentration by a factor of 4-10. Previously reported simulations with both two-moment and size-resolved microphysics have shown biases of a similar nature. The observed ice properties are notably similar to those reported from recent tropical measurements. Based on several lines of evidence, we speculate that updraft microphysical pathways determining outflow properties in the 20 May case are similar to a tropical regime, likely associated with warm-temperature ice multiplication that is not well understood or well represented in models.

Fridlind, Ann M.↗

Derivation of Aerosol Profiles for MC3E Convection Studies and Use in Simulations of the 20 May Squall Line Case

Advancing understanding of deep convection microphysics via mesoscale modeling studies of well-observed case studies requires observation-based aerosol inputs. Here, we derive hygroscopic aerosol size distribution input profiles from ground-based and airborne measurements for six convection case studies observed during the Midlatitude Continental Convective Cloud Experiment (MC3E) over Oklahoma. We demonstrate use of an input profile in simulations of the only well-observed case study that produced extensive stratiform outflow on 20 May 2011. At well-sampled elevations between -11 and -23 degree C over widespread stratiform rain, ice crystal number concentrations are consistently dominated by a single mode near 400 micrometer in randomly oriented maximum dimension (D[superscript max] ). The ice mass at -23 degree C is primarily in a closely collocated mode, whereas a mass mode near D[superscript max] -1000 micrometer becomes dominant with decreasing elevation to the -11 degree C level, consistent with possible aggregation during sedimentation. However, simulations with and without observation-based aerosol inputs systematically overpredict mass peak D[superscript max] by a factor of 3-5 and underpredict ice number concentration by a factor of 4-10. Previously reported simulations with both two-moment and size-resolved microphysics have shown biases of a similar nature. The observed ice properties are notably similar to those reported from recent tropical measurements. Based on several lines of evidence, we speculate that updraft microphysical pathways determining outflow properties in the 20 May case are similar to a tropical regime, likely associated with warm-temperature ice multiplication that is not well understood or well represented in models.

Fridlind, Ann M.↗

P43K-3884 Characterizing Martian Crater Circulations with the NASA Ames Mars GCM

Observations made in Gale Crater by instruments on the MSL Curiosity Rover show that the diurnal amplitude of the surface pressure is increased and the depth of the Convective Boundary Layer (CBL) is decreased relative to other lander locations on flatter regions of Mars (Haberle et al., 2014; Moores et al., 2015). Mesoscale modeling studies of Gale Crater suggest that crater circulations produce these effects. Tyler & Barnes (2013) show that local upslope/downslope flows along the crater rim and Mt. Sharp amplify the diurnal pressure cycle. These same flows are thought to be at least partly responsible for the suppression of the CBL because upward air flow at the rim and in the center (due to Mt. Sharp) forces subsidence over the lowest regions of the crater during the day. Regional flows, largely due to the location of Gale near the dichotomy boundary, may also play a role in shaping the circulation internal to the crater. Whether the behavior of the CBL and the amplified diurnal pressure cycle are phenomena observed in craters morphologically different from Gale (i.e. bowl-shaped, irregular, degraded) is not yet understood. We will explore these questions by characterizing the behavior of these processes as they are shaped by the morphology of craters greater than 100 km in diameter. We use the NASA Ames Mars Global Circulation Model (GCM) that now utilizes the NOAA/GFDL cubed-sphere finite-volume dynamical core to examine ~100 craters of varying size and shape from a database of known Martian craters (Robbins & Hynek, 2014). Run at 7.5 km resolution, the GCM is capable of resolving surface winds, temperature, and pressure inside craters of this size allowing for the analysis of dozens of craters simulated at various seasons and within the context of synoptic and global-scale phenomena.

Batterson, Courtney MaryLou↗

Comprehensive Observations and Geostatistics of Slope Streaks within the Olympus Mons Aureole

Slope streaks (SS) are enigmatic linear features characterized by relatively low-albedo features that appear and fade on high-albedo slopes on Mars. Despite numerous hypotheses proposed to explain their formation, the primary mechanism behind SS remains elusive. Here, we examine 702 SS features using 32 multitemporal Context Imager (CTX) images and mesoscale modeling data obtained from a site (centered at 31.230°N, 216.281°E) in the Olympus Mons Aureole region. Our investigation revealed several key findings that shed light on the dynamics of SS formation and fading. We discovered a significant preference for SS formation on south-facing (equator-facing) slopes compared to north-facing slopes, with SS being over seven times more likely to occur on the former. Furthermore, SS formation was found to be seasonal with significantly enhanced by a factor of ∼6 near the equinoxes (from solar longitude Ls 337°-42° and 136°-227°) compared to other times of the year. Our analysis also revealed a correlation between the rates of SS formation and fading, with scree slopes exhibited the fastest-fading SS also experiencing the highest rates of newly-formed SS. Additionally, we measured the median starting and stopping slopes of SS to be 23.4° and 14.9°, respectively, significantly below the angle of repose of sand. These low slopes suggest the necessity for an energetic trigger mechanism to initiate SS formation. Infrared spectroscopy revealed that the principal distinction between the material inside and outside of a SS lies in the reduced abundance of dust within the streak. Notably, this site demonstrates the highest rates of SS formation (a 29.6% increase in new SS per Mars year) and fading (a 12% fading of SS per Mars year) ever quantified. These elevated rates may be attributed to the site's topography, which facilitates calm surficial nighttime winds throughout the year, leading to widespread dust deposition. Daytime downhill winds near the ridgelines of S-facing slopes may then trigger movement of newly deposited dust aggregates or alternatively, a Knudsen pump phenomenon could serve as a potential trigger for SS. We propose that SS is inhibited during the northern summer solstice season due to the relatively clear aphelion atmosphere, which limits dust deposition. Similarly, the lack of triggering activity during the winter solstice is likely due to more moderate daytime winds. Overall, our findings are consistent with the wind-triggered dry avalanche hypothesis as a plausible explanation for SS formation.

D E Stillman↗