Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stratification”

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 199 records · Page 11

Origin of Stability in Sedimentation

Velocity fluctuations and particle concentrations are studied in a liquid fluidized bed to investigate the origin of steady state sedimentation. Both the velocity fluctuations and the particle concentrations are found to strongly depend on height. A flux balance model shows why the bed is stable: velocity fluctuations drive a downward particle flux that is compensated by an upward particle flux stemming from an excess flow velocity due to the stratification in concentration. Our results show that in steady state the magnitudes of the fluctuations are related to the degree of stratification.

Segre, P. N.↗

Large Eddy Simulation of Aircraft Wake Vortices in a Homogeneous Atmospheric Turbulence: Vortex Decay and Descent

The effects of ambient turbulence on decay and descent of aircraft wake vortices are studied using a validated, three-dimensional: large-eddy simulation model. Numerical simulations are performed in order to isolate the effect of ambient turbulence on the wake vortex decay rate within a neutrally-stratified atmosphere. Simulations are conducted for a range of turbulence intensities, by injecting wake vortex pairs into an approximately homogeneous and isotropic turbulence field. The decay rate of the vortex circulation increases clearly with increasing ambient turbulence level, which is consistent with field observations. Based on the results from the numerical simulations, simple decay models are proposed as functions of dimensionless ambient turbulence intensity (eta) and dimensionless time (T) for the circulation averaged over a range of radial distances. With good agreement with the numerical results, a Gaussian type of vortex decay model is proposed for weak turbulence: while an exponential type of Tortex decay model can be applied for strong turbulence. A relationship for the vortex descent based on above vortex decay model is also proposed. Although the proposed models are based on simulations assuming neutral stratification, the model predictions are compared to Lidar vortex measurements observed during stable, neutral, and unstable atmospheric conditions. In the neutral and unstable atmosphere, the model predictions appear to be in reasonable agreement with the observational data, while in the stably-stratified atmosphere, they largely underestimate the observed circulation decay with consistent overestimation of the observed vortex descent. The underestimation of vortex decay during stably-stratified conditions suggests that stratification has an important influence on vortex decay when ambient levels of turbulence are weak.

Han, Jongil↗

T wave alternans as a predictor of recurrent ventricular tachyarrhythmias in ICD recipients: prospective comparison with conventional risk markers

INTRODUCTION: The current standard for arrhythmic risk stratification is electrophysiologic (EP) testing, which, due to its invasive nature, is limited to patients already known to be at high risk. A number of noninvasive tests, such as determination of left ventricular ejection fraction (LVEF) or heart rate variability, have been evaluated as additional risk stratifiers. Microvolt T wave alternans (TWA) is a promising new risk marker. Prospective evaluation of noninvasive risk markers in low- or moderate-risk populations requires studies involving very large numbers of patients, and in such studies, documentation of the occurrence of ventricular tachyarrhythmias is difficult. In the present study, we identified a high-risk population, recipients of an implantable cardioverter defibrillator (ICD), and prospectively compared microvolt TWA with invasive EP testing and other risk markers with respect to their ability to predict recurrence of ventricular tachyarrhythmias as documented by ICD electrograms. METHODS AND RESULTS: Ninety-five patients with a history of ventricular tachyarrhythmias undergoing implantation of an ICD underwent EP testing, assessment of TWA, as well as determination of LVEF, baroreflex sensitivity, signal-averaged ECG, analysis of 24-hour Holter monitoring, and QT dispersion from the 12-lead surface ECG. The endpoint of the study was first appropriate ICD therapy for electrogram-documented ventricular fibrillation or tachycardia during follow-up. Kaplan-Meier survival analysis revealed that TWA (P < 0.006) and LVEF (P < 0.04) were the only significant univariate risk stratifiers. EP testing was not statistically significant (P < 0.2). Multivariate Cox regression analysis revealed that TWA was the only statistically significant independent risk factor. CONCLUSIONS: Measurement of microvolt TWA compared favorably with both invasive EP testing and other currently used noninvasive risk assessment methods in predicting recurrence of ventricular tachyarrhythmias in ICD recipients. This study suggests that TWA might also be a powerful tool for risk stratification in low- or moderate-risk patients, and needs to be prospectively evaluated in such populations.

