Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance 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,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Relative Importance of Vital Rates to Population Dynamics of Wood Ducks

The wood duck (Aix sponsa) is a common and important cavity-nesting duck in North America; however, we know very little about how changes in vital rates influence population growth rate (λ). We used estimates of fertility and survival of female wood ducks from our nest-box studies in South Carolina, Alabama, and Georgia, USA, to create a stage-based matrix population model. We conducted perturbation analyses and ranked elasticity values to examine the relative importance of 17 component vital rates to λ. Female survival is influenced by nest success, so we recognized this female heterogeneity in our analyses. Four vital rates showed the greatest importance to λ. Analytic elasticities were greatest for breeding season and nonbreeding season survival of females that nested successfully, followed by nest success and female recruitment to the breeding population. Differences in female quality were important to λ. Next, we used process variation of vital rates and conducted life-stage simulation analyses (LSA) followed by variance decomposition to determine the amount of variation in λ explained by each vital rate. Female recruitment to the breeding population explained 57.7% of the variation in λ followed by nest success (11.4%), and breeding and nonbreeding season survival of females that nested successfully (9.3% and 9.4%, respectively). Together these 4 vital rates explained 88% of the variation in λ. Mean asymptotic population growth rate (λ = 0.80 ± 0.08 [SD]) from LSA revealed a declining population. Recruitment of females hatched from nest boxes was insufficient to sustain the nest-box population. Furthermore, including yearling (SY) females that were produced outside of nest boxes (i.e., immigrants) increased recruitment rates 1.5 to 2 times more than when only SY females recruited from nest boxes were included. Future research that examines how emigration and immigration interact with survival and reproduction to influence local population dynamics of wood ducks will be important for identifying the value of nest-box programs to wood duck conservation and management.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Lens covariance effects on likelihood analyses of CMB power spectra

Non-Gaussian correlations induced in CMB power spectra by gravitational lensing needs to be included in likelihood analyses for future CMB experiments. We present a simple but accurate likelihood model which includes these correlations and use it for Markov Chain Monte Carlo parameter estimation from simulated lensed CMB maps in the context of ΛCDM and extensions which include the sum of neutrino masses or the dark energy equation of state w. If lensing-induced covariance is not taken into account for a CMB-S4 type experiment, the errors for one combination of parameters in each case would be underestimated by more then a factor of two and lower limits on w could be misestimated substantially. The frequency of falsely ruling out the true model or finding tension with other data sets would also substantially increase. Our analysis also enables a separation of lens and unlensed information from CMB power spectra, which provides for consistency tests of the model and, if combined with other such measurements, a nearly lens-sample-variance free test for systematics and new physics in the unlensed spectrum. This parameterization also leads to a simple effective likelihood that can be used to assist model building in case consistency tests of ΛCDM fail.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Steady induction effects in geomagnetism. Part 1C: Geomagnetic estimation of steady surficial core motions: Application to the definitive geomagnetic reference field models

In the source-free mantle/frozen-flux core magnetic earth model, the non-linear inverse steady motional induction problem was solved using the method presented in Part 1B. How that method was applied to estimate steady, broad-scale fluid velocity fields near the top of Earth's core that induce the secular change indicated by the Definitive Geomagnetic Reference Field (DGRF) models from 1945 to 1980 are described. Special attention is given to the derivation of weight matrices for the DGRF models because the weights determine the apparent significance of the residual secular change. The derived weight matrices also enable estimation of the secular change signal-to-noise ratio characterizing the DGRF models. Two types of weights were derived in 1987-88: radial field weights for fitting the evolution of the broad-scale portion of the radial geomagnetic field component at Earth's surface implied by the DGRF's, and general weights for fitting the evolution of the broad-scale portion of the scalar potential specified by these models. The difference is non-trivial because not all the geomagnetic data represented by the DGRF's constrain the radial field component. For radial field weights (or general weights), a quantitatively acceptable explication of broad-scale secular change relative to the 1980 Magsat epoch must account for 99.94271 percent (or 99.98784 percent) of the total weighted variance accumulated therein. Tolerable normalized root-mean-square weighted residuals of 2.394 percent (or 1.103 percent) are less than the 7 percent errors expected in the source-free mantle/frozen-flux core approximation.

