Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data fitting”

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 1,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

A Probabilistic Mass Estimation Algorithm for a Novel 7- Channel Capacitive Sample Verification Sensor

A document describes an algorithm created to estimate the mass placed on a sample verification sensor (SVS) designed for lunar or planetary robotic sample return missions. A novel SVS measures the capacitance between a rigid bottom plate and an elastic top membrane in seven locations. As additional sample material (soil and/or small rocks) is placed on the top membrane, the deformation of the membrane increases the capacitance. The mass estimation algorithm addresses both the calibration of each SVS channel, and also addresses how to combine the capacitances read from each of the seven channels into a single mass estimate. The probabilistic approach combines the channels according to the variance observed during the training phase, and provides not only the mass estimate, but also a value for the certainty of the estimate. SVS capacitance data is collected for known masses under a wide variety of possible loading scenarios, though in all cases, the distribution of sample within the canister is expected to be approximately uniform. A capacitance-vs-mass curve is fitted to this data, and is subsequently used to determine the mass estimate for the single channel s capacitance reading during the measurement phase. This results in seven different mass estimates, one for each SVS channel. Moreover, the variance of the calibration data is used to place a Gaussian probability distribution function (pdf) around this mass estimate. To blend these seven estimates, the seven pdfs are combined into a single Gaussian distribution function, providing the final mean and variance of the estimate. This blending technique essentially takes the final estimate as an average of the estimates of the seven channels, weighted by the inverse of the channel s variance.

Wolf, Michael↗
DOE OSTI2021

Re-Analysis of DARHT Axis 2 S4 Magnet Sweep Measurements

The Dual-Axis Radiographic Hydrodynamic Test (DARHT) facility provides flash radiography capabilities using two electron Linear Induction Accelerators (LIA’s). Understanding fundamental properties of the electron beams is essential to optimizing spot size and dose for radiography. This work describes distribution measurements of a single-kicked electron beam (16.5MeV, 1.7kA, 60ns) at a location downstream at imaging station C on DARHT Axis 2. The beam at station C is focused by the S4 solenoid located upstream. The beam distribution is measured by imaging optical transition radiation (OTR) from electrons striking a 51-µm thick titanium foil. This data and analysis were originally published by Ekdahl. The results here are a re-analysis of the asymmetric halo emittance contribution as compared to the emittance resulting from the FWHM of the distribution. This is accomplished using xtr2 fits to the data which include non-linear magnetic field effects in S4. Also examined are the calculated S4 spherical aberration contributions to the emittance using xtr.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

Verification Testing of OLI Systems Mixed Solvent Electrolyte Model for the Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 ) System to High Ionic Strength at 25°C

This technical report summarizes model verification results and summary statistics for 41 evaporite mineral solubility cases evaluated by Savannah River National Laboratory using OLI Systems’ aqueous electrolyte thermodynamic modeling software. The 41 verification cases containing a total of 60 solubility curves comprise mineral solubility data from low to high ionic strength at 25°C for the eight-component system Na-K-Mg-Ca-H-Cl-SO 4 -OH-HCO 3 -CO 3 -CO 2 -H 2 O as reported by Harvie et al. (1984). Thermodynamic calculations were executed using OLI Systems’ Stream Analyzer computation module within the OLI Studio software platform (Ver. 11.0, Rev. 11.0.1.9). The Mixed Solvent Electrolyte (MSE) thermodynamic framework was chosen for this investigation because of its superiority in modeling high ionic-strength inorganic salt solutions and actinide redox chemistry and solubility, both of which are relevant to the geological repository conditions at the Waste Isolation Pilot Plant in Carlsbad, New Mexico. Mineral solubility data in various inorganic salt solutions were digitized and extracted from figures generated by Harvie et al. (1984). For each of the 60 solubility curves, a case-specific chemistry model and input file were generated in OLI Studio using OLI Stream Analyzer and the MSE (H 3 O + ion) public databank provided by OLI Systems. Model simulation results were exported to Microsoft Excel to calculate summary statistics and to generate graphs comparing the OLI model predictions to the solubility data. Summary statistics include residuals (model – data) and concordance (accuracy × precision, where precision is indicated by the Pearson correlation coefficient and accuracy accounts for bias and scale differential). Private databanks were not developed, and activity coefficient model regressions were not performed to improve OLI model fits to the data. Of the 41 model verification plots, 83% have a mean of the percent residuals less than or equal to 25%. Similarly, 75% display a concordance greater than or equal to 0.75. Only seven of the 41 verification plots fail to show good agreement between the model and data. Of these seven, three are relevant to the WIPP repository because they involve the Mg-OH-Cl-SO 4 -CO 3 aqueous system. The remaining four address salt solubilities at the pH extremes (strong acid and strong base). It should be noted that in two of the three Mg-OH-Cl-SO 4 -CO 3 system cases, the regressed Harvie et al. (1984) solubility curve also deviated from the data. Lack of agreement between the OLI model-predicted solubility curves and the data is attributable to one or more of the following: specific solid species are not included in the OLI MSE databank; there is significant variation among the different solubility datasets chosen by Harvie et al. (1984); the OLI MSE model’s thermodynamic parameters were determined using different solubility datasets; and the activity coefficient parameters for certain relevant ion-ion and ion-molecule pairs have not been optimized via data regression. Two recommendations for future work are to (1) evaluate solubility data for the Mg-OH-Cl-SO 4 -CO 3 system at high ionic strength and, if necessary, develop a private OLI MSE database that includes missing species and, where necessary, regressed standard state properties and interaction parameters; (2) perform similar verification testing of the OLI model for actinide solubility data.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

