Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial 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 289 records · Page 16

NASA Tech Briefs, June 2012

Topics covered include: iGlobe Interactive Visualization and Analysis of Spatial Data; Broad-Bandwidth FPGA-Based Digital Polyphase Spectrometer; Small Aircraft Data Distribution System; Earth Science Datacasting v2.0; Algorithm for Compressing Time-Series Data; Onboard Science and Applications Algorithm for Hyperspectral Data Reduction; Sampling Technique for Robust Odorant Detection Based on MIT RealNose Data; Security Data Warehouse Application; Integrated Laser Characterization, Data Acquisition, and Command and Control Test System; Radiation-Hard SpaceWire/Gigabit Ethernet-Compatible Transponder; Hardware Implementation of Lossless Adaptive Compression of Data From a Hyperspectral Imager; High-Voltage, Low-Power BNC Feedthrough Terminator; SpaceCube Mini; Dichroic Filter for Separating W-Band and Ka-Band; Active Mirror Predictive and Requirement Verification Software (AMP-ReVS); Navigation/Prop Software Suite; Personal Computer Transport Analysis Program; Pressure Ratio to Thermal Environments; Probabilistic Fatigue Damage Program (FATIG); ASCENT Program; JPL Genesis and Rapid Intensification Processes (GRIP) Portal; Data::Downloader; Fault Tolerance Middleware for a Multi-Core System; DspaceOgreTerrain 3D Terrain Visualization Tool; Trick Simulation Environment 07; Geometric Reasoning for Automated Planning; Water Detection Based on Color Variation; Single-Layer, All-Metal Patch Antenna Element with Wide Bandwidth; Scanning Laser Infrared Molecular Spectrometer (SLIMS); Next-Generation Microshutter Arrays for Large-Format Imaging and Spectroscopy; Detection of Carbon Monoxide Using Polymer-Composite Films with a Porphyrin-Functionalized Polypyrrole; Enhanced-Adhesion Multiwalled Carbon Nanotubes on Titanium Substrates for Stray Light Control; Three-Dimensional Porous Particles Composed of Curved, Two-Dimensional, Nano-Sized Layers for Li-Ion Batteries 23 Ultra-Lightweight; and Ultra-Lightweight Nanocomposite Foams and Sandwich Structures for Space Structure Applications.

Source record↗

Spatial profiling of the interplay between cell type- and vision-dependent transcriptomic programs in the visual cortex

How early sensory experience during “critical periods” of postnatal life affects the organization of the mammalian neocortex at the resolution of neuronal cell types is poorly understood. We previously reported that the functional and molecular profiles of layer 2/3 (L2/3) cell types in the primary visual cortex (V1) are vision-dependent [S. Chenget al.,Cell185, 311–327.e24 (2022)]. Here, we characterize the spatial organization of L2/3 cell types with and without visual experience. Spatial transcriptomic profiling based on 500 genes recapitulates the zonation of L2/3 cell types along the pial–ventricular axis in V1. By applying multitasking theory, we suggest that the spatial zonation of L2/3 cell types is linked to the continuous nature of their gene expression profiles, which can be represented as a 2D manifold bounded by three archetypal cell types. By comparing normally reared and dark reared L2/3 cells, we show that visual deprivation-induced transcriptomic changes comprise two independent gene programs. The first, induced specifically in the visual cortex, includes immediate-early genes and genes associated with metabolic processes. It manifests as a change in cell state that is orthogonal to cell-type-specific gene expression programs. By contrast, the second program impacts L2/3 cell-type identity, regulating a subset of cell-type-specific genes and shifting the distribution of cells within the L2/3 cell-type manifold. Through an integrated analysis of spatial transcriptomics with single-nucleus RNA-seq data, we describe how vision patterns cortical L2/3 cell types during the critical period.

Science & Technology - Other Topics↗

Application of a data-mining method based on Bayesian networks to lesion-deficit analysis

