Search NASA⌕ Search

SEARCH · Search NASA

Results for “Distributed model predictive control”

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 235 records · Page 13

New York Ecological Forecasting: Utilizing NASA Earth Observations to Map Ash Distribution and Inform Emerald Ash Borer Control

Since their first sightings in the U.S. in 2002, emerald ash borer beetles (Agrilus planipennis; EAB) have killed millions of native ash (Fraxinus spp.) trees across 35 states. Infected ash stands frequently exhibit complete mortality, with the predicted result being the functional extinction of native ash in U.S. forests. In August of 2020, EAB was discovered in the 6.1-million-acre Adirondack Park. The team’s partners at the Adirondack Park Invasive Plant Program (APIPP) desired ash tree distribution and EAB susceptibility information to help improve EAB bio-control efficiency and apply the methodology to future invasive programs. To assist, the team mapped ash tree distribution using NASA Earth observations from Landsat 7 Enhanced Thematic Mapper Plus (ETM+) and Shuttle Radar Topography Mission (SRTM), along with hyperspectral imagery from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). Field data from the Monitoring and Managing Ash (MaMA) project, iMapInvasives and iNaturalist databases, and the New York State Department of Environmental Conservation (NYSDEC) provided ground truthing for mapping and modeling. Results indicate that for ash detection, the team’s Spectral Angle Mapping (SAM) hyperspectral classification is slightly more sensitive but less accurate than multispectral Random Forest (RF) classification, though neither method was above a ~20% detection rate. End products include maps of ash extent derived from both imagery types, a model forecasting future spread scenarios based on current EAB presence, and outreach materials. These products inform APIPP’s management decisions and facilitate public awareness of EAB’s threat to communities within the region.

Liam Megraw↗

Simulations of Hurricane Katrina (2005) with the 0.125 degree finite-volume General Circulation Model on the NASA Columbia Supercomputer

Hurricane Katrina was the sixth most intense hurricane in the Atlantic. Katrina's forecast poses major challenges, the most important of which is its rapid intensification. Hurricane intensity forecast with General Circulation Models (GCMs) is difficult because of their coarse resolution. In this article, six 5-day simulations with the ultra-high resolution finite-volume GCM are conducted on the NASA Columbia supercomputer to show the effects of increased resolution on the intensity predictions of Katrina. It is found that the 0.125 degree runs give comparable tracks to the 0.25 degree, but provide better intensity forecasts, bringing the center pressure much closer to observations with differences of only plus or minus 12 hPa. In the runs initialized at 1200 UTC 25 AUG, the 0.125 degree simulates a more realistic intensification rate and better near-eye wind distributions. Moreover, the first global 0.125 degree simulation without convection parameterization (CP) produces even better intensity evolution and near-eye winds than the control run with CP.

Shen, B.-W.↗

Hypervelocity impact facility for simulating materials exposure to impact by space debris

As a result of man's venturing into space, the local debris contributed by his presence exceeds, at some orbital altitudes, that of the natural component. Man's contribution ranges from fuel residue to large derelect satellites that weigh many kilograms. Current debris models are able to predict the growth of the problem and suggest that spacecraft must employ armor or bumper shields for some orbital altitudes now, and that, the problem will become worse as a function of time. The practical upper limit to the velocity distribution is on the order of 40 km/s and is associated with the natural environment. The maximum velocity of the man-made component is in the 14-16 km/s range. The Long Duration Exposure Facility (LDEF) has verified that the 'high probability of impact' particles are in the microgram to milligram range. These particles can have significant effects on coatings, insulators, and thin metallic layers. The surface of thick materials becomes pitted and the local debris component is enhanced by ejecta from the debris spectrum in a controlled environment. The facility capability is discussed in terms of drive geometry, energetics, velocity distribution, diagnostics, and projectile/debris loading. The facility is currently being used to study impact phenomena on Space Station Freedom's solar array structure, other solar array materials, potential structural materials for use in the station, electrical breakdown in the space environment, and as a means of clarifying or duplicating the impact phenomena on the LDEF surfaces. The results of these experiments are described in terms of the mass/velocity distribution incident on selected samples, crater dynamics, and sample geometry.