RedMaPPer: Evolution and Mass Dependence of the Conditional Luminosity Functions of Red Galaxies in Galaxy Clusters

In this work we characterize the luminosity distribution, halo-mass dependence, and redshift evolution of red galaxies in galaxy clusters using the SDSS Data Release 8 redMaPPer cluster sample. We propose a simple prescription for the relationship between the luminosity of both red central and red satellite galaxies and the mass of their host halos, and show that this model is well fit by the data. Using a larger galaxy cluster sample than previously employed in the literature, we find that the luminosities of red central galaxies scale as $\langle \mathrm{log}L\rangle \propto {A}_{L}\mathrm{log}({M}_{200b})$, with A L = 0.39 ± 0.04, and that the scatter of the red central galaxy luminosity at fixed M 200b (${\sigma }_{\mathrm{log}L| M}$) is ${0.23}_{-0.04}^{+0.05}$ dex, with the error bar including systematics due to mis-centering of the cluster finder, photometry, and photometric redshift estimation. Our data prefers a positive correlation between the luminosity of red central galaxies and the observed richness of clusters at a fixed halo mass, with an effective correlation coefficient ${d}_{\mathrm{eff}}={0.36}_{-0.16}^{+0.17}$. The characteristic luminosity of red satellites becomes dimmer from z = 0.3 to z = 0.1 by ~20% after accounting for passive evolution. We estimate the fraction of galaxy clusters where the brightest red galaxy is not the central to be P BNC ~ 20%. We discuss implications of these findings in the context of galaxy evolution and the galaxy–halo connection.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

Interaction of the equatorial midnight pressure bulge and thermospheric zonal winds

The zonal component of the thermospheric neutral wind in the equatorial region is found to flow eastward all night, having a minimum near midnight. The data, obtained by Spencer et al. on the Dynamics Explorer-2 (DE-2) satellite, and reported more recently by Wharton et al., indicate that the minimum feature is significant. The minimum is illustrated by a 4th order Fourier series fit to the data points reported by Wharton et al. A table gives the amplitudes and phases of the four harmonic components used with the estimated error for each amplitude. It is asked whether this behavior is due to the local passage of the equatorial midnight pressure bulge associated with the midnight temperature anomaly. If so, it is to be expected that other data characterizing the midnight temperature anomaly should show consistency with this observation, especially in view of the large changes observed in the zonal velocity in the midnight sector. Consistency in momentum conservation may be checked using available data from previous independent experiments. Such data, taken under similar conditions of solar activity, was substituted into the momentum equation and found to be consistent using a priori estimates of the effect of viscosity.

Herrero, F. A.↗
NASA NTRS1988

A comparison of mapped and measured total ionospheric electron content using global positioning system and beacon satellite observations

