Search NASA⌕ Search

SEARCH · Search NASA

Results for “map validation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,585 records · Page 88
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Museomics Dissects the Genetic Basis for Adaptive Seasonal Coloration in the Least Weasel

Abstract Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat coloration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage color moults across regions with varying winter snow. Whole-genome sequence data were obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris coloration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200-kb genomic region to coloration morph, which was validated by genotyping museum specimens from intermorph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change cosegregating with coloration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favored winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter color variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

Miranda, Inês↗
DOE OSTI2020

effect of ambient pressure on the behavior of single-component fuels in a gasoline multi-hole injector

The injection characteristics of neat ethanol and pure isooctane are studied under different ambient pressure and temperature conditions. Injection underflash-boiling conditions can enhance liquid atomization and evaporation, providing the possibility of improvement in the fuel/air mixing. These superheated conditions often introduce phenomena that are not taken into account in the standard modeling of sprays for engine applications. The present work proposes a numerical investigation of the behavior of Engine Combustion Network's 8-hole spray-G injector, starting at the subcooled nominal condition and reducing the ambient pressure at constant low temperature to reach the flare flash-boiling condition. To initialize the properties of the injected fuel, the flow in the nozzle is simulated with a Eulerian approach, handling the two phases with a mixture model and the phase change, due to cavitation and flash boiling, with the Homogenous Relaxation Model. A map of the mixture's kinematic and thermal behavior is obtained at the interface between the injector and the chamber to initialize the Lagrangian simulations. A literature-based vaporization model is implemented to obtain the proper description of the characteristic features of a multi-hole spray under super-heated conditions, like plume-plume interaction. The numerical representation of the spray is validated in terms of penetration and radial spreading on DBI images, reproducing the light attenuation profiles caused by the presence of the liquid spray. Simulations show that coupled nozzle flow and spray calculations capture the spray morphology and shape better compared to calculations performed without considering the nozzle flow simulation details, especially under flare flash conditions.

Nocivelli, Lorenzo↗
DOE OSTI2018

HIV-1 Envelope Glycoproteins from Diverse Clades Differentiate Antibody Responses and Durability among Vaccinees

Induction of broadly cross-reactive antiviral humoral responses with the capacity to target globally diverse circulating strains is a key goal for HIV-1 immunogen design. A major gap in the field is the identification of diverse HIV-1 envelope antigens to evaluate vaccine regimens for binding antibody breadth. In this study, we define unique antigen panels to map HIV-1 vaccine-elicited antibody breadth and durability. Diverse HIV-1 envelope glycoproteins were selected based on genetic and geographic diversity to cover the global epidemic, with a focus on sexually acquired transmitted/founder viruses with a tier 2 neutralization phenotype. Unique antigenicity was determined by nonredundancy (Spearman correlation), and antigens were clustered using partitioning around medoids (PAM) to identify antigen diversity. Cross-validation demonstrated that the PAM method was better than selection by reactivity and random selection. Analysis of vaccine-elicited V1V2 binding antibody in longitudinal samples from the RV144 clinical trial revealed the striking heterogeneity among individual vaccinees in maintaining durable responses. These data support the idea that a major goal for vaccine development is to improve antibody levels, breadth, and durability at the population level. Elucidating the level and durability of vaccine-elicited binding antibody breadth needed for protection is critical for the development of a globally efficacious HIV vaccine.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

AlloSHP: deconvoluting single homeologous polymorphism for phylogenetic analysis of allopolyploids