Non-NASA Center↗

Ventless pressure control of two-phase propellant tanks in microgravity

This work studies pressurization and pressure control of a large liquid hydrogen storage tank. A finite element model is developed that couples a lumped thermodynamic formulation for the vapor region with a complete solution of the Navier-Stokes and energy equations for the flow and temperature fields in the liquid. Numerical results show that buoyancy effects are strong, even in microgravity, and can reposition a vapor bubble that is initially at the center of the tank to a region near the tank wall in a relatively short time. Long-term tank pressurization with the vapor bubble at the tank wall shows that after an initial transient lasting about a week, the final rate of pressure increase agrees with a purely thermodynamic analysis of the entire tank. However, the final pressure levels are quite different from thermodynamic predictions. Numerical results also show that there is significant thermal stratification in the liquid due to the effects of natural convection. A subcooled jet is used to provide simultaneous cooling and mixing in order to bring the tank pressure back down to its initial value. Three different jet speeds are examined. Although the lowest jet speed is ineffective at controlling the pressure because of insufficient penetration into the liquid region, the highest jet speed is shown to be quite effective at disrupting thermal stratification and reducing the tank pressure in reasonable time.

Weightlessness↗

Cryogenic Tank Modeling for the Saturn AS-203 Experiment

A computational fluid dynamics (CFD) model is developed for the Saturn S-IVB liquid hydrogen (LH2) tank to simulate the 1966 AS-203 flight experiment. This significant experiment is the only known, adequately-instrumented, low-gravity, cryogenic self pressurization test that is well suited for CFD model validation. A 4000-cell, axisymmetric model predicts motion of the LH2 surface including boil-off and thermal stratification in the liquid and gas phases. The model is based on a modified version of the commercially available FLOW3D software. During the experiment, heat enters the LH2 tank through the tank forward dome, side wall, aft dome, and common bulkhead. In both model and test the liquid and gases thermally stratify in the low-gravity natural convection environment. LH2 boils at the free surface which in turn increases the pressure within the tank during the 5360 second experiment. The Saturn S-IVB tank model is shown to accurately simulate the self pressurization and thermal stratification in the 1966 AS-203 test. The average predicted pressurization rate is within 4% of the pressure rise rate suggested by test data. Ullage temperature results are also in good agreement with the test where the model predicts an ullage temperature rise rate within 6% of the measured data. The model is based on first principles only and includes no adjustments to bring the predictions closer to the test data. Although quantitative model validation is achieved or one specific case, a significant step is taken towards demonstrating general use of CFD for low-gravity cryogenic fluid modeling.

Grayson, Gary D.↗

North Pacific Cloud Feedbacks Inferred from Synoptic-Scale Dynamic and Thermodynamic Relationships

This study analyzed daily satellite cloud observations and reanalysis dynamical parameters to determine how mid-tropospheric vertical velocity and advection over the sea surface temperature gradient control midlatitude North Pacific cloud properties. Optically thick clouds with high tops are generated by synoptic ascent, but two different cloud regimes occur under synoptic descent. When vertical motion is downward during summer, extensive stratocumulus cloudiness is associated with near surface northerly wind, while frequent cloudless pixels occur with southerly wind. Examinations of ship-reported cloud types indicates that midlatitude stratocumulus breaks up as the the boundary level decouples when it is advected equatorward over warmer water. Cumulus is prevalent under conditions of synoptic descent and cold advection during winter. Poleward advection of subtropical air over colder water causes stratification of the near-surface layer that inhibits upward mixing of moisture and suppresses cloudiness until a fog eventually forms. Averaging of cloud and radiation data into intervals of 500-hPa vertical velocity and advection over the SST gradient enables the cloud response to changes in temperature and the stratification of the lower troposphere to be investigated independent of the dynamics.

