Search NASA⌕ Search

SEARCH · Search NASA

Results for “Resolution refinement”

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 505 records · Page 28

The Role of Mesoscale Variability on Plankton Dynamics in the North Atlantic

The intensive field observational phase of JGOFS in the North Atlantic Ocean has shown the importance of oceanic mesoscale variability on biogeochemical cycles and on the strength of the ocean biological pump. Mesoscale physical dynamics govern the major time/space scales of bulk biological variability (biomass, production and export). Mesoscale eddies seem to have a strong impact on the ecosystem structure and functioning, but observational evidence is rather limited. For the signature of the mesoscale features to exist in the ecosystem, the comparison of temporal scales of formation and evolution of mesoscale features and reaction of the ecosystem is a key factor. Biological patterns are driven by active changes in biological source and sink terms rather than simply by passive turbulent mixing. A first modelling assessment of the regional balances between horizontal and vertical eddy-induced nutrient supplies in the euphotic zone shows that the horizontal transport predominates over the vertical route in the subtropical gyre, whereas the reverse holds true for the other biogeochemical provinces of the North Atlantic. Presently. despite some difference in numbers, the net impact of modelled eddies yields an enhancement of the biological productivity in most provinces of the North Atlantic Ocean. Key issues remaining include variation on the mesoscale of subsurface particle and dissolved organic matter remineralization, improved knowledge of the ecological response to patterns of variability, synopticity in mesoscale surveys along with refining measures of biogeochemical time/space variability. Eventual success of assimilation of in situ and satellite data, still in its infancy in coupled physical/biogeochemical models, will be crucial to achieve JGOFS synthesis in answering which data are most informative, standing stocks or rates, and which ones are relevant. Depending on which end of the spectrum quantification of the effect of mesoscale features on production and community structure is required, complementary strategies are offered. Either one may choose to increase resolution of models up to the very fine mesoscale features scale (a few kms) for the high end, or to include a parametric representation of eddies for the low end.

Garcon, Veronique, C.↗

Intelligent Systems: Terrestrial Observation and Prediction Using Remote Sensing Data

NASA has made science and technology investments to better utilize its large space-borne remote sensing data holdings of the Earth. With the launch of Terra, NASA created a data-rich environment where the challenge is to fully utilize the data collected from EOS however, despite unprecedented amounts of observed data, there is a need for increasing the frequency, resolution, and diversity of observations. Current terrestrial models that use remote sensing data were constructed in a relatively data and compute limited era and do not take full advantage of on-line learning methods and assimilation techniques that can exploit these data. NASA has invested in visualization, data mining and knowledge discovery methods which have facilitated data exploitation, but these methods are insufficient for improving Earth science models that have extensive background knowledge nor do these methods refine understanding of complex processes. Investing in interdisciplinary teams that include computational scientists can lead to new models and systems for online operation and analysis of data that can autonomously improve in prediction skill over time.

Coughlan, Joseph C.↗

Radar Reconnaissance of Near-Earth Asteroids

Groundbased radar observations of NEAs can help identify space resources with commercial potential and can dramatically reduce the cost and risk of the initial spacecraft missions to those objects. The near-Earth asteroid (NEA) population is thought to contain approximately 1500 objects as large as a kilometer, approximately 300,000 as large as 100 meters, and more than 100,000,000 as large as 10 meters. More than 10% of the NEAs are more accessible in terms of mission delta-V (i.e., fuel required) than the Moon, Mars, or the moons of Mars. Fewer than 1000 NEAs have been found, but the discovery rate is increasing rapidly. Once an asteroid is discovered (necessarily by wide-field optical telescopes), radar can provide otherwise unavailable information about its size, shape, spin state, and surface properties if it approaches within the range of the Goldstone (California) or Arecibo (Puerto Rico) radar telescopes. Asteroids generally appear as unresolved points through groundbased optical telescopes, but radar measurements of the distribution of echo power in time delay (range) and Doppler frequency (radial velocity) can yield images with resolution as fine as a decameter. Image sequences that furnish adequate orientational coverage can be inverted to construct geologically detailed 3-D models, to define the rotation state precisely, and to constrain the object's internal density distribution. Estimates of radar scattering properties characterize the surface's cm-to-m-scale roughness as well as its bulk density, which for asteroids depends primarily on porosity and metal concentration. A useful spin-off of radar detection is orbit refinement that simplifies navigation of flyby and rendezvous spacecraft.

