Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “frequency estimation”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

fluxfinder: An R Package for Reproducible Calculation and Initial Processing of Greenhouse Gas Fluxes From Static Chamber Measurements

Fluxes of greenhouse gases are a critical component of the earth's natural climate, but anthropogenic emissions have created an imbalance and resulted in global climate change. Quantifying the emission of these gases is vital to our understanding of their sources and sinks, both natural and anthropogenic. The static chamber method, in which a system of interest is enclosed, and gas concentrations are measured over time, is widely used to estimate fluxes of greenhouse gases. With the development of instruments such as infrared gas analyzers (IRGAs) supporting high-frequency concentration data, there is a growing need for open-source workflows to calculate fluxes. Here we present fluxfinder, an R package designed to support reproducible calculations and processing of greenhouse gas fluxes measured with the static chamber method. The package includes raw data file parsing from widely used IRGAs, metadata matching, unit conversion, flux estimations, and initial quality assurance/quality control (QA/QC). Diagnostic graphical plots provide a transparent way to differentiate between measurement issues and nonlinear behavior. The package is also designed to be easily integrated with the gasfluxes package for further fitting of nonlinear concentration-time models, allowing alternative or additional flux QA/QC. The fluxfinder package offers a flexible workflow that is easily adaptable to promote open and reproducible greenhouse gas flux estimations.

Wilson, Stephanie J.↗
NASA NTRS1975

Input design for identification of aircraft stability and control derivatives

An approach for designing inputs to identify stability and control derivatives from flight test data is presented. This approach is based on finding inputs which provide the maximum possible accuracy of derivative estimates. Two techniques of input specification are implemented for this objective - a time domain technique and a frequency domain technique. The time domain technique gives the control input time history and can be used for any allowable duration of test maneuver, including those where data lengths can only be of short duration. The frequency domain technique specifies the input frequency spectrum, and is best applied for tests where extended data lengths, much longer than the time constants of the modes of interest, are possible. These technqiues are used to design inputs to identify parameters in longitudinal and lateral linear models of conventional aircraft. The constraints of aircraft response limits, such as on structural loads, are realized indirectly through a total energy constraint on the input. Tests with simulated data and theoretical predictions show that the new approaches give input signals which can provide more accurate parameter estimates than can conventional inputs of the same total energy. Results obtained indicate that the approach has been brought to the point where it should be used on flight tests for further evaluation.

Gupta, N. K.↗
DOE OSTI2022

Deep Neural Network High Spatiotemporal Resolution Precipitation Estimation (Deep-STEP) Using Passive Microwave and Infrared Data

Recent developments in “headline-making” deep neural networks (DNNs), specifically convolutional neural networks (CNNs), along with advancements in computational power, open great opportunities to integrate massive amounts of real-time observations to characterize spatiotemporal structures of surface precipitation. This study aims to develop a CNN algorithm, named Deep Neural Network High Spatiotemporal Resolution Precipitation Estimation (Deep-STEP), that ingests direct satellite passive microwave (PMW) brightness temperatures (Tbs) at emission and scattering frequencies combined with infrared (IR) Tbs from geostationary satellites and surface information to automatically extract geospatial features related to the precipitable clouds. These features allow the end-to-end Deep-STEP algorithm to instantaneously map surface precipitation intensities with a spatial resolution of 4 km. The main advantages of Deep-STEP, as compared to current state-of-the-art techniques, are 1) it learns and estimates complex precipitation systems directly from raw measurements in near–real time, 2) it uses the automatic spatial neighborhood feature extraction approach, and 3) it fuses coarse-resolution PMW footprints with IR images to reliably retrieve surface precipitation at a high spatial resolution. We anticipate our proposed DNN algorithm to be a starting point for more sophisticated and efficient precipitation retrieval systems in terms of accuracy, fine spatial pattern detection skills, and computational costs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Continuous snow depth, ground interface temperature and shallow soil temperature measurements from 2021-10-1 to 2022-6-14, Seward Peninsula, Alaska

