Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data- limited”

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 163 records · Page 9

The light curve at 10 microns of Algol near secondary minimum

Observations of Algol (Beta Persei) at 10 microns are presented which are of relevance to the system's behavior around the time of its secondary minimum. No evidence is found for a significant excess at 10 microns. The data show a smooth decline by 0.3 mag which starts at about phase 0.42, while the limited data obtained at phases greater than 0.50 are consistent with a symmetrical curve about phase 0.50.

Nadeau, D.↗

HEPA Filter Age Evaluation Report

In 2020, PNNL completed a literature search for high-efficiency particulate air (HEPA) filter age information. HEPA filters are used in many operations to remove particulate matter from effluent exhaust streams. They are thought to degrade over time both during proper storage and during normal operational service; however, the rate at which the filters degrade remains unknown. This brings into question if age is an adequate indicator of HEPA filter performance. Data from six previous reports were obtained from the literature search and combined to create a data set of 1600 operating filters. Filter usage was identified from multiple facilities. The various types of filters (e.g., axial flow, self-contained, and standard 24 x 24 x 11.5 inches; and both separator and separatorless) reported were constructed to the requirements Section FC (HEPA Filters) or Section FK (Special HEPA Filters) of the American Society of Mechanical Engineers AG-1 code. Filters were presumed to have continuously met the operational criteria and in particular passed both annual efficiency tests and annual DP measurements. The environmental conditions within the exhaust system were also assumed adequate for long-term filter operation. The data, by the nature of the reports, excludes rejected filters from quality assurance evaluations, intake, or installation testing. The collective data set shows over half of were operating past the current 10-year Department of Energy (DOE) limit. Data was evaluated for age lifetime using a linear trendline, survival function, probability functions, and failure rate; financial impacts were also addressed. This report supports the notion that HEPA filters can operate safely and efficiently under proper maintenance well past the 10-year lifetime established by DOE. Using the results of the four age evaluation approaches, they collectively point to the reasonableness of an operating HEPA filter lifetime of 20 years. The results are not necessarily definitive, but nevertheless, they are promising. The analysis provides reasonable assurance that when implemented using a graded approach with well-defined performance and operational requirements, extending the service life for HEPA filters beyond 10 years is low risk.

42 ENGINEERING↗

Extrapolation of the Rainflow-Counted Load Ranges for Fatigue Assessment of the Wind Turbine's Blades

Wind turbine design standards recommend the use of statistical modeling coupled with extrapolation of the short-term load data to long-term periods for fatigue reliability assessment. However, statistical error and computational expense can limit the accuracy of such approaches. In the case of wind turbine blades, the errors are more significant because of the high material fatigue exponent that makes the damage estimations more sensitive to variations. In addition, due to different excitation sources, the flapwise load range histogram is not unimodal, and thus its statistical modeling is complex. In the present work, we provide three methods for statistical modeling of the flapwise bending moment ranges including a novel approach based on frequency-based separation of the modes. The first two methods are simplified approaches for modeling the most crucial load ranges using unimodal distributions and the third method involves multimodal distribution fitting. The research is based on 3600 10-minute aeroelastic simulations of DTU 10MW case study wind turbine from which a benchmark damage equivalent load (DEL) is calculated. The DEL calculated by each of the three proposed methods is compared to this reference. The results show that the conventional approach based on using 6 seeds as well as using mixture models fitted on the limited data lead to under-conservative results with errors up to 23%. On the other hand, the simplified unimodal approaches provided in this work can provide conservative estimations of the fatigue damage with mean values 5% and 12% higher than the benchmark. However, the variability of the DEL estimates is higher when using unimodal extrapolation of the load ranges, and the data can be conservative by 17.5%. The proposed unimodal fits suggested for modeling and extrapolation of the blade's load ranges provide less errors relatively and most importantly conservative DEL estimations while maintaining computational efficiency.

blade fatigue↗

A comparative statistical study of long-term agroclimatic conditions affecting the growth of US winter wheat: Distributions of regional monthly average precipitation on the Great Plains and the state of Maryland and the effect of agroclimatic conditions on yield in the state of Kansas