Although lesion-deficit analysis (LDA) has provided extensive information about structure-function associations in the human brain, LDA has suffered from the difficulties inherent to the analysis of spatial data, i.e., there are many more variables than subjects, and data may be difficult to model using standard distributions, such as the normal distribution. We herein describe a Bayesian method for LDA; this method is based on data-mining techniques that employ Bayesian networks to represent structure-function associations. These methods are computationally tractable, and can represent complex, nonlinear structure-function associations. When applied to the evaluation of data obtained from a study of the psychiatric sequelae of traumatic brain injury in children, this method generates a Bayesian network that demonstrates complex, nonlinear associations among lesions in the left caudate, right globus pallidus, right side of the corpus callosum, right caudate, and left thalamus, and subsequent development of attention-deficit hyperactivity disorder, confirming and extending our previous statistical analysis of these data. Furthermore, analysis of simulated data indicates that methods based on Bayesian networks may be more sensitive and specific for detecting associations among categorical variables than methods based on chi-square and Fisher exact statistics.

NASA Discipline Neuroscience↗

The complex structural and chemical nature of monolithic U-10Mo fuel and Zr barrier layer

Nanoscale microstructural characterization by advanced transmission electron microscopy techniques on a U-10Mo/Zr barrier layer monolithic fuel plate was performed to evaluate the microstructural evaluation after high burn-up. Gas bubble superlattice evolution, grain restructuring, and evolution of the Zr interaction layer is investigated through detailed electron microscopy characterization. The use of automated crystallographic orientation mapping to irradiated U-10Mo fuel highlights that the restructured ultra-fine grains are separated by high angle grain boundaries at a burn up of 4.42 × 10 21 fissions/cm 3 . Additionally, advanced chemical analysis and multi-variable statistical analysis shows spatial clustering of solid fission product precipitates. Finally, characterization of a newly observed porous nanocrystalline Zr region in the barrier layer is studied. Finally, this work provides insights into the grain subdivision and restructuring process while using advanced microscopy techniques to analyze fission products in neutron-irradiated U-10Mo fuel.

36 MATERIALS SCIENCE↗

A high resolution spatial and temporal multiple Doppler analysis of a microburst and its application to aircraft flight simulation

Multiple Doppler radar data collected during the Joint Airport Weather Studies (JAWS) Project is used to synthesize the three-dimensional wind in the region of a microburst. The particular microburst used in this study is the strongest one to date for which three-dimensional winds have been recovered. As a diagnostic tool, a six-DOF numerical aircraft model having characteristics similar to 727 series aircraft is used to investigate jet transport aircraft response to observed microburst winds during simulated approaches and departures. Simple pilot control laws are used to adjust thrust, pitch, roll, and yaw so as to maintain given approach or departure parameters. Generally, when horizontal wind shear along the approach or departure path is 0.01/s or greater, the model is unable to maintain the desired approach path and suffers a significant reduction in climb performance during a go-around or departure. Although the mean wind shear along a path gives a good qualitative measure of the wind shear threat to a jet transport, different paths with similar mean shears can yield markedly different results, as do the same paths through the microburst at different times. These findings are a direct consequence of the fine temporal and spatial scale of microburst winds. During any given modeled aircraft traverse through the region of highest horizontal shear, time variations in the microburst wind field are shown to have an insignificant effect on the modeled flight path. This is because the traverse period is short (30 s) compared to the lifetime of a microburst (300-600 s).

Elmore, K. L.↗

Spatially Resolved Outflows in a Seyfert Galaxy at z=2.39

We present the first spatially resolved analysis of rest-frame optical and UV (UltraViolet) imaging and spectroscopy for a lensed galaxy at z equals 2.39 hosting a Seyfert active galactic nucleus (AGN). Proximity to a natural guide star has enabled observations with high signal-to-noise ratio using Very Large Telescope SINFONI (Spectrograph for INtegral Field Observations in the Near Infrared) plus adaptive optics (AO) of rest-frame optical diagnostic emission lines, which exhibit an underlying broad component with full width at half maximum approximately 700 kilometers per second in both the Balmer and forbidden lines. Measured line ratios place the outflow robustly in the region of the ionization diagnostic diagrams associated with AGNs. This unique opportunity - combining gravitational lensing, AO guiding, redshift, and AGN activity - allows for a magnified view of two main tracers of the physical conditions and structure of the interstellar medium in a star-forming galaxy hosting a weak AGN at Cosmic Noon. By analyzing the spatial extent and morphology of the Lyman alpha spectral line and dust-corrected Hydrogen spectral line emission, disentangling the effects of star formation and AGN ionization on each tracer, and comparing the AGN-induced mass outflow rate to the host star formation rate, we find that the AGN does not significantly impact the star formation within its host galaxy.