Background The genomic and evolutionary study of allopolyploid organisms involves multiple copies of homeologous chromosomes, making their assembly, annotation, and phylogenetic analysis challenging. Bioinformatics tools and protocols have been developed to study polyploid genomes, but sometimes require the assembly of their genomes, or at least the genes, limiting their use. Results We have developed AlloSHP, a command-line tool for detecting and extracting single homeologous polymorphisms (SHPs) from the subgenomes of allopolyploid species. This tool integrates three main algorithms, WGA, VCF2ALIGNMENT and VCF2SYNTENY, and allows the detection of SHPs for the study of diploid-polyploid complexes with available diploid progenitor genomes, without assembling and annotating the genomes of the allopolyploids under study. AlloSHP has been validated on three diploid-polyploid plant complexes, Brachypodium, Brassica, and Triticum-Aegilops, and a set of synthetic hybrid yeasts and their progenitors of the genus Saccharomyces. The results and congruent phylogenies obtained from the four datasets demonstrate the potential of AlloSHP for the evolutionary analysis of allopolyploids with a wide range of ploidy and genome sizes. Conclusions AlloSHP combines the strategies of simultaneous mapping against multiple reference genomes and syntenic alignment of these genomes to call SHPs, using as input data a single VCF file and the reference genomes of the known or closest extant diploid progenitor species. This novel approach provides a valuable tool for the evolutionary study of allopolyploid species, both at the interspecific and intraspecific levels, allowing the simultaneous analysis of a large number of accessions and avoiding the complex process of assembling polyploid genomes.

Allopolyploids↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Stress Relaxation Related to Spontaneous Thin Film Buckling: Correlation between Finite Element Calculations and Micro Diffraction Analysis

Compressive residual stresses generated during thin film deposition may lead to undesirable film damage, such as delamination, buckling, and flaking, ultimately leading to the failure of the device employing the film. Understanding the residual stress generation and role in these damage mechanisms is necessary to preserve thin film integrity and optimize its functional properties. Thin shell theory has been used for decades to predict buckling but the results have not yet been correlated with experimental data since the techniques used to measure stress in metallic films were not able to do so at the required micron scale until recently. Micro scanning X-ray diffraction now enables the direct mapping of the local stress of metallic films. In this paper, finite element method based on thin shell theory and synchrotron X-ray micro diffraction have been used to determine stress maps of thin film buckling patterns. Calculations of the stress distribution in the metallic films have been performed taking into account the buckling geometry determined from optical measurements. Stress distributions over gold blisters and tungsten wrinkles obtained with the two techniques are in fair agreement and allow for the accurate determination of the stress relaxation profile from the bottom to the top of the buckling, validating the thin shell theory model.

36 MATERIALS SCIENCE↗
NASA NTRS1991

A hybrid numerical method for three-dimensional spatially-developing free-shear flows

The present novel algorithm for 3D incompressible Navier-Stokes equations treats a domain which is (1) infinite in the vertical direction, using a mapped spectral method; (2) finite in the streamwise direction, using a classical Fourier method; and (3) homogeneous in the spanwise direction, using high-order compact finite differencing. A projection method is employed which ensures the exact conservation of mass, as well as the satisfaction of the boundary conditions at infinity. The novel aspects of these methods are noted, and the code they constitute is validated in light of several test cases. Results are presented for two- and three-dimensional mixing layers.

Buell, Jeffrey C.↗
NASA NTRS2003

Preliminary Optical And Electric Field Pulse Statistics From Storm Overflights During The Altus Cumulus Electrification Study

The Altus Cumulus Electrification Study (ACES) was conducted during the month of August, 2002 in an area near Key West, Florida. One of the goals of this uninhabited aerial vehicle (UAV) study was to collect high resolution optical pulse and electric field data from thunderstorms. During the month long campaign, we acquired 5294 lightning generated optical pulses with associated electric field changes. Most of these observations were made while close to the top of the storms. We found filtered mean and median 10-10% optical pulse widths of 875 and 830 microns respectively while the 50-50% mean and median optical pulse widths are 422 and 365 microns respectively. These values are similar to previous results as are the 10-90% mean and median rise times of 327 and 265 microns. The peak electrical to optical pulse delay mean and median were 209 and 145 microns which is longer than one would expect from theoretical results. The results of the pulse analysis will contribute to further validation of the Optical Transient Detector (OTD) and the Lightning Imaging Sensor (LIS) satellites. Pre-launch estimates of the flash detection efficiency were based on a small sample of optical pulse measurements associated with less than 350 lightning discharges collected by NASA U-2 aircraft in the early 1980s. Preliminary analyses of the ACES measurements show that we have greatly increased the number of optical pulses available for validation of the LIS and other orbital lightning optical sensors. Since the Altus was often close to the cloud tops, many of the optical pulses are from low-energy pulses. From these low-energy pulses, we can determine the fraction of optical lightning pulses below the thresholds of LIS, OTD, and any future satellite-based optical sensors such as the geostationary Lightning Mapping Sensor.