A histogram analysis of average monthly precipitation over 30 and 84 year periods for both Maryland and Kansas was made and the results compared. A second analysis, a statistical assessment of the effect of average monthly precipitation on Kansas winter wheat yield was made. The data sets covered the three periods of 1941-1970, 1887-1970, and 1887-1921. Analyses of the limited data sets used (only the average monthly precipitation and temperature were correlated against yield) indicated that fall precipitation values, especially those of September and October, were more important to winter wheat yield than were spring values, particularly for the period 1941-1970.

Welker, J.↗

Global Carbon Budget 2016

Accurate assessment of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere the global carbon budget is important to better understand the global carbon cycle, support the development of climate policies, and project future climate change. Here we describe data sets and methodology to quantify all major components of the global carbon budget, including their uncertainties, based on the combination of a range of data, algorithms, statistics, and model estimates and their interpretation by a broad scientific community. We discuss changes compared to previous estimates and consistency within and among components, alongside methodology and data limitations. CO2 emissions from fossil fuels and industry (EFF) are based on energy statistics and cement production data, respectively, while emissions from land-use change (ELUC), mainly deforestation, are based on combined evidence from land-cover change data, fire activity associated with deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the annual changes in concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. The variability in SOCEAN is evaluated with data products based on surveys of ocean CO2 measurements. The global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms of the global carbon budget and compared to results of independent dynamic global vegetation models. We compare the mean land and ocean fluxes and their variability to estimates from three atmospheric inverse methods for three broad latitude bands. All uncertainties are reported as +/- 1(sigma), reflecting the current capacity to characterize the annual estimates of each component of the global carbon budget. For the last decade available (2006-2015), EFF was 9.3+/-0.5 GtC/yr, ELUC 1.0+/-0.5 GtC/yr,GATM 4.5+/-0.1 GtC/yr, SOCEAN 2.6+/-0.5 GtC/yr, and SLAND 3.1+/-0.9 GtC/yr. For year 2015 alone, the growth in EFF was approximately zero and emissions remained at 9.9+/-0.5 GtC/yr, showing a slowdown in growth of these emissions compared to the average growth of 1.8/yr that took place during 2006-2015.Also, for 2015, ELUC was 1.3+/-0.5 GtC/yr, GATM was 6.3+/-0.2 GtC/yr, SOCEAN was 3.0+/-0.5 GtC/yr, and SLAND was 1.9+/-0.9 GtC/yr. GATM was higher in 2015 compared to the past decade (2006-2015), reflecting a smaller SLAND for that year. The global atmospheric CO2 concentration reached 399.4+/-0.1 ppm averaged over 2015. For 2016, preliminary data indicate the continuation of low growth in EFF with +0.2% (range of -1.0 to +1.8% ) based on national emissions projections for China and USA, and projections of gross domestic product corrected for recent changes in the carbon intensity of the economy for the rest of the world. In spite of the low growth of EFF in 2016, the growth rate in atmospheric CO2 concentration is expected to be relatively high because of the persistence of the smaller residual terrestrial sink (SLAND) in response to El Nino conditions of 2015-2016. From this projection of EFF and assumed constant ELUC for 2016, cumulative emissions of CO2 will reach 565+/-55 GtC (2075+/-205 GtCO2) for 1870-2016, about 75% from EFF and 25% from ELUC. This living data update documents changes in the methods and data sets used in this new carbon budget compared with previous publications of this data set.

Quéré, Corinne Le↗

Direct Feed High-Level Waste APPS Model Glass Testing (DFHLW APPS) Matrix, Phase 2