Ostro, S. J.↗

Radar Reconnaissance of Near-Earth Asteroids

Groundbased radar observations of NEAs can help identify space resources with commercial potential and can dramatically reduce the cost and risk of the initial spacecraft missions to those objects. The near Earth asteroid (NEA) population is thought to contain approx. 1500 objects as large as a kilometer, approx. 300,000 as large as 100 meters, and more than 100,000,000 as large as 10 meters. More than 10% of the NEAs are more accessible in terms of mission delta-V (i.e., fuel required) than the Moon, Mars, or the moons of Mars. Fewer than 1000 NEAs have been found, but the discovery rate is increasing rapidly. Once an asteroid is discovered (necessarily by wide-field optical telescopes), radar can provide otherwise unavailable information about its size, shape, spin state, and surface properties if it approaches within the range of the Goldstone (California) or Arecibo (Puerto Rico) radar telescopes. Asteroids generally appear as unresolved points through groundbased optical telescopes, but radar measurements of the distribution of echo power in time delay (range) and Doppler frequency (radial velocity) can yield images with resolution as fine as a decameter. Image sequences that furnish adequate orientational coverage can be inverted to construct geologically detailed 3-D models, to define the rotation state precisely, and to constrain the object's internal density distribution. Estimates of radar scattering properties characterize the surface's cm-to-m-scale roughness as well as its bulk density, which for asteroids depends primarily on porosity and metal concentration. A useful spin-off of radar detection is orbit refinement that simplifies navigation of flyby and rendezvous spacecraft.

Ostro, S. J.↗

Analysis of Leaf Area Index and Fraction of PAR Absorbed by Vegetation Products from the Terra MODIS Sensor: 2000-2005

The analysis of two years of Collection 3 and five years of Collection 4 Terra Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FPAR) data sets is presented in this article with the goal of understanding product quality with respect to version (Collection 3 versus 4), algorithm (main versus backup), snow (snow-free versus snow on the ground), and cloud (cloud-free versus cloudy) conditions. Retrievals from the main radiative transfer algorithm increased from 55% in Collection 3 to 67% in Collection 4 due to algorithm refinements and improved inputs. Anomalously high LAI/FPAR values observed in Collection 3 product in some vegetation types were corrected in Collection 4. The problem of reflectance saturation and too few main algorithm retrievals in broadleaf forests persisted in Collection 4. The spurious seasonality in needleleaf LAI/FPAR fields was traced to fewer reliable input data and retrievals during the boreal winter period. About 97% of the snow covered pixels were processed by the backup Normalized Difference Vegetation Index-based algorithm. Similarly, a majority of retrievals under cloudy conditions were obtained from the backup algorithm. For these reasons, the users are advised to consult the quality flags accompanying the LAI and FPAR product.

Yang, Wenze↗

Accelerate microstructure evolution simulation using graph neural networks with adaptive spatiotemporal resolution

Abstract Surrogate models driven by sizeable datasets and scientific machine-learning methods have emerged as an attractive microstructure simulation tool with the potential to deliver predictive microstructure evolution dynamics with huge savings in computational costs. Taking 2D and 3D grain growth simulations as an example, we present a completely overhauled computational framework based on graph neural networks with not only excellent agreement to both the ground truth phase-field methods and theoretical predictions, but enhanced accuracy and efficiency compared to previous works based on convolutional neural networks. These improvements can be attributed to the graph representation, both improved predictive power and a more flexible data structure amenable to adaptive mesh refinement. As the simulated microstructures coarsen, our method can adaptively adopt remeshed grids and larger timesteps to achieve further speedup. The data-to-model pipeline with training procedures together with the source codes are provided.

