Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “temporal computing”

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 883 records · Page 49

Implicit upwind-Euler solution algorithms for unstructured-grid applications

The development of implicit upwind algorithms for the solution of the three-dimensional, time-dependent Euler equations on unstructured tetrahedral meshes is described. The implicit temporal discretization involves either a two-sweep Gauss-Seide relaxation procedure, a two-sweep Point-Jacobi relaxation procedure, or a single-sweep Point-Implicit procedure; the upwind spatial discretization is based on the flux-difference splitting of Roe. Detailed descriptions of the three implicit solution algorithms are given, and calculations for the Boeing 747 transport configuration are presented to demonstrate the algorithms. Advantages and disadvantages of the implicit algorithms are discussed. A steady-state solution for the 747 configuration, obtained at transonic flow conditions using a mesh of over 100,000 cells, required less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the general capability.

Batina, John T.↗

Implicit upwind solution algorithms for three-dimensional unstructured meshes

The development of implicit upwind algorithms for the solution of the three-dimensional, time-dependent Euler equations on unstructured tetrahedral meshes is described. The implicit temporal discretization involves either a two-sweep Gauss-Seide relaxation procedure, a two-sweep Point-Jacobi relaxation procedure, or a single-sweep Point-Implicit procedure; the upwind spatial discretization is based on the flux-difference splitting of Roe. Detailed descriptions of the three implicit solution algorithms are given, and calculations for the Boeing 747 transport configuration are presented to demonstrate the algorithms. Advantages and disadvantages of the implicit algorithms are discussed. A steady-state solution for the 747 configuration, obtained at transonic flow conditions using a mesh of over 100,000 cells, required less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the general capability.

Batina, John T.↗

Implicit upwind-Euler solution algorithms for unstructured-grid applications

The development of implicit upwind algorithms for the solution of the three-dimensional, time-dependent Euler equations on unstructured tetrahdral meshes is described. The implicit temporal discretization involves either a two-sweep Gauss-Seidel relaxation procedure, a two-sweep Point-Jacobi relaxation procedure, or a single-sweep Point-Implicit procedure; the upwind spatial discretization is based on the flux-difference splitting of Roe. Detailed descriptions of the three implicit solution algorithms are given, and calculations for the Boeing 747 transport configuration are presented to demonstrate the algorithms. Advantages and disadvantages of the implicit algorithms are discussed. A steady-state solution for the 747 configuration, obtained at transonic flow conditions using s mesh of over 100,000 cells, required less than one hour of CPU time on a Cray-2 computer, thus demonstrating the speed and robustness of the general capability.

Batina, John T.↗

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗

Gridded Hourly Text Products: A TRMM Data Reduction Approach

The quantity of precipitation data from satellite-based observations is a blessing and a curse. The sheer volume of the data makes it difficult for many researchers to use in targeted applications. This volume increases further as algorithm improvements lead to the reprocessing of mission data. In addition to the overall volume of data, the size and format complexity of orbital granules contribute to the difficulty in using all the available data. Finally, the number of different instruments available to measure rainfall and related parameters further contributes to the volume concerns. In summary, we have an embarrassment of riches. The science team of the Tropical Rainfall Measuring Mission (TRMM) recognized this dilemma and has developed a strategy to address it. The TRMM Science Data and Information System (TSDIS) produces, at the direction of the Joint TRMM Science Team, a number of instantaneous rainfall products. The TRMM Microwave Imager (TMI), the Precipitation Radar and a Combined TMI/PR are the key "instruments" used in this production. Each of these products contains an entire orbit of data. The algorithm code computes not just rain rates but a large number of other physical parameters as well as information needed for monitoring algorithm performance. That makes these products very large. For example, a single orbit of TMI rain rate product is 99 MB, a single orbit of the combined product yields a granule that is 158 MB, while the 80 vertical levels of rain information from the PR yields an orbital product of 253 MB. These are large products that are often difficult for science users to electronically transfer to their sites especially if they want a large period of time. Level 3 gridded products are much smaller, but their 5 or 30 day temporal resolution is insufficient for many researchers. In addition, TRMM standard products are produced in the HDF format. While a large number of user-friendly tools are available to hide the details of the format (including a toolkit developed at TSDIS for the TRMM science team), many potential users shy away

Stocker, Erich↗

Biomass Burning Emissions in the Cerrado of Brazil Computed with Remote Sensing Data and GIS