Norris, Joel R.↗

Situational Lightning Climatologies for Central Florida: Phase IV

The threat of lightning is a daily concern during the warm season in Florida. Research has revealed distinct spatial and temporal distributions of lightning occurrence that are strongly influenced by large-scale atmospheric flow regimes. Previously, the Applied Meteorology Unit (AMU) calculated the gridded lightning climatologies based on seven flow regimes over Florida for 1-, 3- and 6-hr intervals in 5-, 10-,20-, and 30-NM diameter range rings around the Shuttle Landing Facility (SLF) and eight other airfields in the National Weather Service in Melbourne (NWS MLB) county warning area (CWA). In this update to the work, the AMU recalculated the lightning climatologies for using individual lightning strike data to improve the accuracy of the climatologies. The AMU included all data regardless of flow regime as one of the stratifications, added monthly stratifications, added three years of data to the period of record and used modified flow regimes based work from the AMU's Objective Lightning Probability Forecast Tool, Phase II. The AMU made changes so the 5- and 10-NM radius range rings are consistent with the aviation forecast requirements at NWS MLB, while the 20- and 30-NM radius range rings at the SLF assist the Spaceflight Meteorology Group in making forecasts for weather Flight Rule violations during Shuttle landings. The AMU also updated the graphical user interface with the new data.

Bauman, William H., III↗

Estimating Volume, Biomass, and Carbon in Hedmark County, Norway Using a Profiling LiDAR

A profiling airborne LiDAR is used to estimate the forest resources of Hedmark County, Norway, a 27390 square kilometer area in southeastern Norway on the Swedish border. One hundred five profiling flight lines totaling 9166 km were flown over the entire county; east-west. The lines, spaced 3 km apart north-south, duplicate the systematic pattern of the Norwegian Forest Inventory (NFI) ground plot arrangement, enabling the profiler to transit 1290 circular, 250 square meter fixed-area NFI ground plots while collecting the systematic LiDAR sample. Seven hundred sixty-three plots of the 1290 plots were overflown within 17.8 m of plot center. Laser measurements of canopy height and crown density are extracted along fixed-length, 17.8 m segments closest to the center of the ground plot and related to basal area, timber volume and above- and belowground dry biomass. Linear, nonstratified equations that estimate ground-measured total aboveground dry biomass report an R(sup 2) = 0.63, with an regression RMSE = 35.2 t/ha. Nonstratified model results for the other biomass components, volume, and basal area are similar, with R(sup 2) values for all models ranging from 0.58 (belowground biomass, RMSE = 8.6 t/ha) to 0.63. Consistently, the most useful single profiling LiDAR variable is quadratic mean canopy height, h (sup bar)(sub qa). Two-variable models typically include h (sup bar)(sub qa) or mean canopy height, h(sup bar)(sub a), with a canopy density or a canopy height standard deviation measure. Stratification by productivity class did not improve the nonstratified models, nor did stratification by pine/spruce/hardwood. County-wide profiling LiDAR estimates are reported, by land cover type, and compared to NFI estimates.

Nelson, Ross↗

Analytic Modeling of Pressurization and Cryogenic Propellant

An analytic model for pressurization and cryogenic propellant conditions during all mission phases of any liquid rocket based vehicle has been developed and validated. The model assumes the propellant tanks to be divided into five nodes and also implements an empirical correlation for liquid stratification if desired. The five nodes include a tank wall node exposed to ullage gas, an ullage gas node, a saturated propellant vapor node at the liquid-vapor interface, a liquid node, and a tank wall node exposed to liquid. The conservation equations of mass and energy are then applied across all the node boundaries and, with the use of perfect gas assumptions, explicit solutions for ullage and liquid conditions are derived. All fluid properties are updated real time using NIST Refprop.1 Further, mass transfer at the liquid-vapor interface is included in the form of evaporation, bulk boiling of liquid propellant, and condensation given the appropriate conditions for each. Model validation has proven highly successful against previous analytic models and various Saturn era test data and reasonably successful against more recent LH2 tank self pressurization ground test data. Finally, this model has been applied to numerous design iterations for the Altair Lunar Lander, Ares V Core Stage, and Ares V Earth Departure Stage in order to characterize Helium and autogenous pressurant requirements, propellant lost to evaporation and thermodynamic venting to maintain propellant conditions, and non-uniform tank draining in configurations utilizing multiple LH2 or LO2 propellant tanks. In conclusion, this model provides an accurate and efficient means of analyzing multiple design configurations for any cryogenic propellant tank in launch, low-acceleration coast, or in-space maneuvering and supplies the user with pressurization requirements, unusable propellants from evaporation and liquid stratification, and general ullage gas, liquid, and tank wall conditions as functions of time.

