Search NASA⌕ Search

SEARCH · Search NASA

Results for “Smoothing algorithm”

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 469 records · Page 26

Impact of Random and Periodic Surface Roughness on P- and L-band Radiometry

L-band passive microwave remote sensing is currently considered a robust technique for global monitoring of soil moisture. However, soil roughness complicates the relationship between brightness temperature and soil moisture, with current soil moisture retrieval algorithms typically assuming a constant roughness parameter globally, leading to a potential degradation in retrieval accuracy. This current investigation established a tower-based experiment site in Victoria, Australia. P-band (~40-cm wavelength/0.75 GHz) was compared with L-band (~21-cm wavelength/1.41 GHz) over random and periodic soil surfaces to determine if there is an improvement in brightness temperature simulation and soil moisture retrieval accuracy for bare soil conditions, due to reduced roughness impact when using a longer wavelength. The results showed that P-band was less impacted by random and periodic roughness than L-band, evidenced by more comparable statistics across different roughness conditions. The roughness effect from smooth surfaces (e.g., 0.8-cm root-mean-square height and 11.1-cm correlation length) could be potentially ignored at both P- and L-band with satisfactory simulation and retrieval performance. However, for rougher soil (e.g., 1.6-cm root-mean-square height and 6.8-cm correlation length), the roughness impact needed to be accounted for at both P- and L-band, with P-band observations showing less impact than L-band. Moreover, a sinusoidal soil surface with 10-cm amplitude and 80-cm period substantially impacted the brightness temperature simulation and soil moisture retrieval at both P- and L-band, which could not be fully accounted for using the SMOS and SMAP default roughness parameters. However, when retrieving roughness parameters along with soil moisture, the ubRMSE at P-band over periodic soil was improved to a similar level (0.01-0.02 m3/m3) as that of smooth flat soil (0.01 m3/m3), while L-band showed higher ubRMSE over the periodic soil (0.03-0.04 m3/m3) than over smooth flat soil (0.01 m3/m3). Accordingly, periodic roughness effects were reduced by using observations at P-band.

Soil roughness↗

An Optimal Cell Detection Technique for Automated Patch Clamping

While there are several hardware techniques for the automated patch clamping of cells that describe the equipment apparatus used for patch clamping, very few explain the science behind the actual technique of locating the ideal cell for a patch clamping procedure. We present a machine vision approach to patch clamping cell selection by developing an intelligent algorithm technique that gives the user the ability to determine the good cell to patch clamp in an image within one second. This technique will aid the user in determining the best candidates for patch clamping and will ultimately save time, increase efficiency and reduce cost. The ultimate goal is to combine intelligent processing with instrumentation and controls in order to produce a complete turnkey automated patch clamping system capable of accurately and reliably patch clamping cells with a minimum amount of human intervention. We present a unique technique that identifies good patch clamping cell candidates based on feature metrics of a cell's (x, y) position, major axis length, minor axis length, area, elongation, roundness, smoothness, angle of orientation, thinness and whether or not the cell is only particularly in the field of view. A patent is pending for this research.

McDowell, Mark↗

Calculating C-grids with fine and embedded mesh regions

A program for calculating a C-type mesh around airfoil like shapes is described. The Jameson/Caughey approach is used: a parabolic transformation coupled with a shearing transformation. The algebraic algorithm is capable of efficiently generating nearly orthogonal grids. A high degree of grid control is possible. The user may specify grid boundaries, number of grid lines, and location of (and spacing in) trailing edge and leading edge fine mesh areas. The capability of embedding fine mesh regions, for use with new adaptive grid techniques, is being developed. Grids generated by the program were used in Euler flow flow calculatons by W. Usab. Results superior to results calculated on previous O-type grids were obtained. Specifically, calculations converged faster using C-grids rather than 0-grids, total pressure loss spikes at the trailing edge of the airfoil were eliminated, and the Ni method converged with zero artificial smoothing for a subcritical case (resulting overall total pressure loss was then nearly zero). These improvements were attributed to higher grid orthogonality, especially at the trailing edge. The program itself is fairly straightforward. Roughly half of the 800 code lines are comment lines.

Loyd, B.↗

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗

Harmonization of global land use change and management for the period 850–2100 (LUH2) for CMIP6