The dataset contains co-located snow depth, ground interface temperature, and shallow soil temperature measured at 98 discrete locations in a watershed located along the Nome-Teller road at mile marker 27 (referred to as T27), and at 53 discrete locations on a hillslope located along the Kougarok road at mile marker 64 (referred to as K64), in Seward Peninsula, Alaska. The dataset aims to understand the local heterogeneity of snow depth, snow temperature, and soil temperature dynamics and their interactions in a discontinuous permafrost region. The dataset is also valuable to train and evaluate machine learning and physical models to predict snow depth or the impact of snow depth on ground surface temperature. At each location, temperatures above and below the ground surface were measured by a pair of vertically deployed distributed temperature profiling probes designed based on Dafflon et al (2022). The probes have high precision temperature sensors spaced at 5 or 10 cm. The mean daily snow depth was estimated by identifying the pair of consecutive sensors with maximum drop of daily temperature high-frequency fluctuations. The ground interface temperature was measured by the sensor located 1-5 cm above the ground surface at 15-minute intervals. The shallow soil temperature was measured by the sensor located 1-5 cm below the ground surface at 15-minute intervals. The dataset includes a description of the probe locations in the "Probe_locations_*.csv" file and the 3 data files (Snow_depths_*.csv, Ground_interface_temperatures_*.csv, Shallow_soil_temperatures_*.csv). * is either T27 or K64, which are the two study sites. In each data file, each column corresponds to a measurement location. Metadata files include data descriptions (_dd.csv) for tabular data. All included files are listed and described in xxxx_flmd.csv. A more detailed description of data processing, along with an updated dataset incorporating multiple seasons and improved snow depth estimation is available at https://doi.org/10.15485/2480365 (Wang et al., 2025a, Wang et al., 2025b).The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2009

90 GHz and 150 GHz Observations of the Orion M42 Region. A Submillimeter to Radio Analysis

We have used the new 90GHz MUSTANG camera on the Robert C. Green Bank Telescope (GBT)to map the bright Huygens region of the star-forming region M42 with a resolution of 9" and a sensitivity of 2.8 mJy/beam. Ninety GHz is an interesting transition frequency, as MUSTANG detects both the free-free emission characteristic of the H II region created by the Trapezium stars, normally seen at lower frequencies, and thermal dust emission from the background OMCI molecular cloud, normally mapped at higher frequencies. We also present similar data from the 150 GHz GISMO camera taken on the IRAM 30 m telescope. This map has 15" resolution. By combining the MUSTANG data with 1.4, 8. and 31 GHz radio data from the VLA and GBT, we derive a new estimate of the emission measure averaged electron temperature of T(sub e) = 11376+/-1050 K by an original method relating free-free emission intensities at optically thin and optically thick frequencies. Combining Infrared Space Observatory-long wavelength spectrometer (ISO-LWS) data with our data, we derive a new estimate of the dust temperature and spectral emissivity index within the 80" ISO-LWS beam toward Orion KL/BN, T(sub d) = 42+/-3 K and Beta(sub d) = 1.3+/-0.1. We show that both T(sub d) and Beta(sub d) decrease when going from the H II region and excited OMCI interface to the denser UV shielded part OMCI (Orion KL/BN, Orion S). With a model consisting of only free-free and thermal dust emission, we are able to fit data taken at frequencies from 1.5 GHz to 854 GHz (350 micrometers).

Dicker, S. R.↗
DOE OSTI2024

Real-Time Inertia Estimation Tool Implementation Based on Probing Signals

As renewable energy penetration increases and the traditional generators retire in power grids, system inertia decreases and exhibits significant daily fluctuations. Furthermore, fast frequency responses (FFRs) provided by the inverter-based resources (IBRs) begin to playa very critical role and bring new challenges to real-time system inertia monitoring due to the difficulties of quantifying its artificial inertia contribution. Thus, necessitating an accurate real-time inertia estimation tool will not only benefit the secure power grid operations, but also provide insights on assessing the artificial inertia contribution from the IBRs. This paper presents a probing-based real-time inertia estimation tool that has been validated through a power-hardware-in-the-Ioop (PHIL) test system using identical hardware battery energy storage system (BESS) and control in an actual power grid. Preliminary results indicate high estimation accuracy of the developed tool and pave the way for the field test and deployment.

inverter-based resources↗
DOE OSTI2022

In-flight polarization angle calibration for LiteBIRD: blind challenge and cosmological implications

