Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space physics”

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

Incorporating Parallel Computing into the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

Atmospheric data assimilation is a method of combining actual observations with model forecasts to produce a more accurate description of the earth system than the observations or forecast alone can provide. The output of data assimilation, sometimes called the analysis, are regular, gridded datasets of observed and unobserved variables. Analysis plays a key role in numerical weather prediction and is becoming increasingly important for climate research. These applications, and the need for timely validation of scientific enhancements to the data assimilation system pose computational demands that are best met by distributed parallel software. The mission of the NASA Data Assimilation Office (DAO) is to provide datasets for climate research and to support NASA satellite and aircraft missions. The system used to create these datasets is the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The core components of the the GEOS DAS are: the GEOS General Circulation Model (GCM), the Physical-space Statistical Analysis System (PSAS), the Observer, the on-line Quality Control (QC) system, the Coupler (which feeds analysis increments back to the GCM), and an I/O package for processing the large amounts of data the system produces (which will be described in another presentation in this session). The discussion will center on the following issues: the computational complexity for the whole GEOS DAS, assessment of the performance of the individual elements of GEOS DAS, and parallelization strategy for some of the components of the system.

Larson, Jay W.↗

The Structure and Dynamics of the Solar Corona

This report covers technical progress during the first year of the NASA Space Physics Theory contract between NASA and Science Applications International Corporation. Under this contract SAIC, the University of California, Irvine (UCI), and the Jet Propulsion Laboratory (JPL), have conducted research into theoretical modeling of active regions, the solar corona, and the inner heliosphere, using the MHD model. During the period covered by this report we have published 26 articles in the scientific literature. These publications are listed in Section 4 of this report. In the Appendix we have attached reprints of selected articles.

Mikic, Zoran↗

On the Circulation Manifold for Two Adjacent Lifting Sections

The circulation functional relative to two adjacent lifting sections is studied for two cases. In the first case we consider two adjacent circles. The circulation is computed as a function of the displacement of the secondary circle along the axis joining the two centers and of the angle of attack of the secondary circle, The gradient of such functional is computed by deriving a set of elliptic functions with respect both to their argument and to their Period. In the second case studied, we considered a wing-flap configuration. The circulation is computed by some implicit mappings, whose differentials with respect to the variation of the geometrical configuration in the physical space are found by divided differences. Configurations giving rise to local maxima and minima in the circulation manifold are presented.

Zannetti, Luca↗

A SPDS Node to Support the Systematic Interpretation of Cosmic Ray Data

The purpose of this project was to establish and maintain a Space Physics Data System (SPDS) node that supports the analysis and interpretation of current and future galactic cosmic ray (GCR) measurements by (1) providing on-line databases relevant to GCR propagation studies; (2) providing other on-line services, such as anonymous FTP access, mail list service and pointers to e-mail address books, to support the cosmic ray community; (3) providing a mechanism for those in the community who might wish to submit similar contributions for public access; (4) maintaining the node to assure that the databases remain current; and (5) investigating other possibilities, such as CD-ROM, for public dissemination of the data products. Shortly after the original grant to support these activities was established at Louisiana State University a detailed study of alternate choices for the node hardware was initiated. The chosen hardware was an Apple Workgroup Server 9150/120 consisting of a 120 MHz PowerPC 601 processor, 32 MB of memory, two I GB disks and one 2 GB disk. This hardware was ordered and installed and has been operating reliably ever since. A preliminary version of the database server was available during the first year effort and was used as part of the very successful SPDS demonstration during the Rome, Italy International Cosmic Ray Conference. For this server version we were able to establish the html and anonymous FTP server software, develop a Web page structure which can be easily modified to include new items, provide an on-line database of charge changing total cross sections, include the cross section prediction software of Silberberg & Tsao as well as Webber, Kish and Schrier for download access, and provide an on-line bibliography of the cross section measurement references by the Transport Collaboration. The preliminary version of this SPDS Cosmic Ray node was examined by members of the C&H SPDS committee and returned comments were used to refine the implementation.

Source record↗

Orion Nebula and Planetary Nebulae