Voorhies, Coerte V.↗
DOE OSTI2025

HOD-dependent systematics in Emission Line Galaxies for the DESI 2024 BAO analysis

The Dark Energy Spectroscopic Instrument (DESI) will provide precise measurements of Baryon Acoustic Oscillations (BAO) to constrain the expansion history of the Universe and set stringent constraints on dark energy. Therefore, precise control of the global error budget due to various systematic effects is required for the DESI 2024 BAO analysis. In this work, we estimate the level of systematics induced in the DESI BAO analysis due the assumed Halo Occupation Distribution (HOD) model for the Emission Line Galaxy (ELG) tracer. We make use of mock galaxy catalogs constructed by fitting various HOD models to early DESI data, namely the One-Percent survey data. Our analysis includes typical HOD models for the ELG tracer used in the literature as well as extensions to the baseline models. Among the extensions, we consider various recipes for galactic conformity and assembly bias. We use 25 AbacusSummit simulations under the ΛCDM cosmology for each HOD model and perform independent analyses in Fourier space and in configuration space. To recover the BAO signal from our mocks we perform BAO reconstruction and apply the control variates technique to reduce sample variance noise. Our BAO analyses can recover the isotropic BAO parameter α iso within 0.1% and the Alcock Paczynski parameter α AP within 0.3%. Overall, we find that the systematic error due to the HOD dependence is below 0.17%, with the Fourier space analysis being more robust against the HOD systematics. We conclude that our analysis pipeline is robust enough against the HOD systematics for the ELG tracer in the DESI 2024 BAO analysis, for the assumptions made.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2000

Comparison of On-Orbit Performance of Rate Sensing Gyroscopes

This work presents results from the study of a large volume of spacecraft flight data pertaining to gyroscope performance. We have examined long and short term trends of gyroscope biase of first, an exact solution for the time-dependence of the attitude part of the state error covariance, averaged over the three spacecraft axes, is given. This solution is more complete than the usual cubic polynomial for the variance in that it includes coupling among the three gyroscope axes and includes an important term arising from the initial correlations. Second, several continuous 24-hour spans of gyroscope data are examined to verify the short-term statistical model. This analysis demonstrates that in-flight data can be used to determine the strength of the white noise driving the random walk of the gyroscope bias. This may be useful for postlaunch improvement to the noise model and for diagnosing the health of the gyroscope. Third, the long-term trends in gyroscope biases show a nearly linear systematic variation over time scales of years. This has been found on three different missions. While the random walk model is adequate as a basis for onboard Kalman filters or for state estimation using relatively short time spans, these trends indicate that some applications could benefit by accounting for the secular changes in the biases. One example is a new gyroscope calibration method that is under development that allows for multi-epoch bias solutions.

Sedlak, Joseph↗
DOE OSTI2017

Inferring subsurface heterogeneity from push-drift tracer tests

In this study, we consider the late-time tailing in a tracer test performed with a push-drift methodology (i.e., quasi-radial injection followed by drift under natural gradient). Numerical simulations of such tests are performed on 1000 multi-Gaussian 2-D log-hydraulic conductivity field realizations of varying heterogeneity, each under eight distinct mean flow directions. The ensemble pdfs of solute return times are found to exhibit power law tails for each considered variance of the log-hydraulic conductivity field, σ$2\atop{ln K}$. The tail exponent is found to relate straightforwardly to σ$2\atop{ln K}$ and, within the parameter space we explored, to be independent of push-phase pumping rate, pumping duration, and local-scale dispersivity. We conjecture that individual push-drift tracer tests in wells with screened intervals much greater than the vertical correlation length of the aquifer will exhibit quasi-ergodicity and that their tail exponent may be used to infer σ$2\atop{ln K}$. Lastly, we calibrate a predictive relationship of this sort from our Monte Carlo study, and apply it to data from a push-drift test performed at a site of approximately known heterogeneity—closely matching the existing best estimate of heterogeneity.

58 GEOSCIENCES↗
DOE OSTI2022

SCALE 6.2.4 Validation: Radiation Shielding