Fischer, Travis C.↗

Discovery and analysis of low-surface-brightness galaxies in the environment of NGC 1052

The environment of NGC 1052 has recently attracted much attention because of the presence of low-surface-brightness galaxies (LSBGs) with apparently “exotic” properties, making it a region of high interest for the detection of new objects. Here, we used public deep photometric data from the Dark Energy Camera Legacy Survey to carry out a comprehensive search for LSBGs over a wide region of 6 × 6 degrees, equivalent to 2 × 2 Mpc at the distance of NGC 1052. We detected 42 LSBGs with r eff > 5 arcsec and μ g (0) > 24 mag arcsec -2 , of which 20 are previously undetected objects. Among all the newly detected objects, RCP 32 stands out with extreme properties: r eff = 23.0 arcsec and $\langle$μ g $\rangle$ eff = 28.6 mag arcsec -2 . This makes RCP 32 one of the lowest surface brightness galaxies ever detected through integrated photometry, located at just 10 arcmin from the extensively studied NGC 1052-DF2. We explored the presence of globular clusters (GCs) in the LSBGs. We marginally detected a GC system in RCP 32, and argue that this LSBG is of great interest for follow-up observations given its extremely low baryon density. After analyzing the distribution of galaxies with available spectroscopy, we identified a large-scale structure of approximately 1 Mpc that is well isolated in redshift space and centered on NGC 1052. The spatial correlation analysis between the LSBGs and this large-scale structure suggests their association. However, when exploring the distribution of effective radius, we find an overpopulation of large LSBGs (r eff > 15 arcsec) located close to the line of sight of NGC 1052. We argue that this is suggestive of a substructure with similar radial velocity in sight projection, but at a closer distance, to which some of these apparently larger LSBGs could be associated. However, possible effects derived from tidal interactions are worthy of further study. Our work expands the catalog of LSBGs with new interesting objects and provides a detailed environmental context for the study of LSBGs in this region.

79 ASTRONOMY AND ASTROPHYSICS↗

Particle charge distributions in the effluent of a flow-through atmospheric pressure low temperature plasma

Atmospheric pressure low-temperature plasmas are often utilized to perform particle synthesis, treatment, and removal. It is well-known that dust particles are highly negatively charged in these plasmas; however, little is known about dust particle charging behavior as particles leave the plasma volume and pass through the spatial afterglow region. In this work, monodisperse particles of various sizes and work functions were introduced into an atmospheric pressure radiofrequency capacitively coupled flow-through plasma. Dust particle electrical mobility distributions downstream of the flow-through plasma were measured utilizing a differential mobility analyzer in conjunction with a condensation particle counter at various gas flow velocities. Charge distributions were determined from the measured electrical mobility distributions. Experiments confirm that particles become less negatively charged, and even net-positively charged after leaving the plasma volume, with a distribution that follows a shifted Boltzmann charge distribution. Additionally, particle charge in the effluent of the flow-through plasma is negligibly dependent on work function but highly size and flow velocity dependent. Larger particles were shown to have a higher magnitude of charge under all studied conditions; however, particle polarity was switchable by varying gas flow velocity. The charging dynamics were simulated utilizing a constant number Monte Carlo model that accounts for electron temperature decay and the transition from ambipolar to free diffusion of electrons and ions in the spatial afterglow. Simulation results also suggest that, at the same flow velocity, larger particles obtain a greater magnitude of charge, negative or positive. Furthermore, the decrease in electron mobility and the difference between ion and electron convective loss rates create an ion-rich region in the plasma effluent that promotes ion–particle collisions and drives particle charge removal and even reversal of polarity. Larger particles more favorably collide with energetic species in these environments, which results in higher charge states.

atmospheric pressure low temperature plasma↗

MeV–GeV Gamma-Ray Emission from SNR G327.1–1.1 Discovered by the Fermi-LAT

