Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial regression”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

The Viability of Trajectory Analysis for Diagnosing Dynamical and Chemical Influences on Ozone Concentrations in the UTLS

The viability of trajectory analysis for diagnosing the interplay between chemistry and dynamics is investigated by comparing ozone mixing ratios modelled using air-parcel pathways to values observed along flight tracks during ATTREX (Airborne Tropical TRopopause EXperiment). Trajectories are initiated at the locations of ozone observations and tracked backward in time to their sources at termini of backward trajectories. The modelled values of ozone utilize 3-dimensional analysis fields from WACCM (Whole Atmosphere Community Climate Model) (a chemical-climate model with dynamical fields nudged towards MERRA (Modern-Era Retrospective Analysis and Research Applications) reanalysis) and ERA-interim (product of ECMWF - the European Centre for Medium-Range Weather Forecasts) to determine source mixing ratios with chemical production and loss terms derived from the ozone chemistry used in WACCM. A statistical base of modelled ozone is constructed with 6 trajectory platforms (adiabatic, diabatic, and kinematic forced by ERA-interim and MERRA), two chemical models (WACCM chemistry and no chemistry), and 4 trajectory lengths (5, 10, 20, and 30 days). Linear regression is employed to separate systematic errors from random errors and to characterize the impact of source mixing ratios, path length, vertical motion, and chemistry on modelled ozone errors. Errors in the analysis ozone fields are large, if not dominant, contributors to model error. Random errors are particularly large for point-by-point comparisons, however averaging over 800 km (75 minutes) flight segments substantially reduces random error and exposes systematic errors. Of the two analysis ozone data sets, WACCM, which incorporates detailed chemistry, provides the smaller systematic errors while ERA-interim, which has crude chemistry but assimilates observational data, has the smaller random errors. Of the different trajectory platforms, adiabatic calculations produce the smaller random errors (irrespective of the use of chemistry) but both vertical motion and chemistry are required to optimally reduce systematic errors. These results suggest that meaningful analysis of dynamical and chemical interactions that control ozone mixing ratios are viable on spatial scales larger than a few reanalysis grid spaces, that errors in the analyzed ozone data sets are large but not prohibitively so, and that vertical velocities and heating rates from reanalysis data, while problematic, contain useful information [on the ozone concentrations in the UTLS (Upper Troposphere/Lower Stratosphere)].

ozone↗
NASA NTRS2002

Doppler-derived myocardial systolic strain rate is a strong index of left ventricular contractility

BACKGROUND: Myocardial fiber strain is directly related to left ventricular (LV) contractility. Strain rate can be estimated as the spatial derivative of velocities (dV/ds) obtained by tissue Doppler echocardiography (TDE). The purposes of the study were (1) to determine whether TDE-derived strain rate may be used as a noninvasive, quantitative index of contractility and (2) to compare the relative accuracy of systolic strain rate against TDE velocities alone. METHODS AND RESULTS: TDE color M-mode images of the interventricular septum were recorded from the apical 4-chamber view in 7 closed-chest anesthetized mongrel dogs during 5 different inotropic stages. Simultaneous LV volume and pressure were obtained with a combined conductance-high-fidelity pressure catheter. Peak elastance (Emax) was determined as the slope of end-systolic pressure-volume relationships during caval occlusion and was used as the gold standard of LV contractility. Peak systolic TDE myocardial velocities (Sm) and peak (epsilon'(p)) and mean (epsilon'(m)) strain rates obtained at the basal septum were compared against Emax by linear regression. Emax as well as TDE systolic indices increased during inotropic stimulation with dobutamine and decreased with the infusion of esmolol. A stronger association was found between Emax and epsilon'(p) (r=0.94, P<0.01, y=0.29x+0.46) and epsilon'(m) (r=0.88, P<0.01) than for Sm (r=0.75, P<0.01). CONCLUSIONS: TDE-derived epsilon'(p) and epsilon'(m) are strong noninvasive indices of LV contractility. These indices appear to be more reliable than S(m), perhaps by eliminating translational artifact.