This report summarizes the research performed at Rice University related to NASA-Ames University consortium grant NCC2-5199 during the two year period 1996 September 1 through 1998 August 31. The research program, titled Orion Nebula and Planetary Nebulae, involved the analysis of Hubble Space Telescope (HST) imagery and spectroscopy of the Orion Nebula and of the planetary nebulae NGC 6818 and NGC 6210. In addition, we analyzed infrared spectra of the Orion Nebula taken with the Infrared Space Observatory (ISO) The primary collaborators at NASA-Ames were Drs. R. H. Rubin, A. G. C. M. Tielens, S. W. J. Colgan, and S. D. Lord (Tielens & Lord has since changed institutions). Other collaborators include Drs. P. G. Martin (CITA, Toronto), G. J. Ferland (U. KY), J. A. Baldwin (CTIO, Chile), J. J. Hester (ASU), D. K. Walter (SCSU), and P. Harrington (U. MD). In addition to the Principal Investigator, Professor Reginald J. Dufour of the Department of Space Physics & Astronomy, the research also involved two students, Mr. Matthew Browning and Mr. Brent Buckalew. Mr. Browning will be graduating from Rice in 1999 May with a B.A. degree in Physics and Mr. Buckalew continues as a graduate student in our department, having recently received a NASA GSRP research fellowship (sponsored by Ames). The collaboration was very productive, with two refereed papers already appearing in the literature, several others in preparation, numerous meeting presentations and two press releases. Some of our research accomplishments are highlighted below. Attached to the report are copies of the two major publications. Note that this research continues to date and related extensions of it recently has been awarded time with the HST for 1999-2000.

Dufour, Reginald J.↗

Impact of a Merged Precipitation Data on Global Soil Moisture Variability

Accurate soil moisture information has proved to be important to climate simulations and climate and weather forecasts. However, many difficulties exist that limit our understanding of soil moisture distribution and variability. One of them is the lack of accurate precipitation with appropriate spatial and temporal resolution. Precipitation as an input forcing to the land surface greatly influences soil moisture characteristics and variability. To improve precipitation data quality, an algorithm has been developed to generate a spatially and temporally continuous 3-hourly global precipitation data for the period of 1987 to present. This precipitation product is a combination of the precipitation from Special Sensor Microwave Imager (SSMI) with the Goddard Earth Observing System-1 Data Assimilation System (GEOS-1 DAS) employing a Physical-space Statistical Analysis System (PSAS). In this study we investigate the impact of this merged/analyzed precipitation data on the global soil moisture variability using an Off-line Land-surface GEOS Assimilation (OLGA) system. Two OLGA integrations starting from 1987 to 1993 are performed forced with the analyzed and GEOS-1 DAS precipitation respectively. We examine the spatial and temporal characteristics of soil moisture variability in response to the analyzed precipitation. The influence of this merged precipitation on the soil moisture variability and regional hydrological budget is estimated throughout the comparison with the results forced with the GEOS-1 DAS precipitation only. In the OLGA the sut@-grid scale horizontal heterogeneity is explicitly represented on the tile space. This provides a means to assess the role of the surface moisture heterogeneity in the interaction with the surface atmosphere and surface hydrological budget, and to validate OLGA results at tile space with in situ observation. ABRACOS (Anglo-Brazilian Amazonian Climate Observation Study), FIFE (First ISLSCP Field Experiment) I and HAPEX data will be used in the validation.

Yang, Runhua↗

Mercury's Magnetosphere

Among the major discoveries made by the Mariner 10 mission to the inner planets was the existence of an intrinsic magnetic field at Mercury with a dipole moment of approx. 300 nT R(sup 3, sub M). This magnetic field is sufficient to stand off the solar wind at an altitude of about 1 R(sub M) (i.e. approx. 2439 km). Hence, Mercury possesses a 'magnetosphere' from which the so]ar wind plasma is largely excluded and within which the motion of charged particles is controlled by the planetary magnetic field. Despite its small size relative to the magnetospheres of the other planets, a Mercury orbiter mission is a high priority for the space physics community. The primary reason for this great interest is that Mercury unlike all the other planets visited thus far, lacks a significant atmosphere; only a vestigial exosphere is present. This results in a unique situation where the magnetosphere interacts directly with the outer layer of the planetary crust (i.e. the regolith). At all of the other planets the topmost regions of their atmospheres become ionized by solar radiation to form ionospheres. These planetary ionospheres then couple to electrodynamically to their magnetospheres or, in the case of the weakly magnetized Venus and Mars, directly to the solar wind. This magnetosphere-ionosphere coupling is mediated largely through field-aligned currents (FACs) flowing along the magnetic field lines linking the magnetosphere and the high-latitude ionosphere. Mercury is unique in that it is expected that FACS will be very short lived due to the low electrical conductivity of the regolith. Furthermore, at the earth it has been shown that the outflow of neutral atmospheric species to great altitudes is an important source of magnetospheric plasma (following ionization) whose composition may influence subsequent magnetotail dynamics. However, the dominant source of plasma for most of the terrestrial magnetosphere is the 'leakage'of solar wind across the magnetopause and more direct entry through the northern and southern cusps. Although Mariner 10 did not return plasma composition measurements, the Hermean magnetosphere should be ideal for measuring the manner and rate of solar wind plasma entry due to the lack of strong internal atmospheric sources. Finally, the solar wind conditions experienced by Mercury as it orbits the Sun at 0.31 to 0.47 AU are quite different from those typically encountered by the Earth. This may allow for new understanding of the external factors affecting the transfer of mass, momentum and energy from the solar wind to planetary magnetospheres. This article provides a brief overview of what is now known about Mercury's magnetosphere and why it is a priority target for future planetary missions.