Abstract We report the discovery of MeV–GeV γ -ray emission by the Fermi-LAT positionally coincident with the TeV pulsar wind nebula (PWN) HESS J1554–550 within the host supernova remnant (SNR) G327.1–1.1. The γ -ray emission is point-like and faint but significant (>4 σ ) in the 300 MeV–2 TeV energy range. We report here the Fermi-LAT analysis of the observed γ -ray emission followed by a detailed multiwavelength investigation to understand the nature of the emission. The central pulsar powering the PWN within G327.1–1.1 has not been detected in any wave band; however, it is likely embedded within the X-ray nebula, which is displaced from the center of the radio nebula. The γ -ray emission is faint and therefore a pulsation search to determine if the pulsar may be contributing is not feasible. Prior detailed multiwavelength reports revealed an SNR system that is old, τ ∼ 18,000 yr, where the interaction of the reverse shock with the PWN is underway or has recently occurred. We find that the γ -ray emission agrees remarkably well with a detailed broadband model constructed in a prior report based on independent hydrodynamical and semianalytic simulations of an evolved PWN. We further investigate the physical implications of the model for the PWN evolutionary stage incorporating the new Fermi-LAT data and attempt to model the distinct particle components based on a spatial separation analysis of the displaced PWN counterparts.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluation of Satellite-Based Rainfall Estimates in the Lower Mekong River Basin (Southeast Asia)

Satellite-based precipitation is an essential tool for regional water resource applications that requires frequent observations of meteorological forcing, particularly in areas that have sparse rain gauge networks. To fully realize the utility of remotely sensed precipitation products in watershed modeling and decision-making, a thorough evaluation of the accuracy of satellite-based rainfall and regional gauge network estimates is needed. In this study, Tropical Rainfall Measuring Mission (TRMM) Multi-Satellite Precipitation Analysis (TMPA) 3B42 v.7 and Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS) daily rainfall estimates were compared with daily rain gauge observations from 2000 to 2014 in the Lower Mekong River Basin (LMRB) in Southeast Asia. Monthly, seasonal, and annual comparisons were performed, which included the calculations of correlation coefficient, coefficient of determination, bias, root mean square error (RMSE), and mean absolute error (MAE). Our validation test showed TMPA to correctly detect precipitation or no-precipitation 64.9% of all days and CHIRPS 66.8% of all days, compared to daily in-situ rainfall measurements. The accuracy of the satellite-based products varied greatly between the wet and dry seasons. Both TMPA and CHIRPS showed higher correlation with in-situ data during the wet season (June–September) as compared to the dry season (November–January). Additionally, both performed better on a monthly than an annual time-scale when compared to in-situ data. The satellite-based products showed wet biases during months that received higher cumulative precipitation. Based on a spatial correlation analysis, the average r-value of CHIRPS was much higher than TMPA across the basin. CHIRPS correlated better than TMPA at lower elevations and for monthly rainfall accumulation less than 500 mm. While both satellite-based products performed well, as compared to rain gauge measurements, the present research shows that CHIRPS might be better at representing precipitation over the LMRB than TMPA.

Dandridge, Chelsea↗

Assessing heterogeneity of patient and health system delay among TB in a population with internal migrants in China

Backgrounds The diagnostic delay of tuberculosis (TB) contributes to further transmission and impedes the implementation of the End TB Strategy. Therefore, we aimed to describe the characteristics of patient delay, health system delay, and total delay among TB patients in Shanghai, identify areas at high risk for delay, and explore the potential factors of long delay at individual and spatial levels. Method The study included TB patients among migrants and residents in Shanghai between January 2010 and December 2018. Patient and health system delays exceeding 14 days and total delays exceeding 28 days were defined as long delays. Time trends of long delays were evaluated by Joinpoint regression. Multivariable logistic regression analysis was employed to analyze influencing factors of long delays. Spatial analysis of delays was conducted using ArcGIS, and the hierarchical Bayesian spatial model was utilized to explore associated spatial factors. Results Overall, 61,050 TB patients were notified during the study period. Median patient, health system, and total delays were 12 days (IQR: 3–26), 9 days (IQR: 4–18), and 27 days (IQR: 15–43), respectively. Migrants, females, older adults, symptomatic visits to TB-designated facilities, and pathogen-positive were associated with longer patient delays, while pathogen-negative, active case findings and symptomatic visits to non-TB-designated facilities were associated with long health system delays (LHD). Spatial analysis revealed Chongming Island was a hotspot for patient delay, while western areas of Shanghai, with a high proportion of internal migrants and industrial parks, were at high risk for LHD. The application of rapid molecular diagnostic methods was associated with reduced health system delays. Conclusion Despite a relatively shorter diagnostic delay of TB than in the other regions in China, there was vital social-demographic and spatial heterogeneity in the occurrence of long delays in Shanghai. While the active case finding and rapid molecular diagnosis reduced the delay, novel targeted interventions are still required to address the challenges of TB diagnosis among both migrants and residents in this urban setting.