Rose, M. F.↗

Aero-Propulsive Damping Characterization for eVTOL Aircraft Using Free Motion Wind-Tunnel Testing

This paper describes an electric vertical takeoff and landing (eVTOL) aircraft system identification method applied using three degree-of-freedom (3DOF) free motion wind-tunnel testing. The approach, similar to flight-test system identification, allows for efficient mathematical model development of the aero-propulsive moments applied on an eVTOL vehicle, including aerodynamic damping effects. The approach is demonstrated using a subscale tiltrotor eVTOL aircraft mounted on a new 3DOF wind-tunnel apparatus. To execute the test, a model-based 3DOF control system is designed to track attitude commands and transition the aircraft based on the freestream dynamic pressure. While the flight controller is active, orthogonal phase-optimized multisine inputs are injected into the attitude command and control effector command signals to enable collection of informative data for model identification. Aero-propulsive models are then identified at several reference conditions in the transition flight envelope using the equation-error method in the frequency domain. The identified models are shown to have a good fit to the modeling data and good prediction capability of data not used for model identification. The method yields aerodynamic damping estimates using less wind-tunnel test time compared to traditional forced oscillation experiments and supplements static wind-tunnel testing to produce a comprehensive transition aero-propulsive model suitable for use in flight dynamics simulations.

RAVEN↗

Predicting System Accidents with Model Analysis During Hybrid Simulation

Standard discrete event simulation is commonly used to identify system bottlenecks and starving and blocking conditions in resources and services. The CONFIG hybrid discrete/continuous simulation tool can simulate such conditions in combination with inputs external to the simulation. This provides a means for evaluating the vulnerability to system accidents of a system's design, operating procedures, and control software. System accidents are brought about by complex unexpected interactions among multiple system failures , faulty or misleading sensor data, and inappropriate responses of human operators or software. The flows of resource and product materials play a central role in the hazardous situations that may arise in fluid transport and processing systems. We describe the capabilities of CONFIG for simulation-time linear circuit analysis of fluid flows in the context of model-based hazard analysis. We focus on how CONFIG simulates the static stresses in systems of flow. Unlike other flow-related properties, static stresses (or static potentials) cannot be represented by a set of state equations. The distribution of static stresses is dependent on the specific history of operations performed on a system. We discuss the use of this type of information in hazard analysis of system designs.

Malin, Jane T.↗

Development of a Common Research Model for Applied CFD Validation Studies

The development of a wing/body/nacelle/pylon/horizontal-tail configuration for a common research model is presented, with focus on the aerodynamic design of the wing. Here, a contemporary transonic supercritical wing design is developed with aerodynamic characteristics that are well behaved and of high performance for configurations with and without the nacelle/pylon group. The horizontal tail is robustly designed for dive Mach number conditions and is suitably sized for typical stability and control requirements. The fuselage is representative of a wide/body commercial transport aircraft; it includes a wing-body fairing, as well as a scrubbing seal for the horizontal tail. The nacelle is a single-cowl, high by-pass-ratio, flow-through design with an exit area sized to achieve a natural unforced mass-flow-ratio typical of commercial aircraft engines at cruise. The simplicity of this un-bifurcated nacelle geometry will facilitate grid generation efforts of subsequent CFD validation exercises. Detailed aerodynamic performance data has been generated for this model; however, this information is presented in such a manner as to not bias CFD predictions planned for the fourth AIAA CFD Drag Prediction Workshop, which incorporates this common research model into its blind test cases. The CFD results presented include wing pressure distributions with and without the nacelle/pylon, ML/D trend lines, and drag-divergence curves; the design point for the wing/body configuration is within 1% of its max-ML/D. Plans to test the common research model in the National Transonic Facility and the Ames 11-ft wind tunnels are also discussed.

Vassberg, John C.↗

The Sensitivity of African Easterly Waves to Eastern Tropical Atlantic Sea-Surface Temperatures