Human land use activities have resulted in large changes to the biogeochemical and biophysical properties of the Earth's surface, with consequences for climate and other ecosystem services. In the future, land use activities are likely to expand and/or intensify further to meet growing demands for food, fiber, and energy. As part of the World Climate Research Program Coupled Model Intercomparison Project (CMIP6), the international community has developed the next generation of advanced Earth system models (ESMs) to estimate the combined effects of human activities (e.g., land use and fossil fuel emissions) on the carbon–climate system. A new set of historical data based on the History of the Global Environment database (HYDE), and multiple alternative scenarios of the future (2015–2100) from Integrated Assessment Model (IAM) teams, is required as input for these models. With most ESM simulations for CMIP6 now completed, it is important to document the land use patterns used by those simulations. Here we present results from the Land-Use Harmonization 2 (LUH2) project, which smoothly connects updated historical reconstructions of land use with eight new future projections in the format required for ESMs. The harmonization strategy estimates the fractional land use patterns, underlying land use transitions, key agricultural management information, and resulting secondary lands annually, while minimizing the differences between the end of the historical reconstruction and IAM initial conditions and preserving changes depicted by the IAMs in the future. The new approach builds on a similar effort from CMIP5 and is now provided at higher resolution (0.25°×0.25°) over a longer time domain (850–2100, with extensions to 2300) with more detail (including multiple crop and pasture types and associated management practices) using more input datasets (including Landsat remote sensing data) and updated algorithms (wood harvest and shifting cultivation); it is assessed via a new diagnostic package. The new LUH2 products contain > 50 times the information content of the datasets used in CMIP5 and are designed to enable new and improved estimates of the combined effects of land use on the global carbon–climate system.

land use land cover change↗

Convergence Acceleration of the Navier-Stokes Equations Through Time-Derivative Preconditioning

Chorin's method of artificial compressibility is extended to both compressible and incompressible fluids by using physical arguments to define artificial fluid properties that make up a local preconditioning matrix. In particular, perturbation expansions are used to provide appropriate temporal derivatives for the equations of motion at both low speeds and low Reynolds numbers. These limiting forms are then combined into a single function that smoothly merges into the physical time derivatives at high speeds so that the equations are left unchanged at transonic, high Reynolds number conditions. The effectiveness of the resulting preconditioning procedures for the Navier-Stokes equations is demonstrated for a wide speed and Reynolds number ranges by means of stability results and computational solutions. Nevertheless, the preconditioned equations sometimes fail to provide a solution for applications for which the non-preconditioned equations converge. Often this is because the reduced dissipation in the preconditioned equations results in an unsteady solution while the more dissipative non-preconditioned equations result in a steady state. Problems of this type represent a computational challenge; it is important to distinguish between non-convergence of algorithms, and the non-existence of steady state solutions.

Merkle, Charles L.↗

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗

In-Time UAV Flight-Trajectory Estimation and Tracking Using Bayesian Filters

Rapid increase of UAV operation in the next decade in areas of on-demand delivery, medical transportation services, law enforcement, traffic surveillance and several others pose potential risks to the low altitude airspace above densely populated areas. Safety assessment of airspace demands the need for a novel UAV traffic management (UTM) framework for regulation and tracking of the vehicles. Particularly for low-altitude UAV operations, quality of GPS measurements feeding into the UAV is often compromised by loss of communication link caused by presence of trees or tall buildings in proximity to the UAV flight path. Inaccurate GPS locations may yield to unreliable monitoring and inaccurate prognosis of remaining battery life and other safety metrics which rely on future expected trajectory of the UAV. This work therefore proposes a generalized monitoring and prediction methodology for autonomous UAVs using in-time GPS measurements. Firstly, a typical 4D smooth trajectory generation technique from a series of waypoint locations with associated expected times-of-arrival based on B-spline curves is presented. Initial uncertainty in the vehicle's expected cruise velocity is quantified to compute confidence intervals along the entire flight trajectory using error interval propagation approach. Further, the generated planned trajectory is considered as the prior knowledge which is updated during its flight with incoming GPS measurements in order to estimate its current location and corresponding kinematic profiles. Estimation of position is denoted in dicrete state-space representation such that position at a future time step is derived from position and velocity at current time step and expected velocity at the future time step. A linear Bayesian filtering algorithm is employed to efficiently refine position estimation from noisy GPS measurements and update the confidence intervals. Further, a dynamic re-planning strategy is implemented to incorporate unexpected detour or delay scenarios. Finally, critical challenges related to uncertainty quantification in trajectory prognosis for autonomous vehicles are identified, and potential solutions are discussed at the end of the paper. The entire monitoring framework is demonstrated on real UAV flight experiments conducted at the NASA Langley Research Center.