Slavin, J. A.↗

Massively Parallel Assimilation of TOGA/TAO and Topex/Poseidon Measurements into a Quasi Isopycnal Ocean General Circulation Model Using an Ensemble Kalman Filter

A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Keppenne, Christian L.↗

Continuation of the University of Michigan Sounding Rocket Program

This interim report gives the activity on this grant through 31 October 1998, and the planned activities for the next fiscal year. A Project Initiation Conference (PIC) was held at the Wallops Flight Facility (WFF), with planning for the launch of flight 36.174 in Fall 1998. However, following the PIC, the WFF personnel decided that the experiment telemetry system had to be completely rewired, which required substantial design and wiring changes on the experiment side of the payload. These costs had not been anticipated in the FY98 budget, and used up considerable resources of the experimenters. In addition, the telemetry re-design was not completed until Fall 1998, so that there was no opportunity for a launch in 1998 (the Jupiter window closes in late November 1998 and it will reopen in late summer 1999). A Design Review was held by telecon with WFF personnel on 26 October 1998, with plans for the 36.174 launch now in late August 1999 with a launch window which extends into winter 2000. On a more positive note, we have been developing a vacuum UV calibration chamber in Michigan's Space Physics Research Lab, for future measurements of our experiment components, and this facility is now undergoing initial measurements of various calibration lamps in preparation for testing of the telescope UV reflective efficiency. This facility will greatly increase our ability to optimize our experiment for UV efficiency, which has been lacking in the past, and it will be made available for support of other NASA missions at Michigan and for visitors from other institutions. For the next flight 36.174, the following activities will be accomplished: (1) prepare the experiment, including replacement and upgrading of critical components, at University of Michigan; (2) pre-integrate the payload at the Wallops Flight Facility; (3) perform final far-UV calibration of instrument (4) perform final alignment, integration, and electrical checks at WSMR; (5) launch payload from WSMR and check condition of recovered payload; (6) perform post-flight calibration, if necessary and applicable (7) reduce and analyze flight data. The following pages contain information from the experimenter package for the Design Review for the next flight, 36.174, listing the plans, changes, and schedule for this flight.

Clarke, John↗

A Global Data Assimilation System for Atmospheric Aerosol

We will give an overview of an aerosol data assimilation system which combines advances in remote sensing of atmospheric aerosols, aerosol modeling and data assimilation methodology to produce high spatial and temporal resolution 3D aerosol fields. Initially, the Goddard Aerosol Assimilation System (GAAS) will assimilate TOMS, AVHRR and AERONET observations; later we will include MODIS and MISR. This data assimilation capability will allows us to integrate complementing aerosol observations from these platforms, enabling the development of an assimilated aerosol climatology as well as a global aerosol forecasting system in support of field campaigns. Furthermore, this system provides an interactive retrieval framework for each aerosol observing satellites, in particular TOMS and AVHRR. The Goddard Aerosol Assimilation System (GAAS) takes advantage of recent advances in constituent data assimilation at DAO, including flow dependent parameterizations of error covariances and the proper consideration of model bias. For its prognostic transport model, GAAS will utilize the Goddard Ozone, Chemistry, Aerosol, Radiation and Transport (GOCART) model developed at NASA/GSFC Codes 916 and 910.3. GOCART includes the Lin-Rood flux-form, semi-Langrangian transport model with parameterized aerosol chemistry and physical processes for absorbing (dust and black carbon) and non-absorbing aerosols (sulfate and organic carbon). Observations and model fields are combined using a constituent version of DAO's Physical-space Statistical Analysis System (PSAS), including its adaptive quality control system. In this talk we describe the main components of this assimilation system and present preliminary results obtained by assimilating TOMS data.