Biomass burnin is a common force in much of the developing tropical world where it has wide-ranging environmental impacts. Fire is a component of tropical deforestation and is 0 p often used to clear broad expanses of land for shifting agriculture and cattle ranching. Frequent burning in the tropical savannas is a distinct problem from that of primary forest. In Brazil, most of the burning occurs in the cerrado which occupies approximately 1,800,000 km2, primarily on the great plateau in central Brazil. Wildland and agricultural fires are dramatic sources of regional air pollution in central Brazil. Biomass burning is an important source of a large number of trace gases including greenhouse gases and other chemically active species. Knowledge of trace gas emissions from biomass burning in Brazil is limited by a number of factors, most notably relative emission factors for gases from specific fire types/fuels and accurate estimates of temporal and spatial distribution and extent of fire activity. Estimates of trace gas emissions during September 1992 will be presented that incorporates a digital map of vegetation classes, pyrogenic emission factors calculated from ground and aircraft missions, and Instituto Nacional de Pesquisas Espaciais (INPE) fire products derived from Advanced Very High Resolution Radiometer (AVHRR) data. The regional emissions calculated from National Oceanographic and Atmospheric Administration (NOAA) AVHRR estimates of fire activity will provide an independent estimate for comparison with results obtained by the National Aeronautics and Space Administration (NASA) Transport and Atmospheric Chemistry Near the Equator - Atlantic (TRACE-A) experiments.

Guild, Liane S.↗

An Exploration of Discontinuous Time Synchronous Averaging for Helicopter HUMS Using Cruise and Terminal Area Vibration Data

Recent research using NASA Ames AH-1 and OH-58C helicopters, and NASA Glenn test rigs, has shown that in-flight vibration data are typically non-stationary [l-4]. The nature and extent of this non-stationarity is most likely produced by several factors operating simultaneously. The aerodynamic flight environment and pilot commands provide continuously changing inputs, with a complex dynamic response that includes automatic feedback control from the engine regulator. It would appear that the combined effects operate primarily through an induced torque profile, which causes concomitant stress modulation at the individual internal gear meshes in the transmission. This notion is supported by several analyses, which show that upwards of 93% of the vibration signal s variance can be explained by knowledge of torque alone. That this relationship is stronger in an AH-1 than an OH-58, where measured non-stationarity is greater, suggests that the overall mass of the vehicle is an important consideration. In the lighter aircraft, the unsteady aerodynamic influences transmit relatively greater unsteady dynamic forces on the mechanical components, quite possibly contributing to its greater non-stationarity . In a recent paper using OH-58C pinion data [5], the authors have shown that in computing a time synchronous average (TSA) for various single-value metric computations, an effective trade-off can be obtained between sample size and measured stationarity by using data from only a single mesh cycle. A mesh cycle, which is defined as the number of rotations required for the gear teeth to return to their original mating position, has the property of representing all of the discrete phase angles of the opposing gears exactly once in the average. Measured stationarity is probably maximized because a single mesh cycle of the pinion gear occurs over a very short span of time, during which time-dependent non-stationary effects are kept to a minimum. Clearly, the advantage of local stationarity diminishes as the temporal duration of the cycle increases. This is most evident for a planetary mesh cycle, which can take several minutes to complete.

Huff, Edward M.↗

Intrinsic Dimensionality as a Metric for the Impact of Mission Design Parameters

High-resolution space-based spectral imaging of the Earth's surface delivers critical information for monitoring changes in the Earth system as well as resource management and utilization. Orbiting spectrometers are built according to multiple design parameters, including ground sampling distance (GSD), spectral resolution, temporal resolution, and signal-to-noise ratio. Different applications drive divergent instrument designs, so optimization for wide-reaching missions is complex. The Surface Biology and Geology component of NASA's Earth System Observatory addresses science questions and meets applications needs across diverse fields, including terrestrial and aquatic ecosystems, natural disasters, and the cryosphere. The algorithms required to generate the geophysical variables from the observed spectral imagery each have their own inherent dependencies and sensitivities, and weighting these objectively is challenging. Here, we introduce intrinsic dimensionality (ID), a measure of information content, as an applications-agnostic, data-driven metric to quantify performance sensitivity to various design parameters. ID is computed through the analysis of the eigenvalues of the image covariance matrix, and can be thought of as the number of significant principal components. This metric is extremely powerful for quantifying the information content in high-dimensional data, such as spectrally resolved radiances and their changes over space and time. We find that the ID decreases for coarser GSD, decreased spectral resolution and range, less frequent acquisitions, and lower signal-to-noise levels. This decrease in information content has implications for all derived products. ID is simple to compute, providing a single quantitative standard to evaluate combinations of design parameters, irrespective of higher-level algorithms, products, applications, or disciplines.