Banerjee, Portia↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars, and ongoing field programs such as DoE/ARM. These programs are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity.

Cahalan, Robert F.↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.

Cahalan, Robert F.↗

The Acoustic Influence of Cell Depth on the Rotordynamic Characteristics of Smooth-Rotor/Honeycomb-Stator Annular Gas Seals

A two-control volume is employed for honeycomb-stator/smooth-rotor seals, with a conventional control-volume used for the through flow and a 'capacitance accumulator' model for the honeycomb cells. The control volume for the honeycomb cells is shown to cause a dramatic reduction in the effective acoustic velocity of the main flow, dropping the lowest acoustic frequency into the frequency range of interest for rotordynamics. In these circumstances, the impedance functions for the seals can not be modeled with conventional (frequency-independent) stiffness, damping, and mass coefficients. More general transfer functions are required to account for the reaction forces, and calculated here as a lead-lag term for the direct force function and a lag term for the cross-coupled function. These first order functions are simple compared to transfer functions for magnetic bearings or foundations, For synchronous response to imbalance, they can be approximated by running-speed-dependent stiffness and damping coefficients in conventional rotordynamic codes. Correct predictions for stability and transient response will require more general algorithms, pressumably using a state-space format.

Childs, Dara W.↗

Towards Autonomous Lunar Resource Excavation via Reinforcement Learning

To continue on a sustainable and flexible path, NASA needs to address the challenge of collecting and moving large amounts of regolith at the destination. NASA’s Regolith Advanced Surface Systems Operations Robot (RASSOR) is principally designed to mine and deliver regolith for In-Situ Resource Utilization (ISRU) processing. RASSOR’s design enables it to efficiently collect and deposit regolith, return collected material for processing, and myriad related ISRU activities. To reliably perform these operations on the lunar surface, RASSOR software and sensory systems need to be robust and maximize the information extracted from a reduced sensor payload. Herein, we present preliminary findings from the Intelligent Capabilities Enhanced RASSOR project. We created reduced-order simulation environments to develop autonomous trenching controllers via reinforcement learning and prototype state estimation architectures. The goal of reinforcement learning is for an agent to learn a policy (task strategy) through interactions with an environment. When the agent performs an action, a change occurs in environment state and a numerical reward is received which informs the agent whether the action performed was good or not. Since reinforcement learning algorithms learn through trial-and-error, a simulation is a desirable first environment for development and learning. We developed two simulations, the first is a 2D excavation simulation developed to facilitate parameter selection, and a 3D simulation developed using a game physics engine, to simulate simplified soil interactions and increase the fidelity of the dynamic models of the robotic agents. The development of this 3D simulation has enabled the training of additional sensing capabilities and research both at the granular mechanics and operations levels. We experimented with various virtual sensor payloads to identify a combination that enabled efficient excavation operation and learning. Our reward function is based on how much material is excavated per step. A penalty is also received for leaving the dig site and to smooth the acceleration of the drum arms. We implemented pseudo time-of-flight sensors to report distance from each drum to ground and the height above ground which was found to be more efficient than existing solutions. Our findings suggest that reinforcement learning for autonomous operations has learned viable trenching strategies within 3000 training episodes in our simplified 2D environment and helped identify desirable sensing capabilities, arrangements, and considerations such as the positioning of time-of-flight sensors. Future work includes expanding our simulation to more complex environments and scenarios, and transfer learning from simulation to RASSOR 2.0 hardware for deployment in the Regolith Test Bin at NASA's Kennedy Space Center.

rassor↗

Flush Airdata Sensing (FADS) System Calibration Procedures and Results for Blunt Forebodies

Blunt-forebody pressure data are used to study the behavior of the NASA Dryden Flight Research Center flush airdata sensing (FADS) pressure model and solution algorithm. The model relates surface pressure measurements to the airdata state. Spliced from the potential flow solution for uniform flow over a sphere and the modified Newtonian impact theory, the model was shown to apply to a wide range of blunt-forebody shapes and Mach numbers. Calibrations of a sphere, spherical cones, a Rankine half body, and the F-14, F/A-18, X-33, X-34, and X-38 configurations are shown. The three calibration parameters are well-behaved from Mach 0.25 to Mach 5.0, an angle-of-attack range extending to greater than 30 deg, and an angle-of-sideslip range extending to greater than 15 deg. Contrary to the sharp calibration changes found on traditional pitot-static systems at transonic speeds, the FADS calibrations are smooth, monotonic functions of Mach number and effective angles of attack and sideslip. Because the FADS calibration is sensitive to pressure port location, detailed measurements of the actual pressure port locations on the flight vehicle are required and the wind-tunnel calibration model should have pressure ports in similar locations. The procedure for calibrating a FADS system is outlined.

