Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “applied 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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Advanced techniques and technology for efficient data storage, access, and transfer

Advanced techniques for efficiently representing most forms of data are being implemented in practical hardware and software form through the joint efforts of three NASA centers. These techniques adapt to local statistical variations to continually provide near optimum code efficiency when representing data without error. Demonstrated in several earlier space applications, these techniques are the basis of initial NASA data compression standards specifications. Since the techniques clearly apply to most NASA science data, NASA invested in the development of both hardware and software implementations for general use. This investment includes high-speed single-chip very large scale integration (VLSI) coding and decoding modules as well as machine-transferrable software routines. The hardware chips were tested in the laboratory at data rates as high as 700 Mbits/s. A coding module's definition includes a predictive preprocessing stage and a powerful adaptive coding stage. The function of the preprocessor is to optimally process incoming data into a standard form data source that the second stage can handle.The built-in preprocessor of the VLSI coder chips is ideal for high-speed sampled data applications such as imaging and high-quality audio, but additionally, the second stage adaptive coder can be used separately with any source that can be externally preprocessed into the 'standard form'. This generic functionality assures that the applicability of these techniques and their recent high-speed implementations should be equally broad outside of NASA.

Rice, Robert F.↗
NASA NTRS1991

Statistical evaluation of thermal advection and stratification effects in scatterometer observations

The effects of thermal advection and atmospheric stratification are statistically evaluated using Seasat scatterometer observations as a data base. The results indicate that, whenever the surface winds or wind stress are related to the atmospheric pressure field, the appropriate stratification and baroclinic corrections should be applied. Without such corrections, errors of 15-25 percent are likely to arise in the surface fluxes computed from model low-level winds or pressure measurements.

Levy, Gad↗
NASA NTRS1994

A comparison of optical rain gauge and radar data from TOGA/COARE

A comparison between rain gauge data and radar data from the Tropical Oceans Global Atmosphere/Coupled Ocean Atmosphere Response Experiment (TOGA/COARE) was studied. The rain gauge data from an echo that passed over the Xiangyanghong #5 on December 24, 1992 was compared to what the MIT radar saw at this location from the R/V Vickers, 103.4 km to the east. The precipitation measured by the rain gauge peaked at 108 mm/hr 92 seconds into the period before tapering off 11 1/2 minutes later. This sharp gradient was evident in a PPI plot of the radar reflectivities and the percentage area-rainfall for the radar data statistics. The percentage area curve was converted to rain rates using a GATE Z-R and compared to a percentage time curve of rain rates according to the rain gauge. A four minute running average applied to the rain gauge rates improved the comparison of peak rates between the rain gauge and radar. Differences in peaks between rain rates observed by the rain gauge and reflectivities observed by the radar could be due to variations in rainfall rates within a single radar data bin. For example, two measurements of reflectivity such as 37 and 47 dBZ within the same bin would result in a 44 dBZ average. This range in rates from 12 mm/hr to 74 mm/hr is observed in 30 seconds by the rain gauges within the first two minutes of the radar echo passage.

Galusha, Linda↗
NASA NTRS1999

Quasi-Static Probabilistic Structural Analyses Process and Criteria

Current deterministic structural methods are easily applied to substructures and components, and analysts have built great design insights and confidence in them over the years. However, deterministic methods cannot support systems risk analyses, and it was recently reported that deterministic treatment of statistical data is inconsistent with error propagation laws that can result in unevenly conservative structural predictions. Assuming non-nal distributions and using statistical data formats throughout prevailing stress deterministic processes lead to a safety factor in statistical format, which integrated into the safety index, provides a safety factor and first order reliability relationship. The embedded safety factor in the safety index expression allows a historically based risk to be determined and verified over a variety of quasi-static metallic substructures consistent with the traditional safety factor methods and NASA Std. 5001 criteria.

Goldberg, B.↗
NASA NTRS2006

Proteomic Retrieval from Nucleic Acid Depleted Space-Flown Human Cells