Intrinsic dimensionality↗

Satellite data for surface-mine inventory

To determine the feasibility of satellite data for surface-mine inventory, particularly as it applies to coal, a case study was conducted in Maryland. A band-ratio method was developed to measure disturbed surface areas, and it proved to be extendible both temporally and geographically. This method was used to measure area changes in the region over three time periods from September 1972 through July 1974 and to map the entire two-county area for 1973. For mines ranging between 31 and 244 acres (12 to 98 hectares) the measurement accuracy of total affected acreage was determined to be 92%. Mines of 120 acres (50 hectares) and larger were measured with greater accuracy, some within one percent of the actual area. The ability to identify, classify, and measure strip-mine surfaces in a two-county area (1,541 square kilometers - 595 square miles) of western Maryland was demonstrated through the use of computer processing. On the basis of these results the use of LANDSAT satellite data and multilevel sampling of aircraft and field verification inspections, multispectral analysis of digital data is shown to be an effective, rapid, and accurate means of monitoring the surface mining cycle.

Anderson, A. T.↗

Preliminary evaluation of the airborne imaging spectrometer for vegetation analysis in the Klamath National Forest of northeastern California

The experiences and results associated with a project entitled Preliminary Evaluation of the Airborne Imaging Spectrometer for Vegetation Analysis is documented. The primary goal of the project was to provide ground truth, manual interpretation, and computer processing of data from an experimental flight of the Airborne Infrared Spectrometer (AIS) to determine the extent to which high spectral resolution remote sensing could differentiate among plant species, and especially species of conifers, for a naturally vegetated test site. Through the course of the research, JPL acquired AIS imagery of the test areas in the Klamath National Forest, northeastern California, on two overflights of both the Dock Well and Grass Lake transects. Over the next year or so, three generations of data was also received: first overflight, second overflight, and reprocessed second overflight. Two field visits were made: one trip immediately following the first overflight to note snow conditions and temporally-related vegetation states at the time of the sensor overpass; and a second trip about six weeks later, following acquisition of prints of the images from the first AIS overpass.

Strahler, A. H.↗

Model based estimation of image depth and displacement

Passive depth and displacement map determinations have become an important part of computer vision processing. Applications that make use of this type of information include autonomous navigation, robotic assembly, image sequence compression, structure identification, and 3-D motion estimation. With the reliance of such systems on visual image characteristics, a need to overcome image degradations, such as random image-capture noise, motion, and quantization effects, is clearly necessary. Many depth and displacement estimation algorithms also introduce additional distortions due to the gradient operations performed on the noisy intensity images. These degradations can limit the accuracy and reliability of the displacement or depth information extracted from such sequences. Recognizing the previously stated conditions, a new method to model and estimate a restored depth or displacement field is presented. Once a model has been established, the field can be filtered using currently established multidimensional algorithms. In particular, the reduced order model Kalman filter (ROMKF), which has been shown to be an effective tool in the reduction of image intensity distortions, was applied to the computed displacement fields. Results of the application of this model show significant improvements on the restored field. Previous attempts at restoring the depth or displacement fields assumed homogeneous characteristics which resulted in the smoothing of discontinuities. In these situations, edges were lost. An adaptive model parameter selection method is provided that maintains sharp edge boundaries in the restored field. This has been successfully applied to images representative of robotic scenarios. In order to accommodate image sequences, the standard 2-D ROMKF model is extended into 3-D by the incorporation of a deterministic component based on previously restored fields. The inclusion of past depth and displacement fields allows a means of incorporating the temporal information into the restoration process. A summary on the conditions that indicate which type of filtering should be applied to a field is provided.

Damour, Kevin T.↗

A New Time-dependent Model for the Martian Radiation Environment