Corpening, Jeremy H.↗

Three-Phased Wake Vortex Decay

A detailed parametric study is conducted that examines vortex decay within turbulent and stratified atmospheres. The study uses a large eddy simulation model to simulate the out-of-ground effect behavior of wake vortices due to their interaction with atmospheric turbulence and thermal stratification. This paper presents results from a parametric investigation and suggests improvements for existing fast-time wake prediction models. This paper also describes a three-phased decay for wake vortices. The third phase is characterized by a relatively slow rate of circulation decay, and is associated with the ringvortex stage that occurs following vortex linking. The three-phased decay is most prevalent for wakes imbedded within environments having low-turbulence and near-neutral stratification.

Proctor, Fred H.↗

Peak Wind Tool for General Forecasting

The expected peak wind speed of the day is an important forecast element in the 45th Weather Squadron's (45 WS) daily 24-Hour and Weekly Planning Forecasts. The forecasts are used for ground and space launch operations at the Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS). The 45 WS also issues wind advisories for KSC/CCAFS when they expect wind gusts to meet or exceed 25 kt, 35 kt and 50 kt thresholds at any level from the surface to 300 ft. The 45 WS forecasters have indicated peak wind speeds are challenging to forecast, particularly in the cool season months of October - April. In Phase I of this task, the Applied Meteorology Unit (AMU) developed a tool to help the 45 WS forecast non-convective winds at KSC/CCAFS for the 24-hour period of 0800 to 0800 local time. The tool was delivered as a Microsoft Excel graphical user interface (GUI). The GUI displayed the forecast of peak wind speed, 5-minute average wind speed at the time of the peak wind, timing of the peak wind and probability the peak speed would meet or exceed 25 kt, 35 kt and 50 kt. For the current task (Phase II ), the 45 WS requested additional observations be used for the creation of the forecast equations by expanding the period of record (POR). Additional parameters were evaluated as predictors, including wind speeds between 500 ft and 3000 ft, static stability classification, Bulk Richardson Number, mixing depth, vertical wind shear, temperature inversion strength and depth and wind direction. Using a verification data set, the AMU compared the performance of the Phase I and II prediction methods. Just as in Phase I, the tool was delivered as a Microsoft Excel GUI. The 45 WS requested the tool also be available in the Meteorological Interactive Data Display System (MIDDS). The AMU first expanded the POR by two years by adding tower observations, surface observations and CCAFS (XMR) soundings for the cool season months of March 2007 to April 2009. The POR was expanded again by six years, from October 1996 to April 2002, by interpolating 1000-ft sounding data to 100-ft increments. The Phase II developmental data set included observations for the cool season months of October 1996 to February 2007. The AMU calculated 68 candidate predictors from the XMR soundings, to include 19 stability parameters, 48 wind speed parameters and one wind shear parameter. Each day in the data set was stratified by synoptic weather pattern, low-level wind direction, precipitation and Richardson Number, for a total of 60 stratification methods. Linear regression equations, using the 68 predictors and 60 stratification methods, were created for the tool's three forecast parameters: the highest peak wind speed of the day (PWSD), 5-minute average speed at the same time (A WSD), and timing of the PWSD. For PWSD and A WSD, 30 Phase II methods were selected for evaluation in the verification data set. For timing of the PWSD, 12 Phase\I methods were selected for evaluation. The verification data set contained observations for the cool season months of March 2007 to April 2009. The data set was used to compare the Phase I and II forecast methods to climatology, model forecast winds and wind advisories issued by the 45 WS. The model forecast winds were derived from the 0000 and 1200 UTC runs of the 12-km North American Mesoscale (MesoNAM) model. The forecast methods that performed the best in the verification data set were selected for the Phase II version of the tool. For PWSD and A WSD, linear regression equations based on MesoNAM forecasts performed significantly better than the Phase I and II methods. For timing of the PWSD, none of the methods performed significantly bener than climatology. The AMU then developed the Microsoft Excel and MIDDS GUls. The GUIs display the forecasts for PWSD, AWSD and the probability the PWSD will meet or exceed 25 kt, 35 kt and 50 kt. Since none of the prediction methods for timing of the PWSD performed significantly better thanlimatology, the tool no longer displays this predictand. The Excel and MIDDS GUIs display forecasts for Day-I to Day-3 and Day-I to Day-5, respectively. The Excel GUI uses MesoNAM forecasts as input, while the MIDDS GUI uses input from the MesoNAM and Global Forecast System model. Based on feedback from the 45 WS, the AMU added the daily average wind speed from 30 ft to 60 ft to the tool, which is one of the parameters in the 24-Hour and Weekly Planning Forecasts issued by the 45 WS. In addition, the AMU expanded the MIDDS GUI to include forecasts out to Day-7.