36 MATERIALS SCIENCE↗

Regional-scale fault-to-structure earthquake simulations with the EQSIM framework: Workflow maturation and computational performance on GPU-accelerated exascale platforms

Continuous advancements in scientific and engineering understanding of earthquake phenomena, combined with the associated development of representative physics-based models, is providing a foundation for high-performance, fault-to-structure earthquake simulations. However, regional-scale applications of high-performance models have been challenged by the computational requirements at the resolutions required for engineering risk assessments. The EarthQuake SIMulation (EQSIM) framework, a software application development under the US Department of Energy (DOE) Exascale Computing Project, is focused on overcoming the existing computational barriers and enabling routine regional-scale simulations at resolutions relevant to a breadth of engineered systems. This multidisciplinary software development—drawing upon expertise in geophysics, engineering, applied math and computer science—is preparing the advanced computational workflow necessary to fully exploit the DOE’s exaflop computer platforms coming online in the 2023 to 2024 timeframe. Achievement of the computational performance required for high-resolution regional models containing upward of hundreds of billions to trillions of model grid points requires numerical efficiency in every phase of a regional simulation. This includes run time start-up and regional model generation, effective distribution of the computational workload across thousands of computer nodes, efficient coupling of regional geophysics and local engineering models, and application-tailored highly efficient transfer, storage, and interrogation of very large volumes of simulation data. This article summarizes the most recent advancements and refinements incorporated in the workflow design for the EQSIM integrated fault-to-structure framework, which are based on extensive numerical testing across multiple graphics processing unit (GPU)-accelerated platforms, and demonstrates the computational performance achieved on the world’s first exaflop computer platform through representative regional-scale earthquake simulations for the San Francisco Bay Area in California, USA.

58 GEOSCIENCES↗

CHESS 2025: Orthorectified airborne RGB imagery from NEON AOP surveys

This dataset provides Level 1 (L1) and Level 3 (L3) orthorectified Red-Green-Blue (RGB) imagery collected for the 2025 Colorado Headwaters Ecological Spectroscopy Study (CHESS). This high-resolution imagery is a photographic record of red, green, and blue visible light from sunlight reflected off of the Earth’s surface. The data comprise full-color images of the ground surface and are primarily intended to provide context to imaging spectroscopy and light detection and ranging (LiDAR) data. The data were acquired over three study domains in the Upper Gunnison river basin: the upper East River watershed (CRBU); Almont Triangle and Taylor Canyon (ALMO); and Upper Taylor River watershed (UPTA) between 2025-06-13 and 2025-07-15. RGB images were acquired using the PhaseOne IXM-RS150F high-resolution digital camera onboard the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP). The package data include both an L1 product comprising one camera frame per file and an L3 mosaic aligned to the Universal Transverse Mercator (UTM) Zone 13N grid and the World Geodetic System (WGS) 84 projection. Both products are provided in geotif (.tif) format at 0.1 m ground resolution. The bulk of the imagery was collected during the main CHESS field campaign from June 13 to July 15, 2025. Additional images of a portion of the Upper Taylor (UPTA) domain were collected on September 18, 2025, to fill gaps in imagery identified after the main campaign was complete. RGB camera imagery is not radiometrically calibrated, and therefore pixel values should not be exploited for scientific analysis. Pixel values have undergone a manual adjustment to enhance feature identification. The imagery is rigorously geolocated which does allow for reliable geometric information to be retrieved. To generate the orthorectified imagery, the NEON AOP camera captured visible spectrum in red, green, and blue bands. The raw images were then processed using NEON’s camera orthorectification workflow. A boresight calibration flight was made to build a complete camera, distortion, and alignment model. Color balance/white balance and exposure correction were applied to the raw RGB images. The corrected images were orthorectified by ray-tracing image pixels to a lidar-derived digital surface model (DSM) mesh using the refined camera model, outputting orthorectified raster pixels on a regular grid. Flightline-level data were mosaicked by selecting per-pixel contributions from overlapping orthorectified images using line-of-sight (LOS) zenith angle minimization to reduce edge distortions. CHESS Project Description: The Colorado Headwaters Ecological Spectroscopy Study (CHESS) comprised a multi-week airborne remote sensing and field observation campaign in the Upper Gunnison Basin, Colorado, conducted in June and July of 2025. Airborne remote sensing was conducted by the National Ecological Observatory Network Airborne Observation Platform (NEON AOP), concurrent with a field campaign run by the Rocky Mountain Biological Laboratory (RMBL), the Lawrence Berkeley National Laboratory (LBNL) and SLAC National Accelerator Laboratory Watershed Function Science Focus Area (SFA), and NASA-JPL (Jet Propulsion Laboratory) Earth Surface Mineral Dust Source Investigation (EMIT) program. Between June 10 and July 18, 2025, the NEON AOP flight team collected high-resolution aerial imaging spectroscopy and Light Detection and Ranging (LiDAR) data over three domains: the Upper East River (CRBU), Almont Triangle (ALMO), and the Upper Taylor Basin (UPTA). In coordination with the flights, a field campaign acquired ground-truth observations, including observations of vegetation composition, foliar traits, forest demography, and subsurface properties in 18 core sampling areas within the domains. Additional surface water observations were taken at over 380 point locations. All CHESS campaign datasets can be found within the CHESS ESS-DIVE data portal: https://data.ess-dive.lbl.gov/portals/chess. Funding Acknowledgement: Field and remote-sensing data acquisition was performed under a grant from the National Aeronautics and Space Administration (80NSSC24K1005). This work was also supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