The results of two regional atmospheric model simulations are compared to assess the influence of the eastern tropical Atlantic sea-surface temperature maximum on local precipitation, transient easterly waves and the West African summer monsoon. Both model simulations were initialized with reanalysis 2 data (US National Center for Environmental Prediction and Department of Energy) on 15 May 2006 and extended through 6 October 2006, forced by synchronous reanalysis 2 lateral boundary conditions introduced four times daily. One simulation uses 2006 reanalysis 2 sea-surface temperatures, also updated four times daily, while the second simulation considers ocean forcing absent the sea-surface temperature maximum, achieved here by subtracting 3 K at every ocean grid point between 0 and 15 N during the entire simulation. The simulation with 2006 sea-surface temperature forcing produces a realistic distribution of June-September mean precipitation and realistic westward propagating swaths of maximum rainfall, based on validation against Tropical Rainfall Measuring Mission (TRMM) estimates. The simulation without the sea-surface temperature maximum produces only 57% of the control June-September total precipitation over the eastern tropical Atlantic and about 83% of the Sahel precipitation. The simulation with warmer ocean temperatures generates generally stronger circulation, which in turn enhances precipitation by increasing moisture convergence. Some local precipitation enhancement is also attributed to lower vertical thermal stability above the warm water. The study shows that the eastern tropical Atlantic sea-surface temperature maximum enhances the strength of transient easterly waves and broadens the spatial extent of associated precipitation. However, large-scale circulation and its interaction with the African continent, and not sea-surface temperatures, control the timing and trajectories of the waves.

Druyan, Leonard M.↗

Solar Cycle/Seasonal Variations of H, D, H2 and He Distributions and Escape on Mars as Determined by the Mars Thermosphere Global Circulation Model (MTGCM)

The Mars Thermosphere General Circulation Model (MTGCM) was exercised for Ls = 90 (aphelion) solar minimum, and Ls = 270 perihelion) solar maximum conditions. Simulated MTGCM outputs (i.e. helium density distributions) were compared to those previously observed for Earth and Venus. Winter polar night bulges of helium are predicted on Mars, similar to those observed on the nightside of Venus and in the winter polar regions of Earth. A poster on this research was presented at the European Geophysical Society Meeting (EGS) in 2003. This research paves the way for what might be expected in the polar night regions of Mars during upcoming aerobraking and mapping Campaigns. Lastly, Mars thermosphere (approx. 100-130 km) winter polar warming was observed at high Northern latitudes during the perihelion season, but not at high Southern latitudes during the opposite aphelion season. Presumably, the Mars thermospheric circulation is responsible for the dynamically controlled heating needed to warm polar night temperatures above radiative equilibrium values. Again, MTGCM simulations were conducted for Ls = 90 and Ls = 270 conditions; polar temperatures were examined and found to be much warmer at Northern high latitudes (perihelion) than at Southern high latitudes (aphelion), similar to Mars aerobraking datasets. The Mars thermospheric circulation is found to be stronger during perihelion solstice conditions than during aphelion conditions, owing to both stronger seasonal solar and dust heating during Mars perihelion. An invited talk was given at the Spring AGU 2004 on this research. A forthcoming GRL paper was drafted on this same topic, but not submitted before the termination of this 1-year grant.

Bougher, Stephen↗

Internal reversing flow in a tailpipe offtake configuration for SSTOVL aircraft

A generic one-third scale model of a tailpipe offtake system for a supersonic short takeoff vertical landing (SSTOVL) aircraft was tested at LeRC Powered Lift Facility. The model consisted of a tailpipe with twin elbows, offtake ducts, and flow control nozzles, plus a small ventral nozzle and a blind flange to simulate a blocked cruise nozzle. The offtake flow turned through a total angle of 177 degrees relative to the tailpipe inlet axis. The flow split was 45 percent to each offtake and 10 percent to the ventral nozzle. The main test objective was to collect data for comparison to the performance of the same configuration predicted by a computational fluid dynamics (CFD) analysis. Only the experimental results are given - the analytical results are published in a separate paper. Performance tests were made with unheated air at tailpipe-to-ambient pressure ratios up to 5. The total pressure loss through the offtakes was as high as 15.5 percent. All test results are shown as graphs, contour plots, and wall pressure distributions. The complex flow patterns in the tailpipe and elbows at the offtake openings are described with traversing flow angle probe and paint streak flow visualization data.

Mcardle, Jack G.↗