Barrett, Joe H., III↗

Experimental and Numerical Investigation of Reduced Gravity Fluid Slosh Dynamics for the Characterization of Cryogenic Launch and Space Vehicle Propellants

As space programs increasingly investigate various options for long duration space missions the accurate prediction of propellant behavior over long periods of time in microgravity environment has become increasingly imperative. This has driven the development of a detailed, physics-based understanding of slosh behavior of cryogenic propellants over a range of conditions and environments that are relevant for rocket and space storage applications. Recent advancements in computational fluid dynamics (CFD) models and hardware capabilities have enabled the modeling of complex fluid behavior in microgravity environment. Historically, launch vehicles with moderate duration upper stage coast periods have contained very limited instrumentation to quantify propellant stratification and boil-off in these environments, thus the ability to benchmark these complex computational models is of great consequence. To benchmark enhanced CFD models, recent work focuses on establishing an extensive experimental database of liquid slosh under a wide range of relevant conditions. In addition, a mass gauging system specifically designed to provide high fidelity measurements for both liquid stratification and liquid/ullage position in a micro-gravity environment has been developed. This pUblication will summarize the various experimental programs established to produce this comprehensive database and unique flight measurement techniques.

Walls, Laurie K.↗

Enthalpy-Based Thermal Evolution of Loops: II. Improvements to the Model

This paper further develops the zero-dimensional (0D) hydrodynamic coronal loop model "Enthalpy-based Thermal Evolution of Loops" (EBTEL) originally proposed by Klimchuk et al (2008), which studies the plasma response to evolving coronal heating. It has typically been applied to impulsive heating events. The basis of EBTEL is the modelling of mass exchange between the corona and transition region and chromosphere in response to heating variations, with the key parameter being the ratio of transition region to coronal radiation. We develop new models for this parameter that now include gravitational stratification and a physically motivated approach to radiative cooling. A number of examples are presented, including nanoflares in short and long loops, and a small flare. It is found that while the evolution of the loop temperature is rather insensitive to the details of the model, accurate tracking of the density requires the inclusion of our new features. In particular, we are able to now obtain highly over-dense loops in the late cooling phase and decreases to the coronal density arising due to stratification. The 0D results are compared to a 1D hydro code (Hydrad). The agreement is acceptable, with the exception of the flare case where some versions of Hydrad can give significantly lower densities. This is attributed to the method used to model the chromosphere in a flare. EBTEL is suitable for general use as a tool for (a) quick-look results of loop evolution in response to a given heating function and (b) situations where the modelling of hundreds or thousands of elemental loops is needed. A single run takes a few seconds on a contemporary laptop.