Mach, D. A.↗
NASA NTRS2012

Advances in Spectral-Spatial Classification of Hyperspectral Images

Recent advances in spectral-spatial classification of hyperspectral images are presented in this paper. Several techniques are investigated for combining both spatial and spectral information. Spatial information is extracted at the object (set of pixels) level rather than at the conventional pixel level. Mathematical morphology is first used to derive the morphological profile of the image, which includes characteristics about the size, orientation, and contrast of the spatial structures present in the image. Then, the morphological neighborhood is defined and used to derive additional features for classification. Classification is performed with support vector machines (SVMs) using the available spectral information and the extracted spatial information. Spatial postprocessing is next investigated to build more homogeneous and spatially consistent thematic maps. To that end, three presegmentation techniques are applied to define regions that are used to regularize the preliminary pixel-wise thematic map. Finally, a multiple-classifier (MC) system is defined to produce relevant markers that are exploited to segment the hyperspectral image with the minimum spanning forest algorithm. Experimental results conducted on three real hyperspectral images with different spatial and spectral resolutions and corresponding to various contexts are presented. They highlight the importance of spectral–spatial strategies for the accurate classification of hyperspectral images and validate the proposed methods.

hyperspectral image↗
DOE OSTI2026

Design optimization of MAPS-based detectors using a data-driven fast simulation approach

A parametric simulation tool for pixel sensors is presented. A realistic pixel response is simulated purely based on measurement input, without requiring detailed knowledge of the underlying manufacturing process. As such, it provides an efficient alternative to the use of Technology Computer-Aided Design simulations, which typically depend on proprietary process information. Due to its parametric approach, the package is fast and thus particularly useful for larger detector systems and high hit rate environments. This work presents measurements, simulation and its validation for the MALTA2 sensor. It is a small collection electrode monolithic active pixel sensor produced in the Tower 180 nm complementary metal-oxide-semiconductor imaging process. Modifications to the sensor’s periphery, mainly in the hit merger, are studied in order to optimize the performance for tracking and calorimetry. This optimization is of special interest as part of the MALTA3 sensor redesign in the 65 nm Tower Partners Semiconductor Co. process.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

IRMA

IRMA (In)elastic Representation of Materials As S(α,β) evaluations IRMA turns one phonon model into three outputs that usually require three separate tool chains: an evaluated nuclear-data file, predicted neutron-scattering spectra, and scattering kernels for Monte Carlo transport. The three outputs draw on a single, consistent description of the material, so the evaluation, the spectroscopy that can validate it, and the transport that uses it always agree about the physics. Nuclear data. IRMA writes ENDF-6 File 7 thermal scattering evaluations on automatically constructed (α, β) grids. This part reimplements and generalizes NJOY's LEAPR: the classic kernels reproduce freshly generated NJOY2016 tapes digit for digit and published reference tapes to about 1e-4, and the generalized paths add the exact coherent one-phonon term, anisotropic Debye-Waller tensors, coherent elastic for arbitrary crystals, and a per-species partition for polyatomic materials. The tapes feed NJOY, AMPX, FUDGE, and every transport code downstream of them. Neutron spectroscopy. The irma.spectra forward model projects the same physics onto an instrument's kinematics and resolution: INS spectra for VISION and generic indirect geometries, and 2-D S(Q,E) powder maps for direct-geometry spectrometers, from a phonopy model or straight from a phonon DOS. It can be used to predict a proposed measurement before beam time; in analysis, it supplies the calculated single-scattering counterpart of a measured spectrum, from the same material description the evaluation was built from. Monte Carlo transport. The irma.ncrystal exporter writes per-temperature scattering kernels for the companion NCrystal plugin, so McStas, OpenMC, and other NCrystal-aware codes sample the same physics. The exported kernels carry the per-site anisotropic Debye-Waller tensors, keeping directional coherent-elastic physics that NCrystal's standard scalar treatment does not represent. With the same physics inside a transport code, an entire beamline becomes a virtual experiment: IRMA's end-to-end validation ran a custom McStas implementation of the ARCS spectrometer, assembled from the existing McVine and McStas models, against measured data. From a bare crystal structure. The irma mlip front end builds the phonon model itself: a structure file and a choice of potential are enough. Nine pretrained machine-learned interatomic potentials are supported, on a laptop CPU, with no first-principles calculation; an approximate phonon model for a new material costs minutes, not a DFT campaign, and the build emits prefilled inputs for all three outputs. The result is a good starting point rather than a finished evaluation: survey-quality physics with every parameter exposed for review. A converged atomistic calculation enters the same way, as a phonopy model, when higher fidelity is needed.

