Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial Generative Model”

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 685 records · Page 38

Imaging the Complex Geometry of a Magma Reservoir Using Fem-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluid-filled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗

Performance evaluation of the JPL TOPSAR system: An across-track interferometric SAR system for topographic mapping

An across-track interferometric (XTI) SAR has the capability of determining the three dimensional target location from the recorded radar echos thus providing high spatial resolution both horizontally and vertically. The output of an XTI SAR will typically include both a radar backscatter map and a digital elevation model (DEM). The accuracy of the derived DEM is determined by factors such as data acquisition geometry, signal-to-noise ratio, knowledge of platform position and attitude as well as the accuracy of the processing system. In the summer of 1992 the authors acquired TOPSAR data over the Ft. Irwin area in California, a desert area with significant relief (height standard deviation 150 m). Very accurate DEM's have been derived for this area by the Topographic Engineering Center (TEC) using digital correlation methods on 1:20,000 scale digitized photographs. Corner reflectors were deployed in the area, and their locations were determined to cm accuracies by the Defense Mapping Agency (DMA) using differential GPS techniques. DEM's generated from the acquired radar data have been rotated and translated to overlay the reference DEM's provided by TEC. A detailed description of the errors and their characteristics will be given. The standard deviations of the radar data sets we acquired measured over a 5.6 by 7 km area were 1.9 and 2.3 m. The corresponding numbers for flat areas were 1.1 and 2.0 m and for mountain areas 3.3 and 2.2 m.

Madsen, Soren N.↗

A 3-level Model for Schumann-Runge O2 Laser-Induced Fluorescence

A three level model has been developed for the analysis of Schumann-Runge band (B(sup 3)Sigma(sup -)(sub u ) <- X(sup 3)Sigma(sup -)(sub g)) laser-induced fluorescence of molecular oxygen, O2. Such a model is required due to the severe lower state depletion which can occur when transitions having relatively large absorption cross-sections are excited. Such transitions are often utilized via ArF* or KrF* excimer or dye-laser excitation in high temperature environments. The rapid predissociation of the upper state prevents substantial repopulation of the lower state by collisional processes, and the lower state may be largely depleted, even at laser fluences as low as 10-100 mJ/sq cm. The resulting LIF signal in such cases no longer varies linearly with laser pulse energy, and the extent of the sublinear behavior varies with the particular rovibrational transition of interest. Relating the measured signal to the lower state population, then, necessitates the use of exceedingly low laser fluences. These low fluences in turn lead to the need to compromise spatial resolution in order to generate sufficient signal.

Diskin, Glenn S.↗

Global Burned Area and Biomass Burning Emissions from Small Fires

In several biomes, including croplands, wooded savannas, and tropical forests, many small fires occur each year that are well below the detection limit of the current generation of global burned area products derived from moderate resolution surface reflectance imagery. Although these fires often generate thermal anomalies that can be detected by satellites, their contributions to burned area and carbon fluxes have not been systematically quantified across different regions and continents. Here we developed a preliminary method for combining 1-km thermal anomalies (active fires) and 500 m burned area observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) to estimate the influence of these fires. In our approach, we calculated the number of active fires inside and outside of 500 m burn scars derived from reflectance data. We estimated small fire burned area by computing the difference normalized burn ratio (dNBR) for these two sets of active fires and then combining these observations with other information. In a final step, we used the Global Fire Emissions Database version 3 (GFED3) biogeochemical model to estimate the impact of these fires on biomass burning emissions. We found that the spatial distribution of active fires and 500 m burned areas were in close agreement in ecosystems that experience large fires, including savannas across southern Africa and Australia and boreal forests in North America and Eurasia. In other areas, however, we observed many active fires outside of burned area perimeters. Fire radiative power was lower for this class of active fires. Small fires substantially increased burned area in several continental-scale regions, including Equatorial Asia (157%), Central America (143%), and Southeast Asia (90%) during 2001-2010. Globally, accounting for small fires increased total burned area by approximately by 35%, from 345 Mha/yr to 464 Mha/yr. A formal quantification of uncertainties was not possible, but sensitivity analyses of key model parameters caused estimates of global burned area increases from small fires to vary between 24% and 54%. Biomass burning carbon emissions increased by 35% at a global scale when small fires were included in GFED3, from 1.9 Pg C/yr to 2.5 Pg C/yr. The contribution of tropical forest fires to year-to-year variability in carbon fluxes increased because small fires amplified emissions from Central America, South America and Southeast Asia-regions where drought stress and burned area varied considerably from year to year in response to El Nino-Southern Oscillation and other climate modes.

