Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral statistics”

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,135 records · Page 63

Correlation between aircraft MSS and LIDAR remotely sensed data on a forested wetland in South Carolina

Wetlands in a portion of the Savannah River swamp forest, the Steel Creek Delta, were mapped using April 26, 1985 high-resolution aircraft multispectral scanner (MSS) data. Due to the complex spectral characteristics of the wetland vegetation, it was necessary to implement several techniques in the classification of the MSS imagery of the Steel Creek Delta. In particular, when performing unsupervised classification, an iterative cluster busting technique was used which simplified the cluster labeling process. In addition to the MSS data, light detecting and ranging (LIDAR) data were acquired by National Aeronautics and Space Administration (NASA) personnel along two flightlines over the Steel Creek Delta. These data were registered with the wetland classification map and correlated. Statistical analyses demonstrated that the laser derived canopy height information was significantly correlated with the Steel Creek Delta wetland classes encountered along the profiling transect of the LIDAR data.

Jensen, John R.↗

Self-diagnosis of model suitability for continuous measurements of stream-dissolved organic carbon derived from in situ UV–visible spectroscopy

Application of high-frequency monitoring of dissolved organic carbon (DOC) is difficult in instances where training datasets are challenging to develop (e.g., remote locations) and the relationship between optical features and DOC concentration changes due to environmental or landscape shifts (e.g., climate or land-use change). We developed and compared three partial least squares (PLS) models using in situ water level measurements, conductivity, and UV–Vis spectral attenuation to predict DOC. Two site-specific models were developed using data from a hillslope-dominated forest or a low-relief wetland-pond-dominated stream catchment. The third model, using data from both sites, exhibited the best performance (DOC range = 4–15.5 mg C L −1 , mean = 8.38 mg C L −1 , training RMSE = 0.34 mg C L −1 , internal validation RMSE = 0.50 mg C L −1 , external validation RMSE = 2.43 mg C L −1 ). We further demonstrate using PLS model statistics to monitor performance and elucidate when and how models should be updated. These statistics, Hotelling's T 2 and squared prediction errors, are useful consistency checks for the predictions made and detect underlying inconsistencies that, if undetected, can reduce the robustness of DOC prediction. For example, via the T 2 statistic, we identified the summer–autumn transition as a period when DOC composition differed from what was represented in the training dataset. We also determined that elevated SUVA 254 values contributed to the overall bias observed in predictions made during the subsequent year as part of the external validation. This enabled the application of a bias correction that reduced the RMSE from 2.43 to 0.89 mg C L −1 . The method presented here could be applied to future monitoring programs enabling model updates to monitor DOC fluxes accurately from optical datasets (e.g., attenuance or fluorescence) in the face of developing datasets in remote locations or environmental change. In conclusion, implementation of this approach may also identify possible regime shifts or landscape and hydrologic change associated with climate and other environmental changes relevant to terrestrial to aquatic fluxes.

UV-VIS spectroscopy↗

Fractal applications to complex crustal problems

Complex scale-invariant problems obey fractal statistics. The basic definition of a fractal distribution is that the number of objects with a characteristic linear dimension greater than r satisfies the relation N = about r exp -D where D is the fractal dimension. Fragmentation often satisfies this relation. The distribution of earthquakes satisfies this relation. The classic relationship between the length of a rocky coast line and the step length can be derived from this relation. Power law relations for spectra can also be related to fractal dimensions. Topography and gravity are examples. Spectral techniques can be used to obtain maps of fractal dimension and roughness amplitude. These provide a quantitative measure of texture analysis. It is argued that the distribution of stress and strength in a complex crustal region, such as the Alps, is fractal. Based on this assumption, the observed frequency-magnitude relation for the seismicity in the region can be derived.

Turcotte, Donald L.↗

Optical and Gamma Ray Space Observations

Results of the first year of data acquired from several Earth observatories concerning the optical counterparts gamma ray bursts (GRB) are presented. From the present statistics, it seems to be obvious that typical GRB's have optical emission at the time of the burst at a level at least below 1/(F(sub gamma)/F(sub opt)) approximately equal to 1/0.5 and optical emission a few hours after the burst is lower by a factor of 10 to 200 than the simultaneous emission. Given the fact that GRB spectra are rather broad over the observed energy range of say 20 keV up to 100 MeV, the observations indicate that the broad spectral shape may not continue into the optical range. After the confirmation of the isotropic distribution of GRB's by the BATSE experiment the interpretation now tends to put the sources at cosmological distances.