Manned space activities have been until present time limited to the near-Earth environment, most of them to low Earth orbit (LEO) scenarios, with only some of the Apollo missions targeted to the Moon. In current times most human exploration and development of space (HEDS) activities are related to the development of the International Space Station (ISS), and therefore take place in the LEO environment. A natural extension of HEDS activities will be going beyond LEO, and reach asteroids, Mars, Jupiter, Saturn, the Kuiper belt and the outskirts of the Solar System. Such long journeys onboard spacecraft outside the protective umbrella of the geomagnetic field will require higher levels of protection from the radiation environment found in the deep space for both astronauts and equipment. So, it is important to have available a tool for radiation shielding which takes into account the radiation environments found all along the interplanetary space and at the different bodies encountered in the Solar System. Moreover, the radiation protection is one of the two NASA highest concerns and priorities. A tool integrating different radiation environments with shielding computation techniques especially tailored for deep space mission scenario is instrumental in view of this exigency. In view of manned missions targeted to Mars, for which radiation exposure is one of the greatest problems and challenges to be tackled, it is of fundamental importance to have available a tool which allows to know which are the particle flux and spectra at any time at any point of the Martian surface. With this goal in mind, a new model for the radiation environment to be found on the planet Mars due to Galactic Cosmic Rays (GCR) has been developed. Solar modulated primary particles rescaled for Mars conditions are transported within the Martian atmosphere, with temporal properties modeled with variable timescales, down to the surface, with altitude and backscattering patterns taken into account. The tool allows analysis for manned Mars landing missions, as well as planetary science studies, e.g. subsurface water and volatile inventory studies. This Mars environmental model is available through the SIREST website, a project of NASA Langley Research Center.

G De Angelis↗

Presentation Summary: State of the AGN: Progress Toward Understanding Black Hole Accretion Processes

Accretion of plasma onto black holes power some of the most powerfulsystems in the cosmos. Supermassive black holes at the centers ofgalaxies represent the high-mass limit of these objects and so accountfor the most luminous accretors. As a result, their influence spansvast spatial and temporal scales of cosmic phenomena: intraclusterheating, intergalactic media, galactic feedback and star formation,kiloparsec-scale jets/outflows, variability over time scales of minutesto centuries, and luminous multi-wavelength electromagnetic emissionextending all the way down to its event horizon. Their intrigue isheightened by the fact that they lie at the intersection of variousphysical laws---e.g., general relativistic gravity,magnetohydrodynamics, radiation, high-energy particle physics,thermodynamics, and photo-ionization---which all must be reconciled toarrive at a fundamental understanding and probe for new physics, liketests of general relativity. These physics ingredients must beincorporated into simulations performed somehow on dynamical scalesranging from that of the event horizon to parsec-scales. Fortunately,new computational and theoretical techniques---such as GPU computing,radiation transport, and and novel gridding techniques---are enablingprogress to larger scales, more degrees of freedom, and even to binarysystems. Some of the topics we will survey include recent progress onsimulating the relationship between the disk-jet interaction, howtilted black holes behave, radiation-dominated flow, and how binaryAGN affect the standard picture of black hole accretion. Along theway, we will highlight how new technologies have enabled thesescientific rewards. Future directions and open questions will beprovided to inspire discussion and interaction during the session.

Black Holes↗

Buoyancy Effects on Flow Structure and Instability of Low-Density Gas Jets