Total ionospheric electron contents (TEC) were measured by global positioning system (GPS) dual-frequency receivers developed by the Jet Propulsion Laboratory. The measurements included P-code (precise ranging code) and carrier phase data for six GPS satellites during multiple five-hour observing sessions. A set of these GPS TEC measurements were mapped from the GPS lines of sight to the line of sight of a Faraday beacon satellite by statistically fitting the TEC data to a simple model of the ionosphere. The mapped GPS TEC values were compared with the Faraday rotation measurements. Because GPS transmitter offsets are different for each satellite and because some GPS receiver offsets were uncalibrated, the sums of the satellite and receiver offsets were estimated simultaneously with the TEC in a least squares procedure. The accuracy of this estimation procedure is evaluated indicating that the error of the GPS-determined line of sight TEC can be at or below 1 x 10 to the 16th el/sq cm. Consequently, the current level of accuracy is comparable to the Faraday rotation technique; however, GPS provides superior sky coverage.

Lanyi, Gabor E.↗
NASA NTRS1998

Stellar Occultation Probe of Triton's Atmosphere

The goals of this research were (i) to better characterize Triton's atmospheric structure by probing a region not well investigated by Voyager and (ii) to begin acquiring baseline data for an investigation of the time evolution of the atmosphere which will set limits on the thermal conductivity of the surface and the total mass of N2 in the atmosphere. Our approach was to use observations (with the Kuiper Airborne Observatory) of a stellar occultation by Triton that was predicted to occur on 1993 July 10. As described in the attached reprint, we achieved these objectives through observation of this occultation and a subsequent one with the KAO in 1995. We found new results about Triton's atmospheric structure from the analysis of the two occultations observed with the KAO and ground-based data. These stellar occultation observations made both in the visible and infrared, have good spatial coverage of Triton including the first Triton central-flash observations, and are the first data to probe the 20-100 km altitude level on Triton. The small-planet light curve model of Elliot and Young (AJ 103, 991-1015) was generalized to include stellar flux refracted by the far limb, and then fitted to the data. Values of the pressure, derived from separate immersion and emersion chords, show no significant trends with latitude indicating that Triton's atmosphere is spherically symmetric at approximately 50 km altitude to within the error of the measurements. However, asymmetry observed in the central flash indicates the atmosphere is not homogeneous at the lowest levels probed (approximately 20 km altitude). From the average of the 1995 occultation data, the equivalent-isothermal temperature of the atmosphere is 47 +/- 1 K and the atmospheric pressure at 1400 km radius (approximately 50 km altitude) is 1.4 +/- 0.1 microbar. Both of these are not consistent with a model based on Voyager UVS and RSS observations in 1989 (Strobel et al, Icarus 120, 266-289). The atmospheric temperature from the occultation is 5 K colder than that predicted by the model and the observed pressure is a factor of 1.8 greater than the model.

Elliot, James L.↗
NASA NTRS1998

The Radio luminosity Function of Radio-Loud Quasars from the 7C Redshift Survey

We present a complete sample of 24 radio-loud quasars (RLQs) from the new 7C Redshift Survey. Every quasar with a low-frequency (151 MHz) radio flux-density S(sub 151) > 0.5 Jy in two regions of the sky covering 0.013 sr is included; 23 of these have sufficient extended flux to meet the selection criteria, 18 of these have steep radio spectra (hereafter denoted as SSQs). The key advantage of this sample over most samples of RLQs is the lack of an optical magnitude limit. By combining the 7C and 3CRR samples, we have investigated the properties of RLQs as a function of redshift z and radio luminosity L(sub 151). We derive the radio luminosity function (RLF) of RLQs and find that the data are well fitted by a single power-law with slope alpha(sub 1) = 1.9 +/- 0.1 (for H(sub 0) = 50 km/s.Mpc, OMEGA(sub M) = 1, OMEGA(sub DELTA) = 0). We find that there must be a break in the RLQ RLF at log(sub 10)(L(sub 151)/W Hz.sr) approximately < or = 27, in order for the models to be consistent with the 7C and 6C source counts. The z-dependence of the RLF follows a one-tailed gaussian which peaks at z = 1.7 +/- 0.2. We find no evidence for a decline in the co-moving space density of RLQs at higher redshifts. A positive correlation between the radio and optical luminosities of SSQs is observed, confirming a result of Serjeant. We are able to rule out this correlation being due to selection effects or biases in our combined sample. The radio-optical correlation and best-fit model RLF enable us to estimate the distribution of optical magnitudes of quasars in samples selected at low radio frequencies, We con- clude that for samples with S(sub 151) approximately < or = 1 Jy one must use optical data significantly deeper than the POSS-I limit (R approximately equal 20), in order to avoid severe incompleteness.