Compared to experiments utilizing humans in microgravity, cell-based approaches to questions about subsystems of the human system afford multiple advantages, such as crew safety and the ability to achieve statistical significance. To maximize the science return from flight samples, an optimized method was developed to recover protein from samples depleted of nucleic acid. This technique allows multiple analyses on a single cellular sample and when applied to future cellular investigations could accelerate solutions to significant biomedical barriers to human space exploration. Cell cultures grown in American Fluoroseal bags were treated with an RNA stabilizing agent (RNAlater - Ambion), which enabled both RNA and immunoreactive protein analyses. RNA was purified using an RNAqueous(registered TradeMark) kit (Ambion) and the remaining RNA free supernatant was precipitated with 5% trichloroacetic acid. The precipitate was dissolved in SDS running buffer and tested for protein content using a bicinchoninic acid assay (1) (Sigma). Equal loads of protein were placed on SDS-PAGE gels and either stained with CyproOrange (Amersham) or transferred using Western Blotting techniques (2,3,4). Protein recovered from RNAlater-treated cells and stained with protein stain, was measured using Imagequant volume measurements for rectangles of equal size. BSA treated in this way gave quantitative data over the protein range used (Fig 1). Human renal cortical epithelial (HRCE) cells (5,6,7) grown onboard the International Space Station (ISS) during Increment 3 and in ground control cultures exhibited similar immunoreactivity profiles for antibodies to the Vitamin D receptor (VDR) (Fig 2), the beta isoform of protein kinase C (PKC ) (Fig 3), and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Fig 4). Parallel immunohistochemical studies on formalin-fixed flight and ground control cultures also showed positive immunostaining for VDR and other biomarkers (Fig 5). These results are consistent with data from additional antigenic recovery experiments performed on human Mullerian tumor cells cultured in microgravity (8).

Hammond, D. K.↗
NASA NTRS2008

A New Concurrent Multiscale Methodology for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Yamakov, Vesselin↗
NASA NTRS2008

An Embedded Statistical Method for Coupling Molecular Dynamics and Finite Element Analyses

The coupling of molecular dynamics (MD) simulations with finite element methods (FEM) yields computationally efficient models that link fundamental material processes at the atomistic level with continuum field responses at higher length scales. The theoretical challenge involves developing a seamless connection along an interface between two inherently different simulation frameworks. Various specialized methods have been developed to solve particular classes of problems. Many of these methods link the kinematics of individual MD atoms with FEM nodes at their common interface, necessarily requiring that the finite element mesh be refined to atomic resolution. Some of these coupling approaches also require simulations to be carried out at 0 K and restrict modeling to two-dimensional material domains due to difficulties in simulating full three-dimensional material processes. In the present work, a new approach to MD-FEM coupling is developed based on a restatement of the standard boundary value problem used to define a coupled domain. The method replaces a direct linkage of individual MD atoms and finite element (FE) nodes with a statistical averaging of atomistic displacements in local atomic volumes associated with each FE node in an interface region. The FEM and MD computational systems are effectively independent and communicate only through an iterative update of their boundary conditions. With the use of statistical averages of the atomistic quantities to couple the two computational schemes, the developed approach is referred to as an embedded statistical coupling method (ESCM). ESCM provides an enhanced coupling methodology that is inherently applicable to three-dimensional domains, avoids discretization of the continuum model to atomic scale resolution, and permits finite temperature states to be applied.

Saether, E.↗
NASA NTRS2006

Analysis of Rapidly Developing Low Cloud Ceilings in a Stable Environment

Forecasters at the Space Meteorology Group (SMG) issue 30 to 90 minute forecasts for low cloud ceilings at the Space Shuttle Landing Facility (TTS) to support Space Shuttle landings. Mission verification statistics have shown ceilings to be the number one forecast challenge for SMG. More specifically, forecasters at SMG are concerned with any rapidly developing clouds/ceilings below 8000 ft in a stable, capped thermodynamic environment. Therefore, the Applied Meteorology Unit (AMU) was tasked to examine archived events of rapid stable cloud formation resulting in ceilings below 8000 ft, and document the atmospheric regimes favoring this type of cloud development. The AMU examined the cool season months of November to March during the years of 1993-2003 for days that had low-level inversions and rapid, stable low cloud formation that resulted in ceilings violating the Space Shuttle Flight Rules. The AMU wrote and modified existing code to identify inversions from the morning (-10 UTC) Cape Canaveral, FL rawinsonde (XMR) during the cool season and output pertinent sounding information. They parsed all days with cloud ceilings below 8000 ft at TTS, forming a database of possible rapidly-developing low ceiling events. Days with precipitation or noticeable fog burn-off situations were excluded from the database. In the first phase of this work, only the daytime hours were examined for possible ceiling development events since low clouds are easier to diagnose with visible satellite imagery. Phase II of this work includes expanding the database to include nighttime cases which is underway as this abstract is being written. For the nighttime cases, the AMU will analyze both the 00 UTC soundings and the 10 UTC soundings to examine those data for the presence of a low-level inversion. The 00 UTC soundings will probably not have a surface-based inversion, but the presence of inversions or "neutral" layers aloft and below 8,000 ft will most likely help define the stable regime, being a thermodynamically "capped" environment. Occurrences of elevated low-level inversions or stable layers will be highlighted in conjunction with nights that experienced a possible development or onset of cloud ceilings below 8,000 ft. Using these criteria to narrow down the database, the AMU will then use archived IR satellite imagery for these possible events. This presentation summarizes the composite meteorological conditions for 20 daytime event days with rapid low cloud ceiling formation and 48 non-events days consisting of advection or widespread low cloud ceilings and describes two sample cases of daytime rapidly-developing low cloud ceilings. The authors will also summarize the work from the nighttime cases and describe a representative sample case from this data set.