Cobleigh, Brent R.↗

Flush Airdata Sensing (FADS) System Calibration Procedures and Results for Blunt Forebodies

Blunt-forebody pressure data are used to study the behavior of the NASA Dryden Flight Research Center flush airdata sensing (FADS) pressure model and solution algorithm. The model relates surface pressure measurements to the airdata state. Spliced from the potential flow solution for uniform flow over a sphere and the modified Newtonian impact theory, the model was shown to apply to a wide range of blunt-forebody shapes and Mach numbers. Calibrations of a sphere, spherical cones, a Rankine half body, and the F-14, F/A-18, X-33, X-34, and X-38 configurations are shown. The three calibration parameters are well-behaved from Mach 0.25 to Mach 5.0, an angle-of-attack range extending to greater than 30 deg, and an angle-of-sideslip range extending to greater than 15 deg. Contrary to the sharp calibration changes found on traditional pitot-static systems at transonic speeds, the FADS calibrations are smooth, monotonic functions of Mach number and effective angles of attack and sideslip. Because the FADS calibration is sensitive to pressure port location, detailed measurements of the actual pressure port locations on the flight vehicle are required and the wind-tunnel calibration model should have pressure ports in similar locations. The procedure for calibrating a FADS system is outlined.

Cobleigh, Brent R.↗

An Improved Retrieval of Tropospheric Nitrogen Dioxide from GOME

We present a retrieval of tropospheric nitrogen dioxide (NO2) columns from the Global Ozone Monitoring Experiment (GOME) satellite instrument that improves in several ways over previous retrievals, especially in the accounting of Rayleigh and cloud scattering. Slant columns, which are directly fitted without low-pass filtering or spectral smoothing, are corrected for an artificial offset likely induced by spectral structure on the diffuser plate of the GOME instrument. The stratospheric column is determined from NO2 columns over the remote Pacific Ocean to minimize contamination from tropospheric NO2. The air mass factor (AMF) used to convert slant columns to vertical columns is calculated from the integral of the relative vertical NO2 distribution from a global 3-D model of tropospheric chemistry driven by assimilated meteorological data (Global Earth Observing System (GEOS)-CHEM), weighted by altitude dependent scattering weights computed with a radiative transfer model (Linearized Discrete Ordinate Radiative Transfer), using local surface albedos determined from GOME observations at NO2 wavelengths. The AMF calculation accounts for cloud scattering using cloud fraction, cloud top pressure, and cloud optical thickness from a cloud retrieval algorithm (GOME Cloud Retrieval Algorithm). Over continental regions with high surface emissions, clouds decrease the AMT by 20- 30% relative to clear sky. GOME is almost twice as sensitive to tropospheric NO2 columns over ocean than over land. Comparison of the retrieved tropospheric NO2 columns for July 1996 with GEOS-CHEM values tests both the retrieval and the nitrogen oxide radical (NOx) emissions inventories used in GEOS-CHEM. Retrieved tropospheric NO2 columns over the United States, where NOx emissions are particularly well known, are within 18% of GEOS-CHEM columns and are strongly spatially correlated (r = 0.78, n = 288, p less than 0.005). Retrieved columns show more NO2 than GEOS-CHEM columns over the Transvaal region of South Africa and industrial regions of the northeast United States and Europe. They are lower over Houston, India, eastern Asia, and the biomass burning region of central Africa, possibly because of biases from absorbing aerosols.

Martin, Randall V.↗

A New Hybrid Spatio-temporal Model for Estimating Daily Multi-year PM2.5 Concentrations Across Northeastern USA Using High Resolution Aerosol Optical Depth Data