Source record↗

Spectral Proper Orthogonal Decomposition of uPSP Measurements in Recent NASA Ames Wind Tunnel Test

This paper discusses Spectral Proper Orthogonal Decomposition (SPOD) of the Unsteady Pressure-Sensitive Paint (uPSP) measurements in recent NASA Ames wind tunnel test. The uPSP measurements were collected using Innovative Scientific Solutions, Inc. (ISSI) porous, fast-response pressure-sensitive paint, 40 ISSI four-inch air-cooled Light-Emitting Diodes, and 8 Phantom v2512 high-speed cameras at 10,000 frames per second in the uPSP Launch Vehicle Demonstration Test (LVDT) of the Space Launch System (SLS) vehicle in the 11-by 11-foot transonic test section of the Unitary Plan Wind Tunnel at NASA Ames Research Center in April 2024. SPOD is derived from a space-time proper orthogonal decomposition problem for statistically stationary flows. SPOD modes are determined in the frequency domain. Each SPOD mode oscillates at a single frequency. SPOD can be viewed as an extension of the Discrete Fourier Transform composition and the Dynamic Mode Decomposition. In this paper, the outputs of SPOD of the uPSP measurements in the uPSP LVDT are presented and the effectiveness of SPOD in the identification, diagnosis and analysis of the aerodynamic and aeroacoustic phenomena is demonstrated. The unsteady and dynamic property of the pressure field on the surface of the SLS Block 1B crew vehicle is presented with the visualization of the SPOD modes of the uPSP measurements in the tests of a Mach sweep run of the uPSP LVDT. The SPOD outputs were generated with the execution in parallel of a code in Python, with the library of Message Passing Interface for parallel processing, on the NASA Pleiades supercomputer. The work described in this paper is a part of NASA’s development of a new state-of-the-art uPSP capability in production wind tunnels. Funding was provided by the NASA Aerosciences Evaluation and Test Capabilities Portfolio Office.

Aeroacoustics↗

Atmospheric Correction of Satellite Ocean-Color Imagery During the PACE Era

The Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission will carry into space the Ocean Color Instrument (OCI), a spectrometer measuring at 5 nm spectral resolution in the ultraviolet (UV) to near infrared (NIR) with additional spectral bands in the shortwave infrared (SWIR), and two multi-angle polarimeters that will overlap the OCI spectral range and spatial coverage, i. e., the Spectrometer for Planetary Exploration (SPEXone) and the Hyper-Angular Rainbow Polarimeter (HARP2). These instruments, especially when used in synergy, have great potential for improving estimates of water reflectance in the post Earth Observing System (EOS) era. Extending the top-of-atmosphere (TOA) observations to the UV, where aerosol absorption is effective, adding spectral bands in the SWIR, where even the most turbid waters are black and sensitivity to the aerosol coarse mode is higher than at shorter wavelengths, and measuring in the oxygen A-band to estimate aerosol altitude will enable greater accuracy in atmospheric correction for ocean color science. The multi-angular and polarized measurements, sensitive to aerosol properties (e.g., size distribution, index of refraction), can further help to identify or constrain the aerosol model, or to retrieve directly water reflectance. Algorithms that exploit the new capabilities are presented, and their ability to improve accuracy is discussed. They embrace a modern, adapted heritage two-step algorithm and alternative schemes (deterministic, statistical) that aim at inverting the TOA signal in a single step. These schemes, by the nature of their construction, their robustness, their generalization properties, and their ability to associate uncertainties, are expected to become the new standard in the future. A strategy for atmospheric correction is presented that ensures continuity and consistency with past and present ocean-color missions while enabling full exploitation of the new dimensions and possibilities. Despite the major improvements anticipated with the PACE instruments, gaps/issues remain to be filled/tackled. They include dealing properly with whitecaps, taking into account Earth-curvature effects, correcting for adjacency effects, accounting for the coupling between scattering and absorption, modeling accurately water reflectance, and acquiring a sufficiently representative dataset of water reflectance in the UV to SWIR. Dedicated efforts, experimental and theoretical, are in order to gather the necessary information and rectify inadequacies. Ideas and solutions are put forward to address the unresolved issues. Thanks to its design and characteristics, the PACE mission will mark the beginning of a new era of unprecedented accuracy in ocean-color radiometry from space.