2018 NEON and 2025 CHESS Campaigns↗

Grid refinement studies of turbine rotor-stator interaction

Results from a three-dimensional, time-accurate Navier-Stokes simulation of rotor-stator interaction in an axial turbine stage are presented. The present study uses a fine grid in the spanwise direction to better resolve endwall and tip cllearance effects and complements coarse-grid calculations that were reported earlier. A realistic turbine stage with 22 stator vanes and 28 rotor blades is simulated as a single-stator, single-rotor airfoil combination with the stator geometry modified to properly account for blockage effects. This is in contrast to the earlier coarse-grid calculations where the rotor geometry was modified. The improved grid resolution and the unmodified rotor geometry result in a more accurate simulation of the flow field, particularly in the rotor channel where the interaction effects are more severe. The numerical results are compared to experimental data wherever possible and to earlier calculations.

Madavan, N. K.↗

Flow Characterization of the NASA Langley Unitary Plan Wind Tunnel, Test Section 2: Computational Results

This is an abstract for an invited paper at the AIAA Aviation Conference, June 2021. The work described here is part of an effort of coordinated experiments in the Unitary Plan Wind Tunnel (UPWT) facility at the NASA Langley Research Center (LaRC) and matching CFD simulations. The primary goal of the work is to assess the productivity and true predictive accuracy of CFD, absent any guidance from experiment, in the high supersonic speed range as compared to experiments performed at the NASA LaRC’s UPWT facility. This report concerns CFD simulation of the primary flow-path in the empty wind tunnel, including the settling chamber, nozzle, test section, and some of the tunnel downstream of the test section. The Mach number in the test section ranges from M~2.4 to M~4.6, and the required area ratio variation is achieved by translation of a nozzle block which constricts the area of a loosely S-shaped throat. Flow past protuberances in the settling chamber and into this S-bend throat are predicted by CFD to generate streamwise vorticity that may, or may not, persist through the throat and into the test section as coherent vortices. Some flow conditions are notably unsteady at frequencies well below those of turbulence, due to unsteady separated flow ahead of the nozzle block. The bulk flow moves at velocities ranging from 'walking speed' in the settling chamber to M~4.6 in the test section. Heat transfer to the settling chamber walls and buoyancy are significant at high Mach number. Subtle variations in surface curvature in the nozzle generate Mach waves that propagate into the test section. CFD of the empty tunnel serves two purposes. Firstly, the full-tunnel simulations are used to provide upstream boundary conditions for CFD of vehicle aerodynamics which are generally performed in a domain confined to the wind tunnel test section; these companion studies are addressed in other papers. Secondly, it is a challenging test for CFD to resolve all of the empty tunnel flow phenomena relevant to flow in the test section. It requires a more complete definition of geometry than was originally anticipated. In addition, it requires good spatial and temporal accuracy, and turbulence modeling that performs well on specific phenomena such as corner flows. The boundary conditions and solution algorithms must perform well from incompressible to almost hypersonic speeds. The final state of the pre-test CFD was a product of an iterative self-improvement process. The initial simulations of the empty tunnel were deficient in many respects, but hints to those deficiencies were recognized in the solutions, and remedies were implemented. Possible further improvements will be studied in the post-test phase when comparisons with experimental data are possible. The CFD was performed by five separate collaborative teams using four different flow solvers: FUN3D, Overflow, Star-CCM+ and USM3D. The level of effort of these teams varied significantly, but each made important contributions to the goals of the work. A concerted effort to use uncertainty quantification methods (UQ) in CFD is also a goal of this work. To this end, variations in CFD results due to grid refinement, turbulence modeling, and boundary conditions have been characterized. Code-to-code variation is another means of assessing CFD uncertainty. All CFD solvers predict similarity among the primary flow features; these include the variation of Mach number due to changes in Reynolds number, and the bulk flow angularity due to tunnel-wall curvature. All CFD solvers also predict similar trends in secondary flows, such as the downwash in the side-wall boundary layers. Three of the high-spatial resolution simulations give similar predictions of a complex secondary flow phenomena, streamwise vortices generated in the S-bend throat that persist into the side-wall boundary layers of the test section. Two of the highest-resolution simulations, run with the same turbulence model in different CFD solvers, gave encouragingly similar predictions of a complex tertiary flow phenomenon, small transient "sprites" of upwash flows, resulting from vortices that presumably originate in the separated flow near the leading edge of the nozzle block. While the CFD was done in a "blind pre-test" mode, requests from the experimental team for CFD results pertaining to unsteadiness and total temperature variations in the test section resulted in CFD runs that included a cooled wall in settling chamber. This then led to a change in the standard practice for running the Overflow results, which would not have occurred without this "release" of this experimental information. CFD was also used to guide some measurements. The paper will focus on establishing the consensus among CFD results and understanding differences among those results. Initial findings from the efforts to characterize CFD uncertainty have been done and will be included in the paper.