Ramic, Kemal [Oak Ridge National Laboratory (ORNL)↗
DOE OSTI2022

Genotype–environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence

Abstract Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush ( Catharus guttatus ), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Estimating snow cover from high-resolution satellite imagery by thresholding blue wavelengths

We report the extent and duration of snow cover, a critical component of the hydrologic cycle and the global climate system, is expected to shift dramatically under climate change. Therefore, developing high-resolution assessments of snow cover change is crucial for estimating the impact of changing snow cover on watershed and ecosystems processes in cold regions. Remote sensing tools provide a powerful method for mapping snow-covered area (SCA) across a landscape. The most common method for estimating SCA utilizes the normalized difference snow index (NDSI), which relies on spectral measurements in the shortwave-infrared wavelengths (SWIR). NDSI can effectively estimate catchment- to regional-scale SCA, but it cannot be used to assess fine-scale SCA because of current limitations on the spatial resolution of satellite-derived SWIR measurements. Here, we map SCA using a threshold of blue wavelengths and high-resolution satellite imagery. The thresholding method, which we call the Blue Snow Threshold algorithm (BST), has previously been used with digital camera imagery. We refine and automate the algorithm for use with cloud-free high-resolution satellite imagery and find that the BST can be used to assess fine-scale SCA. For validation, we compared BST-derived estimates of SCA to a) airborne lidar surveys, b) Landsat fractional SCA, and c) snow disappearance dates from Snow Telemetry (SNOTEL) stations. When compared to airborne lidar surveys of SCA, the BST predicted SCA had a range of F-scores between 0.81 and 0.94 in four study areas in California and Colorado. We also found general agreement between SCA and snow disappearance at multiple SNOTEL sites across the western United States. Given the relatively recent availability of high-resolution satellite imagery with spectral measurements in the visible wavelengths but lacking in SWIR, the BST offers a reliable and easy-to-apply tool for examining fine-scale snow-related processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Robust Neural Network for Extracting Dynamics from Electrostatic Force Microscopy Data

Advances in scanning probe microscopy (SPM) methods such as time-resolved electrostatic force microscopy (trEFM) now permit the mapping of fast local dynamic processes with high resolution in both space and time, but such methods can be time-consuming to analyze and calibrate. Here, we design and train a regression neural network (NN) that accelerates and simplifies the extraction of local dynamics from SPM data directly in a cantilever-independent manner, allowing the network to process data taken with different cantilevers. We validate the NN’s ability to recover local dynamics with a fidelity equal to or surpassing conventional, more time-consuming, calibrations using both simulated and real microscopy data. We apply this method to extract accurate photoinduced carrier dynamics on n = 1 butylammonium lead iodide, a halide perovskite semiconductor film that is of interest for applications in both solar photovoltaics and quantum light sources. Lastly, we use SHapley Additive exPlanations to evaluate the robustness of the trained model, confirm its cantilever-independence, and explore which parts of the trEFM signal are important to the network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Lunar Prospector Orbit Determination Uncertainties Using the High Resolution Lunar Gravity Models