ocean color↗

Stacked reverberation mapping of high-redshift quasars in DESI. I. Feasibility analysis

The broad-line region of quasars has long been probed by reverberation mapping techniques that measure time lags between continuum and broad emission-line variations. Stacked reverberation mapping has been proposed as a less observationally expensive alternative to traditional methods. This ensemble approach also reduces biases from small-number statistics. The Dark Energy Spectroscopic Instrument (DESI) is conducting the most extensive spectroscopic survey of quasars to date. We create mock light curves emulating expected DESI quasar observations at redshifts $1.48\lt z\lt 5.2$ and luminosities $44.68 \le \log \lambda L_{1350 \mathring{\rm A}{}} / \mathrm{erg\, s^{-1}} \le 45.99$ to test stacked reverberation mapping feasibility using sparse spectroscopic data paired with well-sampled photometric data. The pipeline, using the lag estimation code JAVELIN (Just Another Vehicle for Estimating Lags In Nuclei), successfully recovers the simulated C IV lags within 1σ of the true values using spectroscopic light curves composed of only a few spectral epochs (2–10) with irregular cadences. We investigate how observational factors, including C IV flux error magnitude, number of stacked quasars, and spectral epoch count, affect performance. This work motivates a pathway for future stacked reverberation mapping projects with large-scale spectroscopic surveys of quasars having $\ge 2$ spectroscopic observations. Our results suggest an economical alternative for constraining and extending the radius–luminosity relation to higher redshifts and luminosities. Subsequently, this relation can be employed more reliably in single-epoch black hole mass measurements and quasar cosmology in these distant regimes.

quasars: general, quasars: supermassive black hole↗

Numerical wave modeling for operational and survival analyses of wave energy converters at the US Navy Wave Energy Test Site in Hawaii

The US Navy Wave Energy Test Site (WETS) is currently operating three grid-connected berths off Marine Corps Base Hawaii in support of technological development through scaled-model testing and pre-commercial prototyping. The location on the east shore of Oahu is open to year-round wind waves and seasonal north swells as well as severe seas from subtropical systems and hurricanes. We have assembled a spectral wave model system comprising WAVEWATCH III and SWAN on a hierarchy of global, regional, and nearshore computational grids. With wind forcing from global forecast and rean-alysis datasets as well as their regional downscaling, the system produces operational 7-day wave forecasts and a long-term hindcast. The daily forecasts, validated with real-time buoy measurements, facilitate safe deployment, operation, and retrieval of wave energy converters (WECs). The wave hindcast from 1979 to 2017 defines the intra- and inter-annual variations for statistical analysis of sea states and energy resources. Due to limited occurrences of hurricanes in Hawaii waters, their wind forcing is determined parametrically with simulated tracks from climate model downscaling. Finally, systematic analysis of the hindcast and simulation datasets provides significant wave heights from trade winds, swells, subtropical systems, and hurricanes with return periods up to 100 years for WECs survival analysis.

16 TIDAL AND WAVE POWER↗

Study of recreational land and open space using Skylab imagery

The author has identified the following significant results. An analysis of the statistical uniqueness of each of the signatures of the Gratiot-Saginaw State Game Area was made by computing a matrix of probabilities of misclassification for all possible signature pairs. Within each data set, the 35 signatures were then aggregated into a smaller set of composite signatures by combining groups of signatures having high probabilities of misclassification. Computer separation of forest denisty classes was poor with multispectral scanner data collected on 5 August 1973. Signatures from the scanner data were further analyzed to determine the ranking of spectral channels for computer separation of the scene classes. Probabilities of misclassification were computed for composite signatures using four separate combinations of data source and channel selection.

Sattinger, I. J.↗

Regional Characteristics of the Relationship Between Columnar AOD and Surface PM2.5: Application of Lidar Aerosol Extinction Profiles over Baltimore-Washington Corridor During DISCOVER-AQ