Willott, Chris J.↗
NASA NTRS1997

Model of visual contrast gain control and pattern masking

We have implemented a model of contrast gain and control in human vision that incorporates a number of key features, including a contrast sensitivity function, multiple oriented bandpass channels, accelerating nonlinearities, and a devisive inhibitory gain control pool. The parameters of this model have been optimized through a fit to the recent data that describe masking of a Gabor function by cosine and Gabor masks [J. M. Foley, "Human luminance pattern mechanisms: masking experiments require a new model," J. Opt. Soc. Am. A 11, 1710 (1994)]. The model achieves a good fit to the data. We also demonstrate how the concept of recruitment may accommodate a variant of this model in which excitatory and inhibitory paths have a common accelerating nonlinearity, but which include multiple channels tuned to different levels of contrast.

NASA Discipline Space Human Factors↗
NASA NTRS1989

Dynamic sea surface topography, gravity and improved orbit accuracies from the direct evaluation of SEASAT altimeter data

A method for the simultaneous solution of dynamic ocean topography, gravity and orbits using satellite altimeter data is described. A GEM-T1 based gravitational model called PGS-3337 that incorporates Seasat altimetry, surface gravimetry and satellite tracking data has been determined complete to degree and order 50. The altimeter data is utilized as a dynamic observation of the satellite's height above the sea surface with a degree 10 model of dynamic topography being recovered simultaneously with the orbit parameters, gravity and tidal terms in this model. PGS-3337 has a geoid uncertainty of 60 cm root-mean-square (RMS) globally, with the uncertainty over the altimeter tracked ocean being in the 25 cm range. Doppler determined orbits for Seasat, show large improvements, with the sub-30 cm radial accuracies being achieved. When altimeter data is used in orbit determination, radial orbital accuracies of 20 cm are achieved. The RMS of fit to the altimeter data directly gives 30 cm fits for Seasat when using PGS-3337 and its geoid and dynamic topography model. This performance level is two to three times better than that achieved with earlier Goddard earth models (GEM) using the dynamic topography from long-term oceanographic averages. The recovered dynamic topography reveals the global long wavelength circulation of the oceans with a resolution of 1500 km. The power in the dynamic topography recovery is now found to be closer to that of oceanographic studies than for previous satellite solutions. This is attributed primarily to the improved modeling of the geoid which has occurred. Study of the altimeter residuals reveals regions where tidal models are poor and sea state effects are major limitations.

Marsh, J. G.↗
DOE OSTI2024

Forward Modeling of 3-D Ion Properties in Jupiter’s Magnetosphere Using Juno/JADE-I Data

The Jovian Auroral Distributions Experiment Ion sensor (JADE-I) on NASA’s Juno mission provides in-situ measurements of ions from 0.1 to 46.2 keV/q inside Jupiter’s magnetosphere. JADE-I is used to study the plasma with two types of datasets from the same measurement: Time-of-flight (TOF) and SPECIES. The TOF dataset provides mass-per-charge measurements with a range of 1–64 amu/q but oversamples particles over 6π steradian viewing per spacecraft spin and has little directional information. On the other hand, the SPECIES dataset can provide a good measurement of the flow direction but does not provide mass-per-charge information due to the telemetry limit. In this study, we developed a 2-step forward modeling method that combines the advantages and avoids the disadvantages of TOF and SPECIES data to derive the 3-D properties of heavy ions. Assuming that the ion velocity distribution can be described with the kappa distribution, we first perform the forward model fit of the TOF data to calculate the relative abundance of heavy ion species. Then we fix the relative abundance and perform the second forward model fit on the SPECIES data. Here, using this method, we obtain the densities of different heavy ions, the shared temperature and kappa value, and the 3-D flow velocity vector. Some data examples of the equatorial plasma disk before Perijove 24 are included to demonstrate the method. Plasma properties can then be mapped to explore spatial and temporal variabilities in Jupiter’s magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Empirical estimation of weather-driven yield shocks using biophysical characteristics for U.S. rainfed and irrigated maize, soybeans, and winter wheat