We present a demonstration of the in-flight polarization angle calibration for the JAXA/ISAS second strategic large class mission, LiteBIRD, and estimate its impact on the measurement of the tensor-to-scalar ratio parameter, r, using simulated data. We generate a set of simulated sky maps with CMB and polarized foreground emission, and inject instrumental noise and polarization angle offsets to the 22 (partially overlapping) LiteBIRD frequency channels. Our in-flight angle calibration relies on nulling the EB cross correlation of the polarized signal in each channel. This calibration step has been carried out by two independent groups with a blind analysis, allowing an accuracy of the order of a few arc-minutes to be reached on the estimate of the angle offsets. Both the corrected and uncorrected multi-frequency maps are propagated through the foreground cleaning step, with the goal of computing clean CMB maps. We employ two component separation algorithms, the Bayesian-Separation of Components and Residuals Estimate Tool (B-SeCRET), and the Needlet Internal Linear Combination (NILC). We find that the recovered CMB maps obtained with algorithms that do not make any assumptions about the foreground properties, such as NILC, are only mildly affected by the angle miscalibration. However, polarization angle offsets strongly bias results obtained with the parametric fitting method. Once the miscalibration angles are corrected by EB nulling prior to the component separation, both component separation algorithms result in an unbiased estimation of the r parameter. While this work is motivated by the conceptual design study for LiteBIRD, its framework can be broadly applied to any CMB polarization experiment. In particular, the combination of simulation plus blind analysis provides a robust forecast by taking into account not only detector sensitivity but also systematic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2021

The NANOGrav 12.5 Year Data Set: Monitoring Interstellar Scattering Delays

We extract interstellar scintillation parameters for pulsars observed by the NANOGrav radio pulsar timing program. Dynamic spectra for the observing epochs of each pulsar were used to obtain estimates of scintillation timescales, scintillation bandwidths, and the corresponding scattering delays using a stretching algorithm to account for frequency-dependent scaling. We were able to measure scintillation bandwidths for 28 pulsars at1500 MHz and 15 pulsars at 820 MHz. We examine scaling behavior for 17 pulsars and find power-law indices ranging from−0.7 to−3.6, though these may be biased shallow due to insufficient frequency resolution at lower frequencies. We were also able to measure scintillation timescales for six pulsars at 1500 MHz and seven pulsars at820 MHz. There is fair agreement between our scattering delay measurements and electron-density model predictions for most pulsars. We derive interstellar scattering-based transverse velocities assuming isotropic scattering and a scattering screen halfway between the pulsar and Earth. We also estimate the location of the scattering screens assuming proper motion and interstellar scattering-derived transverse velocities are equal. We find no correlations between variations in scattering delay and either variations in dispersion measure or flux density. For most pulsars for which scattering delays are measurable, we find that time-of-arrival uncertainties for a given epoch are larger than our scattering delay measurements, indicating that variable scattering delays are currently subdominant in our overall noise budget but are important for achieving precisions of tens of nanoseconds or less.

Jacob E Turner↗
NASA NTRS1986

Modified LTR robust control for flexible structures

A procedure is developed for dealing with performance and robustness issues in the design of multi-input multi-output compensators for lightly damped flexible structures. The procedure is based upon representing errors in the plant design model as structured uncertainties, and applying a modified version of the Loop Transfer Recovery (LTR) design method. Real parameter errors, such as frequency errors, damping errors or modal displacement errors can be treated. The approach involves adjusting the cost function in the regulator problem and the process noise model in the estimator problem in a particular manner which reflects the assumed structure of the modeling errors. Numerical examples dealing with the control of a large flexible space antenna with uncertain frequencies demonstrate a considerable improvement over standard LTR methods. Convenient design parameters can be varied until a satisfactory compromise is achieved between performance and robustness.

Blelloch, P. A.↗
NASA NTRS2015

Numerical Simulations of the XR-5 Hall Thruster for Life Assessment at Different Operating Conditions

NASA's Jet Propulsion Laboratory has been investigating the applicability of Aerojet Rocketdyne's XR-5 thruster, a 4.5 kW class Hall thruster, for deep-space missions. Major considerations for qualifying the XR-5 for deep-space missions are demonstration of a wide throttling envelope and a usable life capability in excess of 10,000 h. Numerical simulations with the 2-D axisymmetric code Hall2De are employed to inform the qualification process by assessing erosion rates at the thruster surfaces in a wide range of throttling conditions without the need for conducting costly endurance testing. In previous work at JPL by Jorns et al., the anomalous collision frequency distribution for 11 different throttling conditions of the XR-5 spanning 0.3-4.5 kW were identified based on probe measurements of the electron temperature in the near plume region. In this paper, we provide estimates for the erosion rates at the channel walls and pole covers for the same 11 conditions. Uncertainties in the plasma measurements and in the anomalous collision frequency distribution are addressed by determining upper and lower bounds of the erosion rates. Results suggest that erosion of the walls only occurs in the last 5% of the acceleration channel and the rate of such erosion decreases as the geometry of the thruster changes in time due to magnetic shielding. A quasi-zero-erosion state is eventually achieved in all the examined throttling conditions. Examination of the results for pole surface erosion and estimated cathode life indicates that the XR-5 propellant throughput capability will exceed 700 kg, which provides 50% margin over the usable throughput capability of 466 kg as already demonstrated in wear testing.