Internal reversing flow in a tailpipe offtake configuration for SSTOVL aircraft

A generic one-third scale model of a tailpipe offtake system for a supersonic short takeoff vertical landing (SSTOVL) aircraft was tested at LeRC Powered Lift Facility. The model consisted of a tailpipe with twin elbows, offtake ducts, and flow control nozzles, plus a small ventral nozzle and a blind flange to simulate a blocked cruise nozzle. The offtake flow turned through a total angle of 177 degrees relative to the tailpipe inlet axis. The flow split was 45 percent to each offtake and 10 percent to the ventral nozzle. The main test objective was to collect data for comparison to the performance of the same configuration predicted by a computational fluid dynamics (CFD) analysis. Only the experimental results are given - the analytical results are published in a separate paper. Performance tests were made with unheated air at tailpipe-to-ambient pressure ratios up to 5. The total pressure loss through the offtakes was as high as 15.5 percent. All test results are shown as graphs, contour plots, and wall pressure distributions. The complex flow patterns in the tailpipe and elbows at the offtake openings are described with traversing flow angle probe and paint streak flow visualization data.

Mcardle, Jack G.↗

Integrated Vehicle Ground Vibration Testing in Support of NASA Launch Vehicle Loads and Controls Analysis

The National Aeronautics and Space Administration (NASA) Ares Projects Office (APO) is continuing to make progress toward the final design of the Ares I crew launch vehicle and Ares V cargo launch vehicle. Ares I and V will form the space launch capabilities necessary to fulfill NASA's exploration strategy of sending human beings to the Moon, Mars, and beyond. As with all new space vehicles there will be a number of tests to ensure the design can be Human Rated. One of these is the Integrated Vehicle Ground Vibration Test (IVGVT) that will be measuring responses of the Ares I as a system. All structural systems possess a basic set of physical characteristics unique to that system. These unique characteristics include items such as mass distribution, frequency and damping. When specified, they allow engineers to understand and predict how a structural system like the Ares I launch vehicle behaves under given loading conditions. These physical properties of launch vehicles may be predicted by analysis or measured through certain types of tests. Generally, these properties are predicted by analysis during the design phase of a launch vehicle and then verified through testing before the vehicle is Human Rated. The IVGVT is intended to measure by test the fundamental dynamic characteristics of Ares I during various phases of operational/flight. This testing includes excitations of the vehicle in lateral, longitudinal, and torsional directions at vehicle configurations representing different trajectory points. During the series of tests, properties such as natural frequencies, mode shapes, and transfer functions are measured directly. These data will then be used to calibrate loads and Guidance, Navigation, and Controls (GN&C) analysis models for verifying analyses of Ares I. NASA launch vehicles from Saturn to Shuttle have undergone Ground Vibration Tests (GVTs) leading to successful launch vehicles. A GVT was not performed on the unmanned Delta III. This vehicle was lost during launch. Subsequent analyses indicated that had a GVT been conducted on the vehicle, problems with vehicle modes and control may have been discovered and corrected, avoiding loss of the vehicle/mission. This paper will address GVT planning, set-up, conduction and analyses, for the Saturn and Shuttle programs, and also focus on the current and on-going planning for the Ares I and V IVGVT.

Tuma, Margaret L.↗

Case Study of Using Flocad to Model a Ground Test Station Ln2 Heat Exchanger

This paper presents a case study of using FLOCAD to model the heat transfer and boiling flow behavior of a Liquid Nitrogen (LN2) ground station equipment heat exchanger. The unique aspect of this work is that is one of very few studies documenting the use of FLOCAD in the spacecraft thermal control community. The results of the simulations described herein were used to guide and plan the execution of thermal vacuum testing spacecraft hardware using a ground test heat exchanger such as the one modeled herein. The paper will review the theory of two-phase boiling flow and correlations used by SINDA/FLUINT as well as the terminology and nomenclature used within Thermal Desktop / FLOCAD regarding two-phase boiling flow heat transfer modeling and simulation. Results for using the ground station heat exchanger as a heat sink in a thermal vacuum test set up are shown in order to demonstrate the set-up and execution of a typical FLOCAD model. Results indicate that properly modeling of the two-phase behavior is critical in order to ascertain the time constant associated with the heat exchanger, as well as understanding the temperature distribution across the test equipment and prediction of required LN2 flowrates to be used during testing. The paper is meant to augment the FLOCAD user’s manual and serve as a tutorial for thermal engineers and analysts wishing to learn and apply FLOCAD to industrial problems.