Bauman, William H., III↗
NASA NTRS2016

Uncertainty Assessment of the NASA Earth Exchange Global Daily Downscaled Climate Projections (NEX-GDDP) Dataset

The NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset is comprised of downscaled climate projections that are derived from 21 General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5) and across two of the four greenhouse gas emissions scenarios (RCP4.5 and RCP8.5). Each of the climate projections includes daily maximum temperature, minimum temperature, and precipitation for the periods from 1950 through 2100 and the spatial resolution is 0.25 degrees (approximately 25 km x 25 km). The GDDP dataset has received warm welcome from the science community in conducting studies of climate change impacts at local to regional scales, but a comprehensive evaluation of its uncertainties is still missing. In this study, we apply the Perfect Model Experiment framework (Dixon et al. 2016) to quantify the key sources of uncertainties from the observational baseline dataset, the downscaling algorithm, and some intrinsic assumptions (e.g., the stationary assumption) inherent to the statistical downscaling techniques. We developed a set of metrics to evaluate downscaling errors resulted from bias-correction ("quantile-mapping"), spatial disaggregation, as well as the temporal-spatial non-stationarity of climate variability. Our results highlight the spatial disaggregation (or interpolation) errors, which dominate the overall uncertainties of the GDDP dataset, especially over heterogeneous and complex terrains (e.g., mountains and coastal area). In comparison, the temporal errors in the GDDP dataset tend to be more constrained. Our results also indicate that the downscaled daily precipitation also has relatively larger uncertainties than the temperature fields, reflecting the rather stochastic nature of precipitation in space. Therefore, our results provide insights in improving statistical downscaling algorithms and products in the future.

climate projection↗
NASA NTRS2016

Uncertainty Assessment of the NASA Earth Exchange Global Daily Downscaled Climate Projections (NEX-GDDP) Dataset

The NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset is comprised of downscaled climate projections that are derived from 21 General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5) and across two of the four greenhouse gas emissions scenarios (RCP4.5 and RCP8.5). Each of the climate projections includes daily maximum temperature, minimum temperature, and precipitation for the periods from 1950 through 2100 and the spatial resolution is 0.25 degrees (approximately 25 km by 25 km). The GDDP dataset has received warm welcome from the science community in conducting studies of climate change impacts at local to regional scales, but a comprehensive evaluation of its uncertainties is still missing. In this study, we apply the Perfect Model Experiment framework (Dixon et al. 2016) to quantify the key sources of uncertainties from the observational baseline dataset, the downscaling algorithm, and some intrinsic assumptions (e.g., the stationary assumption) inherent to the statistical downscaling techniques. We developed a set of metrics to evaluate downscaling errors resulted from bias-correction ("quantile-mapping"), spatial disaggregation, as well as the temporal-spatial non-stationarity of climate variability. Our results highlight the spatial disaggregation (or interpolation) errors, which dominate the overall uncertainties of the GDDP dataset, especially over heterogeneous and complex terrains (e.g., mountains and coastal area). In comparison, the temporal errors in the GDDP dataset tend to be more constrained. Our results also indicate that the downscaled daily precipitation also has relatively larger uncertainties than the temperature fields, reflecting the rather stochastic nature of precipitation in space. Therefore, our results provide insights in improving statistical downscaling algorithms and products in the future.

general circulation model (GCM)↗
NASA NTRS2019

The Cumulus and Stratocumulus CloudSat-CALIPSO Dataset (CASCCAD)