Robert Edward Childs↗

Vacuum deposition of iridium on large astronomical mirrors for use in the far UV

An iridium coating has been deposited by electron-beam evaporation on a 0.91-m mirror which serves as the telescope primary of a sounding rocket instrument for far-UV spectrometry. The evaporation was carried out by applying 8 kV at 400 mA to the electron gun. Zone refined Ir of 99.99% purity was used, and the electron beam was electromagnetically swept over the surface of the evaporant. Under these conditions, deposition rates of 0.55 A/sec were achieved. The reflectance distribution achieved at a wavelength of 584 A was extremely uniform; the mean reflectance was 21.2% with a standard deviation of only 0.3%. This represents a substantial improvement over Al + MgF2 and Al + LiF coatings for applications involving multiple reflections and weak signals, as might be expected in a high-resolution spectrograph studying distant celestial objects.

Herzig, H.↗

The Unit of Lightning

For the past century, scientists have made quantitative measurements of lightning discharges. In the process, they refined the definition of a lightning unit, or basic quantum of lightning, in order to base it on observable parameters. Although many components of a lightning discharge have been identified, lightning usually occur in groups of discharges or pulses that, although complex, can be organized into units of flashes. This unit definition is based mainly on measurements of lightning from electric field, video, and ground flash lightning locating networks. More recent instrumentation with various combinations of high sensitivity, high temporal, or high spatial resolution often measure signals produced by lightning that do not cleanly divide into flashes. The data from these systems indicate the need or a more fundamental unit for lightning. Such a unit would be of benefit for both basic understanding of lightning and comparing lightning information between instruments. Without a common lightning unit definition, intercomparisons are difficult. For an example, the Lightning Detection And Ranging system (LDAR) at Kennedy Space Center (KSC) have detected ,flash' rates as high as 600 per minute while analysis based on the Advanced Ground Based Field Mill network (AGBFM) detect only 33 "flashes" per minute in the same area and time periods. The satellite based Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) sometimes see single "flashes" that last longer than 10 seconds. Such long duration flashes are not found in electric field records or ground strike location system data sets. The unit of lightning should be based on the fundamental components of the lightning discharge. This should make the unit as generally applicable as possible. For example, studies of NO(x) production by lightning depend on parameters of the individual lightning channels and not the summary flash characteristics. For such studies, the best unit of lightning may be closer to the classical stroke in the electric field data sets, flash in the LDAR data sets, and group in the OTD/LIS data sets. For this type of study, however, more than just the stroke count information is needed. In this paper, we will outline the current definitions of the unit of lightning and propose new basic definitions for the unit of lightning. Since different instruments detect lightning in different ways and different studies may need almost mutually exclusive aspects of the lightning discharge, there may not be a single definition of the unit of lightning that can be applied to all systems and all studies. However, common components can be found since all instrumentation are detecting aspects of the same phenomenon.