The use of satellite-based aerosol optical depth (AOD) to estimate fine particulate matter PM(sub 2.5) for epidemiology studies has increased substantially over the past few years. These recent studies often report moderate predictive power, which can generate downward bias in effect estimates. In addition, AOD measurements have only moderate spatial resolution, and have substantial missing data. We make use of recent advances in MODIS satellite data processing algorithms (Multi-Angle Implementation of Atmospheric Correction (MAIAC), which allow us to use 1 km (versus currently available 10 km) resolution AOD data.We developed and cross validated models to predict daily PM(sub 2.5) at a 1X 1 km resolution across the northeastern USA (New England, New York and New Jersey) for the years 2003-2011, allowing us to better differentiate daily and long term exposure between urban, suburban, and rural areas. Additionally, we developed an approach that allows us to generate daily high-resolution 200 m localized predictions representing deviations from the area 1 X 1 km grid predictions. We used mixed models regressing PM(sub 2.5) measurements against day-specific random intercepts, and fixed and random AOD and temperature slopes. We then use generalized additive mixed models with spatial smoothing to generate grid cell predictions when AOD was missing. Finally, to get 200 m localized predictions, we regressed the residuals from the final model for each monitor against the local spatial and temporal variables at each monitoring site. Our model performance was excellent (mean out-of-sample R(sup 2) = 0.88). The spatial and temporal components of the out-of-sample results also presented very good fits to the withheld data (R(sup 2) = 0.87, R(sup)2 = 0.87). In addition, our results revealed very little bias in the predicted concentrations (Slope of predictions versus withheld observations = 0.99). Our daily model results show high predictive accuracy at high spatial resolutions and will be useful in reconstructing exposure histories for epidemiological studies across this region.

Air pollution↗

Cloud Spatial Structure and 3D Radiative Transfer

Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape, spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars and lidars. Ongoing field programs such as DoE/ARM are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS (Earth Observing System) instruments which are beginning to provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", (G). G depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple estimates of cloud radiative properties that are based on the "Independent Pixel Approximation" or IPA. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is underway to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.gsfc.nasa.gov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate the effects of cloud inhomogeneity and 3D radiation.

Cahalan, Robert F.↗

Strawman Philosophical Guide for Developing International Network of GPM GV Sites

The creation of an international network of ground validation (GV) sites that will support the Global Precipitation Measurement (GPM) Mission's international science programme will require detailed planning of mechanisms for exchanging technical information, GV data products, and scientific results. An important component of the planning will be the philosophical guide under which the network will grow and emerge as a successful element of the GPM Mission. This philosophical guide should be able to serve the mission in developing scientific pathways for ground validation research which will ensure the highest possible quality measurement record of global precipitation products. The philosophical issues, in this regard, partly stem from the financial architecture under which the GV network will be developed, i.e., each participating country will provide its own financial support through committed institutions -- regardless of whether a national or international space agency is involved.At the 1st International GPM Ground Validation Workshop held in Abingdon, UK in November-2003, most of the basic tenants behind the development of the international GV network were identified and discussed. Therefore, with this progress in mind, this presentation is intended to put forth a strawman philosophical guide supporting the development of the international network of GPM GV sites, noting that the initial progress has been reported in the Proceedings of the 1st International GPM GV Workshop -- available online. The central philosophical issues themselves, all flow from the fact that each participating institution can only bring to the table, GV facilities and scientific personnel that are affordable to the sanctioning (funding) national agency (be that a research, research-support, or operational agency). This situation imposes on the network, heterogeneity in the measuring sensors, data collection periods, data collection procedures, data latencies, and data reporting capabilities. Therefore, in order for the network to be effective in supporting the central scientific goals of the GPM mission, there must be a basic agreed upon doctrine under which the network participants function vis-a-vis: (1) an overriding set of general scientific requirements, (2) a minimal set of policies governing the free flow of GV data between the scientific participants, (3) a few basic definitions concerning the prioritization of measurements and their respective value to the mission, (4) a few basic procedures concerning data formats, data reporting procedures, data access, and data archiving, and (5) a simple means to differentiate GV sites according to their level of effort and ability to perform near real-time data acquisition - data reporting tasks. Most important, in case they choose to operate as a near real-time data collection-data distribution site, they would be expected to operate under a fairly narrowly defined protocol needed to ensure smooth GV support operations. This presentation will suggest measures responsive to items (1) - (5) from which to proceed,. In addition, this presentation will seek to stimulate discussion and debate concerning how much heterogeneity is tolerable within the eventual GV site network, given that the any individual GV site can only be considered scientifically useful if it supports the achievement of the central GPM Mission goals. Only ground validation research that has a direct connection to the space mission should be considered justifiable given the overarching scientific goals of the mission. Therefore each site will have to seek some level of accommodation to what the GPM Mission requires in the way of retrieval error characterization, retrieval error detection and reporting, and generation of GV data products that support assessment and improvement of the mission's standard precipitation retrieval algorithms. These are all important scientific issues that will be best resolved in open scientific debate.

Smith, Eric A.↗