Low clouds continue to contribute greatly to the uncertainty in cloud feedback estimates. Depending on whether a region is dominated by cumulus (Cu) or stratocumulus (Sc) clouds, the interannual low-cloud feedback is somewhat different in both spaceborne and large-eddy simulation studies. Therefore, simulating the correct amount and variation of the Cu and Sc cloud distributions could be crucial to predict future cloud feedbacks. Here we document spatial distributions and profiles of Sc and Cu clouds derived from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) and CloudSat measurements. For this purpose, we create a new dataset called the Cumulus And Stratocumulus CloudSat-CALIPSO Dataset (CASCCAD), which identifies Sc, broken Sc, Cu under Sc, Cu with stratiform outflow and Cu. To separate the Cu from Sc, we design an original method based on the cloud height, horizontal extent, vertical variability and horizontal continuity, which is separately applied to both CALIPSO and combined CloudSat–CALIPSO observations. First, the choice of parameters used in the discrimination algorithm is investigated and validated in selected Cu, Sc and Sc–Cu transition case studies. Then, the global statistics are compared against those from existing passive- and active-sensor satellite observations. Our results indicate that the cloud optical thickness – as used in passive-sensor observations – is not a sufficient parameter to discriminate Cu from Sc clouds, in agreement with previous literature. Using clustering-derived datasets shows better results although one cannot completely separate cloud types with such an approach. On the contrary, classifying Cu and Sc clouds and the transition between them based on their geometrical shape and spatial heterogeneity leads to spatial distributions consistent with prior knowledge of these clouds, from ground-based, ship-based and field campaigns. Furthermore, we show that our method improves existing Sc–Cu classifications by using additional information on cloud height and vertical cloud fraction variation. Finally, the CASCCAD datasets provide a basis to evaluate shallow convection and stratocumulus clouds on a global scale in climate models and potentially improve our understanding of low-level cloud feedbacks. The CASCCAD dataset (Cesana, 2019, https://doi.org/10.5281/zenodo.2667637) is available on the Goddard Institute for Space Studies (GISS) website at https://data.giss.nasa.gov/clouds/casccad/ (last access: 5 November 2019) and on the zenodo website at https://zenodo.org/record/2667637 (last access: 5 November 2019).

Cesana, Gregory V.↗
NASA NTRS2022

Estimating pixel-level uncertainty in ocean color retrievals from MODIS

The spectral distribution of marine remote sensing reflectance, R(rs), is the fundamental measurement of ocean color science, from which a host of bio-optical and biogeochemical properties of the water column can be derived. Estimation of uncertainty in these derived properties is thus dependent on knowledge of the uncertainty in satellite-retrieved R(rs) (u(c)(R(rs))) at each pixel. Uncertainty in R(rs), in turn, is dependent on the propagation of various uncertainty sources through the R(rs) retrieval process, namely the atmospheric correction (AC). A derivative-based method for uncertainty propagation is established here to calculate the pixel-level uncertainty in R(rs), as retrieved using NASA’s multiple-scattering epsilon (MSEPS) AC algorithm and verified using Monte Carlo (MC) analysis. The approach is then applied to measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite, with uncertainty sources including instrument random noise, instrument systematic uncertainty, and forward model uncertainty. The uc(Rrs) is verified by comparison with statistical analysis of coincident retrievals from MODIS and in situ Rrs measurements, and our approach performs well in most cases. Based on analysis of an example 8-day global products, we also show that relative uncertainty in R(rs) at blue bands has a similar spatial pattern to the derived concentration of the phytoplankton pigment chlorophyll-a (chl-a), and around 7.3%, 17.0%, and 35.2% of all clear water pixels (chl-a ≤ 0.1 mg/cu.m) with valid u(c)(R(rs)) have a relative uncertainty ≤ 5% at bands 412 nm, 443 nm, and 488 nm respectively, which is a common goal of ocean color retrievals for clear waters. While the analysis shows that u(c)(R(rs)) calculated from our derivative-based method is reasonable, some issues need further investigation, including improved knowledge of forward model uncertainty and systematic uncertainty in instrument calibration.

Pixel-level uncertainty↗
NASA NTRSDate not supplied

Inter-Calibrating the CERES SW Instrument Fluxes Using Deep Convective Clouds(DCC) Within the Same Sun-Synchronous Orbit

The NASA CERES EBAF product provides the scientific community observed TOA fluxes to monitor the Earth’s energy imbalance and to validate climate models. To provide a seamless EBAF 24-year record, the CERES instrument calibration must be stable and consistent across satellite records. Once the Terra and Aqua spacecraft are decommissioned at the end of 2025, the CERES project will rely on CERES instrument observations on the SNPP and NOAA-20 1:30 PM sun-synchronous orbits. The SNPP and NOAA-20 satellite orbit placement will prevent any time matched observations for inter-calibration efforts. The future Libera instrument, which will continue the CERES record and scheduled for launch in 2028 onboard the NOAA-22 satellite, will also need to be inter-calibrated without the aid of time-matched observations. Deep convective clouds (DCC) are the most Lambertian, brightest, tropical Earth invariant targets located at the tropopause making them ideal to radiometrically scale the CERES SW observed radiance to a common reference. The SNPP, NOAA-20 and future NOAA-22 satellites are in the same 16-day repeatable orbits, allowing the DCC targets to be observed with the same angular configuration. The empirically derived SNPP CERES SW channel Bidirectional Reflectance Distribution Function (BRDF) can be easily be applied to the NOAA-20 CERES SW channel radiances since they observe the nearly the same DCC systems. The DCC BRDF corrected radiances are analyzed collectively into radiance probability distribution functions (PDF). By comparing the PDF statistics, the SNPP and NOAA-20 CERES SW channel can be inter-calibrated. The inter-calibration coefficients are validated with the CERES instrument calibration team coefficients to optimize the methodology.

David Doelling↗
NASA NTRSDate not supplied

PCB Quality Metrics that Drive Reliability

Risk-based technology infusion is a deliberate and systematic process that defines the analysis and communication methodology by which new technology is applied and integrated into existing and new designs, identifies technology development needs based on trends analysis and facilitates the identification of shortfalls against performance objectives. This presentation at IPC Works Asia Aerospace 2021 Events provides the audience a snapshot of quality variations in printed-wiring-board quality, as assessed, using experiences in processing and risk analysis of PWB structural integrity coupons. The presentation will focus on printed-wiring-board quality metrics used, the relative type and number of non-conformances observed and trend analysis using statistical methods. Trend analysis shows the top five non-conformances observed across PWB suppliers, the root cause(s) behind these non-conformance and suggestions of mitigation plans. The trends are then matched with the current state of the PWB supplier base and its challenges and opportunities. The presentation further discusses the risk-based SMA approaches and methods being applied at GSFC for evaluating candidate printed wiring board technologies that promote the adoption of higher throughput and faster processing technology for GSFC missions.

Printed circuit boards↗
NASA NTRS1981

Statistical estimation of correlation for nonstationary aircraft noise

Based on the acoustic radiation theory for a moving source, the estimation of the correlation function of a nonstationary acoustic signal is considered. An estimator is obtained for the acoustic correlation, which is shown to be unbiased and consistent. The estimation procedure can be applied to an aircraft noise measurement, where the reflected wave can be eliminated.

Chow, P. L.↗
DOE OSTI2024

Localization of infrasonic sources via Bayesian back projection

SUMMARY A Bayesian framework is investigated for event-specific localization of infrasonic sources using back projection ray tracing. Direction-of-arrival information from array-based detection analysis is used to initialize a back projection ray path originating from the detecting array location and quantifying propagation characteristics from hypothetical source locations. The Fisher statistic, computed from the array’s beam coherence, is mapped into uncertainty in the launch angles of the ray path. Auxiliary parameters previously introduced for solving the Transport equation to compute geometric spreading along ray paths are used to map uncertainty in the ray launch angles into spatial and temporal uncertainties in the ray path. An atmospheric ensemble approach is applied to account for atmospheric uncertainty, and the relation between uncertainties in the atmospheric state and confidence in estimated localization are evaluated using several ensembles with specified variances. The method is evaluated using a synthetic event in the western United States constructed via forward propagation simulations as well as a single-station, multi-arrival detection from a surface explosion in the western United States. Localization results using this event-specific approach are more accurate and exhibit improved precision than existing Bayesian localization methods that leverage generalized, pre-computed propagation statistics.

58 GEOSCIENCES↗
DOE OSTI2025

Breaking the curse of dimensionality: Solving configurational integrals for crystalline solids by tensor networks

Accurately evaluating configurational integrals for dense solids remains a central and difficult challenge in the statistical mechanics of condensed systems. Here, we present a tensor network approach that reformulates the high-dimensional configurational integral for identical-particle crystals into a sequence of computationally efficient summations. We represent the integrand as a high-dimensional tensor and apply tensor-train (TT) decomposition together with a custom TT-cross interpolation. This approach circumvents the need to explicitly construct the full tensor. We introduce tailored rank-1 and rank-2 schemes optimized for sharply peaked Boltzmann probability densities, typical for identical-particle crystals. When applied to the calculation of internal energy and pressure-temperature curves for crystalline Cu and Ar at high (GPa) pressures, as well as the alpha-to-beta phase transition diagram of Sn, our method accurately reproduces molecular dynamics simulation results using tight-binding, machine learning, hierarchical interacting particle–neural network, and modified embedded atom method potentials,all within seconds of computation time.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1999

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