Mach, Douglas M.↗

Non-Linear Analysis of the NASA Super Pressure Balloons: Whole Flight Simulations

Tensys have a long-established background in the shape generation and load analysis of architectural stressed membrane structures. Founded upon their inTENS finite element analysis suite, these activities have broadened to encompass 'lighter than air' structures such as aerostats, airships, hybrid air vehicles and stratospheric balloons. Tensys have acted as consultants to the NASA Super Pressure Balloon (SPB) Program since 2004. Previous papers have focussed upon the application of inTENS to the overall structural and stability analysis of pumpkin type balloons as used in the SPB Program. Particular emphasis has been placed upon the ability to study both stress and stability at all stages of a flight. As the program has developed, increasing modelling fidelity has been introduced as design refinement has moved emphasis from the overall shape to the performance of individual details. Examples include the introduction of a contact capability into inTENS to better represent the separate cap layer of film. Localised investigations have considered the consequence of local geometric anomalies introduced during fabrication and the effects of debonding of the PBO tendon within its sleeve. Analysis to date has used a material model for the polyethylene shell film developed by Dr Rand of Winzen Engineering, supported by a program of fine resolution material tests at the Balloon Research and Development Laboratory (BRDL) at the GSFC Wallops Island Facility. Based upon work first reported by Schapery, this model has been incorporated into so-called 'snapshot' analyses within inTENS. For a given elapsed time point and temperature distribution, non-linear material properties are iteratively updated for the current stress state in each individual element until changes in those material properties are insignificant. This process has to be iterative as the film stresses and material properties are interdependent. Attempts to derive a time-stepping incremental viscoelastic capability using the Schapery Rand model encountered problems when dealing with the varying temperatures associated with balloon deployment, pressurisation and diurnal behaviour. An alternative approach has been developed by Pellegrino et al at Caltech, again with support from the NASA Balloon Program Office. This is a large strain non-linear viscoelastic model that includes film out-of-plane mechanical and thermal effects. Working from a unified base, this model can be utilised for time-stepping analyses in both a modulus or compliance mode, or a combination of both. This new capability is intended to enhance the current NASA balloon design process. This paper presents the implementation of the Caltech model into a specialist finite element suite, inTENS and its application to whole flight simulations.

Wakefield, David↗

Trajectory Specification Applied to Terminal Airspace