This report summarizes the data collected during the batching and melting of a second matrix of Direct Feed High-Level Waste (DFHLW) glasses generated using the preliminary enhanced waste glass models (EWG2.5) and the Britton and Anderson (2024) preliminary DFHLW feed vector. The purpose of these glasses is two-fold: 1. Validate EWG2.5 glass calculations being used in the Aspen Process Performance Simulation (APPS) model. 2. Evaluate and ultimately improve the glass property models and formulation methods used for design of DFHLW glasses as part of an iterative process of data collection and model refinement. Some of the 16 APPS2 glasses tested did not satisfy all target property constraints due to the limited data on DFHLW glass supporting the EWG2.5 models. • One glass, APPS2-10, formed nepheline on canister centerline cooling (CCC) heat-treatment and failed the product consistency test (PCT) response limits. This glass also had high B and Cr release rates for the toxicity characteristic leaching procedure (TCLP). All other glasses were found to satisfy the PCT and TCLP constraints for both quenched and CCC samples. • One glass, APPS2-08, had higher than acceptable viscosity due to magnetite crystallization. • One glass, APPS2-09, formed greater than 2 vol% crystals at 950 °C. As the glass design criterion was that the temperature at 2 vol% crystal (T 2% ) be less than 950 °C, only one glass failed the criteria. However, this criterion is being reevaluated. Four additional glasses formed crystal fractions between 1 and 2 vol% at 950 °C (APPS2-03, -08, -12, and -14). • Four glasses – APPS2-01, -02, -04, and -16 – failed the Monofrax K-3 refractory neck corrosion (k neck ) design limit of 0.04 in. at 1208 °C for 6 d. This is another criterion being reevaluated. Four additional glasses (APPS2-05, -06, -11, and -13) exhibited 0.025 = k neck = 0.04 in. • All 16 glasses passed the sulfur solubility and TCLP constraints. The measured property values were compared to predicted values using EWG2.5 and a selection of other existing models. A few models (e.g., electrical conductivity, TCLP) were found to be adequate for designing DFHLW glasses in the near future, while others require refits or offsets. It is recommended that new property models be developed for EWG3.0, as a large amount of DFHLW glass property data (> 14 × existing data) is expected to be collected in the compositional spaces where no data was previously available. To enable near-term calculations and formulations for designing DFHLW glasses and processing rate estimations, a formulation algorithm with minor modifications will be developed, EWG2.6.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A versatile high-speed x-ray microscope for sub-10 nm imaging

We have developed a next-generation scanning x-ray microscope RASMI (RApid Scanning Microscopy Instrument) for high-throughput tomographic imaging. RASMI is installed at the hard x-ray nanoprobe beamline at NSLS-II and is capable of manipulating 1D multilayer Laue lenses (MLLs) and 2D optics (both zone plates and monolithically assembled 2D MLLs). The sample scanning stage utilizes line-focusing interferometry as an encoder while performing fly-scanning data acquisition. The system can be configured for both position- and time-triggering modes during fly-scanning. The microscope demonstrated a detector-limited data acquisition rate of 1.25 kHz during ptychography measurements. The initial x-ray results yielded a sample-limited resolution of ∼6 nm in 2D. RASMI can be adopted for in-vacuum applications and is a foundation for the next-generation scanning microscopy systems to be developed and commissioned at NSLS-II.

36 MATERIALS SCIENCE↗

Optical Flow for Intermediate Frame Interpolation of Multispectral Geostationary Satellite Data

Applications in areas such as weather tracking and modeling, ecosystem monitoring, wildfire detection, and land-cover change are heavily dependent on spatial and temporal resolutions of satellite observations. However, there are typically trade-offs between spatial and temporal resolutions in dataset selection. For instance, geostationary weather tracking satellites are designed to take snapshots many times throughout the day but sensor hardware limits data collection. In this work we tackle this limitation, developing a method for temporal upsampling of multi-spectral satellite imagery using optical flow video interpolation deep convolutional neural networks. The presented model, extends Super SloMo (SSM) from single optical flow estimates to multichannel where flows are computed per band. We apply this technique on 8 multi-spectral bands of NOAA/NASA's GOES-16 mesoscale dataset to temporally enhance full disk hemispheric snapshots from 15 minutes to 1 minute. Through extensive experimentation, we show SSM vastly outperforms the linear interpolation baseline and that multichannel optical flows improves performance on GOES-16. A visual analysis of optical flow vectors clearly identifies hurricanes and large-scale atmospheric dynamics. Furthermore, we discuss challenges and open questions related to optical flow and temporal interpolation of multispectral geostationary satellite imagery.

Optical Flow↗

Conceptual Designs for Irradiation Creep Testing of SiC in HFIR