Sun, Ruoyao↗

Effect of molecular permanent dipole moment on guest aggregation and exciton quenching in phosphorescent organic light emitting diodes

This study explores the effect of molecular permanent dipole moment (PDM) on aggregation of guest molecules in phosphorescent host–guest organic light-emitting diodes (OLEDs). Through a combination of photoluminescence measurements, high-angle annular dark-field scanning transmission electron microscopy analysis, and an Ising model based physical vapor-deposition simulation, we show that higher PDM of tris[2-phenylpyridinato-C2,N]iridium(III) guest can actually lead to a reduced aggregation relative to tris[bis[2-(2-pyridinyl-N)phenyl-C] (acetylacetonato)iridium(III) when doped into a non-polar host 1,3,5-tris(carbazol-9-yl)benzene. This study further explores the effect of host polarity by using a polar host 3′,5′-di(carbazol-9-yl)-[1,1′-biphenyl]-3,5-dicarbonitrile, and it is shown that the polar host leads to reduced guest aggregation. In conclusion, this study provides a comprehensive understanding of the impact of molecular PDM on OLED material efficiency and stability, providing insights for optimizing phosphorescent OLED materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Upscale integration of normalized difference vegetation index - The problem of spatial heterogeneity

An analysis is conducted of the correspondence between the normalized difference vegetation index (NDVI) calculated from average reflectances, or M(NDVI), and the I(NDVI) that is integrated from individual NDVIs, by simulating AVHRR data from high spatial resolution SPOT 1 radiometer and Landsat TM data. The West African and French sites analyzed at 300-1000 m scale show a strong correlation between the two types of index; the relationship is almost perfectly linear, with a slope that is somewhat dependent on vegetation cover. Effecting the scale change using M(NDVI) instead of I(NDVI) does not introduce significant errors.

Aman, Angora↗

Nth-plant supply: corn stover supplies and costs in a fleet of biorefineries

Feedstock cost and cost variability is expected to increase with the number of biorefineries. To quantify this effect, this spatial-economic analysis simulates feedstock cost and cost variability of an industry based on corn stover as a function of the number of biorefineries. Results are reported for nine scenarios (a base case and sensitivity analysis of four variables – harvest efficiency, sustainability constraints, opportunity cost, and corn grain yield) under deterministic and stochastic simulations, assuming biorefineries using 658 000 Mg (725 000 tons) year- 1 of corn stover in 2019. The resulting supply curves are highly elastic (i.e. little change in cost) for the first 50 of the 121 biorefineries, with price increases in subsequent biorefineries depending on scenario. In the base-case deterministic scenario, weighted-average stover costs are $66 Mg -1 ($60 ton- 1 ), $69 Mg -1 ($62 ton -1 ), and $156 Mg -1 ($142 ton -1 ), at the first, 60th, and 121st biorefineries, respectively. The stochastic simulations, subject to observed 30-year corn yield variability, follow a similar pattern, with price distributions that vary by scenario. The base-case stochastic simulations illustrate minimal cost variability for the first 60 biorefineries, but rapid increases in cost variability in the second half of potential biorefineries, with similar patterns observed in the other scenarios. Of the four variables explored, price was most sensitive to harvest efficiency, followed by sustainability constraints, corn yield, and opportunity cost. Results suggest that, under conventional logistics, about half of the US corn stover resource is reliably available with minimum cost increase and variability. Interactive visualization is available at https://doi.org/10.11578/1828779.

09 BIOMASS FUELS↗

Uncertainty quantification of bank vegetation impacts on the flood flow field in the American River, California, using large‐eddy simulations

Bank vegetation plays a key role in both hydrodynamics and morphodynamics of natural rivers; however, these effects are often unaccounted for in the computational flow dynamics of natural waterways. Recent studies using the large‐eddy simulation (LES), however, have attempted to gain insights into the impacts of bank vegetation on the mean flow field of the natural rivers using a vegetation model, which applies a sink term to the momentum equations of motion. This approach accounts for the effects of the vegetation and provides a practical approach to account for the complex patches of bank vegetation in large‐scale rivers. To implement the vegetation model, a drag coefficient reflecting the overall resistance of vegetal structures to the flow is needed, but due to the lack of calibrated data and range of size, density and type of vegetation, this parameter can be a significant source of uncertainty in the model results. Here, in this study, we use uncertainty quantification (UQ) to investigate the hydrodynamics and bed shear results when a bank vegetation is incorporated in an LES model. To this end, we used the polynomial chaos expansion and Monte Carlo sampling techniques to determine the uncertainties associated with the drag coefficient in the vegetation model and from uncertainties in the bed roughness and inflow discharge. The UQ analysis provided spatially varying confidence levels for the spanwise and vertical distribution of velocity magnitude and for the bed shear stress distributions. In addition, Sobol indices were computed to indicate the relative influence that each parameter had on the overall uncertainty. In general, it was found that uncertainty in flow discharge was the dominant source of uncertainty; however, the drag coefficient in the vegetation model and the bed roughness parameter also made significant contribution to the uncertainty near the banks and bed, respectively.

54 ENVIRONMENTAL SCIENCES↗

Single-Cell Metabolomics with Rapid Determination of Chemical Formulas from Isotopic Fine Structures

Metabolomic measurements can provide functional readouts of cellular states and phenotypes. Here, we present a protocol for single-cell metabolomics that permits direct untargeted detection of a broad number of metabolites under ambient conditions, without the need for sample processing, and with high confidence in the discovery and identification of the molecular formulas for detected metabolites. This protocol describes combining fiber-based laser ablation electrospray ionization (f-LAESI) with a 21 Tesla Fourier transform ion cyclotron resonance mass spectrometer (21T-FTICR-MS) to obtain high confidence molecular formula information about detected metabolites. The f-LAESI source utilizes mid-infrared laser ablation through a sharp optical fiber tip, affording direct ambient analysis of cells without the need for sample processing. Using the 21T-FTICR-MS as a mass analyzer enabled measurement of the isotopic fine structure (IFS) for numerous metabolites simultaneously from single cells, and the IFSs were in turn computationally processed to rapidly determine the corresponding elemental compositions. This metabolomics technique complements other single cell omics measurement methods, helping to resolve complex molecular interactions that take place within cells unattainable from single cell transcriptomic and proteomics methods.

mass spectrometry, FTICR, FTMS, ultrahigh mass res↗

Exploring the Spatial Relationship Between Demographic Indicators and the Built Environment of a City

In addition to global and regional drivers of urbanization, neighborhood development in urban areas across the United States has been shown to be influenced by various local socio-economic factors. These factors, despite varying across socio-economic groups, have large implications regarding a population’s vulnerability to extreme climate events, including heat waves resulting in adverse health impacts. Additionally, the demographics of an urban area can shape its infrastructural characteristics, causing different populations groups to face varying levels of risks and benefits. As a result, the urban morphology and socio-economic characteristics of a city are deeply intertwined; however, their interactions on a finer scale are not yet fully understood. This research aims to better understand the relationships between various socio-economic factors and the built environment of a city, considering variability in building types, and temperature patterns. This research focuses on the city of Las Vegas, NV, and uses spatial data analysis to understand the correlation between of socio-economic characteristics, building morphology, building characteristics, and temperature data to understand the correlation between these various factors. Results of these research shows there is a distinct pattern of clustering of socio-economic characteristics with the city and there is a distinct correlation between age and cost, socio-economic characteristics, and locations of high heat distribution within the city.

Singh, Ridhima↗

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