Non-NASA Center↗
NASA NTRS2011

Modeling the Height of Young Forests Regenerating from Recent Disturbances in Mississippi using Landsat and ICESat data

Many forestry and earth science applications require spatially detailed forest height data sets. Among the various remote sensing technologies, lidar offers the most potential for obtaining reliable height measurement. However, existing and planned spaceborne lidar systems do not have the capability to produce spatially contiguous, fine resolution forest height maps over large areas. This paper describes a Landsat-lidar fusion approach for modeling the height of young forests by integrating historical Landsat observations with lidar data acquired by the Geoscience Laser Altimeter System (GLAS) instrument onboard the Ice, Cloud, and land Elevation (ICESat) satellite. In this approach, "young" forests refer to forests reestablished following recent disturbances mapped using Landsat time-series stacks (LTSS) and a vegetation change tracker (VCT) algorithm. The GLAS lidar data is used to retrieve forest height at sample locations represented by the footprints of the lidar data. These samples are used to establish relationships between lidar-based forest height measurements and LTSS-VCT disturbance products. The height of "young" forest is then mapped based on the derived relationships and the LTSS-VCT disturbance products. This approach was developed and tested over the state of Mississippi. Of the various models evaluated, a regression tree model predicting forest height from age since disturbance and three cumulative indices produced by the LTSS-VCT method yielded the lowest cross validation error. The R(exp 2) and root mean square difference (RMSD) between predicted and GLAS-based height measurements were 0.91 and 1.97 m, respectively. Predictions of this model had much higher errors than indicated by cross validation analysis when evaluated using field plot data collected through the Forest Inventory and Analysis Program of USDA Forest Service. Much of these errors were due to a lack of separation between stand clearing and non-stand clearing disturbances in current LTSS-VCT products and difficulty in deriving reliable forest height measurements using GLAS samples when terrain relief was present within their footprints. In addition, a systematic underestimation of about 5 m by the developed model was also observed, half of which could be explained by forest growth that occurred between field measurement year and model target year. The remaining difference suggests that tree height measurements derived using waveform lidar data could be significantly underestimated, especially for young pine forests. Options for improving the height modeling approach developed in this study were discussed.

Li, Ainong↗
DOE OSTI2020

Distribution Development for the RDX Regional Model at Los Alamos National Laboratory - 20385

Representing uncertainty in model inputs often means finding a balance between uncertainty and physical reality. Developing wide distributions may seem conservative in principle, but this approach may lead to unrealistic model results. Characterizing the current state of knowledge of stochastic inputs presents many challenges, especially if the data available are limited or have limited relevance to the site. Relationships among these inputs may also be important to represent but are typically complex or difficult to define. Often special adjustments must be made to account for reduced credibility in particular data. If parameters are strongly related to one another, a correlation structure may be developed for input to the model. Other techniques such as regression models may be used to incorporate relationships between the information available and the desired parameters. The process of developing distributions must consider details of the model in terms of what the distribution is meant to represent. This paper uses the example of a probabilistic fate and transport model for hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) in the regional aquifer at Los Alamos National Laboratory (LANL). For many parameters, a single draw is applied to all space and time over which the model is run, for a single iteration. This simplification is often made for many reasons, and can often be beneficial, but also adds additional complexity in the distribution development process. Defining the distributional goals as they relate to the modeling process is an important step, which should take place prior to evaluation of the data. Usually, the distributions developed are meant to characterize the average value of the parameter over the spatial and temporal domain of the model. Distribution development requires consideration of many sources of information on the parameter where available, ideally from multiple references. Examples of different sources include data from different references but also from different conditions, measurement methods, or experimental types. Depending on these conditions and the reliability or relevance of particular references, different sources of data may each contribute valuable information but have varying relevance to the site. In these cases, weighting data unequally is a useful way to incorporate this information. As an example, aqueous dispersivity data are available for a variety of rock types. Only a few values are available for the desired rock type, and this is not enough to develop a distribution. Therefore, dispersivity values from other rock types are included in distribution development but are down-weighted such that the best data have the most influence on the distribution developed. In another case, K{sub d} distributions in the model are meant to represent a known composition of multiple soil types. Data from these materials are weighted accordingly to develop a distribution for the weighted average K{sub d} across soil types. Other cases include varying reliability of different sources, and weighting data according to the confidence in these sources. In some cases, input parameters are correlated with one another. An example is advective porosity, which is positively correlated with total porosity and must be less than total porosity. Paired data with both parameters must exist to discern the relationship between the two parameters and if it is necessary to build a correlation structure into the model. In general, correlated parameters may be represented in the model by a multivariate distribution, or perhaps more desirably, capturing the correlation within the developed distributions which may be treated as independent from one another. In the example of porosity, this can be done by transforming advective porosity into a proportion of total porosity which may be drawn independently from the distribution of total porosity. This paper explores how complex data and correlations can be incorporated to meet the distributional goals of the model, using the RDX regional model as a detailed example. (authors)

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