Lopez Ortega, Alejandro↗
DOE OSTI2023

Asteroids seen by JWST-MIRI: Radiometric size, distance, and orbit constraints

Infrared measurements of asteroids are crucial for the determination of physical and thermal properties of individual objects, and for understanding the small-body populations in the solar system as a whole. However, standard radiometric methods can only be applied if the orbit of an object is known, hence its position at the time of the observation. With JWST-MIRI observations the situation will change and many unknown, often very small, solar system objects will be detected. Later orbit determinations are difficult due to the faintness of the objects and the lack of dedicated follow-up concepts. We present MIRI observations of the outer-belt asteroid (10920) 1998 BC1 and an unknown object, detected in all nine MIRI bands in close apparent proximity to (10920). We developed a new method called STM-ORBIT to interpret the multi-band measurements without knowing the object’s true location. The power of the new technique is that it determines the most-likely heliocentric and observer-centric distance and phase angle ranges, allowing us to make a radiometric size estimate. The application to the MIRI fluxes of (10920) was used to validate the method. It leads to a confirmation of the known radiometric size-albedo solution, and puts constraints on the asteroid’s location and orbit in agreement with its true orbit. To back up the validation of the method, we obtained additional ground-based light curve observations of (10920), combined with Gaia data, which indicate a very elongated object (a/b ≥ 1.5), with a spin-pole at (λ, β) ecl = (178°, +81°), with an estimated error of about 20°, and a rotation period of 4.861191 ± 0.000015 h. A thermophysical study of all available JWST-MIRI and WISE measurements leads to a size of 14.5–16.5 km (diameter of an equal-volume sphere), a geometric albedo pV between 0.05 and 0.10, and a thermal inertia in the range 9–35 (best value 15) J m -2 s -0.5 K -1 . For the newly discovered MIRI object, the STM-ORBIT method revealed a size of 100–230 m. The new asteroid must be on a low-inclination orbit (0.7° < i < 2.0°) and it was located in the inner main-belt region during JWST observations. A beaming parameter η larger than 1.0 would push the size even below 100 meters, a main-belt regime that has escaped IR detections so far. This kind of MIRI observations can therefore contribute to formation and evolution studies via classical size-frequency studies, which are currently limited to objects larger than about one kilometer in size. We estimate that MIRI frames with pointings close to the ecliptic and short integration times of only a few seconds will always include a few asteroids; most of them will be unknown objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Forced Component Estimation Statistical Method Intercomparison Project (ForceSMIP)

Anthropogenic climate change is unfolding rapidly, yet its regional manifestation can be obscured by internal variability. A primary goal of climate science is to identify the externally forced climate response from among the noise of internal variability. Separating the forced response from internal variability can be addressed in climate models by using a large ensemble to average over different possible realizations of internal variability. However, with only one realization of the real world, it is a major challenge to isolate the forced response directly in observations. In the Forced Component Estimation Statistical Method Intercomparison Project (ForceSMIP), contributors used existing and newly developed statistical and machine learning methods to estimate the forced response over 1950–2022 within individual realizations of the climate system. Participants used neural networks, linear inverse models, fingerprinting methods, and low-frequency component analysis, among other approaches. These methods were trained using large ensembles from multiple climate models and then applied to observations. Here, we evaluate method performance within large ensembles and investigate the estimates of the forced response in observations. Our results show that many different types of methods are skillful for estimating the forced response in climate models, though the relative skill of individual methods varies depending on the variable and evaluation metric. Methods with comparable skill in models can give a wide range of estimates of the forced response pattern in observations, illustrating the epistemic uncertainty in forced response estimates. ForceSMIP gives new insights into the forced response in observations, its uncertainty, and methods for its estimation.

Climate attribution↗
NASA NTRS2019

A Frequency Diversity Algorithm for Extending the Radar Doppler Velocity Nyquist Interval

Compact millimeter wavelength radars have been widely used for applications such as remote sensing of clouds, guidance avionics and recently, automotive navigation. However, the short wavelength of these radars limit their maximum unambiguous Doppler velocity. A common solution to this problemis to subsequently unfold the Doppler velocity estimate with the staggered pulse repetition time (PRT) algorithm, which requires two different PRTs to be employed in sequence. This work investigates a novel and potentially more rapid method to extend the Doppler velocity Nyquist interval. The investigated algorithm uses a pair of frequency diverse pulses separated by a short time lag for Doppler velocity estimation while the unambiguous range still corresponds to the pulse repetition time employed.