Randerson, J. T.↗

A Thermodynamic Paradigm for Studying Disease Vector’s Habitats & Life Cycles Using NASA’s NextGen Remote Sensing Instruments

Global public health is entering a new informational age through the use of spatial models of disease vector/host ecologies driven by the use of remotely sensed data to measure environmental and structural factors critical in determining disease vector habitats, distributions, life cycles, and host interactions. The vector habitat microclimates can be quantified in terms of the surface energy budget measured by satellites. The epidemiological equations (processes) can be adapted and modified to explicitly incorporate environmental factors and interfaces required by a specific disease and its vector/host cycle. Remote sensing can be used to measure or evaluate or estimate both environment (state functions) and interface (process functions). Remote sensing can be used to measure or evaluate or estimate both environment (state functions) and interface (process functions) defining vector habitats. The products of remote sensing can be integrated directly into the epidemiological equations to significantly enhance our understanding of disease vector’s life cycles and habitats. The next generation of NASA’s remote sensing instruments recently operational will provide a significant enhancement in our ability to study disease vector’s life cycles and habitats. Three new instruments have become operational on the International Space Station (ISS) In June 2018, ECOSTRESS (The ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station), a five channel, thermal IR instrument with 70 m resolution and approximately 4-5 day repeat cycle of day/night pairs along with DESISa hyperspectral sensor system 235 channels (between 400 and 1000 nanometers), 30 m resolution, and data only being acquired on demand. In November 2018. GEDI (Global Ecosystem Dynamics Investigation). GEDI is a high-resolution laser ranging of Earth’s forests and topography. These instruments set the stage for NASA’s new global designated observables set of remote sensing measurements collecting hyperspectral and multispectral thermal data called SBG (Surface Biology Geology). The SBG is currently in the planning stage. The availability and use of level 1-4 data products generated from these ISS instruments will be discussed.

thermodynamics↗

Paths and patterns: the biology and physics of swimming bacterial populations

The velocity distribution of swimming micro-organisms depends on directional cues supplied by the environment. Directional swimming within a bounded space results in the accumulation of organisms near one or more surfaces. Gravity, gradients of chemical concentration and illumination affect the motile behaviour of individual swimmers. Concentrated populations of organisms scatter and absorb light or consume molecules, such as oxygen. When supply is one-sided, consumption creates gradients; the presence of the population alters the intensity and the symmetry of the environmental cues. Patterns of cues interact dynamically with patterns of the consumer population. In suspensions, spatial variations in the concentration of organisms are equivalent to variations of mean mass density of the fluid. When organisms accumulate in one region whilst moving away from another region, the force of gravity causes convection that translocates both organisms and dissolved substances. The geometry of the resulting concentration-convection patterns has features that are remarkably reproducible. Of interest for biology are (1) the long-range organisation achieved by organisms that do not communicate, and (2) that the entire system, consisting of fluid, cells, directional supply of consumables, boundaries and gravity, generates a dynamic that improves the organisms' habitat by enhancing transport and mixing. Velocity distributions of the bacterium Bacillus subtilis have been measured within the milieu of the spatially and temporally varying oxygen concentration which they themselves create. These distributions of swimming speed and direction are the fundamental ingredients required for a quantitative mathematical treatment of the patterns. The quantitative measurement of swimming behaviour also contributes to our understanding of aerotaxis of individual cells.