EFIT-Prime: Probabilistic and physics-constrained reduced-order neural network model for equilibrium reconstruction in DIII-D

We introduce EFIT-Prime, a novel machine learning surrogate model for EFIT (Equilibrium FIT) that integrates probabilistic and physics-informed methodologies to overcome typical limitations associated with deterministic and ad hoc neural network architectures. EFIT-Prime utilizes a neural architecture search-based deep ensemble for robust uncertainty quantification, providing scalable and efficient neural architectures that comprehensively quantify both data and model uncertainties. Physically informed by the Grad–Shafranov equation, EFIT-Prime applies a constraint on the current density J tor and a smoothness constraint on the first derivative of the poloidal flux, ensuring physically plausible solutions. Furthermore, the spatial location of the diagnostics is explicitly incorporated in the inputs to account for their spatial correlation. Extensive evaluations demonstrate EFIT-Prime's accuracy and robustness across diverse scenarios, most notably showing good generalization on negative-triangularity discharges that were excluded from training. Timing studies indicate an ensemble inference time of 15 ms for predicting a new equilibrium, offering the possibility of plasma control in real-time, if the model is optimized for speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine Learning Correlation of Electron Micrographs and ToF-SIMS for the Analysis of Organic Biomarkers in Mudstone

The spatial distribution of organics in geological samples can be used to determine when and how these organics were incorporated into the host rock. Mass spectrometry (MS) imaging can rapidly collect a large amount of data, but ions produced are mixed without discrimination, resulting in complex mass spectra that can be difficult to interpret. Here, we apply unsupervised and supervised machine learning (ML) to help interpret spectra from time-of-flight-secondary ion mass spectrometry (ToF-SIMS) of an organic-carbon-rich mudstone of the Middle Jurassic of England (UK). It was previously shown that the presence of sterane molecular biomarkers in this sample can be detected via ToF-SIMS (Pasterski, M. J. et al., Astrobiology 2023, 23, 936). We use unsupervised ML on scanning electron microscopy–electron dispersive spectroscopy (SEM-EDS) measurements to define compositional categories based on differences in elemental abundances. We then test the ability of four ML algorithms─k-nearest neighbors (KNN), recursive partitioning and regressive trees (RPART), eXtreme gradient boost (XGBoost), and random forest (RF)─to classify the ToF-SIM spectra using (1) the categories assigned via SEM-EDS, (2) organic and inorganic labels assigned via SEM-EDS, and (3) the presence or absence of detectable steranes in ToF-SIMS spectra. In terms of predictive accuracy and balanced accuracy, KNN was the best performing model and RPART the worst. The feature importance, or the specific features of the ToF-SIM spectra used by the models to make classifications, cannot be determined for KNN, preventing posthoc model interpretation. Nevertheless, the feature importance extracted from the other models was useful for interpreting spectra. In conclusion, we determined that some of the organic ions used to classify biomarker containing spectra may be fragment ions derived from kerogen which is abundant in this mudstone sample.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Mitigating contamination in LSS surveys: a comparison of methods