A low-density gas jet injected into a high-density ambient gas is known to exhibit self-excited global oscillations accompanied by large vortical structures interacting with the flow field. The primary objective of the proposed research is to study buoyancy effects on the origin and nature of the flow instability and structure in the near-field of low-density gas jets. Quantitative rainbow schlieren deflectometry, Computational fluid dynamics (CFD) and Linear stability analysis were the techniques employed to scale the buoyancy effects. The formation and evolution of vortices and scalar structure of the flow field are investigated in buoyant helium jets discharged from a vertical tube into quiescent air. Oscillations at identical frequency were observed throughout the flow field. The evolving flow structure is described by helium mole percentage contours during an oscillation cycle. Instantaneous, mean, and RMS concentration profiles are presented to describe interactions of the vortex with the jet flow. Oscillations in a narrow wake region near the jet exit are shown to spread through the jet core near the downstream location of the vortex formation. The effects of jet Richardson number on characteristics of vortex and flow field are investigated and discussed. The laminar, axisymmetric, unsteady jet flow of helium injected into air was simulated using CFD. Global oscillations were observed in the flow field. The computed oscillation frequency agreed qualitatively with the experimentally measured frequency. Contours of helium concentration, vorticity and velocity provided information about the evolution and propagation of vortices in the oscillating flow field. Buoyancy effects on the instability mode were evaluated by rainbow schlieren flow visualization and concentration measurements in the near-field of self-excited helium jets undergoing gravitational change in the microgravity environment of 2.2s drop tower at NASA John H. Glenn Research Center. The jet Reynolds number was varied from 200 to 1500 and jet Richardson number was varied from 0.72 to 0.002. Power spectra plots generated from Fast Fourier Transform (FFT) analysis of angular deflection data acquired at a temporal resolution of 1000Hz reveal substantial damping of the oscillation amplitude in microgravity at low Richardson numbers (~0.002). Quantitative concentration data in the form of spatial and temporal evolutions of the instability data in Earth gravity and microgravity reveal significant variations in the jet flow structure upon removal of buoyancy forces. Radial variation of the frequency spectra and time traces of helium concentration revealed the importance of gravitational effects in the jet shear layer region. Linear temporal and spatio-temporal stability analyses of a low-density round gas jet injected into a high-density ambient gas were performed by assuming hyper-tan mean velocity and density profiles. The flow was assumed to be non parallel. Viscous and diffusive effects were ignored. The mean flow parameters were represented as the sum of the mean value and a small normal-mode fluctuation. A second order differential equation governing the pressure disturbance amplitude was derived from the basic conservation equations. The effects of the inhomogeneous shear layer and the Froude number (signifying the effects of gravity) on the temporal and spatio-temporal results were delineated. A decrease in the density ratio (ratio of the density of the jet to the density of the ambient gas) resulted in an increase in the temporal amplification rate of the disturbances. The temporal growth rate of the disturbances increased as the Froude number was reduced. The spatio-temporal analysis performed to determine the absolute instability characteristics of the jet yield positive absolute temporal growth rates at all Fr and different axial locations. As buoyancy was removed (Fr . 8), the previously existing absolute instability disappeared at all locations establhing buoyancy as the primary instability mechanism in self-excited low-density jets.

Pasumarthi, Kasyap Sriramachandra↗

Toward an Integrated Product of Higher Accuracy and Longer Time Span: Improvement of the Upcoming NASA GEWEX SRB Data Release 4 over Release 3 Based on Validation Against Surface-Based Measurements

The NASA Global Energy and Water Exchanges (GEWEX) Surface Radiation Budget (SRB) is moving toward its next generation, the Release 4 Integrated Product (Rel. 4-IP), and its shortwave algorithm preliminary results as of this writing span 34 years continuously from July 1983 to June 2017 on a quasi-equal-area grid system with 44,016 gridboxes ranging from 1 longitude times 1 latitude around the Equator to 120 longitude times 1 latitude. The data are in the form of 3-hourly, 3-hourly-monthly, daily and monthly means. The Baseline Surface Radiation Network (BSRN) is a project of the Data and Assessments Panel from the GEWEX under the umbrella of the World Climate Research Programme (WCRP). Currently in the archive there are data from 67 sites located in all seven continents spanning from 1992 to the very present. The original BSRN records are at 1-, 2-, 3- or 5-minute intervals, and sophisticated quality-control procedures have been developed to eliminate erroneous records before hourly, daily and monthly means are computed. The Pacific Marine Environmental Laboratory (PMEL) arrays of buoys in the tropical oceans, TAO/TRITON (Pacific Ocean), RAMA (Indian Ocean) and PITARA (Atlantic Ocean), measure shortwave downward fluxes using pyranometers. These buoys transmit measured data to shore in real-time via the Argos satellite system. And these buoys are major components of the CLIVAR climate analysis project and the GOOS, GCOS, and GEOSS observing systems. So far there are data from 64 buoys and the earliest data start from 2000. We have extensively validated the Rel. 4-IP preliminary results against the BSRN and PMEL 3-hourly, daily and monthly means. In this paper, we will show that the new results not only temporally span 10 years longer than its predecessor, Release 3, but agree significantly better with the BSRN data over the land and PMEL data over the ocean.

Taiping Zhang↗

Dynamical Approach Study of Spurious Numerics in Nonlinear Computations