Agricultural yields are highly susceptible to changes in weather system patterns, including both annual and sub-annual changes in temperature and precipitation. Understanding the impacts of future changes in meteorological variables on crop yields have been widely studied in both empirical and process-based models. This work presents a structural econometric approach that combines historical weather data with information about the biophysical growth cycles of maize, winter wheat, and soy to predict the year-to-year yield responses for rainfed crops. Modeled soil moisture and temperature are taken as key predictors, which allows testing the fitted rainfed model’s ability to predict irrigated yields by assuming irrigation produces the numerically optimal level of soil moisture. This approach is grounded in known biophysical processes, which improves confidence in the model predictions. It produces estimates of interannual variability in yields, which enable study of the effects of extreme events on agricultural productivity. Finally, it enables prediction of the potential impacts of changing weather patterns on irrigated crops in areas that are currently primarily rainfed. We present the results of the empirical model, fitted with rainfed data; out-of-sample validation on irrigated crops; and projections of yield shocks under multiple future climate scenarios. Under a bias-corrected GFDL RCP8.5 scenario, this approach predicts average yield change, relative to 2006-2020 average yields, across U.S. counties in the 2040-2060 period of -16.9% and -13.4% for rainfed and irrigated maize, of -19.1% and -18.4% for rainfed and irrigated soy, and of -4.0% and -2.0% for rainfed and irrigated winter wheat. And in the 2070-2090 period of -30.3% and -28.8 % for rainfed and irrigated maize, of -33.1% and -32.8% for rainfed and irrigated soy, and of -6.5% and -5.2%for rainfed and irrigated winter wheat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Improved Neutron Lifetime Measurement with UCN$\tau$

In this work, we report an improved measurement of the free neutron lifetime $\tau_n$ using the UCN$\tau$ apparatus at the Los Alamos Neutron Science Center. We count a total of approximately 38×10 6 surviving ultracold neutrons (UCNs) after storing in UCN$\tau$’s magnetogravitational trap over two data acquisition campaigns in 2017 and 2018. We extract $\tau$ n from three blinded, independent analyses by both pairing long and short storage time runs to find a set of replicate $\tau$ n measurements and by performing a global likelihood fit to all data while self-consistently incorporating the β-decay lifetime. Both techniques achieve consistent results and find a value $\tau$ n =877.75±0.28 stat +0.22/–0.16 syst s. With this sensitivity, neutron lifetime experiments now directly address the impact of recent refinements in our understanding of the standard model for neutron decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1987

Deconvolution of spectra for intimate mixtures

Visible to near infrared reflectance spectra of macroscopic mixtures have been shown to be linear combinations of the reflections of the pure mineral components in the mixture. However, for microscopic mixtures the mixing systematics are in general nonlinear. The systematics may be linearized by conversion of reflectance to single scattering albedo (w), where the equations which relate reflectance to w depend on the method of data collection. Several proposed mixing models may be used to estimate mineral abundances from the reflectance spectra of intimate mixtures. These models are summarized and a revised model is presented. A noniterative (linear) least squares approach was used for curve fitting and the data, measured as bi-directional reflectance with incidence and emergence angles of 30 and 0 deg were converted to w by a simplified version of Hapke's equation for bi-directional reflectance. This model was tested with two mixture series composed of 45 to 75 micron particles: an anorthite-enstatite series and an olivine-magnetite series. The data indicate that the simplified Hapke's equation may be used to convolve reflectance spectra into mineral abundances if appropriate endmembers are known or derived from other techniques. For surfaces that contain a significant component of very low albedo material, a somewhat modified version of this technique will need to be developed. Since the abundances are calculated using a noniterative approach, the application of this method is especially efficient for large spectral data sets, such as those produced by mapping spectrometers.