daSilva, Arlindo↗

Computer Modeling of Sand Transport on Mars Using a Compart-Mentalized Fluids Algorithm (CFA)

It has been postulated that aeolian transport on Mars may be significantly different from that on Earth. From laboratory experiments simulating martian grain transport [2], it has been observed that (saltating) grains striking the bed can cause hundreds of secondary reptation trajectories when impact occurs at speeds postulated for Mars. Some of the ballistically induced trajectories "die ouf' and effectively join the ranks on the creep population that is merely nudged along by impact. Many of the induced reptation trajectories, however, are sufficiently high for the grains to become part of the saltation load (it is irrelevant to the boundary layer how a grain attained its initial lift force). When these grains, in turn, strike the surface, they too are capable of inducing more reptating grains. This cascading effect has been discussed in connection with terrestrial aeolian transport in an attempt to dispel the notion that sand motion is divisible only into creep and saltation loads. On Earth, only a few grains are splashed by impact. On Mars, it may be hundreds. We developed a computer model to address this phenomenon because there are some important ramifications: First, this ratio may mean that martian aeolian transport is dominated by reptation flux rather than saltation. On Earth, the flux would be a roughly balanced mixture between reptation/creep and saltation. On Venus, there would be no transport other than by saltation. In other words, an understanding of planetary aeolian processes may not be necessarily understood by extrapolating from the "Earth case", with only gravity and atmospheric density/viscosity being considered as variables. Second, the reptation flux on Mars may be self sustaining, so that little input is required by the wind once transport has been initiated. The number of grains saturating the boundary layer near the bed may mean that average grain speed on Mars might conceivably be less than that on Earth. This would say much for models of sand comminution on Mars. A multiple-grain transport model using just the equations of grain motion describing lift and drag is impossible to develop owing to stochastic effects --the very effects we wish to model. Also, unless we were to employ supercomputing techniques and extremely complex computer codes that could deal with millions of grains simultaneously, it would also be difficult to model grain transport if we attempted to consider every grain in motion. No existing computer models were found that satisfactorily used the equations of motion to arrive at transport flux numbers for the different populations of saltation and reptation. Modeling all the grains in a transport system was an intractable problem within our resources, and thus we developed what we believe to be a new modeling approach to simulating grain transport. The CFA deals with grain populations, but considers them to belong to various compartmentalized fluid units in the boundary layer. In this way, the model circumvents the multigrain problem by dealing primarily with the consequences of grain transport --momentum transfer between air and grains, which is the physical essence of a dynamic grain-fluid mixture. We thus chose to model the aeolian transport process as a superposition of fluids. These fluids include the air as well as particle populations of various properties. The prime property distinguishing these fluids is upward and downward grain motion. In a normal saltation trajectory, a grain's downwind velocity increases with time, so a rising grain will have a smaller downwind velocity than a failing grain. Because of this disparity in rising and falling grain proper-ties, it seemed appropriate to track these as two separate grain populations within the same physical space. The air itself can be considered a separate fluid superimposed within and interacting with the various grain-cloud "fluids". Additional informaiton is contained in the original.

Marshall, J.↗

The Structure and Dynamics of the Solar Corona

This report covers technical progress during the third year of the NASA Space Physics Theory contract "The Structure and Dynamics of the Solar Corona," between NASA and Science Applications International Corporation, and covers the period June 16, 1998 to August 15, 1999. This is also the final report for this contract. Under this contract SAIC, the University of California, Irvine (UCI), and the Jet Propulsion Laboratory (JPL), have conducted research into theoretical modeling of active regions, the solar corona, and the inner heliosphere, using the MHD model. During the three-year duration of this contract we have published 49 articles in the scientific literature. These publications are listed in Section 3 of this report. In the Appendix we have attached reprints of selected articles. We summarize our progress during the third year of the contract. Full descriptions of our work can be found in the cited publications, a few of which are attached to this report.

Mikic, Zoran↗

CUSP Energetic Particles: Confinement, Acceleration and Implications