The last two decades have been an era when computation is ahead of analysis and when very large scale practical computations are increasingly used in poorly understood multiscale complex nonlinear physical problems and non-traditional fields. Ensuring a higher level of confidence in the predictability and reliability (PAR) of these numerical simulations could play a major role in furthering the design, understanding, affordability and safety of our next generation air and space transportation systems, and systems for planetary and atmospheric sciences, and in understanding the evolution and origin of life. The need to guarantee PAR becomes acute when computations offer the ONLY way of solving these types of data limited problems. Employing theory from nonlinear dynamical systems, some building blocks to ensure a higher level of confidence in PAR of numerical simulations have been revealed by the author and world expert collaborators in relevant fields. Five building blocks with supporting numerical examples were discussed. The next step is to utilize knowledge gained by including nonlinear dynamics, bifurcation and chaos theories as an integral part of the numerical process. The third step is to design integrated criteria for reliable and accurate algorithms that cater to the different multiscale nonlinear physics. This includes but is not limited to the construction of appropriate adaptive spatial and temporal discretizations that are suitable for the underlying governing equations. In addition, a multiresolution wavelets approach for adaptive numerical dissipation/filter controls for high speed turbulence, acoustics and combustion simulations will be sought. These steps are corner stones for guarding against spurious numerical solutions that are solutions of the discretized counterparts but are not solutions of the underlying governing equations.

Yee, H. C.↗

Bird Migration Under Climate Change - A Mechanistic Approach Using Remote Sensing

The broad-scale reductions and shifts that may be expected under climate change in the availability and quality of stopover habitat for long-distance migrants is an area of increasing concern for conservation biologists. Researchers generally have taken two broad approaches to the modeling of migration behaviour to understand the impact of these changes on migratory bird populations. These include models based on causal processes and their response to environmental stimulation, "mechanistic models", or models that primarily are based on observed animal distribution patterns and the correlation of these patterns with environmental variables, i.e. "data driven" models. Investigators have applied the latter technique to forecast changes in migration patterns with changes in the environment, for example, as might be expected under climate change, by forecasting how the underlying environmental data layers upon which the relationships are built will change over time. The learned geostatstical correlations are then applied to the modified data layers.. However, this is problematic. Even if the projections of how the underlying data layers will change are correct, it is not evident that the statistical relationships will remain the same, i.e. that the animal organism may not adapt its' behaviour to the changing conditions. Mechanistic models that explicitly take into account the physical, biological, and behaviour responses of an organism as well as the underlying changes in the landscape offer an alternative to address these shortcomings. The availability of satellite remote sensing observations at multiple spatial and temporal scales, coupled with advances in climate modeling and information technologies enable the application of the mechanistic models to predict how continental bird migration patterns may change in response to environmental change. In earlier work, we simulated the impact of effects of wetland loss and inter-annual variability on the fitness of migratory shorebirds in the central fly ways of North America. We demonstrated the phenotypic plasticity of a migratory population of Pectoral sandpipers consisting of an ensemble of 10,000 individual birds in response to changes in stopover locations using an individual based migration model driven by remotely sensed land surface data, climate data and biological field data. With the advent of new computing capabilities enabled hy recent GPU-GP computing paradigms and commodity hardware, it now is possible to simulate both larger ensemble populations and to incorporate more realistic mechanistic factors into migration models. Here, we take our first steps use these tools to study the impact of long-term drought variability on shorebird survival.

Smith, James A.↗

Structural and functional analyses of human cerebral cortex using a surface-based atlas

We have analyzed the geometry, geography, and functional organization of human cerebral cortex using surface reconstructions and cortical flat maps of the left and right hemispheres generated from a digital atlas (the Visible Man). The total surface area of the reconstructed Visible Man neocortex is 1570 cm2 (both hemispheres), approximately 70% of which is buried in sulci. By linking the Visible Man cerebrum to the Talairach stereotaxic coordinate space, the locations of activation foci reported in neuroimaging studies can be readily visualized in relation to the cortical surface. The associated spatial uncertainty was empirically shown to have a radius in three dimensions of approximately 10 mm. Application of this approach to studies of visual cortex reveals the overall patterns of activation associated with different aspects of visual function and the relationship of these patterns to topographically organized visual areas. Our analysis supports a distinction between an anterior region in ventral occipito-temporal cortex that is selectively involved in form processing and a more posterior region (in or near areas VP and V4v) involved in both form and color processing. Foci associated with motion processing are mainly concentrated in a region along the occipito-temporal junction, the ventral portion of which overlaps with foci also implicated in form processing. Comparisons between flat maps of human and macaque monkey cerebral cortex indicate significant differences as well as many similarities in the relative sizes and positions of cortical regions known or suspected to be homologous in the two species.

NASA Discipline Neuroscience↗