For safe and reliable use of computer codes by the community, accuracy must be clearly evaluated. In particular, the nuclear reactor engineering and licensing field needs accurate tools for radiation shielding modeling. Monaco with Automated Variance Reduction using Importance Calculations (MAVRIC) is one such tool, with built-in variance reduction methods distributed within the SCALE code, and its validity is demonstrated in this report for the released version 6.2.4. Representative benchmarks corresponding to shielding analysis are selected for the validation study. Typical experimental results analyzed from those benchmarks include neutron fluxes, detector count rates, detector energy response functions, neutron and gamma doses, foil neutron activation rates and activities, neutron leakage fluxes, and skyshine dose rates. Thousands of points of comparison between experiment and calculation are presented in this work. Other than rare outliers typically explained by either a lack of information or large uncertainties in the experiment conditions, material, or dimensions, MAVRIC agrees well with the experiment results. MAVRIC is also compared to Monte Carlo N-Particle (MCNP) calculations when available, and both codes generally produce good agreements within estimated uncertainties. The selected benchmarks are obtained from reliable sources such as the International Criticality Safety Benchmark Evaluation Project Handbook (ICSBEP Handbook), the Shielding Integral Benchmark Archive & Database (SINBAD), and other shielding validation work found in the literature. Additional datapoints and benchmarks will be added to future versions of this report to incrementally expand the shielding validation suite incrementally.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1992

Estimation of surface heat and moisture fluxes over a prairie grassland. II - Two-dimensional time filtering and site variability

The behavior of in situ measurements of surface fluxes obtained during FIFE 1987 is examined by using correlative and spectral techniques in order to assess the significance of fluctuations on various time scales, from subdiurnal up to synoptic, intraseasonal, and annual scales. The objectives of this analysis are: (1) to determine which temporal scales have a significant impact on areal averaged fluxes and (2) to design a procedure for filtering an extended flux time series that preserves the basic diurnal features and longer time scales while removing high frequency noise that cannot be attributed to site-induced variation. These objectives are accomplished through the use of a two-dimensional cross-time Fourier transform, which serves to separate processes inherently related to diurnal and subdiurnal variability from those which impact flux variations on the longer time scales. A filtering procedure is desirable before the measurements are utilized as input with an experimental biosphere model, to insure that model based intercomparisons at multiple sites are uncontaminated by input variance not related to true site behavior. Analysis of the spectral decomposition indicates that subdiurnal time scales having periods shorter than 6 hours have little site-to-site consistency and therefore little impact on areal integrated fluxes.

Crosson, William L.↗
NASA NTRS2003

Flexor bias of joint position in humans during spaceflight

The ability to estimate ankle and elbow joint position was tested before, during, and after a 17-day spaceflight. Subjects estimated targeted joint angles during isovelocity (IsoV) joint movements with agonist muscle groups either active or relaxed. These movements included elbow extension (EE) and elbow flexion (EF), and plantarflexion (PF) and dorsiflexion (DF) of the ankle. Subjects also estimated these joint positions while moving the dynamometer at their chosen (variable) velocity (VarV) during EE and PF. For IsoV tests, no differences were observed between active and passive movements for either the ankle or elbow. Compared with those of pre-flight test days, estimates of targeted elbow joint angles were approximately 5 degrees to 15 degrees more flexed in-flight, and returned toward the pre-flight values during recovery. The spaceflight effects for the ankle were inconsistent and less prevalent than those for the elbow. The VarV PF test condition for the 120 degrees target angle at the ankle exhibited approximately 5 degrees to 7 degrees more DF target angle estimates in-flight compared with those pre- or post-flight. In contrast, during IsoV PF there was a tendency for ankle estimates to be approximately 2 degrees to 3 degrees more PF after 2-3 days exposure to spaceflight. These data indicate that during spaceflight the perception of elbow extension is greater than actuality, and are consistent with the interpretation that microgravity induced a flexor bias in the estimation of the actual elbow joint position. Moreover, these effects in joint proprioception during spaceflight were observed in individual isolated single-joint movements during tasks in which vestibular function in maintaining posture were minimal.

STS-78 Shuttle Project↗
DOE OSTI2022

Multi-variance replica exchange SGMCMC for inverse and forward problems via Bayesian PINN

