Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear regression analysis”

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 307 records · Page 17

The OCEAN ICE mooring compilation: a standardised, pan-Antarctic database of ocean hydrography and current time series

Continuous moored time series of temperature, salinity, pressure and current speed and direction are of great importance for understanding the continental shelf and under-ice-shelf dynamics and thermodynamics that govern water mass transformations and ice melting in and around Antarctic marginal seas. In these regions, icebergs and sea ice make ship-based mooring deployment and recovery challenging. Nevertheless, over decades, expeditions around the fringe of Antarctica sporadically deployed and recovered hundreds of moored instruments, including those facilitated through ice shelves boreholes. These datasets tend to be archived in a wide range of data centres, with, to our knowledge, no clear format standardisation. As a result, systematic analysis of historical mooring time series in the marginal seas is often challenging. Here we present the first version of a standardised pan-Antarctic moored hydrography and current time series compilation, with broad international contributions from data centres, research institutes and individual data owners. The mooring records in this compilation span over five decades, from the 1970s to the 2020s, providing an opportunity for a systematic study of the pan-Antarctic water mass transport and shelf connectivity. As a demonstration of the utility of this compilation, we present spectral analysis of the compiled current velocity time series, which unsurprisingly shows the dominating presence of tidal variability within most records. This component of the variability is fitted using multi-linear regression to tidal frequencies, and the tidal fit is removed from the original time series to leave de-tided variability. Given the limited record durations to months to years, de-tided variability is dominated by synoptic (3–10 d period), intraseasonal (10–80 d) and seasonal (∼6 months–1 year) signals. The spatial distribution of the kinetic energy integrated within frequency bands is presented and discussed within respective regional contexts, and future avenues of research are proposed. This data compilation is assembled under the endorsement of Ocean-Cryosphere Exchanges in ANtarctica: Impacts on Climate and the Earth System (OCEAN ICE) project (https://ocean-ice.eu/, last access: 23 October 2025) funded by the European Commission and UK Research and Innovation. It is available and regularly updated in NetCDF format with the SEANOE database at https://doi.org/10.17882/99922 (Zhou et al., 2024a).

54 ENVIRONMENTAL SCIENCES↗

Welch Method and Bootstrapping Applied to Subcritical Gamma Noise

We measured the prompt neutron decay constant 𝛼 of the CROCUS zero-power reactor at the Swiss Federal Institute of Technology Lausanne using cross-power spectral density (CPSD) analysis of gamma-gamma correlations from two trans-stilbene organic scintillators positioned near the reactor core. We measured critical and subcritical states, with water levels ranging from 960 mm (critical) to 800 mm (𝜌=−1.4 $ subcritical). Our analysis used the Welch method, dividing signal segments for fast Fourier transform (FFT) frequency analysis and applying bootstrapping uncertainty quantification that uses Welch-defined segments. Results demonstrated a clear increase in the measured 𝛼 as reactor reactivity decreased, distinguishing critical from subcritical conditions. At the 960-mm critical level, 𝛼 was estimated at 155.9 ± 0.7 s −1 , and for the 800-mm subcritical level, 𝛼 increased significantly to 367.3 ± 6.9 s –1 . A linear regression of subcritical states yielded a critical estimate of 154.0 ± 3.1 s –1 , aligning with the static 𝛼 estimate at critical. The bootstrapping method produced normally distributed 𝛼 estimates, confirming data consistency. The gamma CPSD 𝛼 estimates clearly distinguish reactor states and improve monitoring of zero-power reactors. The future deployment of modular and microreactors as potential candidates for noise analysis is demonstrated in CROCUS, particularly zero-power mock-ups of new designs. The improvement of noise analysis in the subcritical domain from this work will support experimental data for reactor deployment and procedure.

CROCUS↗

Analysis of Forest Foliage Using a Multivariate Mixture Model

Data with wet chemical measurements and near infrared spectra of ground leaf samples were analyzed to test a multivariate regression technique for estimating component spectra which is based on a linear mixture model for absorbance. The resulting unmixed spectra for carbohydrates, lignin, and protein resemble the spectra of extracted plant starches, cellulose, lignin, and protein. The unmixed protein spectrum has prominent absorption spectra at wavelengths which have been associated with nitrogen bonds.

Hlavka, C. A.↗

Air separation via a two-step solar thermochemical cycle based on (Ba,La) x Sr 1-x FeO 3-δ : Thermodynamic analysis

A two-step solar thermochemical cycle was considered for air separation to produce N2 based on (Ba,La) x Sr 1-x FeO 3-δ perovskite reduction/oxidation (redox) reactions for A-site fractions of 0 ≤ x ≤ 0.2. The cycle steps encompassed (1) thermal reduction and O 2 release via concentrated solar input and (2) re-oxidation with air to uptake O 2 and produce high-purity N 2 . Thermogravimetry at temperatures between 400 and 1100 °C in atmospheres of 0.005 to 90% O 2 /Ar at 1 bar was performed to measure equilibrium nonstoichiometries. The compound energy formalism was applied to model redox thermodynamics for both Ba 2+ and La 3+ substitution. Non-linear regression was used to determine the empirical parameters based on the thermogravimetric measurements. The model was used to define partial molar reaction enthalpies and entropies and predicted equilibrium oxygen nonstoichiometry as functions of oxide stoichiometry, site fraction, temperature, and O 2 partial pressure. The thermodynamic analysis showed the materials are appealing for air separation at temperatures below 800 °C.

42 ENGINEERING↗

Reduced-Order Models of Static Power Grids based on Spectral Clustering

For large-scale interconnected power systems that cover large geographical areas, certain electrical studies are required so that appropriate decisions ensure system reliability and low cost. For such studies, it is often neither practical nor necessary to model in detail the entire power system, which is increasingly complex due to a more diverse range of grid assets to choose from in both short and long-term planning. The goal of this paper is to present a methodology to reduce the order of large-scale power networks based on spectral graph theory given that current methods for static network reduction are not scalable. A brief analysis of some spectral clustering properties to determine which graph Laplacian matrix should be used and why is included. The analysis shows that the utilization of the normalized graph Laplacian is more advantageous for clustering purposes. Techniques are proposed to approximate cost functions for the aggregated generators. This is done via linear regression. The reduced-order model obtained with the proposed methodology has an accuracy above 94% and solves the scalability issue commonly present in other reduction methods. If the utilization of the reduced-order model is either constrained to load levels above mid-peak demand, or cost functions of aggregated units are approximated via a piecewise quadratic approach, then the error distribution is in the order of 10^-3. .

Baquedano-Aguilar, Mario D.↗

Application of LANDSAT to the surveillance of lake eutrophication in the Great Lakes basin

The author has identified the following significant results. A step-by-step procedure for establishing and monitoring the trophic status of inland lakes with the use of LANDSAT data, surface sampling, laboratory analysis, and aerial observations were demonstrated. The biomass was related to chlorophyll-a concentrations, water clarity, and trophic state. A procedure was developed for using surface sampling, LANDSAT data, and linear regression equations to produce a color-coded image of large lakes showing the distribution and concentrations of water quality parameters, causing eutrophication as well as parameters which indicate its effects. Cover categories readily derived from LANDSAT were those for which loading rates were available and were known to have major effects on the quality and quantity of runoff and lake eutrophication. Urban, barren land, cropland, grassland, forest, wetlands, and water were included.

Rogers, R. H.↗

Empirical Estimation of Shortest Route Length along U.S. Interstate Highways Based on Great Circle Distance

In this study, 98 regression models were specified for easily estimating shortest distances based on great circle distances along the U.S. interstate highways nationwide and for each of the continental 48 states. This allows transportation professionals to quickly generate distance, or even distance matrix, without expending significant efforts on complicated shortest path calculations. For simple usage by all professionals, all models are present in the simple linear regression form. Only one explanatory variable, the great circle distance, is considered to calculate the route distance. For each geographic scope (i.e., the national or one of the states), two different models were considered, with and without the intercept. Based on the adjusted R-squared, it was observed that models without intercepts generally have better fitness. Additionally, all these models generally have good fitness with the linear regression relationship between the great circle distance and route distance. At the state level, significant variations in the slope coefficients between the state-level models were also observed. Furthermore, a preliminary analysis of the effect of highway density on this variation was conducted.

33 ADVANCED PROPULSION SYSTEMS↗

COVID-19 prevention at institutions of higher education, United States, 2020–2021: implementation of nonpharmaceutical interventions

Background, In early 2020, following the start of the coronavirus disease 2019 (COVID-19) pandemic, institutions of higher education (IHEs) across the United States rapidly pivoted to online learning to reduce the risk of on-campus virus transmission. We explored IHEs’ use of this and other nonpharmaceutical interventions (NPIs) during the subsequent pandemic-affected academic year 2020–2021. Methods, From December 2020 to June 2021, we collected publicly available data from official webpages of 847 IHEs, including all public (n = 547) and a stratified random sample of private four-year institutions (n = 300). Abstracted data included NPIs deployed during the academic year such as changes to the calendar, learning environment, housing, common areas, and dining; COVID-19 testing; and facemask protocols. We performed weighted analysis to assess congruence with the October 29, 2020, US Centers for Disease Control and Prevention (CDC) guidance for IHEs. For IHEs offering ≥50% of courses in person, we used weighted multivariable linear regression to explore the association between IHE characteristics and the summated number of implemented NPIs. Results, Overall, 20% of IHEs implemented all CDC-recommended NPIs. The most frequently utilized NPI was learning environment changes (91%), practiced as one or more of the following modalities: distance or hybrid learning opportunities (98%), 6-ft spacing (60%), and reduced class sizes (51%). Additionally, 88% of IHEs specified facemask protocols, 78% physically changed common areas, and 67% offered COVID-19 testing. Among the 33% of IHEs offering ≥50% of courses in person, having < 1000 students was associated with having implemented fewer NPIs than IHEs with ≥ 1000 students. Conclusions, Only 1 in 5 IHEs implemented all CDC recommendations, while a majority implemented a subset, most commonly changes to the classroom, facemask protocols, and COVID-19 testing. IHE enrollment size and location were associated with degree of NPI implementation. Additional research is needed to assess adherence to NPI implementation in IHE settings.

59 BASIC BIOLOGICAL SCIENCES↗

Extreme Weather Events and Their Relationship to Low Frequency Teleconnection Patterns

A new method for identifying the structure and other characteristics of extreme weather events is introduced and applied to both model simulations and observations. The approach is based on a linear regression model that links daily extreme precipitation amounts for a particular point on the globe to precipitation and related quantities at all other points. We present here some initial results of our analysis of extreme precipitation events over the United States, including how they are influenced by ENSO and various large-scale teleconnection patterns such as the PNA. The results are based on simulations made with the NASA/NCAR AGCM (Lin and Rood 1996). The quality of the simulated climate for the NASA/NCAR AGCM forced with observed SSTs is described in Chang et al. (2001). The runs analyzed here consist of three 20-year runs forced with idealized cold, neutral and warm ENSO SST anomalies (superimposed on the mean seasonal cycle of SST). The idealized warm or cold SST anomalies are fixed throughout each 20- year simulation and consist of the first EOF (+/- 3 standard deviations) of monthly SST data. Comparisons are made with the results obtained from a similar analysis that uses daily NOAA precipitation observations (Higgins et al. 1996) over the United States and NCEP/NCAR reanalysis data for the period 1949-1998.

Chang, Yehui↗

Trivalent f-Element Squarates, Squarate-Oxalates, and Cationic Materials, and the Determination of the Nine-Coordinate Ionic Radius of Cf(III)

The synthesis, structure, and solid-state UV–vis–NIR spectroscopy of four new f-element squarates, M 2 (C 4 O 4 ) 3 (H 2 O) 4 (M = Eu, Am, Cf) and Sm(C 4 O 4 )(C 4 O 3 OH)(H 2 O) 2 ·0.5H 2 O, four new cationic lanthanide squarate chlorides, [M 4 (C 4 O 4 ) 5 (H 2 O) 12 ]Cl 2 ·5H 2 O (M = Eu, Dy, Ho Er), and two new actinide squarate oxalates, M 2 (C 4 O 4 ) 2 (C 2 O 4 )(H 2 O) 4 (M = Am, Cf), are presented. All of the metal centers are trivalent. Single-crystal X-ray diffraction analysis reveals that M 2 (C 4 O 4 ) 3 (H 2 O) 4 and Sm(C 4 O 4 )(C 4 O 3 OH)(H 2 O) 2 ·0.5H 2 O have a two-dimensional sheet structure constructed from MO 7 (H 2 O) 2 monocapped square-antiprismatic (coordination number (CN) = 9) metal centers and SmO 6 (H 2 O) 2 square-antiprismatic (CN = 8) metal centers, respectively, whereas M 2 (C 4 O 4 ) 2 (C 2 O 4 )(H 2 O) 4 have a three-dimensional (3D) structure constructed from MO 7 (H 2 O) 2 monocapped square-antiprismatic (CN = 9) metal centers. Additionally, the cationic framework materials [M 4 (C 4 O 4 ) 5 (H 2 O) 12 ]Cl 2 ·5H 2 O have a 3D structure constructed from two crystallographically unique MO 5 (H 2 O) 3 square-antiprismatic (CN = 8) metal centers. In these structures, the squarate ligands bind to the metal centers with varying coordination modes and denticities. The results of this study provide another example of the nonparallel chemistry between the lanthanides and transplutonium elements. From the crystallographic data for the isotypic series M 2 (C 4 O 4 ) 3 (H 2 O) 4 (M = La–Nd, Sm, Eu) and the linear regression fit to a plot of the unit cell volume as a function of the cube of the ionic radius, the nine-coordinate ionic radius of Cf 3+ was determined to be 1.127 ± 0.003 Å. Lastly, computational analysis of the americium and californium complexes M 2 (C 4 O 4 ) 3 (H 2 O) 4 and M 2 (C 4 O 4 ) 2 (C 2 O 4 )(H 2 O) 4 reveals three important attributes: (i) the 5f orbitals are nonbonding in all cases, with the bonding differences occurring with the empty 6d orbitals; (ii) the Cf complexes exhibit more covalent character than their Am counterparts; and (iii) there is more covalent character in the squarate-oxalate complexes than in the squarate complexes.

Crystallography↗

Evaluation of the CEAS model for barley yields in North Dakota and Minnesota

The CEAS yield model is based upon multiple regression analysis at the CRD and state levels. For the historical time series, yield is regressed on a set of variables derived from monthly mean temperature and monthly precipitation. Technological trend is represented by piecewise linear and/or quadriatic functions of year. Indicators of yield reliability obtained from a ten-year bootstrap test (1970-79) demonstrated that biases are small and performance as indicated by the root mean square errors are acceptable for intended application, however, model response for individual years particularly unusual years, is not very reliable and shows some large errors. The model is objective, adequate, timely, simple and not costly. It considers scientific knowledge on a broad scale but not in detail, and does not provide a good current measure of modeled yield reliability.

Barnett, T. L.↗

Surface radiation budget imagery

Linear regressions which summarize the relationships between the turbulent fluxes of heat and latent heat of vaporization leaving the surface and surface albedos are given. This observation is based on a semi-empirical analysis of the data collected from sets of Landsat/TM (Thematic Mapper) imagery as the top of the atmosphere images have been converted into maps of surface radiation budget. A preliminary analysis shows that in spite of the land surface's temporal and spatial variability, the transfer of latent and sensible heat into the atmosphere would appear to be governed by the land surface albedo.

Kim, Hongsuk H.↗

Validation of MODIS Aerosol Optical Depth Retrievals over a Tropical Urban Site, Pune, India

In the present paper, MODIS (Terra and Aqua; level 2, collection 5) derived aerosoloptical depths (AODs) are compared with the ground-based measurements obtained from AERONET (level 2.0) and Microtops - II sun-photometer over a tropical urban station, Pune (18 deg 32'N; 73 deg 49'E, 559 m amsl). This is the first ever systematic validation of the MODIS aerosol products over Pune. Analysis of the data indicates that the Terra and Aqua MODIS AOD retrievals at 550 nm have good correlations with the AERONET and Microtops - II sun-photometer AOD measurements. During winter the linear regression correlation coefficients for MODIS products against AERONET measurements are 0.79 for Terra and 0.62 for Aqua; however for premonsoon, the corresponding coefficients are 0.78 and 0.74. Similarly, the linear regression correlation coefficients for Microtops measurements against MODIS products are 0.72 and 0.93 for Terra and Aqua data respectively during winter and are 0.78 and 0.75 during pre-monsoon. On yearly basis in 2008-2009, correlation coefficients for MODIS products against AERONET measurements are 0.80 and 0.78 for Terra and Aqua respectively while the corresponding coefficients are 0.70 and 0.73 during 2009-2010. The regressed intercepts with MODIS vs. AERONET are 0.09 for Terra and 0.05 for Aqua during winter whereas their values are 0.04 and 0.07 during pre-monsoon. However, MODIS AODs are found to underestimate during winter and overestimate during pre-monsoon with respect to AERONET and Microtops measurements having slopes 0.63 (Terra) and 0.74 (Aqua) during winter and 0.97 (Terra) and 0.94 (Aqua) during pre-monsoon. Wavelength dependency of Single Scattering Albedo (SSA) shows presence of absorbing and scattering aerosol particles. For winter, SSA decreases with wavelength with the values 0.86 +/- 0.03 at 440 nm and 0.82 +/- 0.04 at 1020nm. In pre-monsoon, it increases with wavelength (SSA is 0.87 +/- 0.02 at 440nm; and 0.88 +/-0.04 at 1020 nm).

More, Sanjay↗

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗

Expanding the application of soil moisture monitoring systems through regression-based transformation

Relative to other geophysical variables, soil moisture (SM) estimates derived from land surface models (LSMs) and land data assimilation systems (LDAS) are difficult to transfer between platforms and applications. This difficulty stems from the highly model-dependent nature of LSM SM estimates and differences in the vertical support of discretized SM values. As a result, operational SM estimates generated by one LSM (or LDAS) cannot generally be directly applied to a hydrologic monitoring or forecast system designed around a second LSM. This lack of transferability is particularly problematic for LDAS applications, where the time, expertise, and computational resources required to generate an operational LDAS analysis cannot be practically duplicated for every LSM-specific application. Here, we develop a set of simple regression tools for translating SM estimates between LSMs and multiple LDAS analyses. Results demonstrate that simple multivariate linear regression - utilizing independent variables based on multi-layer and temporally lagged SM estimates - can significantly improve upon baseline transformation approaches using direct percentile matching. The proposed regression approaches are effective for both the LSM-to-LSM and LDAS-to-LDAS transformation of multi-layer SM percentiles. Application of this approach will expand the utility of existing, high-quality (but LSM-specific) operational sources of SM information like the NASA Soil Moisture Active Passive Level-4 Soil Moisture product.

Soil Moisture↗

Long-wavelength Magnetic and Gravity Anomaly Correlations of Africa and Europe

Preliminary MAGSAT scalar magnetic anomaly data were compiled for comparison with long-wavelength-pass filtered free-air gravity anomalies and regional heat-flow and tectonic data. To facilitate the correlation analysis at satellite elevations over a spherical-Earth, equivalent point source inversion was used to differentially reduce the magnetic satellite anomalies to the radial pole at 350 km elevation, and to upward continue the first radial derivative of the free-air gravity anomalies. Correlation patterns between these regional geopotential anomaly fields are quantitatively established by moving window linear regression based on Poisson's theorem. Prominent correlations include direct correspondences for the Baltic Shield, where both anomalies are negative, and the central Mediterranean and Zaire Basin where both anomalies are positive. Inverse relationships are generally common over the Precambrian Shield in northwest Africa, the Basins and Shields in southern Africa, and the Alpine Orogenic Belt. Inverse correlations also presist over the North Sea Rifts, the Benue Rift, and more generally over the East African Rifts. The results of this quantitative correlation analysis support the general inverse relationships of gravity and magnetic anomalies observed for North American continental terrain which may be broadly related to magnetic crustal thickness variations.

Vonfrese, R. R. B.↗

Long-wavelength magnetic and gravity anomaly correlations on Africa and Europe

Preliminary MAGSAT scalar magnetic anomaly data were compiled for comparison with long-wavelength-pass filtered free-air gravity anomalies and regional heat-flow and tectonic data. To facilitate the correlation analysis at satellite elevations over a spherical-Earth, equivalent point source inversion was used to differentially reduce the magnetic satellite anomalies to the radial pole at 350 km elevation, and to upward continue the first radial derivative of the free-air gravity anomalies. Correlation patterns between these regional geopotential anomaly fields are quantitatively established by moving window linear regression based on Poisson's theorem. Prominent correlations include direct correspondences for the Baltic shield, where both anomalies are negative, and the central Mediterranean and Zaire Basin where both anomalies are positive. Inverse relationships are generally common over the Precambrian Shield in northwest Africa, the Basins and Shields in southern Africa, and the Alpine Orogenic Belt. Inverse correlations also presist over the North Sea Rifts, the Benue Rift, and more generally over the East African Rifts. The results of this quantitative correlation analysis support the general inverse relationships of gravity and magnetic anomalies observed for North American continental terrain which may be broadly related to magnetic crustal thickness variations.

Vonfrese, R. R. B.↗