Understanding irradiation creep of nuclear fuel cladding is important to properly size the initial fuel-cladding gap and understand when pellet-cladding contact is expected to occur due to a combination of fuel swelling and cladding creep-down. Irradiation creep also plays a role in relaxing stresses that develop in-pile. Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are the leading long-term accident-tolerant fuel cladding concept for light-water reactors (LWRs). Although some limited data are available regarding irradiation creep of the individual constituents (fibers, matrix), data regarding irradiation creep of SiC/SiC composites are currently insufficient. Additional data regarding irradiation creep compliance and the rupture lifetime (combination of creep and slow crack growth) are needed to understand material limitations. This work describes the design and development of two irradiation vehicles that are being pursued for testing SiC/SiC concepts in the High Flux Isotope Reactor (HFIR). The first is a passive experiment, referred to as the PRECISE experiment, that leverages the constant coolant pressure of HFIR to compress a metallic bellows and provide a well-characterized load to drive creep in a SiC/SiC dog bone specimen. The total creep strain would be quantified post-irradiation by measuring dimensional changes of the specimen length as well as local dimensional changes within the gauge region. Non-stressed specimens would also be irradiated under the same conditions to provide an indication of dimensional changes due to radiation-induced swelling in the absence of creep. A second, more complex experiment, referred to as the INSITE experiment, is being designed in parallel that would use pneumatics to pressurize a metal bellows and linear variable differential transformers (LVDTs) to measure the specimen displacement in situ during irradiation. Such an experiment would provide significantly more data regarding the evolution of the creep compliance as a function of dose and applied stress within a single experiment but would require significantly more development time and cost to execute. The primary concern with the INSITE experiment is the accuracy, reliability, and expected lifetime of the LVDTs during irradiation at elevated temperatures. Efforts are being made to adjust the experiment design and operating procedure to limit LVDT temperatures and mitigate or otherwise compensate for uncertainties due to factors such as temperature fluctuations, creep in the surrounding structural materials, and drift of the LVDTs. This work describes the experiment designs, thermal and structural analysis that were performed to ensure that the desired temperature and stress conditions can be achieved, some initial sensitivity analyses to predict the evolution of the radiation-induced specimen displacements, and potential sources of uncertainty in the measurements. Out-of-pile testing is being performed in parallel to confirm that the test trains achieve the expected stress states in the specimens and do not result in prohibitive stress concentrators (e.g., in the grip regions) that might risk pre-mature failure. The PRECISE experiments are proceeding toward fabrication and assembly with HFIR insertion planned during fiscal year 2026. The INSITE experiment is progressing toward out-of-pile demonstrations, which will provide more conclusive evidence regarding the feasibility of executing these tests in HFIR or whether alternative displacement monitoring techniques may need to be considered.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Dynamic Boundary Microgrids Under Privatization Considerations

Microgrids have physical, electrical, and logical (data, network, and ownership) boundaries. To power unserved customer loads during an outage, microgrids can extend the traditional operational boundaries. This can become complex when considering microgrid-to-microgrid (M2M) interactions where sensitive information such as competitive microgrid operational data is not shared. This work proposes an optimization method coordinated between microgrid controllers and distribution management systems that limits data sharing. The method involves a competitive bidding strategy that maximizes unserved load coverage while minimizing resource utilization and sensitive operational data sharing among entities. The work is validated on a two-microgrid system with photovoltaic and energy storage systems and curves of load derived from real world residential buildings datasets. Results show that the proposed method, when applied for three distinct use cases of energy storage sufficiency to cover the predefined boundary and/or the expanded boundary, can successfully select and bid the available load coverage.

Starke, Michael [ORNL] (ORCID:0000000221211195)↗

Baseline Cost Model for Hydropower: Documentation (2025)