Anderson, Kevin R.↗

Boreal Winter Predictions with the GEOS-2 GCM: The Role of Boundary Forcing and Initial Conditions

Ensembles of atmospheric General Circulation Model (GCM) seasonal forecasts and long-term simulations (1980-94) are analyzed to assess the controlling influences of boundary forcing and memory of the initial conditions. Both the forecasts and simulations are carried out with version 2 of the Goddard Earth Observing System (GEOS-2) GCM forced with observed sea surface temperatures (SSTs). While much of the focus is on the seasonal time scale (January- March) and the Pacific North American (PNA) region, we also present results for other regions, shorter time scales, and other known modes of variability in the northern hemisphere extratropics. Forecasts of indices of some of the key large-scale modes of variability show that there is considerable variability in skill between different regions of the Northern Hemisphere. The eastern North Atlantic region has the poorest long lead forecast skill showing no skill beyond about 10 days. Skillful seasonal forecasts are primarily confined to the wave-like ENSO response emanating from the tropical Pacific. In the Northern Hemisphere, this is associated with the well-known Pacific/North American (PNA) pattern. Memory of the initial conditions is the major factor leading to skillful extratropical forecasts of lead time less than one month, while SST forcing is the only factor at the seasonal time scale. SST forcing contributes to skillful forecasts at sub- seasonal time scales only over the PNA region. The GEOS-2 GCM produces average (1980-94) signal to noise ratios which are less than one everywhere in the extratropics, except for the subtropical Pacific where they approach 1.5. When confined to the ENSO years, the maximum signal to noise ratios occur in the PNA region where they exceed three. An assessment of the sampling distribution of the forecasts suggests the model's ENSO response is very likely too weak. These results show some sensitivity to the uncertainties in the estimates of the SST forcing fields. In the North Pacific region, the sensitivity to SST forcing manifests itself primarily as changes in the variability of the PNA response, underscoring the need for an ensemble approach to the seasonal prediction problem.

Chang, Yehui↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

Development of A Maximum Anti-Slosh Baffle Pressure Load Model

The sloshing of propellants can affect the stability of a spacecraft and the integrity of the tank structure. Undesirable sloshing can be controlled by the addition of anti-slosh baffles, and the spacing and configuration of baffles are driven by damping requirements. The structural design of the baffle is determined after consideration of many factors, such as the strength and rigidity needed to support the baffle for its lifetime. Therefore, knowledge of distributed pressure loading is important for detailed structural design. In addition, the resultant force and moment produced by the distributed pressure are of direct importance to the design of a vehicle’s control system. Previous experimental investigations have been conducted to determine the liquid pressure loads and slosh damping associated with a rigid ring baffle. The results suggested that when the nondimensional velocity parameter is larger than 1.0, the theories agree with the test. However, when the velocity parameter is less than 1.0, all theories are nonconservative and under-predict the pressure loads. The present study has derived a maximum pressure load on the slosh baffle based on the energy conservation principle. It is verified from CFD that pressure in the slosh flow field can be decomposed into static and transient components. The CFD results confirm that there is a phase shift in pressure across the baffle, which depends on the fluid damping. Higher damping leads to a higher phase shift. The CFD investigation further verifies the proposed theory: the maximum pressure load occurs when the phase shift is 90 degrees. A comparison of the present computational results to the previous comprehensive experimental data validates the maximum pressure theory. When the baffle is submerged, the maximum pressure theory envelopes all the experimental data points.

CFD↗

Equatorial Kelvin waves: A UARS MLS view