Cargill, P. J.↗

Mixed Layer Sub-Mesoscale Parameterization - Part 1: Derivation and Assessment

Several studies have shown that sub-mesoscales (SM 1km horizontal scale) play an important role in mixed layer dynamics. In particular, high resolution simulations have shown that in the case of strong down-front wind, the re-stratification induced by the SM is of the same order of the de-stratification induced by small scale turbulence, as well as of that induced by the Ekman velocity. These studies have further concluded that it has become necessary to include SM in ocean global circulation models (OGCMs), especially those used in climate studies. The goal of our work is to derive and assess an analytic parameterization of the vertical tracer flux under baroclinic instabilities and wind of arbitrary directions and strength. To achieve this goal, we have divided the problem into two parts: first, in this work we derive and assess a parameterization of the SM vertical flux of an arbitrary tracer for ocean codes that resolve mesoscales, M, but not sub-mesoscales, SM. In Part 2, presented elsewhere, we have used the results of this work to derive a parameterization of SM fluxes for ocean codes that do not resolve either M or SM. To carry out the first part of our work, we solve the SM dynamic equations including the non-linear terms for which we employ a closure developed and assessed in previous work. We present a detailed analysis for down-front and up-front winds with the following results: (a) down-front wind (blowing in the direction of the surface geostrophic velocity) is the most favorable condition for generating vigorous SM eddies; the de-stratifying effect of the mean flow and re-stratifying effect of SM almost cancel each other out,

Canuto, V. M.↗

Turbulence in Supercritical O2/H2 and C7H16/N2 Mixing Layers

This report presents a study of numerical simulations of mixing layers developing between opposing flows of paired fluids under supercritical conditions, the purpose of the study being to elucidate chemical-species-specific aspects of turbulence. The simulations were performed for two different fluid pairs O2/H2 and C7H16/N2 at similar reduced initial pressures (reduced pressure is defined as pressure divided by critical pressure). Thermodynamically, O2/H2 behaves more nearly like an ideal mixture and has greater solubility, relative to C7H16/N2, which departs strongly from ideality. Because of a specified smaller initial density stratification, the C7H16/N2 layers exhibited greater levels of growth, global molecular mixing, and turbulence. However, smaller density gradients at the transitional state for the O2/H2 system were interpreted as indicating that locally, this system exhibits enhanced mixing as a consequence of its greater solubility and closer approach to ideality. These thermodynamic features were shown to affect entropy dissipation, which was found to be larger for O2/H2 and concentrated in high-density-gradient-magnitude regions that are distortions of the initial density-stratification boundary. In C7H16/N2, the regions of largest dissipation were found to lie in high-density-gradient-magnitude regions that result from mixing of the two fluids.

Bellan, Josette↗

Using Cloud-to-Ground Lightning Climatologies to Initialize Gridded Lightning Threat Forecasts for East Central Florida

Each morning, the forecasters at the National Weather Service in Melbourn, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in central Florida, especially during the warm season months of May-September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC (0700 AM EST) each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to increase consistency between forecasters while enabling them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) for which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities from National Lightning Detection Network (NLDN) data. The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at the two organizations provided this data and supporting software for the work performed by the AMU. The densities were first stratified by flow regime, then by time in 1-, 3-, 6-, 12-, and 24-hour increments while maintaining the 2.5 km x 2.5 km grid resolution. A CG frequency of occurrence was calculated for each stratification and grid box by counting the number of days with lightning and dividing by the total number of days in the data set. New CG strike densities were calculated for each stratification and grid box by summing the strike number values over all warm seasons, then normalized by dividing the summed values by the number of lightning days. This makes the densities conditional on whether lightning occurred. The frequency climatology values will be used by forecasters as proxy inputs for lightning prObability, while the density climatology values will be used for CG amount. In addition to the benefits outlined above, these climatologies will provide improved temporal and spatial resolution, expansion of the lightning threat area to include adjacent coastal waters, and potential to extend the forecast to include the day-2 period. This presentation will describe the lightning threat index map, discuss the work done to create the maps initialized with climatological guidance, and show examples of the climatological CG lightning densities and frequencies of occurren based on flow regime.