Hydropower currently contributes about 80 GW of conventional and 23 GW of pumped storage capacity to the United States (US) power grid. Previous studies have estimated a considerable amount of remaining US hydropower resources, including non-powered dams (NPD) (Hadjerioua et al., 2012), new stream-reach developments (NSD) (Kao et al., 2014), pumped storage hydropower (PSH) and canal/conduit (Kao et al., 2022). The combined theoretical capacity potential of these various hydropower resources is comparable to the existing US hydropower capacity. There is a continuing interest in developing this hydropower potential, particularly to help meet the increasing demand for electricity. However, available data (Sasthav and Oladosu, 2022) show that the rate of new hydropower development has slowed considerably over time despite the interest of industry stakeholders. This is partly due to the competition from other energy resources and from the highly dispersed nature of remaining hydropower resources, which lead to high information requirements for evaluating the feasibility of potential projects. Cost information provides the most succinct summary of the feasibility of a potential hydropower project required by stakeholders, including developers, investors, policymakers, consumer groups, etc., considering investment options. The best estimates of hydropower costs can be obtained through detailed engineering design and cost assessments of individual projects. However, this approach has high data and resource (time, funds, cross-disciplinary expertise) requirements that render it inapplicable for rapid cost estimation with limited data. Although innovative approaches can overcome some of these impediments (see Oladosu and Ma, 2024 for such an application to potential NPD projects), the development of such approaches still requires significant amounts of resources and are not generally applicable to all hydropower project types. Therefore, statistical and parametric methods using simpler cost specifications remain of significant utility to hydropower stakeholders and are, at the least, complementary to more detailed approaches, particularly when evaluating many potential projects.

13 HYDRO ENERGY↗

Tank Pressure Control Experiment/thermal Phenomena (TPCE/TP)

The 'Tank Pressure Control Experiment/Thermal Phenomena (TPCE/TP)' is a reflight of the tank pressure control experiment (TPCE), flown on STS-43 in a standard Get-Away Special (GAS) container in August 1991. The TPCE obtained extensive video and digital data of the jet induced mixing process in a partially filled tank in low gravity environments. It also provided limited data on the thermal processes involved. The primary objective of the reflight of TPCE is to investigate experimentally the phenomena of liquid superheating and pool nucleate boiling at very low heat fluxes in a long duration low gravity environment. The findings of this experiment will be of direct relevance to space based subcritical cryogenic fluid system design and operation. Experiment hardware and results from the first TPCE are described in outline and graphic form.

Hasan, M. M.↗

Operational Experiences of General Aviation Aircraft

An investigation is currently under way to determine the operational practices and load experiences of general aviation aircraft performing five basic types of operations: twin-engine executive, single-engine executive, personal, instructional, and commercial survey. Limited data obtained to date from aircraft engaged in these operations indicate that aircraft are generally being operated within the limits to which they were designed.

Jewel, Joseph W., Jr.↗

Data Production on Past and Future NASA Missions

Data return is a metric that is commonly publicized for all space science missions. In the early days of the Space Program, this figure was small, and could be described in bits or maybe even megabits. But now, missions are capable of returning data volumes two or three orders of magnitude larger. For example, Voyager 1 and 2 combined produced a little over 5 Terabits of data in 39 years of operation. In contrast, the Cassini mission, launched two decades after Voyager, produced about one and a half times those data volumes in half the time. NISAR, an Earth Science Mission currently in implementation, plans to produce over 28 Petabits of raw data in just 3 years. This means that NISAR will produce about as many data in 30 days as the combined data production of nearly all planetary missions to date. These increases in capability are a result of technology enhancements in two main areas: telecommunications architecture (both space and ground segments) and data storage technology. This paper describes the progression of these two technologies over the course of more than three decades of space missions and provides additional insight into the design of the end-to-end NISAR Data System Architecture. Trends in the data are briefly explored and compared to Moore’s Law which provides only a qualitative model for memory growth but not for data production. In summary, early missions are found to be driven by unrefined processes while later missions, having utilized earlier lessons learned, focus more on improvements to flight and ground capabilities. Data return seems to fall into three categories. First, deep space missions are driven by the large distances that limit data return to the Earth. Next, the orbiter infrastructure around Mars helps these missions generate more data than other deep space spacecraft. Finally, near-Earth missions have the greatest capabilities for the studied metrics due to their close proximity to Earth and the ground network availability.

Xaypraseuth, Peter↗

Integrated Framework of Multisource Data Fusion for Outage Location in Looped Distribution Systems