ABSTRACT Future large-scale structure surveys will measure the locations and shapes of billions of galaxies. The precision of such catalogues will require meticulous treatment of systematic contamination of the observed fields. We compare several existing methods for removing such systematics from galaxy clustering measurements. We show how all the methods, including the popular pseudo-Cℓ Mode Projection and Template Subtraction methods, can be interpreted under a common regression framework and use this to suggest improved estimators. We show how methods designed to mitigate systematics in the power spectrum can be used to produce clean maps, which are necessary for cosmological analyses beyond the power spectrum, and we extend current methods to treat the next-order multiplicative contamination in observed maps and power spectra, which reduced power spectrum errors from $\Delta \chi ^2_{\rm C_\ell }\simeq 10$ to ≃ 1 in simulated analyses. Two new mitigation methods are proposed, which incorporate desirable features of current state-of-the-art methods while being simpler to implement. Investigating the performance of all the methods on a common set of simulated measurements from Year 5 of the Dark Energy Survey, we test their robustness to various analysis cases. Our proposed methods produce improved maps and power spectra when compared to current methods, while requiring almost no user tuning. We end with recommendations for systematics mitigation in future surveys, and note that the methods presented are generally applicable beyond the galaxy distribution to any field with spatial systematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

HostSub_GP: Precise Galaxy Background Subtraction in Transient Long-slit Spectroscopy with Gaussian Processes

We present a novel host galaxy subtraction technique in long-slit spectroscopy for extragalactic transients. Unlike classic methods which generally estimate the background using simple interpolation of local galaxy flux in the 2D spectrum, our approach leverages multi-band archival images of the host galaxies to model the background emission from the galaxy in the 2D spectrum. Such imaging encodes the wavelength-dependent galaxy profile along the slit, and is readily accessible through wide-field imaging surveys. We construct a smooth prior for the 2D galaxy profile with a Gaussian process (GP) based on these reference images, and use another GP to model the correlated deviations from the prior in the observed spectrum. This enables accurate inference of the galaxy flux blended with the transient. On synthetic long-slit data of a spiral galaxy extracted from a Multi Unit Spectroscopic Explorer hyper-spectral cube, the GP method remains robust as long as the host galaxy is spatially resolved and consistently outperforms classic methods. We apply the method to archival Keck spectra of two real transients, SN 2019eix and AT 2019qiz, to further demonstrate how the method uniquely recovers weak spectral features amid strong galaxy contamination, enabling refined constraints on the properties of both transients. We have released the software implementation, HostSub_GP, a scalable toolkit that leverages JAX, with an MIT license.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1999

Spatial distribution of venous gas emboli in the lungs

The distribution of gaseous pulmonary emboli is presumed to be determined by their buoyancy. We hypothesized that regional pulmonary blood flow may also influence their distribution. Therefore, pulmonary blood flow was measured in supine, anesthetized dogs with use of 15-microm fluorescent microspheres at baseline and during N(2) embolism. The animals were killed, and the lungs were excised, air-dried, and diced into approximately 2-cm(3) pieces with weights and spatial coordinates recorded. Embolism was defined as a >10% flow decrease relative to baseline. Vertically, the incidence of embolism increased substantially by 6 +/- 1% per additional centimeter in height compared with baseline (P = 0.0003). Embolism also increased radially by 3 +/- 1%/cm from the hilum (P = 0.002). There was a weaker but statistically significant increase in embolism to pieces with greater baseline flow, 9 +/- 2% for every 1. 0 increase in relative baseline flow (P = 0.008). We conclude that the distribution of gaseous emboli is influenced by buoyancy and flow dynamics within the pulmonary vasculature.

NASA Discipline Cardiopulmonary↗
DOE OSTI2023