Lambert, Winnie↗

Improving Post-Hurricane Katrina Forest Management with MODIS Time Series Products

Hurricane damage to forests can be severe, causing millions of dollars of timber damage and loss. To help mitigate loss, state agencies require information on location, intensity, and extent of damaged forests. NASA's MODerate Resolution Imaging Spectroradiometer (MODIS) Normalized Difference Vegetation Index (NDVI) time series data products offers a potential means for state agencies to monitor hurricane-induced forest damage and recovery across a broad region. In response, a project was conducted to produce and assess 250 meter forest disturbance and recovery maps for areas in southern Mississippi impacted by Hurricane Katrina. The products and capabilities from the project were compiled to aid work of the Mississippi Institute for Forest Inventory (MIFI). A series of NDVI change detection products were computed to assess hurricane induced damage and recovery. Hurricane-induced forest damage maps were derived by computing percent change between MODIS MOD13 16-day composited NDVI pre-hurricane "baseline" products (2003 and 2004) and post-hurricane NDVI products (2005). Recovery products were then computed in which post storm 2006, 2007, 2008 and 2009 NDVI data was each singularly compared to the historical baseline NDVI. All percent NDVI change considered the 16-day composite period of August 29 to September 13 for each year in the study. This provided percent change in the maximum NDVI for the 2 week period just after the hurricane event and for each subsequent anniversary through 2009, resulting in forest disturbance products for 2005 and recovery products for the following 4 years. These disturbance and recovery products were produced for the Mississippi Institute for Forest Inventory's (MIFI) Southeast Inventory District and also for the entire hurricane impact zone. MIFI forest inventory products were used as ground truth information for the project. Each NDVI percent change product was classified into 6 categories of forest disturbance intensity. Stand age and stand type raster data, also provided by MIFI, were used along with the forest disturbance/recovery products to create forest damage stratification products integrating 3 stand type classes, 6 stand age classes, and 6 forest disturbance intensity classes. This stratification product will be used to aid MIFI timber inventory planning and to prepare for damage assessments due to future hurricane events. Validation of MODIS percent NDVI change products was performed by comparing the MODIS percent NDVI change products to those from Landsat data for the same time and MIFI inventory district area.

Lewis, Mark David↗

Exploring the Influence of Impervious Surface Density and Shape on Urban Heat Islands in the Northeast USA Using MODIS and Landsat

Impervious surface area (ISA) from the National Land Cover Database (NLCD) 2001 and land surface temperature (LST) from MODIS averaged over three annual cycles (2003-2005) are used in a spatial analysis to assess the urban heat island (UHI) signature and its relationship to settlement size and shape, development intensity distribution, and land cover composition for 42 urban settlements embedded in forest biomes in the Northeastern United States. Development intensity zones, based on percent ISA, are defined for each urban area emanating outward from the urban core to nearby rural areas and are used to stratify land surface temperature. The stratification is further constrained by biome type and elevation to insure objective intercomparisons between urban zones within an urban settlement and between settlements. Stratification based on ISA allows the definition of hierarchically ordered urban zones that are consistent across urban settlements and scales. In addition to the surrounding ecological context, we find that the settlement size and shape as well as the development intensity distribution significantly influence the amplitude of summer daytime UHI. Within the Northeastern US temperate broadleaf mixed forest, UHI magnitude is positively related to the logarithm of the urban area size. Our study indicates that for similar urban area sizes, the development intensity distribution is one of the major drivers of UHI. In addition to urban area size and development intensity distribution, this analysis shows that both the shape of the urban area and the land cover composition in the surrounding rural area play an important role in modulating the UHI magnitude in different urban settlements. Our results indicate that remotely sensed urban area size and shape as well as the development intensity distribution influence UHI amplitude across regional scales.

Zhang, Ping↗