Non-NASA Center↗

The 2011 Draconid Shower Risk to Earth-Orbiting Satellites

Current meteor shower forecast models project a strong Draconid outburst, possibly a storm, on October 8, 2011, with a duration of approximately 7 hours and peaking between 19 and 21 hours UT. Predicted rates span an order of magnitude, with maximum Zenithal Hourly Rates (ZHRs) ranging from a few tens to several hundred. Calibration of the NASA MSFC Meteoroid Stream Model 1 to radar and optical observations of past apparitions, particularly the 2005 Draconid outburst 2, suggest that the maximum rate will be several hundreds per hour. Given the high spatial density of the Draconid stream, this implies a maximum meteoroid flux of 5-10 Draconids km(exp -2)/hr (to a limiting diameter of 1 mm), some 25-50 times greater than the normal sporadic flux of 0.2 km(exp -2)/ hr for particles of this size. Total outburst fluence, assuming a maximum ZHR of 750, is 15.5 Draconids km(exp -2), resulting in an overall 10x risk increase to spacecraft surfaces vulnerable to hypervelocity impacts by 1 mm particles. It is now established that a significant fraction of spacecraft anomalies produced by shower meteoroids (e.g. OLYMPUS and LandSat 5) are caused by electrostatic discharges produced by meteoroid impacts. In these cases, the charge generated is roughly proportional to v(exp 3.5(4)), giving a Draconid moving at 20 km/s approximately 1/80th the electrical damage potential of a Leonid of the same mass. In other words, a Draconid outburst with a maximum ZHR of 800 presents the same electrical risk as a normal Leonid shower with a ZHR of 15, assuming the mass indices and shower durations are the same. This is supported by the fact that no spacecraft electrical anomalies were reported during the strong Draconid outbursts of 1985 and 1998. However, the lack of past anomalies should not be taken as carte blanche for satellite operators to ignore the 2011 Draconids, as the upcoming outburst will constitute a period of enhanced risk for vehicles in near-Earth space. Each spacecrft is unique, and components have differing damage thresholds; programs are encouraged to conduct analyses to determine whether or not mitigation strategies are necessary for their vehicles.

Cooke, William J.↗

Modeling and Prediction of the Noise from Non-Axisymmetric Jets