The cusp energetic particle (CEP) event is a new magnetospheric phenomenon. The events were detected in the dayside cusp for hours, in which the measured helium ions had energies up to 8 MeV. All of these events were associated with a dramatic decrease and large fluctuations in the local magnetic field strength. During January 1999 - December 1999 covered by this report, I have studied the CEP events by analyzing the POLAR, GEOTAIL, and WIND particle and magnetic field data measured during the geomagnetic quiet periods in 1996 and one geomagnetic storm period in 1998. The simultaneous observations indicated that the ion fluxes in the CEP events were higher than that in both the upstream and the downstream from the bow shock. The pitch angle distribution of the helium ions in the CEP events was found to peak around 90 deg. It was found that the mirror parameter, defined as the ratio of the square root of the integration of the parallel turbulent power spectral component over the ultra-low frequency (ULF) ranges to the mean field in the cusp, is correlated with the intensity of the cusp MeV helium flux, which is a measure of the influence of mirroring interactions and an indication of local effect. It was also found that the turbulent power of the local magnetic field in the ultra-low frequency (ULF) ranges is correlated with the intensity of the cusp energetic helium ions. Such ULF ranges correspond to periods of about 0.33-500 seconds that cover the gyroperiods, the bounce periods, and the drift periods of the tens keV to MeV charged particles when they are temporarily confined in the high-altitude dayside cusp. These observations represent a discovery that the high-altitude dayside cusp is a new acceleration and dynamic trapping region of the magnetosphere. The cusp geometry is connected via gradient and curvature drift of these energized ions to the equatorial plasma sheet as close as the geostationary orbit at local midnight. It implies that the dayside cusp is potentially an important source of magnetospheric particles. The discovery of the CEP events has been recognized as one of the most significant results from POLAR. I was invited to give a talk at 1999 IUGG meeting to interpret the CEP events. This discovery has also been written into the web-based Space Physics Text Book (http://www.oulu.fi/- spaceweb/textbook/cusp.html).

Chen, Jiasheng↗

On the Use of CAD-Native Predicates and Geometry in Surface Meshing

Several paradigms for accessing CAD geometry during surface meshing for CFD are discussed. File translation, inconsistent geometry engines and non-native point construction are all identified as sources of non-robustness. The paper argues in favor of accessing CAD parts and assemblies in their native format, without translation, and for the use of CAD-native predicates and constructors in surface mesh generation. The discussion also emphasizes the importance of examining the computational requirements for exact evaluation of triangulation predicates during surface meshing. The native approach is demonstrated through an algorithm for the generation of closed manifold surface triangulations from CAD geometry. CAD parts and assemblies are used in their native format, and a part's native geometry engine is accessed through a modeler-independent application programming interface (API). In seeking a robust and fully automated procedure, the algorithm is based on a new physical space manifold triangulation technique specially developed to avoid robustness issues associated with poorly conditioned mappings. In addition, this approach avoids the usual ambiguities associated with floating-point predicate evaluation on constructed coordinate geometry in a mapped space. The technique is incremental, so that each new site improves the triangulation by some well defined quality measure. The algorithm terminates after achieving a prespecified measure of mesh quality and produces a triangulation such that no angle is less than a given angle bound, a or greater than pi - 2alpha. This result also sets bounds on the maximum vertex degree, triangle aspect-ratio and maximum stretching rate for the triangulation. In addition to the output triangulations for a variety of CAD parts, the discussion presents related theoretical results which assert the existence of such an angle bound, and demonstrate that maximum bounds of between 25 deg and 30 deg may be achieved in practice.

Aftosmis, M. J.↗

Preliminary Results from an Assimilation of TOMS Aerosol Observations Into the GOCART Model

At NASA Goddard we are developing a global aerosol data assimilation system that combines advances in remote sensing and modeling of atmospheric aerosols. The goal is to provide high resolution, 3-D aerosol distributions to the research community. Our first step is to develop a simple assimilation system for Saharan mineral aerosol. The Goddard Chemistry and Aerosol Radiation model (GOCART) provides accurate 3-D mineral aerosol size distributions that compare well with TOMS satellite observations. Surface, mobilization, wet and dry deposition, convective and long-range transport are all driven by assimilated fields from the Goddard Earth Observing System Data Assimilation System, GEOS-DAS. Our version of GOCART transports sizes from.08-10 microns and only simulates Saharan dust. TOMS radiance observations in the ultra violet provide information on the mineral and carbonaceous aerosol fields. We use two main observables in this study: the TOMS aerosol index (AI) which is directly related to the ratio of the 340 and 380 radiances and the 380 radiance. These are sensitive to the aerosol optical thickness, the single scattering albedo and the height of the aerosol layer. The Goddard Aerosol Assimilation System (GAAS) uses the Data Assimilation Office's Physical-space Statistical Analysis System (PSAS) to combine TOMS observations and GOCART model first guess fields. At this initial phase we only assimilate observations into the the GOCART model over regions of Africa and the Atlantic where mineral aerosols dominant and carbonaceous aerosols are minimal, Our preliminary results during summer show that the assimilation with TOMS data modifies both the aerosol mass loading and the single scattering albedo. Assimilated aerosol fields will be compared with assimilated aerosol fields from GOCART and AERONET observations over Cape Verde.

daSilva, Arlindo↗

The DAO Retrospective Data Assimilation System

The fixed-lag Kalman smoother of Cohn et al. (1994), or a version of a more computationally feasible approximation of it developed by Todling et al. (1998), is under implementation at the Data Assimilation Office (DAO) to build a retrospective data assimilation system (RDAS). The initial version of the RDAS uses the physical-space statistical analysis system and a modification of the incremental analysis update (IAU) procedure of Bloom et al. (1996); a follow up version of the RDAS Nill require the use of the adjoint of the DAO general circulation model (developed by Y. Yang and M Navon). The retrospective procedure is designed to produce improved analyses as well as improved assimilated fields consequently providing an improved climate representation through data assimilation.

Todling, Richardo↗

DAO's Finite-volume/Physical-space Data Assimilation System: Stratospheric Applications

In this talk we describe the next-generation data assimilation system being developed at NASA's Data Assimilation Office (DAO), with emphasis on the applications to stratospheric forecasts and stratospheric constituent transport. This data assimilation system includes the General Circulation Model jointly developed by the DAO and the Climate and Global Dynamics Division (CGDD) at NCAR. This model is based on the finite-volume dynamical core) developed at DAO with physical parameterizations from the NCAR Climate Community Model. The Physical-space Statistical Analysis System (PSAS) is used to combine a first guess from the NASA-NCAR GCM with observational data to provide an updated estimate of the state of the atmosphere. Case studies for the (northern) Winter of 2000 will be discussed.

daSilva, Arlindo↗

Preliminary Results from an Assimilation of Saharan Dust Using TOMS Radiances and the GOCART Model

At NASA Goddard we are developing a global aerosol data assimilation system that combines advances in remote sensing and modeling of atmospheric aerosols. The goal is to provide high resolution, 3-D aerosol distributions to the research community. Our first step is to develop a simple assimilation system for Saharan mineral aerosol. The Goddard Chemistry and Aerosol Radiation model (GOCART) provides accurate 3-D mineral aerosol size distributions. Surface mobilization, wet and dry deposition, convective and long-range transport are all driven by assimilated fields from the Goddard Earth Observing System Data Assimilation System, GEOS-DAS. Our version of GOCART transports sizes from .08-10 microns and only simulates Saharan dust. We draw the assimilation to two observables in this study: the TOMS aerosol index (Al) which is directly related to the ratio of the 340 and 380 radiances and the 380 radiance alone. The forward model that simulates the observables requires the aerosol optical thickness, the single scattering albedo and the height of the aerosol layer from the GOCART fields. The forward model also requires a refractive index for the dust. We test three index values to see which best fits the TOMS observables. These are 1) for Saharan dust reported by Patterson, 2) for a mixture of Saharan dust and a highly reflective material (sea salt or sulfate) and 3) for pure illite. The assimilation works best assuming either pure illite or the dust mixture. Our assimilation cycle first determines values of the aerosol index (Al) and the radiance at 380 nm based on the GOCART aerosol fields. Differences between the observed and GOCART model calculated Al and 380 nm radiance are first analyzed horizontally using the Physical-space Statistical Analysis System (PSAS). A quasi-Newton iteration is then performed to produce analyzed 3D aerosol fields according to parameterized background and observation error covariances. We only assimilate observations into the the GOCART model over regions of Africa and the Atlantic where mineral aerosols are dominant and carbonaceous aerosols are minimal.

Weaver, C. J.↗