Physics-informed neural network (PINN) has been successfully applied in solving a variety of nonlinear non-convex forward and inverse problems. However, the training is challenging because of the non-convex loss functions and the multiple optima in the Bayesian inverse problem. In this work, we propose a multi-variance replica exchange stochastic gradient Langevin dynamics method to tackle the challenge of the multiple local optima in the optimization and the challenge of the multiple modal posterior distribution in the inverse problem. Replica exchange methods are capable of escaping from the local traps and accelerating the convergence; two chains with different temperatures are designed where the low temperature chain aims for the local convergence, and the target of the high temperature chain is to travel globally and explore the whole loss function entropy landscape. However, it may not be efficient to solve mathematical inversion problems by using the vanilla replica method directly since the method doubles the computational cost in evaluating the forward solvers (likelihood functions) in the two chains. To address this issue, we propose to make different assumptions on the energy function estimation and this facilities one to use solvers of different fidelities in the likelihood function evaluation. More precisely, one can use a solver with low fidelity in the high temperature chain while using a solver with high fidelity in the low temperature chain. Our proposed method significantly lowers the computational cost in the high temperature chain, meanwhile preserving the accuracy and converging very fast. Here we give an unbiased estimate of the swapping rate and give an estimation of the discretization error of the scheme. To verify our idea, we design and solve four inverse problems which have multiple modes. The proposed method is also employed to train the Bayesian PINN to solve the forward and inverse problems; faster and more accurate convergence has been observed when compared to the stochastic gradient Langevin dynamics (SGLD) method and vanilla replica exchange methods.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1999

A Multiphase Model for the Intracluster Medium

Constraints on the clustered mass density of the universe derived from the observed population mean intracluster gas fraction of x-ray clusters may be biased by reliance on a single-phase assumption for the thermodynamic structure of the intracluster medium (ICM). We propose a descriptive model for multiphase structure in which a spherically symmetric ICM contains isobaric density perturbations with a radially dependent variance. Fixing the x-ray emission and emission weighted temperature, we explore two independently observable signatures of the model in the parameter space. For bremsstrahlung dominated emission, the central Sunyaev-Zel'dovich (SZ) decrement in the multiphase case is increased over the single-phase case and multiphase x-ray spectra in the range 0.1-20 keV are flatter in the continuum and exhibit stronger low energy emission lines than their single-phase counterpart. We quantify these effects for a fiducial 10e8 K cluster and demonstrate how the combination of SZ and x-ray spectroscopy can be used to identify a preferred location in the plane of the model parameter space. From these parameters the correct value of mean intracluster gas fraction in the multiphase model results, allowing an unbiased estimate of clustered mass density to he recovered.

Nagai, Daisuke↗
DOE OSTI2014

SU-E-T-108: Development of a Novel Clinical Neutron Dose Monitor for Proton Therapy Based On Twin TLD500 Chips in a Small PE Moderator

Purpose: In proton therapy, it could be desirable to measure out-of-field fast neutron doses at critical locations near and outside the patient body. Methods: The working principle of a novel clinical neutron dose monitor is verified by MCNPX simulation. The device is based on a small PE moderator of just 5.5cm side length for easy handling covered with a thermal neutron suppression layer. In the simulation, a polystyrene phantom is bombarded with a standard proton beam. The secondary thermal neutron flux produced inside the moderator by the impinging fast neutrons from the treatment volume is estimated by pairs of α-Al2O3:C (TLD500) chips which are evaluated offline after the treatment either by TL or OSL methods. The first chip is wrapped with 0.5mm natural Gadolinium foil converting the thermal neutrons to gammas via (n,γ) reaction. The second chip is wrapped with a dummy material. The chip centers have a distance of 2cm from each other. Results: The simulation shows that the difference of gamma doses in the TLD500 chips is correlated to the mean fast neutron dose delivered to the moderator material. Different outer shielding materials have been studied. 0.5mm Cadmium shielding is preferred for cost reasons and convenience. Replacement of PE moderator material by other materials like lead or iron at any place is unfavorable. The spatial orientation of the moderator cube is uncritical. Using variance reduction techniques like splitting/Russian roulette, the TLD500 gamma dose simulation give positive differences up to distances of 0.5m from the treatment volume. Conclusion: Applicability and basic layout of a novel clinical neutron dose monitor are demonstrated. The monitor measures PE neutron doses at locations outside the patient body up to distances of 0.5m from the treatment volume. Tissue neutron doses may be calculated using neutron kerma factors.