Despite major efforts to automate air traffic control (ATC), it is still performed by humans today. The complexity and safety-criticality of ATC makes it very difficult to safely automate, but it must be automated to increase airspace capacity (the density of traffic that can be safely managed) and airport throughput (the number of arrivals and departures that an airport can safely handle in a given period of time) beyond what is possible with human controllers. This paper presents the Trajectory Specification (TS) concept, which can help to safely automate ATC. TS is a method of specifying aircraft trajectories such that the position at any given time in flight is restricted to a precisely defined bounding space, removing all ambiguity as to where the flight is allowed to be. The bounding space or volume is determined by tolerances relative to a reference trajectory (position as a function of time). The tolerances are dynamic and are based on the aircraft navigation capabilities and the traffic situation. The tolerances can be a piecewise linear function of time or distance along the route, allowing the tolerances to vary as needed, typically increasing with time for departures and decreasing for arrivals. A Trajectory Specification Language (TSL) is proposed for communicating trajectories from aircraft to ATC as requests and from ATC to aircraft as assignments. The TS concept requires a new generation of airborne Flight Management Systems (FMS) that understand the TSL and can fly the assigned trajectories, but this paper focuses on the ATC functions and the prototype ATC algorithms and software that were developed to test the TS concept. Assuming conformance, TS can guarantee safe separation for an arbitrary length of time even in the event of an ATC system or communication outage. It can help to achieve the high level of safety and reliability needed for ATC automation, and it can also reduce the reliance on ATC backup systems for tactical conflict detection and resolution during normal operation. TS can be applied to any controlled airspace, including enroute, terminal, and urban airspace, but this paper presents algorithms and software for arrival spacing and conflict detection and resolution in the terminal airspace serving a major airport. In a fast-time simulation of a full day of traffic in a major terminal airspace, all conflicts were resolved in near real time, demonstrating the computational feasibility and the preliminary operational feasibility of the TS concept. This paper is a compilation of previous papers, and it adds significant information that was omitted from those papers due to length limitations. It also updates some of the results of those earlier papers due to algorithm refinements and corrections of minor software errors.

air traffic control, trajectory↗

Topography and stratigraphy of Martian polar layered deposits

The first samples of high resolution Viking Orbiter topographic and stratigraphic data for the layered polar deposits of Mars are presented, showing that these deposits are with respect to both slopes and angular relief similar to those in the south. It is also demonstrated that, in conjunction with stereophotogrammetry, photoclinometry holds promise as a tool for detailed layered deposit studies. The spring season photography, which lends itself to photoclinometric analysis, covers the entire area of the north polar deposits. Detailed tests of layered terrain evolution hypotheses will be made, upon refinement of the data by comparison with stereo data. A more promising refining technique will make use of averaging perpendicular to selected sections to enhance SNR. Local reliefs of 200-800 m, and slopes of 1-8 deg, lead to initial calculations of average layer thickness which yields results of 14-46 m, linearly correlated with slope.

Blasius, K. R.↗

Venus small volcano classification and description

The high resolution and global coverage of the Magellan radar image data set allows detailed study of the smallest volcanoes on the planet. A modified classification scheme for volcanoes less than 20 km in diameter is shown and described. It is based on observations of all members of the 556 significant clusters or fields of small volcanoes located and described by this author during data collection for the Magellan Volcanic and Magmatic Feature Catalog. This global study of approximately 10 exp 4 volcanoes provides new information for refining small volcano classification based on individual characteristics. Total number of these volcanoes was estimated to be 10 exp 5 to 10 exp 6 planetwide based on pre-Magellan analysis of Venera 15/16, and during preparation of the global catalog, small volcanoes were identified individually or in clusters in every C1-MIDR mosaic of the Magellan data set. Basal diameter (based on 1000 measured edifices) generally ranges from 2 to 12 km with a mode of 34 km, and follows an exponential distribution similar to the size frequency distribution of seamounts as measured from GLORIA sonar images. This is a typical distribution for most size-limited natural phenomena unlike impact craters which follow a power law distribution and continue to infinitely increase in number with decreasing size. Using an exponential distribution calculated from measured small volcanoes selected globally at random, we can calculate total number possible given a minimum size. The paucity of edifice diameters less than 2 km may be due to inability to identify very small volcanic edifices in this data set; however, summit pits are recognizable at smaller diameters, and 2 km may represent a significant minimum diameter related to style of volcanic eruption. Guest, et al, discussed four general types of small volcanic edifices on Venus: (1) small lava shields; (2) small volcanic cones; (3) small volcanic domes; and (4) scalloped margin domes ('ticks'). Steep-sided domes or 'pancake domes', larger than 20 km in diameter, were included with the small volcanic domes. For the purposes of this study, only volcanic edifices less than 20 km in diameter are discussed. This forms a convenient cutoff since most of the steep-sided domes ('pancake domes') and scalloped margin domes ('ticks') are 20 to 100 km in diameter, are much less numerous globally than are the smaller diameter volcanic edifices (2 to 3 orders of magnitude lower in total global number), and do not commonly occur in large clusters or fields of large numbers of edifices.