The first field campaign of DISCOVER-AQ (Deriving Information on Surface conditions from COlumn and VERtically resolved observations relevant to Air Quality) took place in July 2011 over Baltimore-Washington Corridor (BWC). A suite of airborne remote sensing and in-situ sensors was deployed along with ground networks for mapping vertical and horizontal distribution of aerosols. Previous researches were based on a single lidar station because of the lack of regional coverage. This study uses the unique airborne HSRL (High Spectral Resolution Lidar) data to baseline PM2.5 (particulate matter of aerodynamic diameter less than 2.5 μm) estimates and applies to regional air quality with satellite AOD (Aerosol Optical Depth) retrievals over BWC (∼6500 sq. km). The linear approximation takes into account aerosols aloft above AML (Aerosol Mixing Layer) by normalizing AOD with haze layer height (i.e., AOD/HLH). The estimated PM2.5 mass concentrations by HSRL AOD/HLH are shown within 2 RMSE (Root Mean Square Error ∼9.6 μg/cu. m) with correlation ∼0.88 with the observed over BWC. Similar statistics are shown when applying HLH data from a single location over the distance of 100 km. In other words, a single lidar is feasible to cover the range of 100 km with expected uncertainties. The employment of MPLNET-AERONET (MicroPulse Lidar NETwork - AErosol RObotic NETwork) measurements at NASA GSFC produces similar statistics of PM2.5 estimates as those derived by HSRL. The synergy of active and passive remote sensing aerosol measurements provides the foundation for satellite application of air quality on a daily basis. For the optimal range of 10 km, the MODIS-estimated PM2.5 values are found satisfactory at 27 (out of 36) sunphotometer locations with mean RMSE of 1.6-3.3 μg/cu. m relative to PM2.5 estimated by sunphotometers. The remaining 6 of 8 marginal sites are found in the coastal zone, for which associated large RMSE values ∼4.5-7.8 μg/cu. m are most likely due to overestimated AOD because of water-contaminated pixels.

DISCOVER-AQ↗

TPSAS-NF1676L-13037-DND

Mineral dust has a significant and uncertain role in the direct aerosol radiative forcing of climate. Spaceborne lidars such as CALIOP help reduce these uncertainties through vertical profile measurements of aerosol optical properties. One current limitation to the accurate retrieval of aerosol extinction and optical depth from CALIOP is the assumed relationship between the aerosol extinction to aerosol backscatter (i.e. the extinction-to-backscatter ratio, also referred to here as the lidar ratio or Sa). This problem is especially acute at 1064 nm, where few estimates of the lidar ratio exist. This study uses a dataset of eight underflights of CALIOP during August 2010 by the NASA Langley Research Center airborne High Spectral Resolution Lidar (HSRL) to study Saharan dust transported across the Atlantic Ocean. The standard HSRL profile products include aerosol backscatter coefficients and depolarization ratios at both 532 nm and 1064 nm, and aerosol extinction coefficients (and therefore also lidar ratios) at 532 nm only. In this work, we further derive estimates of aerosol lidar ratios and extinction coefficients at 1064 nm via application of a two-wavelength technique that uses the 532 nm aerosol backscatter coefficients and the 1064 nm attenuated total backscatter profile. Summary statistics of the dust lidar ratio and depolarization at 532 nm and 1064 nm from these eight flights are presented. Implications for the CALIOP dust and polluted dust aerosol types and lidar ratio selection are discussed. In addition to the two-wavelength retrievals of lidar ratio at 1064 nm, a demonstration case of a 1064 nm lidar ratio retrieval over the ocean from CALIOP using the CloudSat measurement of surface scattering cross section as a constraint is presented (Josset et al, 2010).

Raymond R Rogers↗

NuSTAR Observation of GROJ1744−28 at Low Mass Accretion Rate