07 ISOTOPE AND RADIATION SOURCES↗
NASA NTRS2023

The Statistical Emulators of GGCMI Phase 2: Responses of Year-to-Year Variation of Crop Yield to CO 2 , Temperature, Water and Nitrogen Perturbations

Understanding the impact of climate change on year-to-year variation of crop yield is critical to global food stability and security. While crop model emulators are believed to be lightweight tools to replace the models, few emulators have been developed to capture such interannual variation of crop yield in response to climate variability. In this study, we developed a statistical emulator with a machine learning algorithm to reproduce the response of year-to-year variation of four crop yields to CO 2 (C), temperature (T), water (W), and nitrogen (N) perturbations defined in the Global Gridded Crop Model Intercomparison Project (GGCMI) phase 2. The emulators were able to explain more than 52 % of the variance of simulated yield and performed well in capturing the year-to-year variation of global average and gridded crop yield over current croplands in the baseline. With the changes in CO 2 –temperature–water–nitrogen (CTWN) perturbations, the emulators could reproduce the year-to-year variation of crop yield well over most current cropland. The variation of R and the mean absolute error was small under the single CTWN perturbations and dual-factor perturbations. These emulators thus provide statistical response surfaces of yield, including both its mean and interannual variability, to climate factors. They could facilitate spatiotemporal downscaling of crop model simulation, projecting the changes in crop yield variability in the future and serving as a lightweight tool for multi-model ensemble simulation. The emulators enhanced the flexibility of crop yield estimates and expanded the application of large-ensemble simulations of crop yield under climate change.

Climate change↗
DOE OSTI2024

Time series methods for the analysis of soundscapes and other cyclical ecological data

Biodiversity monitoring has entered an era of ‘big data’, exemplified by a near-continuous collection of sounds, images, chemical and other signals from organisms in diverse ecosystems. Such data streams have the potential to help identify new threats, assess the effectiveness of conservation interventions, as well as generate new ecological insights. However, appropriate analytical methods are often still missing, particularly with respect to characterizing cyclical temporal patterns. Here, we present a framework for characterizing and analysing ecological responses that represent nonstationary, complex temporal patterns and demonstrate the value of using Fourier transforms to decorrelate continuous data points. In our example, we use a framework based on three approaches (spectral analysis, magnitude squared coherence, and principal component analysis) to characterize differences in tropical forest soundscapes within and across sites and seasons in Gabon. By reconstructing the underlying, cyclic behaviour of the soundscape for each site, we show how one can identify circadian patterns in acoustic activity. Soundscapes in the dry season had a complex diel cycle, requiring multiple harmonics to represent daily variation, while in the wet season there was less variance attributable to the daily cyclic patterns. Our framework can be applied to most continuous, or near-continuous ecological data collected at a fine temporal resolution, allowing ecologists to explore patterns of temporal autocorrelation at multiple levels for biologically meaningful trends. Such methods will become indispensable as biological big data are used to understand the impact of anthropogenic pressures on biodiversity and to inform efforts to mitigate them.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Measuring the largest angular scale CMB B-mode polarization with galactic foregrounds on a cut sky

We consider the effectiveness of foreground cleaning in the recovery of Cosmic Microwave Background (CMB) polarization sourced by gravitational waves for tensor-to-scalar ratios in the range 0<r<0.1. Using the planned survey area, frequency bands, and sensitivity of the Cosmology Large Angular Scale Surveyor (CLASS), we simulate maps of Stokes Q and U parameters at 40, 90, 150, and 220 GHz, including realistic models of the CMB, diffuse Galactic thermal dust and synchrotron foregrounds, and Gaussian white noise. We use linear combinations (LCs) of the simulated multifrequency data to obtain maximum likelihood estimates of r, the relative scalar amplitude s, and LC coefficients. We find that for 10,000 simulations of a CLASS-like experiment using only measurements of the reionization peak (ℓ⩽23), there is a 95% C.L. upper limit of r<0.017 in the case of no primordial gravitational waves. For simulations with r=0.01, we recover at 68% C.L. r=0.012{sub −0.006}{sup +0.011}. The reionization peak corresponds to a fraction of the multipole moments probed by CLASS, and simulations including 30⩽ℓ⩽100 further improve our upper limits to r<0.008 at 95% C.L. (r=0.010{sub −0.004}{sup +0.004} for primordial gravitational waves with r = 0.01). In addition to decreasing the current upper bound on r by an order of magnitude, these foreground-cleaned low multipole data will achieve a cosmic variance limited measurement of the E-mode polarization’s reionization peak.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Exploring the non-Gaussianity of the cosmic infrared background and its weak gravitational lensing