Aubele, J. C.↗

Distributed Assimilation of Satellite-based Snow Extent for Improving Simulated Streamflow in Mountainous, Dense Forests: An Example Over the DMIP2 Western Basins

Snow cover area affects snowmelt, soil moisture, evapotranspiration, and ultimately streamflow. For the Distributed Model Intercomparison Project - Phase 2 Western basins, we assimilate satellite-based fractional snow cover area (fSCA) from the Moderate Resolution Imaging Spectroradiometer, or MODIS, into the National Weather Service (NWS) SNOW-17 model. This model is coupled with the NWS Sacramento Heat Transfer (SAC-HT) model inside the National Aeronautics and Space Administration's (NASA) Land Information System. SNOW-17 computes fSCA from snow water equivalent (SWE) values using an areal depletion curve. Using a direct insertion, we assimilate fSCAs in two fully distributed ways: 1) we update the curve by attempting SWE preservation, and 2) we reconstruct SWEs using the curve. The preceding are refinements of an existing simple, conceptually-guided NWS algorithm. Satellite fSCA over dense forests inadequately accounts for below-canopy snow, degrading simulated streamflow upon assimilation during snowmelt. Accordingly, we implement a below-canopy allowance during assimilation. This simplistic allowance and direct insertion are found to be inadequate for improving calibrated results, still degrading them as mentioned above. However, for streamflow volume for the uncalibrated runs, we obtain: (1) substantial to major improvements (64-81 %) as a percentage of the control run residuals (or distance from observations), and (2) minor improvements (16-22 %) as a percentage of observed values. We highlight the need for detailed representations of canopy-snow optical radiative transfer processes in mountainous, dense forest regions if assimilation-based improvements are to be seen in calibrated runs over these areas.

Yatheendradas, Soni↗

A Large-Eddy Simulation Of Turbulent Flow Over A Backward Facing Step

Large-eddy simulation (LES) was used to study the turbulent flow over a backward facing step. The simulations examine the details of the turbulent flow’s separated and reattaching regions created by the step geometry. An explicit high-order of accuracy/high resolution finite-difference computational fluid dynamics solver was used for both LES and Reynolds-averaged Navier-Stokes (RANS) simulations. The experimental data from the backward facing step experiment of Driver and Seegmiller was used for comparison. Data presented are derived from a RANS simulation using the Spalart-Allmaras turbulence model and LES simulations using the implicit LES assumption on three different grids of varying resolution. The RANS solution agrees with past work. It does an adequate job of predicting the flowfield and reattachment point, but misses details such as the strength of the recirculating region and magnitudes of the Reynolds stresses. The synthetic eddy method was used to provide a turbulent inflow to the quasi-2D LES domains. For the LES, the coarse grid was of insufficient resolution and under-predicted the skin friction coefficient and over-predicted the normal stresses. The finer grid LES results accurately predict the skin friction coefficient and provide more accurate Reynolds stresses. Triple product velocity correlations are also computed from the LES and the finer grid results did a good job predicting the magnitudes of the curves and the shapes were reasonably represented. It is recommended that future work include investigating the turbulent inflow domain length, width of the quasi-2D domain, and additional grid refinement

turbulence↗