The new source model was combined with the original sound propagation model developed for rectangular jets to produce a new version of the rectangular jet noise prediction code. This code was validated using a set of rectangular nozzles whose geometries were specified by NASA. Nozzles of aspect ratios two, four and eight were studied at jet exit Mach numbers of 0.5, 0.7 and 0.9, for a total of nine cases. Reynolds-averaged Navier-Stokes solutions for these jets were provided to the contactor for use as input to the code. Quantitative comparisons of the predicted azimuthal and polar directivity of the acoustic spectrum were made with experimental data provided by NASA. The results of these comparisons, along with a documentation of the propagation and source models, were reported in a journal article publication (Ref. 4). The complete set of computer codes and computational modules that make up the prediction scheme, along with a user's guide describing their use and example test cases, was provided to NASA as a deliverable of this task. The use of conformal mapping, along with simplified modeling of the mean flow field, for noise propagation modeling was explored for other nozzle geometries, to support the task milestone of developing methods which are applicable to other geometries and flow conditions of interest to NASA. A model to represent twin round jets using this approach was formulated and implemented. A general approach to solving the equations governing sound propagation in a locally parallel nonaxisymmetric jet was developed and implemented, in aid of the tasks and milestones charged with selecting more exact numerical methods for modeling sound propagation, and developing methods that have application to other nozzle geometries. The method is based on expansion of both the mean-flowdependent coefficients in the governing equation and the Green's function in series of orthogonal functions. The method was coded and tested on two analytically prescribed mean flows which were meant to represent noise reduction concepts being considered by NASA. Testing (Ref. 5) showed that the method was feasible for the types of mean flows of interest in jet noise applications. Subsequently, this method was further developed to allow use of mean flow profiles obtained from a Reynolds-averaged Navier-Stokes (RANS) solution of the flow. Preliminary testing of the generalized code was among the last tasks completed under this contract. The stringent noise-reduction goals of NASA's Fundamental Aeronautics Program suggest that, in addition to potentially complex exhaust nozzle geometries, next generation aircraft will also involve tighter integration of the engine with the airframe. Therefore, noise generated and propagated by jet flows in the vicinity of solid surfaces is expected to be quite significant, and reduced-order noise prediction tools will be needed that can deal with such geometries. One important source of noise is that generated by the interaction of a turbulent jet with the edge of a solid surface (edge noise). Such noise is generated, for example, by the passing of the engine exhaust over a shielding surface, such as a wing. Work under this task supported an effort to develop a RANS-based prediction code for edge noise based on an extension of the classical Rapid Distortion Theory (RDT) to transversely sheared base flows (Refs. 6 and 7). The RDT-based theoretical analysis was applied to the generic problem of a turbulent jet interacting with the trailing edge of a flat plate. A code was written to evaluate the formula derived for the spectrum of the noise produced by this interaction and results were compared with data taken at NASA Glenn for a variety of jet/plate configurations and flow conditions (Ref. 8). A longer-term goal of this task was to work toward the development of a high-fidelity model of sound propagation in spatially developing non-axisymmetric jets using direct numerical methods for solving the relevant equations. Working with NASA Glenn Acoustics Branch personnel, numerical methods and boundary conditions appropriate for use in a high-resolution calculation of the full equations governing sound propagation in a steady base flow were identified. Computer codes were then written (by NASA) and tested (by OAI) for an increasingly complex set of flow conditions to validate the methods. The NASA-supplied codes were ported to the High-End Computing resources of the NASA Advanced Supercomputing facility for testing and validation against analytical (where possible) and independent numerical solutions. The cases which were completed during the course of this contract were solutions of the two-dimensional linearized Euler equations with no mean flow, a uniform mean flow and a nonuniform mean flow representative of a parallel flow jet.

Leib, Stewart J.↗

Absolute Calibration of Optical Satellite Sensors Using Libya 4 Pseudo Invariant Calibration Site

The objective of this paper is to report the improvements in an empirical absolute calibration model developed at South Dakota State University using Libya 4 (+28.55 deg, +23.39 deg) pseudo invariant calibration site (PICS). The approach was based on use of the Terra MODIS as the radiometer to develop an absolute calibration model for the spectral channels covered by this instrument from visible to shortwave infrared. Earth Observing One (EO-1) Hyperion, with a spectral resolution of 10 nm, was used to extend the model to cover visible and near-infrared regions. A simple Bidirectional Reflectance Distribution function (BRDF) model was generated using Terra Moderate Resolution Imaging Spectroradiometer (MODIS) observations over Libya 4 and the resulting model was validated with nadir data acquired from satellite sensors such as Aqua MODIS and Landsat 7 (L7) Enhanced Thematic Mapper (ETM+). The improvements in the absolute calibration model to account for the BRDF due to off-nadir measurements and annual variations in the atmosphere are summarized. BRDF models due to off-nadir viewing angles have been derived using the measurements from EO-1 Hyperion. In addition to L7 ETM+, measurements from other sensors such as Aqua MODIS, UK-2 Disaster Monitoring Constellation (DMC), ENVISAT Medium Resolution Imaging Spectrometer (MERIS) and Operational Land Imager (OLI) onboard Landsat 8 (L8), which was launched in February 2013, were employed to validate the model. These satellite sensors differ in terms of the width of their spectral bandpasses, overpass time, off-nadir-viewing capabilities, spatial resolution and temporal revisit time, etc. The results demonstrate that the proposed empirical calibration model has accuracy of the order of 3% with an uncertainty of about 2% for the sensors used in the study.