ClimSim: A large multi-scale dataset for hybrid physics-ML climate emulation

Modern climate projections lack adequate spatial and temporal resolution due to computational constraints. A consequence is inaccurate and imprecise predictions of critical processes such as storms. Hybrid methods that combine physics with machine learning (ML) have introduced a new generation of higher fidelity climate simulators that can sidestep Moore’s Law by outsourcing compute-hungry, short, high-resolution simulations to ML emulators. However, this hybrid ML-physics simulation approach requires domain-specific treatment and has been inaccessible to MLexperts because of lack of training data and relevant, easy-to-use workflows. Wepresent ClimSim, the largest-ever dataset designed for hybrid ML-physics research. It comprises multi-scale climate simulations, developed by a consortium of climate scientists and ML researchers. It consists of 5.7 billion pairs of multivariate input and output vectors that isolate the influence of locally-nested, high-resolution, high-fidelity physics on a host climate simulator’s macro-scale physical state. The dataset is global in coverage, spans multiple years at high sampling frequency, and is designed such that resulting emulators are compatible with downstream coupling into operational climate simulators. We implement a range of deterministic and stochastic regression baselines to highlight the ML challenges and their scoring. The data (https://huggingface.co/datasets/LEAP/ClimSim_high-res2) and code(https://leap-stc.github.io/ClimSim)arereleasedopenlytosupport the development of hybrid ML-physics and high-fidelity climate simulations for the benefit of science and society.

artificial intelligence, machine learning↗
NASA NTRS2009

Watershed and Hydrodynamic Modeling for Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation and Seagrasses in Mobile Bay

There is a continued need to understand how human activities along the northern Gulf of Mexico coast are impacting the natural ecosystems. The gulf coast is experiencing rapid population growth and associated land cover/land use change. Mobile Bay, AL is a designated pilot region of the Gulf of Mexico Alliance (GOMA) and is the focus area of many current NASA and NOAA studies, for example. This is a critical region, both ecologically and economically to the entire United States because it has the fourth largest freshwater inflow in the continental USA, is a vital nursery habitat for commercially and recreational important fisheries, and houses a working waterfront and port that is expanding. Watershed and hydrodynamic modeling has been performed for Mobile Bay to evaluate the impact of land use change in Mobile and Baldwin counties on the aquatic ecosystem. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use Scenarios in 1948, 1992, 2001, and 2030. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on observed trends. All land use scenarios were developed to a common land classification system developed by merging the 1992 and 2001 National Land Cover Data (NLCD). The LSPC model output provides changes in flow, temperature, sediments and general water quality for 22 discharge points into the Bay. These results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment concentrations on a grid with four vertical profiles throughout the Bay s aquatic ecosystems. The models were calibrated using in-situ data collected at sampling stations in and around Mobile bay. This phase of the project has focused on sediment modeling because of its significant influence on light attenuation which is a critical factor in the health of submerged aquatic vegetation. The impact of land use change on sediment concentrations was evaluated by analyzing the LSPC and EFDC sediment simulations for the four land use scenarios. Such analysis was also performed for storm and non-storm periods. In- situ data of total suspended sediments (TSS) and light attenuation were used to develop a regression model to estimate light attenuation from TSS. This regression model was used to derive marine light attenuation estimates throughout Mobile bay using the EFDC TSS outputs. The changes in sediment concentrations and associated impact on light attenuation in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on seagreasses and Submerged Aquatic Vegetation (SAV) habitat. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of sediments due to land use driven flow changes with the restoration potential of SAVs.

Estes, Maurice G.↗
DOE OSTI2023

Biogeophysical Effects of Land-Use and Land-Cover Change Not Detectable in Warmest Month

Land-use and land-cover changes (hereafter simply “land use”) alter climates biogeophysically by affecting surface fluxes of energy and water. Yet, near-surface temperature responses to land use across observational versus model-based studies and spatial-temporal scales can be inconsistent. Here we assess the prevalence of the historical land use signal of daily maximum temperatures averaged over the warmest month of the year (t LU ) using regularized optimal fingerprinting for detection and attribution. We use observations from the Climatic Research Unit and Berkeley Earth alongside historical simulations with and without land use from phase 6 of the Coupled Model Intercomparison Project to reconstruct an experiment representing the effects of land use on climate. To assess the signal of land use at spatially resolved continental and global scales, we aggregate all input data across reference regions and continents, respectively. At both scales, land use does not comprise a significantly detectable set of forcings for two of four Earth system models and their multimodel mean. Furthermore, using a principal component analysis, we find that t LU is mostly composed of the nonlocal effects of land use rather than its local effects. These findings show that, at scales relevant for climate attribution, uncertainties in Earth system model representations of land use are too high relative to the effects of internal variability to confidently assess land use.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Microstrip Ring Resonator for Soil Moisture Measurements

Accurate determination of spatial soil moisture distribution and monitoring its temporal variation have a significant impact on the outcomes of hydrologic, ecologic, and climatic models. Development of a successful remote sensing instrument for soil moisture relies on the accurate knowledge of the soil dielectric constant (epsilon(sub soil)) to its moisture content. Two existing methods for measurement of dielectric constant of soil at low and high frequencies are, respectively, the time domain reflectometry and the reflection coefficient measurement using an open-ended coaxial probe. The major shortcoming of these methods is the lack of accurate determination of the imaginary part of epsilon(sub soil). In this paper a microstrip ring resonator is proposed for the accurate measurement of soil dielectric constant. In this technique the microstrip ring resonator is placed in contact with soil medium and the real and imaginary parts of epsilon(sub soil) are determined from the changes in the resonant frequency and the quality factor of the resonator respectively. The solution of the electromagnetic problem is obtained using a hybrid approach based on the method of moments solution of the quasi-static formulation in conjunction with experimental data obtained from reference dielectric samples. Also a simple inversion algorithm for epsilon(sub soil) = epsilon'(sub r) + j(epsilon"(sub r)) based on regression analysis is obtained. It is shown that the wide dynamic range of the measured quantities provides excellent accuracy in the dielectric constant measurement. A prototype microstrip ring resonator at L-band is designed and measurements of soil with different moisture contents are presented and compared with other approaches.

Sarabandi, Kamal↗
NASA NTRS2002

Smooth Pursuit Saccade Amplitude Modulation During Exposure to Microgravity

Russian investigators have reported changes in pursuit tracking of a vertically moving point stimulus during space flight. Early in microgravity, changes were manifested by decreased eye movement amplitude (undershooting) and the appearance of correction saccades. As the flight progressed, pursuit of the moving point stimulus deteriorated while associated saccadic movements were unchanged. Immediately postflight there was an improved execution of active head movements indicating that the deficiencies in pursuit function noted in microgravity may be of central origin. In contrast, tests of two cosmonauts showed that horizontal and vertical smooth pursuit were unchanged inflight. However, results of corresponding saccadic tasks showed a tendency toward the overshooting of a horizontal target early inflight with high accuracy developing later inflight, accompanied by an increased saccade velocity and a trend toward decreased saccade latency. Based on these equivocal results, we have further investigated the effects of space flight on the smooth pursuit mechanism during and after short duration flight, and postflight on returning MIR crewmembers. Sinusoidal target movement was presented horizontally at frequencies of 0.33 and 1.0 Hz. Subjects were asked to perform two trials for each stimulus combination: (1) moving eyes-only (EO) and (2) moving eyes and head (EH) with the target motion. Peak amplitude was 30 deg for 0.33 Hz trials and 15 deg for the 1.0 Hz trials. The relationship between saccade amplitude and peak velocity were plotted as a main sequence for each phase of flight, and linear regression analysis allowed us to determine the slope of each main sequence plot. The linear slopes were then combined for each flight phase for each individual subject. The main sequence for both EO and EH trials at both the 0.33 and 1.0 Hz frequencies during flight for the short duration flyers showed a reduction in saccade velocity and amplitude when compared to the preflight main sequence . This difference in the regression slopes between flight phase, head/eye condition (EO or EH), and pursuit target frequency was observed across all subjects (statistically significant at the p<0.02, df= 2). It is interesting to note that postflight for the short duration flyers there was an immediate recovery to the preflight main sequence across all trials. There were no significant differences observed between the preflight slopes for either head movement condition (EO vs. EH). When the immediate postflight (R+O) regression slopes were compared with the preflight slopes, there was a tendency (not significant) for both saccade amplitude and peak velocity to increase during the postflight testing. This tendency had vanished by R+ 1. Of particular interest was the redistribution of saccades during the latter stages of the flight and immediately postflight in the EO condition. At the 1.0 Hz frequency the saccades tended to be clustered near the lowest target velocity. It was also interesting to note that gaze performance (eye in skull + head in space) was consistently better during the EH condition; a finding also observed by our Russian colleagues. As the results of the long duration flight become available we expect that they will not only show that postflight effects will be similar to those observed during the short duration flights, but will also last for a greater period of time following flight. It is not clear what mechanism is responsible for the decreased peak saccadic velocity during flight unless the change is related to the control of retinal slip. For example, it is possible that saccades will tend to initially undershoot their targets by a small percentage and these saccades are then followed, if vision is available, by a small augmenting corrective saccade. It has been postulated that the functional significance of this undershooting tendency is to maintain the spatial representation of the target on the same side of the fovea (as opposedo racing across the fovea) and hence in the same cerebral hemisphere that initiated the primary saccade thus minimizing delays caused by an intra-hemispheric transfer of information . One could also speculate that with saccade velocities greater than normal, additional corrective saccades would be required to bring the target back on the fovea. A less plausible explanation of our findings could be fatigue. Yet it seems unlikely that our subjects would show lower velocities on all inflight test days while showing increased saccade velocities immediately following space flight where fatigue is usually the greatest. Finally, the redistribution effect noted late inflight is likely caused by adaptive changes. Overall, corrective saccades appeared to be used in maintaining gaze on target; reducing retinal slip and assisting space travelers in maintaining clear vision throughout the different phases of the space flight.

Reschke, M. F.↗
DOE OSTI2021

Understanding the Drivers of Mobility during the COVID-19 Pandemic in Florida, USA Using a Machine Learning Approach

As of March 2021, the State of Florida, U.S.A. had accounted for approximately 6.67% of total COVID-19 (SARS-CoV-2 coronavirus disease) cases in the U.S. The main objective of this research is to analyze mobility patterns during a three month period in summer 2020, when COVID-19 case numbers were very high for three Florida counties, Miami-Dade, Broward, and Palm Beach counties. To investigate patterns, as well as drivers, related to changes in mobility across the tri-county region, a random forest regression model was built using sociodemographic, travel, and built environment factors, as well as COVID-19 positive case data. Mobility patterns declined in each county when new COVID-19 infections began to rise, beginning in mid-June 2020. While the mean number of bar and restaurant visits was lower overall due to closures, analysis showed that these visits remained a top factor that impacted mobility for all three counties, even with a rise in cases. Our modeling results suggest that there were mobility pattern differences between counties with respect to factors relating, for example, to race and ethnicity (different population groups factored differently in each county), as well as social distancing or travel-related factors (e.g., staying at home behaviors) over the two time periods prior to and after the spike of COVID-19 cases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Bayesian time-varying occupancy model for West Nile virus in Ontario, Canada

Occupancy models determine the true presence or absence of a species by adjusting for imperfect detection in surveys. They often assume that species presences can be detected only if sites are occupied during a sampling season. We extended these models to estimate occupancy rates that vary throughout a sampling season as well as account for spatial dependence among sites. For these methods, we constructed a fast Gibbs sampler with the Pólya-Gamma augmentation strategy to conduct inference on covariate effects. We applied these methods to evaluate how environmental conditions and surveillance practices are associated with the presence of West Nile virus in mosquito traps across Ontario, Canada from 2002 to 2017. We found that urban land cover and warm temperatures drove viral occupancy, whereas viral testing on pools with higher proportions of Culex mosquitoes was more likely to result in a positive test for West Nile virus. Models with time-varying occupancy effects achieved much lower Watanabe-Akaike information criteria than models without such effects. Our final model had strong predictive performance on test data that included some of the most extreme seasons, demonstrating the promise of these methods in the study of pathogens spread by mosquito vectors.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1994

Calibration of TOPEX/POSEIDON at Platform Harvest

We present estimates for the mean bias of the TOPEX/POSEIDON NASA altimeter (ALT) and the Centre National d'Etudes Spatiales altimeter (SSALT) using in-situ data gathered at Platform Harvest during the first 36 cycles of the mission. Data for 21 overflights of the ALT and six overflights of the SSALT have been analyzed. The analysis includes an independent assessment of in-situ measurements of sea level, the radial component of the orbit, wet tropospheric path delay, and ionospheric path delay. (The sign convention used is such that, to correct the geophysical data record values for sea level, add the bias algebraically. Unless otherwise stated, the uncertainty in a given parameter is depicted by +/- sigma(sub x), where sigma(sub x) is the sample standard deviation of x about the mean.) Tide gauges at Harvest provide estimates of sea level with an uncertainty of +/- 1.5 cm. The uncertainty in the radial component of the orbit is estimated to be +/- 1.3 cm. In-situ measurements of tropopsheric path delay at Harvest compare to within +/- 1.3 cm of the TOPEX/POSEIDON microwave radiometer, and in-situ measurements of the ionospheric path delay compare to within -0.4 +/- 0.7 cm of the dual-frequency ALT and 1.1 +/- 0.6 cm of Doppler orbitography and radiopositioning integrated by satellite. We obtain mean bias estimates of -14.5 +/- 2.9 cm for the ALT and +0.9 +/- 3.1 cm for the SSALT (where the uncertainties are based on the standard deviation of the estimated mean (sigma(sub bar x/y), which is derived from sample statistics and estimates for errors that cannot be observed). These results are consistent with independent estimates for the relative bias between the two altimeters. A linear regression applied to the complete set of data shows that there is a discernable secular trend in the time series for the ALT bias estimates. A preliminary analysis of data obtained through cycle 48 suggests that the apparent secular drift may be the result of a poorly sampled annual signal.

Christensen, E. J.↗
DOE OSTI2024

Studying single-electron traps in newly fabricated Skipper-CCDs for the Oscura experiment using the pocket-pumping technique

Understanding and characterizing very low-energy (⁠ ∼eV) background sources is a must in rare-event searches. Oscura, an experiment aiming to probe electron recoils from sub-GeV dark matter using a 10 kg skipper-CCD detector, has recently fabricated its first two batches of sensors. In this work, we present the characterization of defects/contaminants identified in the buried-channel region of these newly fabricated skipper-CCDs. These defects/contaminants produce deferred charge from trap emission in the images next to particle tracks, which can be spatially resolved due to the sub-electron resolution achieved with these sensors. Using the trap-pumping technique, we measured the energy and cross section associated with these traps in three Oscura prototype sensors from different fabrication batches which underwent different gettering methods during fabrication. Results suggest that the type of defects/contaminants is more closely linked to the fabrication batch rather than to the gettering method used. The exposure-dependent single-electron rate (SER) of one of these sensors was measured ∼100 m underground, yielding (1.8±0.3)×10−3e−/pix/day at 131 K. The impact of the identified traps on the measured exposure-dependent SER is evaluated via a Monte Carlo simulation. Results suggest that the exposure-dependent SER of Oscura prototype sensors would be lower in lower background environments as expected.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