ABSTRACT Gravitational lensing deflects the paths of photons, altering the statistics of cosmic backgrounds and distorting their information content. We take the cosmic infrared background (CIB), which provides plentiful information about galaxy formation and evolution, as an example to probe the effect of lensing on non-Gaussian statistics. Using the Websky simulations, we first quantify the non-Gaussianity of the CIB, revealing additional detail on top of its well-measured power spectrum. To achieve this, we use needlet-like multipole-band filters to calculate the variance and higher-point correlations. Using our simulations, we show the two-, three- and four-point spectra, and compare our calculated power spectra and bispectra to Planck values. We then lens the CIB, shell-by-shell with corresponding convergence maps, to capture the broad redshift extent of both the CIB and its lensing convergence. The lensing of the CIB changes the three- and four-point functions by a few tens of per cent at large scales, unlike with the power spectrum, which changes by less than two per cent. We expand our analyses to encompass the full intensity probability distribution functions (PDFs) involving all n-point correlations as a function of scale. In particular, we use the relative entropy between lensed and unlensed PDFs to create a spectrum of templates that can allow estimation of lensing. The underlying CIB model is missing the important role of star bursting, which we test by adding a stochastic lognormal term to the intensity distributions. The novel aspects of our filtering and lensing pipeline should prove useful for any radiant background, including line intensity maps.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1984

Thematic Mapper Data Quality and Performance Assessment in Renewable Resources/agriculture Remote Sensing

It is believed that the increased spatial resolution will provide solutions to proportion estimation error due to mixed pixels, and the increased spectral resolution will provide for the identification of important agricultural features such as crop stage, and condition. The results of analyses conducted relative to these hypothesis from sample segments extracted from the 4-band Detroit scene and the 7-band Mississippi County, Arkansas engineering test scene are described. Several studies were conducted to evaluate the geometric and radiometric performance of the TM to determine data viability for the more pertinent investigations of TM utility. In most cases this requirement was more than sufficiently satisfied. This allowed the opportunity to take advantage of detailed ground observations for several of the sample segments to assess class separability and detection of other important features with TM. The results presented regarding these TM characteristics show that not only is the increased definition of the within scene variance captured by the increased spatial and spectral resolution, but that the mid-IR bands (5 and 7) are necessary for optimum crop type classification. Both qualitative and quantitative results are presented that describe the improvements gained with the TM both relative to the MSS and on its own merit.

Bizzell, R. M.↗
DOE OSTI2022

Photometric redshift estimation with convolutional neural networks and galaxy images: Case study of resolving biases in data-driven methods

Deep-learning models have been increasingly exploited in astrophysical studies, but these data-driven algorithms are prone to producing biased outputs that are detrimental for subsequent analyses. In this work, we investigate two main forms of biases: class-dependent residuals, and mode collapse. We do this in a case study, in which we estimate photometric redshift as a classification problem using convolutional neural networks (CNNs) trained with galaxy images and associated spectroscopic redshifts. We focus on point estimates and propose a set of consecutive steps for resolving the two biases based on CNN models, involving representation learning with multichannel outputs, balancing the training data, and leveraging soft labels. The residuals can be viewed as a function of spectroscopic redshift or photometric redshift, and the biases with respect to these two definitions are incompatible and should be treated individually. We suggest that a prerequisite for resolving biases in photometric space is resolving biases in spectroscopic space. Experiments show that our methods can better control biases than benchmark methods, and they are robust in various implementing and training conditions with high-quality data. Our methods hold promises for future cosmological surveys that require a good constraint of biases, and they may be applied to regression problems and other studies that make use of data-driven models. Nonetheless, the bias-variance tradeoff and the requirement of sufficient statistics suggest that we need better methods and optimized data usage strategies.

79 ASTRONOMY AND ASTROPHYSICS↗