radiometric calibration↗

A Decision Support Information System for Urban Landscape Management Using Thermal Infrared Data

In this paper, we describe efforts to use remote sensing data within the purview of an information support system, to assess urban thermal landscape characteristics as a means for developing more robust models of the Urban Heat Island (UHI) effect. We also present a rationale on how we have successfully translated the results from the study of urban thermal heating and cooling regimes as identified from remote sensing data, to decision-makers, planners, government officials, and the public at large in several US cities to facilitate better understanding of how the UHI affects air quality. Additionally, through the assessment of the spatial distribution of urban thermal landscape characteristics using remote sensing data, it is possible to develop strategies to mitigate the UHI that hopefully will in turn, drive down ozone levels and improve overall urban air quality. Four US cities have been the foci for intensive analysis as part of our studies: Atlanta, GA, Baton Rouge, LA, Salt Lake City, UT, and Sacramento, CA. The remote sensing data for each of these cities has been used to generate a number of products for use by "stakeholder" working groups to convey information on what the effects are of the UHI and what measures can be taken to mitigate it. In turn, these data products are used to both educate and inform policy-makers, planners, and the general public about what kinds of UHI mitigation strategies are available.

Quattrochi, Dale A.↗

Holographic Gratings for Optical Processing

Investigation of astronomical objects and tracking of man-made space objects lead to generation of huge amount of information for optical processing. Traditional big-size optical elements (such as optical telescopes) have a tendency for increasing aperture size in order to improve sensitivity. This tendency leads to increasing of weight and costs of optical systems and stimulate search for the new, more adequate technologies. One approach to meet these demands is based on developing of holographic optical elements using new polymeric materials. We have investigated possibility to use new material PQ-PMMA (phenantrenequinone-doped PMMA (Polymethyl Methacrylate)) for fabrication of highly selective optical filters and fast spatial-temporal light modulators. This material was originally developed in Russia and later was tested in CalTech as a candidate material for optical storage. Our theoretical investigation predicts the possibility of realization of fast spatial and temporal light modulation, using volume reflection-type spectral filter. We have developed also model of holographic-grating recording in PQ-PMMA material, based on diffusional amplification. This mechanism of recording allow to receive high diffraction efficiency during recording of reflection-type volume holographic grating (holographic mirror). We also investigated recording of dynamic gratings in the photorefractive crystals LiNbO3 (LN) for space-based spectroscopy and for adaptive correction of aberrations in the telescope's mirrors. We have shown, that specific 'photogalvanic' mechanism of holographic grating recording in LN allow to realize recording of blazed gratings for volume and surface gratings. Possible applications of dynamic gratings in LN for amplification of images, transmitted through an imaging fiber guide was also demonstrated.

Nickolai Kukhtarev↗

Polarization Evolution of the Afterglow of GRB 030329

The association of a supernova with GRB 030329l strongly supports the collapsar model of gamma-ray bursts (GRBs), where a relativistic jet forms after the progenitor star collapses. Such jets cannot be spatially resolved because of their cosmological distances. Their existence is conjectured based on breaks in GRB afterglow light curves and the theoretical desire to reduce the GRB energy requirements. Temporal evolution of polarization may provide independent evidence for the jet structure of the relativistic outflow. Small-level polarization (approx. 1-3%) has been reported for a few bursts, but the temporal evolution of polarization properties could not be established. Here, we report polarimetric observations of the afterglow of GRB 030329 with high signal-to-noise and high sampling frequency. We establish the polarization light curve, detect sustained polarization at the percent level, and find significant variability. The data imply that the afterglow magnetic field has small coherence length and is mostly random, probably generated by turbulence, in contrast with the high polarization detected in the prompt gamma-rays from GRB 02120618. Our results suggest a different structure and origin of the magnetic field in the prompt vs. afterglow emission regions.