Accurate outage location is essential for expediting post-outage power restoration, minimizing outage duration, and enhancing the resilience of distribution networks. With the advent of advanced metering infrastructure, data-driven outage location methods have significantly advanced beyond traditional approaches that rely on manual inspections. However, existing methods still face critical challenges, like reliance on single-source data, limited ability to handle partially observable systems or difficulties with loop networks. To the best of our knowledge, no single approach has comprehensively addressed all of these challenges at once. To this end, this paper proposes a comprehensive multisource data fusion framework for outage locations via probabilistic graph networks. The framework consists of three key phases. First, a novel method for reconstituting distribution networks with loops is developed, transforming looped networks into multiple radial subnetworks that retain all outage causalities of the original network. Second, Bayesian network (BN) models are established for each subnetwork, integrating multiple data sources and network structures. Finally, a joint Gibbs sampling mechanism, featuring forward and backward information flow, is designed to merge data from separate BN models and maximize the utilization of limited evidence, ensuring accurate outage location identification. In conclusion, the framework was validated on two modified public test systems, and comparative studies confirmed its effectiveness.

24 POWER TRANSMISSION AND DISTRIBUTION↗

The Design and Implementation of a Secure Datastore Based on Ethereum Smart Contract

In this paper, we present a secure datastore based on an Ethereum smart contract. Our research is guided by three research questions. First, we will explore to what extend a smart-contract-based datastore should resemble a traditional database system. Second, we will investigate how to store the data in a smart-contract-based datastore for maximum flexibility while minimizing the gas consumption. Third, we seek answers regarding whether or not a smart-contract-based datastore should incorporate complex processing such as data encryption and data analytic algorithms. The proposed smart-contract-based datastore aims to strike a good balance between several constraints: (1) smart contracts are publicly visible, which may create a confidentiality concern for the data stored in the datastore; (2) unlike traditional database systems, the Ethereum smart contract programming language (i.e., Solidity) offers very limited data structures for data management; (3) all operations that mutate the blockchain state would incur financial costs and the developers for smart contracts must make sure sufficient gas is provisioned for every smart contract call, and ideally, the gas consumption should be minimized. Our investigation shows that although it is essential for a smart-contract-based datastore to offer some basic data query functionality, it is impractical to offer query flexibility that resembles that of a traditional database system. Furthermore, we propose that data should be structured as tag-value pairs, where the tag serves as a non-unique key that describes the nature of the value. We also conclude that complex processing should not be allowed in the smart contract due to the financial burden and security concerns. The tag-based secure datastore designed this way also defines its applicative perimeter, i.e., only applications that align with our strategy would find the proposed datastore a good fit. Those that would rather incur higher financial cost for more data query flexibility and/or less user burden on data pre- and post-processing would find the proposed database too restrictive.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Jupiter radiation belt engineering model

On the basis of earth observations of the HF and UHF radio emission generated near Jupiter, the presence of energetic charged particles trapped in the planet's dipole magnetic field has been inferred. For electrons, energies of the order of 10 MeV and peak fluxes of the order of 10 to the 7th power per square centimeter per second can be derived from the data for equatorial regions about two planetary radii from the dipole. Energetic protons and lower-energy electrons and protons are also expected, but the limited data require that their fluxes be based on theory or earth analogy. Because descriptions available in the literature suggest large associated uncertainties, both nominal and limiting models for the charged-particle populations of Jupiter's belts are derived. These new engineering models describe electron and proton fluxes and their distributions in energy and position in forms suitable as space vehicle design criteria.

Divine, N.↗

Atmospheric transmission computer program CP

A computer program is described which allows for calculation of the effects of carbon dioxide, water vapor, methane, ozone, carbon monoxide, and nitrous oxide on earth resources remote sensing techniques. A flow chart of the program and operating instructions are provided. Comparisons are made between the atmospheric transmission obtained from laboratory and spacecraft spectrometer data and that obtained from a computer prediction using a model atmosphere and radiosonde data. Limitations of the model atmosphere are discussed. The computer program listings, input card formats, and sample runs for both radiosonde data and laboratory data are included.

Pitts, D. E.↗