Mustard, John F.↗
NASA NTRS1992

Modelling atmospheric scatterers using spacecraft observations

Voyager images of Triton indicate considerable spatial variability in the concentration of at least two different scattering components in the atmosphere. Data from high phase angle limb scans were fit to Mie scattering models to derive mean particle sizes, number densities, and vertical extent for both types of scattering material at ten different locations between 15 deg S and 70 deg S. These fits reveal a thin haze at latitudes equatorward of 25-30 deg S. The imaging data can be fit reasonably well by both conservatively scattering and absorbing hazes with particle sizes near 0.18 micron and optical depths of order 0.001-0.01. Rayleigh scattering haze fits the imaging data somewhat less well, and can be totally ruled out by combining the imaging and UVS measurements. At high southern latitudes, Triton displays clouds below an altitude of approximately 8 km, as well as the haze at higher altitudes. The clouds have particle sizes which may range from 0.7-2.0 microns, or may be near 0.25 micron. The atmospheric optical depth poleward of 30 deg S must be generally greater than 0.1, but need not be more than 0.3. Horizontal inhomogeneities are quite noticeable, especially at longitudes east of (i.e., higher than) 180 deg.

Rages, Kathy A.↗
NASA NTRS2008

Visual Data Analysis for Satellites

The Visual Data Analysis Package is a collection of programs and scripts that facilitate visual analysis of data available from NASA and NOAA satellites, as well as dropsonde, buoy, and conventional in-situ observations. The package features utilities for data extraction, data quality control, statistical analysis, and data visualization. The Hierarchical Data Format (HDF) satellite data extraction routines from NASA's Jet Propulsion Laboratory were customized for specific spatial coverage and file input/output. Statistical analysis includes the calculation of the relative error, the absolute error, and the root mean square error. Other capabilities include curve fitting through the data points to fill in missing data points between satellite passes or where clouds obscure satellite data. For data visualization, the software provides customizable Generic Mapping Tool (GMT) scripts to generate difference maps, scatter plots, line plots, vector plots, histograms, timeseries, and color fill images.

Lau, Yee↗
NASA NTRS2007

Remote and Ground Truth Spectral Measurement Comparisons

FORMOSAT III are a set of six research satellites from Taiwan that were launched in April 2006. The satellites are in 800 km, 71 degree inclination orbits and separated by 24 degrees in ascending node. Ground truth spectral measurements were taken of outer surface materials on FORMOSAT III. From those measurements, a computer model was built to predict the spectral reflectance, which included phase angle and orientation of the spacecraft relative to the observer. However, materials exposed to the space environment have exhibited spectral changes including a darkening and a reddening of the spectra. This reddening was seen as an increase in slope of the reflectance as the wavelength increases. Therefore, the model of pristine materials was augmented to include the space weathering effects. Remote data were collected on two of the six FORMOSAT satellites using the 1.6 meter telescope at AMOS (Air Force Maui Optical and Supercomputing) site with the Spica spectrometer. Due to the separation in ascending node, observations were made on whichever one of the six satellites was visible on that specific night. Three nights of data were collected using the red (6000 9500 angstroms) filter and two nights of data were collected using the blue (3200 -6600 angstroms) filter. A comparison of the data showed a good match to the pristine model for the blue filter region. The absorption feature near 5500 angstroms due to the copper colored Kapton multi-layer insulation (MLI) was very apparent in the remote samples and a good fit to the data was seen in both satellites observed. The features in the red filter regime agreed with the pristine model up through 7000 angstroms where the reddening begins and the slope of the remote sample increases. A comparison of the two satellites showed similar features in the red and blue filter regions, i.e. the satellites were aging at the same rate. A comparison of the pristine model to the first month of remote measurements showed the amount by which the satellite had reddened. This information was used in the space weathered model to predict the second month s reflectance measurement. The results showed that the reddening did not increase from the first month of observations to the second. A third month of data will be necessary to determine if the reddening occurs instantly and then stabilizes with time or if the effect continues to increase with time.

Abercromby, Kira Jorgensen↗
NASA NTRS2012

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↗