V Venkatesh↗
NASA NTRS1994

Evaluation of the TOPEX/POSEIDON dual-frequency ionosphere correction

The NASA altimeter on board TOPEX/POSEIDON exploits the difference in the delays of the Ku and C band radar pulses to estimate an ionosphere correction to the range measurement. The dependence of the ionosphere correction on ocean and satellite parameters is less than 1 cm. The standard deviation of the 1-s averaged ionosphere correction depends on the height of the ocean waves and ranges from 5 to 14 mm. The accuracy of the ionosphere correction is better than 1 cm at the 1 sigma confidence level. The ionosphere correction should be averaged over 140 km (20 s) along track in order to minimize its noise without sacrificing its accuracy. Ionosphere models must achieve an independent sample spacing of 900 km or less in order to allow a single-frequency altimeter to have an ionosphere correction comparable in accuracy to that of the NASA dual-frequency altimeter.

Imel, David A.↗
DOE OSTI2026

Data for A Hybrid Biophysical-Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81–0.94) and H (R2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

AI/ML↗
NASA NTRS1995

Measurements of Holographic Properties of Bacteriorhodopsin Films

Several different bacteriorhodopsin (BR) films are characterized with respect to general holographic properties. Experimental measurements include diffraction efficiency and sensitivity as functions of writing intensity and grating frequency, hologram thermal decay behavior, diffraction efficiency as a function of grating tilt within the film and modulation depth, and estimates of the refractive index change from the diffraction efficiency data. The films studied include those made from wildtype BR and the genetic variants D96N and D96N/T46V. The maximum diffraction efficiency measured was just over 3% for a hydrated wildtype BR film. The films' holographic properties were found to be relatively insensitive to grating frequency and grating tilt angle. The diffraction efficiency dropped off more sharply as a function of modulation depth than a purely linear medium, and only the hydrated wildtype film exhibited significant behavior variation with different writing intensities.

Downie, John D.↗
DOE OSTI2024

The Atacama Cosmology Telescope: Mitigating the Impact of Extragalactic Foregrounds for the DR6 Cosmic Microwave Background Lensing Analysis

We investigate the impact and mitigation of extragalactic foregrounds for the cosmic microwave background (CMB) lensing power spectrum analysis of Atacama Cosmology Telescope (ACT) data release 6 (DR6) data. Two independent microwave sky simulations are used to test a range of mitigation strategies. We demonstrate that finding and then subtracting point sources, finding and then subtracting models of clusters, and using a profile bias-hardened lensing estimator together reduce the fractional biases to well below statistical uncertainties, with the inferred lensing amplitude, A lens , biased by less than 0.2σ. We also show that another method where a model for the cosmic infrared background (CIB) contribution is deprojected and high-frequency data from Planck is included has similar performance. Other frequency-cleaned options do not perform as well, either incurring a large noise cost or resulting in biased recovery of the lensing spectrum. In addition to these simulation-based tests, we also present null tests on the ACT DR6 data for sensitivity of our lensing spectrum estimation to differences in foreground levels between the two ACT frequencies used, while nulling the CMB lensing signal. These tests pass whether the nulling is performed at the map or bandpower level. The CIB-deprojected measurement performed on the DR6 data is consistent with our baseline measurement, implying that contamination from the CIB is unlikely to significantly bias the DR6 lensing spectrum. This collection of tests gives confidence that the ACT DR6 lensing measurements and cosmological constraints presented in companion papers to this work are robust to extragalactic foregrounds.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1987

Comparison of measured and modeled radiation, heat and water vapor fluxes: FIFE pilot study

The feasibility of using radio frequency receivers to collect data from automated weather stations to model fluxes of latent heat, sensible heat, and radiation using routine weather data collected by automated weather stations was tested and the estimated fluxes were compared with fluxes measured over wheat. The model Cupid was used to model the fluxes. Two or more automated weather stations, interrogated by radio frequency and other means, were utilized to examine some of the climatic variability of the First ISLSCP (International Satellite Land-Surface Climatology Project) Field Experiment (FIFE) site, to measure and model reflected and emitted radiation streams from various locations at the site and to compare modeled latent and sensible heat fluxes with measured values. Some bidirectional reflected and emitted radiation data were collected from 23 locations throughout the FIFE site. Analysis of these data along with analysis of the measured sensible and latent heat fluxes is just beginning.

Blad, Blaine L.↗