Data from the Microwave Limb Sounder (MLS) instrument on the Upper Atmosphere Research Satellite (UARS) are used to compare two periods of Kelvin wave activity during different stages of the equatorial quasi-biennial oscillation. The analysis is carried out using an asynoptic mapping technique. A wide bandpass filter is used to isolate the frequency bands where Kelvin waves have been identified in previous studies. Time-height and time-latitude plots of the bandpassed data are used to identify Kelvin wave activity in the temperature and ozone fields. Frequency spectra of temperature and ozone amplitudes are constructed to further analyze the latitudinal and meridional distribution of Kelvin wave activity in zonal wavenumbers 1 and 2. The characteristics identified in these plots agree well with theoretical predictions and previous observations of middle atmosphere Kelvin waves. The time-height and time-latitude plots support the existence of Kelvin waves in discrete frequency bands; the slow, fast, and ultrafast Kelvin modes are all identified in the data. The characteristics of these modes do not vary much despite different mean flow conditions in the two periods examined. For the Kelvin wave-induced perturbations in ozone, the change from a transport-dominated regime below 10 hPa to a photochemically controlled regime above 10 hPa is clearly apparent in the height dependence of the phase difference between temperature and ozone. The ratios of the ozone perturbation amplitude to the temperature perturbation amplitude for the various observed Kelvin wave modes are in agreement with model estimates and LIMS (Limb Infrared Monitor of the Stratosphere) observations in the lower half of the region sampled but appear to be too large in the upper stratosphere and lower mesosphere.

Canziani, Pablo O.↗

Magnetization of the Lunar Crust

Magnetic fields measured by the satellite Lunar Prospector show large scale features resulting from remanently magnetized crust. Vector data synthesized at satellite altitude from a spherical harmonic model of the lunar crustal field, and the radial component of the magnetometer data, have been used to produce spatially continuous global magnetization models for the lunar crust. The magnetization is expressed in terms of localized basis functions, with a magnetization solution selected having the smallest root-mean square magnetization for a given fit to the data, controlled by a damping parameter. Suites of magnetization models for layers with thicknesses between 10 and 50 km are able to reproduce much of the input data, with global misfits of less than 0.5 nT (within the uncertainties of the data), and some surface field estimates. The magnetization distributions show robust magnitudes for a range of model thicknesses and damping parameters, however the magnetization direction is unconstrained. These global models suggest that magnetized sources of the lunar crust can be represented by a 30 km thick magnetized layer. Average magnetization values in magnetized regions are 30-40 mA/m, similar to the measured magnetizations of the Apollo samples and significantly weaker than crustal magnetizations for Mars and the Earth. These are the first global magnetization models for the Moon, providing lower bounds on the magnitude of lunar crustal magnetization in the absence of multiple sample returns, and can be used to predict the crustal contribution to the lunar magnetic field at a particular location.

Lunar↗

Biogeochemical modeling at mass extinction boundaries

The causes of major mass extinctions is a subject of considerable interest to those concerned with the history and evolution of life on earth. The primary objectives of the proposed plan of research are: (1) to develop quantitative time-dependent biogeochemical cycle models, coupled with an ocean atmosphere in order to improve the understanding of global scale physical, chemical, and biological processes that control the distribution of elements important for life at times of mass extinctions; and (2) to develop a comprehensive data base of the best available geochemical, isotopic, and other relevant geologic data from sections across mass extinction boundaries. These data will be used to constrain and test the biogeochemical model. These modeling experiments should prove useful in: (1) determining the possible cause(s) of the environmental changes seen at bio-event boundaries; (2) identifying and quantifying little-known feedbacks among the oceans, atmosphere, and biosphere; and (3) providing additional insights into the possible responses of the earth system to perturbations of various timescales. One of the best known mass extinction events marks the Cretaceous/Tertiary (K/T) boundary (66 Myr ago). Data from the K/T boundary are used here to constrain a newly developed time-dependent biogeochemical cycle model that is designed to study transient behavior of the earth system. Model results predict significant fluctuations in ocean alkalinity, atmospheric CO2, and global temperatures caused by extinction of calcareous plankton and reduction in the sedimentation rates of pelagic carbonates and organic carbon. Oxygen-isotome and other paleoclimatic data from K/T time provide some evidence that such climatic fluctuations may have occurred, but stabilizing feedbacks may have acted to reduce the ocean alkalinity and carbon dioxide fluctuations.

Rampino, M. R.↗