Greiner, Jochen↗

Intraseasonal Variations in Tropical Energy Balance: Relevance to Climate Sensitivity?

Intraseasonal variability of deep convection represents a fundamental mode of organization for tropical convection. While most studies of intraseasonal oscillations (ISOs) have focused on the spatial propagation and dynamics of convectively coupled circulations, here we examine the projection of ISOs on the tropically-averaged heat and moisture budget. One unresolved question concerns the degree to which observable variations in the "fast" processes (e.g. convection, radiative / turbulent fluxes) can inform our understanding of feedback mechanisms operable in the context of climate change. Our analysis use daily data from satellite observations, the Modern Era analysis for Research and Applications (MERRA), and other model integrations to address these questions: (i) How are tropospheric temperature variations related to that tropical deep convection and the associated ice cloud fractional amount (ICF), ice water path (IWP), and properties of warmer liquid clouds? (ii) What role does moisture transport play vis-a-vis ocean latent heat flux in enabling the evolution of deep convection to mediate PBL - free atmospheric temperature equilibration? (iii) What affect do convectively generated upper-tropospheric clouds have on the TOA radiation budget? Our methodology is similar to that of Spencer et al., (2007 GRL ) whereby a composite time series of various quantities over 60+ ISO events is built using tropical mean tropospheric temperature signal as a reference to which the variables are related at various lag times (from -30 to +30 days). The area of interest encompasses the global oceans between 20oN/S. The increase of convective precipitation cannot be sustained by evaporation within the domain, implying strong moisture transports into the tropical ocean area. The decrease in net TOA radiation that develops after the peak in deep convective rainfall, is part of the response that constitutes a "discharge" / "recharge" mechanism that facilitates tropical heat balance maintenance on these time scales. However, water vapor and hydrologic scaling relationships for this mode of variability cast doubt on the utility of ISO variations as proxies for climate sensitivity response to external radiatively forced (e.g. greenhouse gas-induced) climate change.

Robertson, Franklin R.↗

Review of rigorous coupled-wave analysis and of homogeneous effective medium approximations for high spatial-frequency surface-relief gratings

A review of the rigorous coupled-wave analysis as applied to the diffraction of electro-magnetic waves by gratings is presented. The analysis is valid for any polarization, angle of incidence, and conical diffraction. Cascaded and/or multiplexed gratings as well as material anisotropy can be incorporated under the same formalism. Small period rectangular groove gratings can also be modeled using approximately equivalent uniaxial homogeneous layers (effective media). The ordinary and extraordinary refractive indices of these layers depend on the gratings filling factor, the refractive indices of the substrate and superstrate, and the ratio of the freespace wavelength to grating period. Comparisons of the homogeneous effective medium approximations with the rigorous coupled-wave analysis are presented. Antireflection designs (single-layer or multilayer) using the effective medium models are presented and compared. These ultra-short period antireflection gratings can also be used to produce soft x-rays. Comparisons of the rigorous coupled-wave analysis with experimental results on soft x-ray generation by gratings are also included.

Glytsis, Elias N.↗

Decadal Prediction Skill in the GEOS-5 Forecast System

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office?s GEOS-5 Atmosphere-Ocean General Circulation Model (AOGCM). The hindcasts are initialized every December from 1959 to 2010 following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multi-variate ensemble optimal interpolation ocean and sea-ice reanalysis, and from the atmospheric reanalysis (MERRA, the Modern-Era Retrospective Analysis for Research and Applications) generated using the GEOS-5 atmospheric model. The forecast skill of a three-member-ensemble mean is compared to that of an experiment without initialization but forced with observed CO2. The results show that initialization acts to increase the forecast skill of Northern Atlantic SST compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. The annual-mean Atlantic Meridional Overturning Circulation (AMOC) index is predictable up to a 5-year lead time, consistent with the predictable signal in upper ocean heat content over the Northern Atlantic. While the skill measured by Mean Squared Skill Score (MSSS) shows 50% improvement up to 10-year lead forecast over the subtropical and mid-latitude Atlantic, however, prediction skill is relatively low in the subpolar gyre, due in part to the fact that the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region appears to be unrealistic. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.