Context.Neutron stars in low-mass X-ray binaries (LMXBs) are important systems to study the physics of accretion onto compact ob-jects. The system GRO J1744−28 is particularly interesting, as it usually shows clear pulsations as well as X-ray bursts. Additionally,there are claims for a magnetic field of 5×1011G through the detection of a Cyclotron Resonant Scattering Feature (CRSF).Aims.We present the spectral analysis of GRO J1744−28 using∼29 ks of NuSTAR data taken in 2017 February at a low luminosityof 3.2×1036erg s−1(3–50 keV). Our goal is to study the variability of the source spectrum with pulse phase and to search for the claimed CRSF. Methods.The continuum spectrum is modeled with an absorbed power-law with exponential cut-off, and an additional iron line component. We find no obvious indications for a CRSF and therefore perform a detailed cyclotron line search using statistical methods.We perform this search on pulse phase-averaged as well as phase-resolved spectra. Results.GRO J1744−28 was observed at one of the lowest luminosities studied so far. The previously detected Type II X-ray bursts re absent. Clear pulsations at a period of 2.141124(9) Hz are detected. The pulse profile shows an indication of a secondary peak,which was not seen at higher flux. The 4σupper limit for the strength of a CRSF in the 3–20 keV band is 0.07 keV, lower than the strength of the line found at higher luminosity. Conclusions.The detection of pulsations shows that the source did not enter the “propeller” regime, even though the source flux of4.15×10−10erg cm−2s−1was almost one order of magnitude below the threshold for the propeller regime claimed in previous studies on this source. The transition into the propeller regime in GRO J1744−28 must therefore be below a luminosity of 3.2×1036erg s−1(3–50 keV), which implies a surface magnetic field.2.9×1011G and mass accretion rate.1.7×1016g s−1. A change of the CRSF depth as function of luminosity is not unexpected and has been observed in other sources. This result possibly implies a change in emission geometry as function of mass accretion rate to reduce the depth of the line below our detection limit.

Ole König↗

Probing and Tuning Strain‐Localized Exciton Emission in 2D Material Bubbles at Room Temperature

In monolayer transition metal dichalcogenides bubbles-nanoscale deformations typically exhibiting a dome-like shape-Excitons are confined by the strain effect, which exhibits extraordinary emission properties, such as single photon generation, enhanced light emission, and spectrally tunable excitonic states. While the strain profiles of these bubbles are extensively studied, this work provides an approach 1) to directly visualize the associated exciton properties in bubbles formed in WSe2 monolayer, revealing an intrinsic emission wavelength shift of ≈40 nm, and 2) actively modify local strain, enabling further exciton emission tuning over a range of 50 nm. These are achieved by emission mapping and nanoindentation using a dielectric near-field probe, which enables the detection of local emission spectra and emission lifetimes within individual bubbles. Statistical analysis of 67 bubbles uncovers an emission wavelength distribution centered around 780 nm. Furthermore, saturation behavior in the power-dependent studies and the associated lifetime change reveal the localized nature of the strain-induced states. These findings provide direct insights into the strain-localized emission dynamics in bubbles and establish a robust framework for non-destructive, reversible, and predictable nanoscale emission control, presenting a potential avenue for developing next-generation tunable quantum optical sources.

2D materials bubbles↗

Novel application of neutrinos to evaluate U.S. nuclear weapons performance

There is a growing realization that neutrinos can be used as a diagnostic tool to better understand the inner workings of a nuclear weapon. Robust estimates demonstrate that an Inverse Beta Decay (IBD) neutrino scintillation detector built at the Nevada Test Site with a 1000-ton active target mass at a standoff distance of 500 m would detect thousands of antineutrino events per nuclear test. This would provide less than 4% statistical error on the measured antineutrino rate and 5% error on antineutrino energy. Extrapolating this to an error on the test device explosive yield requires knowledge from evaluated nuclear databases, non-equilibrium fission rates, and assumptions on internal neutron fluxes. Initial calculations demonstrate that the total number of neutrinos emitted per fission in the first 10 3 s after a short pulse of 239 Pu fission is about a factor of two less than that from Pu fissioning under steady state conditions. Furthermore, there are significant energy spectral differences as a function of time after the pulse that must be considered. These and other model dependencies will be discussed in the paper. In the absence of nuclear weapons testing, many of the technical and theoretical challenges of a full nuclear test could be mitigated with a low cost smaller scale 20 ton fiducial mass IBD demonstration detector placed near a pulsed reactor. Potential reactors include the Texas A&M University TRIGA 1 GW–10 ms pulsed facility or the Sandia Annular Core Research Reactor. The short duty cycle and repeatability of pulses would provide critical real environment testing and measurements, which would be valuable for planning a possible real test shot in the future. Furthermore, the antineutrino rate as a function of time data would provide unique constraints on fission databases and model assumptions. Finally, there are impactful science drivers such as sensitive searches for ∼1 eV 2 sterile neutrinos and ∼MeV scale axions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