The Lunar Prospector (LP) mission began on January 6, 1998, when the LP spacecraft was launched from Cape Canaveral, Florida. The objectives of the mission were to determine whether water ice exists at the lunar poles, generate a global compositional map of the lunar surface, detect lunar outgassing, and improve knowledge of the lunar magnetic and gravity fields. Orbit determination of LP performed at the Jet Propulsion Laboratory (JPL) is conducted as part of the principal science investigation of the lunar gravity field. This paper will describe the JPL effort in support of the LP Gravity Investigation. This support includes high precision orbit determination, gravity model validation, and data editing. A description of the mission and its trajectory will be provided first, followed by a discussion of the orbit determination estimation procedure and models. Accuracies will be examined in terms of orbit-to-orbit solution differences, as a function of oblateness model truncation, and inclination in the plane-of-sky. Long term predictions for several gravity fields will be compared to the reconstructed orbits to demonstrate the accuracy of the orbit determination and oblateness fields developed by the Principal Gravity Investigator.

Carranza, Eric↗
NASA NTRS2001

TRMM Validation: Integration of TRMM Satellite and Ground-Validation Data to Determine Tropical Heating

Steiner and Houze showed from ground validation data that the Tropical Rain Measuring Mission (TRMM) satellite Precipitation Radar's (PR's) twice daily only sampling should lead to an uncertainty of approximately 20% in rain estimates. They further showed that the uncertainties are smallest at the 5-7.5 km level. Schumacher and Houze used Kwajalein ground validation data to show that the TRMM PR misses only 2.3% of the near surface rainfall but does not see 46% of the area where rain occurs, because of the 17 dBZ PR reflectivity threshold. Houze discusses how the TRMM data extend earlier tropical convective studies to global coverage of the vertical profile of latent heating via the TRMM PR''s ability to distinguish and globally map convective and stratiform precipitation. Process studies carried out under this TRMM grant Yuter and Houze and Yuter et al. studied ship-based radar observations in the tropical eastern Pacific ITCZ. The eastern Pacific precipitation process is different from the western Pacific (the COARE area); rain is heavier but the clouds are not as deep. These process differences may affect the ability to remotely sense precipitation accurately in the two regions. Satellite microwave data were able to detect the precipitation as long as the rain areas exceeded 10 km in dimension. However, the microwave algorithms had difficulty distinguishing light and heavy rain. Satellite IR algorithms only partially detected the rain because the tops of the smaller and more short-lived rain clouds were sometimes not cold enough for the IR algorithms to detect them. Houze et al. focused on the west Pacific precipitating mesoscale convective systems and showed how their precipitation and internal dynamics vary in relation to the slowly varying large-scale heating-driven circulation, which has a structure described by a combination of Kelvin and Rossby wave response to the near-equatorial convective heating constituted by the mesoscale convective systems. Ship and aircraft radar data were used in this study.

Houze, Robert A., Jr.↗
NASA NTRS2006

Experimental Validation of a Closed Brayton Cycle System Transient Simulation

The Brayton Power Conversion Unit (BPCU) is a closed cycle system with an inert gas working fluid. It is located in Vacuum Facility 6 at NASA Glenn Research Center. Was used in previous solar dynamic technology efforts (SDGTD). Modified to its present configuration by replacing the solar receiver with an electrical resistance heater. The first closed-Brayton-cycle to be coupled with an ion propulsion system. Used to examine mechanical dynamic characteristics and responses. The focus of this work was the validation of a computer model of the BPCU. Model was built using the Closed Cycle System Simulation (CCSS) design and analysis tool. Test conditions were then duplicated in CCSS. Various steady-state points. Transients involving changes in shaft rotational speed and heat input. Testing to date has shown that the BPCU is able to generate meaningful, repeatable data that can be used for computer model validation. Results generated by CCSS demonstrated that the model sufficiently reproduced the thermal transients exhibited by the BPCU system. CCSS was also used to match BPCU steady-state operating points. Cycle temperatures were within 4.1% of the data (most were within 1%). Cycle pressures were all within 3.2%. Error in alternator power (as much as 13.5%) was attributed to uncertainties in the compressor and turbine maps and alternator and bearing loss models. The acquired understanding of the BPCU behavior gives useful insight for improvements to be made to the CCSS model as well as ideas for future testing and possible system modifications.

Johnson, Paul K.↗