Ham, Yoo-Geun↗

An innovations approach to decoupling of multibody dynamics and control

The problem of hinged multibody dynamics is solved using an extension of the innovations approach of linear filtering and prediction theory to the problem of mechanical system modeling and control. This approach has been used quite effectively to diagonalize the equations for filtering and prediction for linear state space systems. It has similar advantages in the study of dynamics and control of multibody systems. The innovations approach advanced here consists of expressing the equations of motion in terms of two closely related processes: (1) the innovations process e, a sequence of moments, obtained from the applied moments T by means of a spatially recursive Kalman filter that goes from the tip of the manipulator to its base; (2) a residual process, a sequence of velocities, obtained from the joint-angle velocities by means of an outward smoothing operations. The innovations e and the applied moments T are related by means of the relationships e = (I - L)T and T = (I + K)e. The operation (I - L) is a causal lower triangular matrix which is generated by a spatially recursive Kalman filter and the corresponding discrete-step Riccati equation. Hence, the innovations and the applied moments can be obtained from each other by means of a causal operation which is itself casually invertible.

Rodriguez, G.↗

Observations of polar ice using enhanced resolution microwave scatterometer data

Increasing interest in the role of ice sheets in regulating global climate has created a need for synoptic interannual recording and monitoring of the earth's major ice sheets. Based on their ability to 'see' into ice, microwave remote sensing instruments are well-suited for monitoring the polar regions. Although microwave radar has increased capabilities over some aspects of visible wavelength systems, there are tradeoffs between spatial and temporal resolution. Synthetic aperture radar (SAR) can make high resolution images, but is limited in ground coverage over polar regions because of high data rates and a lack of receiving stations. Microwave scatterometers, originally designed for measuring winds over the ocean, are capable of greater ground coverage in polar regions, but are inherently low resolution instruments. However, a recently developed resolution enhancement technique can generate enhanced resolution scatterometer images. In this paper the algorithm is applied to create enhanced resolution images from Seasat and ERS-1 scatterometer data over. Time series of medium-scale radar images for Greenland and Antarctica and the surrounding seas are presented. The growth and shrinkage of the southern polar sea-ice throughout a full annual cycle is illustrated with the seasonal cycle of patterns on the Antarctic and Greenland Ice Sheets clearly evident. Using the time series and simple scattering models, the locations of key ice facies in Greenland are determined.

Long, David G.↗

High-Fidelity Simulation of Turbulent Flow Past a Gaussian Bump

A spanwise-periodic computation of a turbulent flow past a Gaussian bump is performed in the form of a hybrid direct numerical simulation and wall-resolved large-eddy simulation. A fourth-order spatially-accurate flow solver is employed to perform the simulation, using 10.2 billion grid points for a Reynolds number of 170000 based on the bump height. The key findings from the simulation are reported in the acceleration and deceleration flow regions associated with the bump shape. Significant anisotropy in the normal Reynolds stresses, along both the wall-normal and streamwise directions, is observed within the acceleration region. The ratio between the Reynolds shear stress and turbulent kinetic energy in that region also experiences significant deviations from the norms of a zero pressure gradient turbulent boundary layer. The chosen Reynolds number generates strong flow separation in the adverse pressure gradient region, which is in contrast with a previous simulation at half the Reynolds number that only indicated incipient separation. An internal layer generated in the acceleration region evolves into a free shear layer that develops in the deceleration region and separates. Proper modeling of this inner layer appears crucial to predict the flow separation. Surface curvature effects on the attached flow development are also discussed.

Turbulent Boundary Layer↗