H i asymmetries in LVHIS, VIVA, and HALOGAS galaxies

ABSTRACT We present an analysis of morphological, kinematic, and spectral asymmetries in observations of atomic neutral hydrogen (H i) gas from the Local Volume H i Survey (LVHIS), the VLA Imaging of Virgo in Atomic Gas (VIVA) survey, and the Hydrogen Accretion in Local Galaxies Survey. With the aim of investigating the impact of the local environment density and stellar mass on the measured H i asymmetries in future large H i surveys, we provide recommendations for the most meaningful measures of asymmetry for use in future analysis. After controlling for stellar mass, we find signs of statistically significant trends of increasing asymmetries with local density. The most significant trend we measure is for the normalized flipped spectrum residual (Aspec), with mean LVHIS and VIVA values of 0.204 ± 0.011 and 0.615 ± 0.068 at average weighted 10th nearest-neighbour galaxy number densities of log (ρ10/Mpc−3) = −1.64 and 0.88, respectively. Looking ahead to the Widefield ASKAP L-band Legacy All-sky Blind survey on the Australian Square Kilometre Array Pathfinder, we estimate that the number of detections will be sufficient to provide coverage over 5 orders of magnitude in both local density and stellar mass increasing the dynamic range and accuracy with which we can probe the effect of these properties on the asymmetry in the distribution of atomic gas in galaxies.

Reynolds, T. N.↗

Gamma-Ray Burst Lines

The evidence for spectral features in gamma-ray bursts is summarized. As a guide for evaluating the evidence, the properties of gamma-ray detectors and the methods of analyzing gamma-ray spectra are reviewed. In the 1980's, observations indicated that absorption features below 100 kev were present in a large fraction of bright gamma-ray bursts. There were also reports of emission features around 400 kev. During the 1990's the situation has become much less clear. A small fraction of bursts observed with BATSE have statistically significant low-energy features, but the reality of the features is suspect because in several cases the data of the BATSE detectors appear to be inconsistent. Furthermore, most of the possible features appear in emission rather than the expected absorption. Analysis of data from other instruments has either not been finalized or has not detected lines.

Briggs, Michael S.↗

First Ground-Based Infrared Solar Absorption Measurements of Free Tropospheric Methanol (CH3OH): Multidecade Infrared Time Series from Kitt Peak (31.9 deg N 111.6 deg W): Trend, Seasonal Cycle, and Comparison with Previous Measurements

Atmospheric CH3OH (methanol) free tropospheric (2.09-14-km altitude) time series spanning 22 years has been analyzed on the basis of high-spectral resolution infrared solar absorption spectra of the strong n8 band recorded from the U.S. National Solar Observatory on Kitt Peak (latitude 31.9degN, 111.6degW, 2.09-km altitude) with a 1-m Fourier transform spectrometer (FTS). The measurements span October 1981 to December 2003 and are the first long time series of CH3OH measurements obtained from the ground. The results were analyzed with SFIT2 version 3.93 and show a factor of three variations with season, a maximum at the beginning of July, a winter minimum, and no statistically significant long-term trend over the measurement time span.

Rinsland, Curtis P.↗

First Ground-Based Infrared Solar Absorption Measurements of Free Tropospheric Methanol (CH3OH): Multidecade Infrared Time Series from Kitt Peak (31.9 deg N 111.6 deg W): Trend, Seasonal Cycle, and Comparison with Previous Measurements

Atmospheric CH3OH (methanol) free tropospheric (2.09-14-km altitude) time series spanning 22 years has been analyzed on the basis of high-spectral resolution infrared solar absorption spectra of the strong vs band recorded from the U.S. National Solar Observatory on Kitt Peak (latitude 31.9degN, 111.6degW, 2.09-km altitude) with a 1-m Fourier transform spectrometer (FTS). The measurements span October 1981 to December 2003 and are the first long time series of CH3OH measurements obtained from the ground. The results were analyzed with SFIT2 version 3.93 and show a factor of three variations with season, a maximum at the beginning of July, a winter minimum, and no statistically significant long-term trend over the measurement time span.

Rinsland, Curtis P.↗