Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “binary uncertainty”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

The luminosities of the binary Cepheids SU Cyg, SU Cas, and W Sgr

Absolute magnitudes for binary classical Cepheids were derived by ratioing IUE low dispersion spectra in regions where the blue companion dominates with spectra of main sequence spectral type standards. The IUE spectra are also used to determine the magnitude difference between the Cepheid and the blue companion. Absolute magnitudes determined for SU Cyg and W Sgr agree with the absolute magnitudes of Sandage and Tammann and Caldwell within the + or - Om2 estimated uncertainty, but not with the fainter absolute magnitudes of Schmidt. The absolute magnitude of SU Cas is in better agreement with the PLC relations if it is pulsating in the first overtone mode.

Evans, Nancy Remage↗

Unveiling the chemical fingerprint of phosphorus-rich stars: I. In the infrared region of APOGEE-2

The origin of phosphorus, one of the essential elements for life on Earth, is currently unknown. Prevalent models of Galactic chemical evolution (GCE) systematically underestimate the amount of P compared to observations, especially at low metallicities. The recently discovered P-rich ([P/Fe] ≳ 1.2 dex) and metal-poor ([Fe/H] ≃ –1.0 dex) giants further challenge the GCE models, calling current theories on stellar nucleosynthesis into question. Since the observed low-mass giants are not expected to produce their high P contents themselves, our primary goal is to find clues on their progenitor or polluter. By increasing the number of known P-rich stars, we aim to narrow down a statistically reliable chemical abundance pattern that defines these peculiar stars. In this way, we place more robust constraints on the nucleosynthetic mechanism that causes the unusually high P abundances. In the long term, identifying the progenitor of the P-rich stars may contribute to the search for the source of P in our Galaxy. We performed a detailed chemical abundance analysis based on the high-resolution near-infrared (H band) spectra from the latest data release (DR17) of the APOGEE-2 survey. Employing the BACCHUS code, we measured the abundances of 13 elements in the inspected sample, which is mainly composed of a recent collection of Si-enhanced giants. We also analyzed the orbital motions and compared the abundance results to possible nucleosynthetic formation scenarios, and also to detailed GCE models. These models were produced with the OMEGA+ chemical evolution code, using four different massive star yield sets to investigate different scenarios for massive star evolution. We enlarged the sample of confirmed P-rich stars from 16 to a group of 78 giants, which represents the largest sample of P-rich stars to date. The sample includes the first detection of a P-rich star in a Galactic globular cluster. Significant enhancements in O, Al, Si, and Ce, as well as systematic correlations among the studied elements, unveil the unique chemical fingerprint of the P-rich stars. In contrast, the high [Mg/Fe] and [(C+N)/Fe] found in some of the P-rich stars with respect to P-normal stars is not confirmed over the full sample because of the current uncertainties. Strikingly, the strong overabundance in the α-element Si is accompanied by normal Ca and S abundances. This is at odds with current stellar nucleosynthesis models of massive stars. Our analysis of the orbital motion showed that the P-rich stars do not belong to a locally specific population in the Galaxy. In addition, we confirm that the majority of the sample stars are not part of binary systems.

79 ASTRONOMY AND ASTROPHYSICS↗

From Existing and New Nuclear and Astrophysical Constraints to Stringent Limits on the Equation of State of Neutron-Rich Dense Matter

Through continuous progress in nuclear theory and experiment and an increasing number of neutron-star (NS) observations, a multitude of information about the equation of state (EOS) for matter at extreme densities is available. To constrain the EOS across its entire density range, this information needs to be combined consistently. However, the impact and model dependency of individual observations vary. Given their growing number, assessing the various methods is crucial to compare the respective effects on the EOS and discover potential biases. For this purpose, we present a broad compendium of different constraints and apply them individually to a large set of EOS candidates within a Bayesian framework. Specifically, we explore different ways of how chiral effective field theory and perturbative quantum chromodynamics can be used to place a likelihood on EOS candidates. We also investigate the impact of nuclear experimental constraints, as well as different radio and x-ray observations of NS masses and radii. This is augmented by reanalyses of the existing data from binary neutron star coalescences, in particular of GW170817, with improved models for the tidal waveform and kilonova light curves, which we also utilize to construct a tight upper limit of 2.39 M ⊙ on the TOV mass based on GW170817’s remnant. Our diverse set of constraints is eventually combined to obtain stringent limits on NS properties. We organize the combination in a way to distinguish between constraints where the systematic uncertainties are deemed small and those that rely on less conservative assumptions. For the former, we find the radius of the canonical 1.4 M ⊙ neutron star to be R 1.4 = 12.2 6 − 0.91 + 0.80 km and the TOV mass at M TOV = 2.2 5 − 0.22 + 0.42 M ⊙ (95% credibility). Including all the presented constraints yields R 1.4 = 12.2 0 − 0.48 + 0.50 km and M TOV = 2.3 0 − 0.20 + 0.07 M ⊙ . When comparing these limits to individual data points, we find that the quoted radius of HESS J1731-347 displays noticeable tension with other constraints. Constraining microphysical properties of the EOS proves more challenging. For instance, the symmetry energy slope is restricted to L sym = 48 − 25 + 21 MeV , where this constraint is mainly dominated by our reanalysis of the PREX-II and CREX experiment. Published by the American Physical Society 2025

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Zwicky Transient Facility: Data Processing, Products, and Archive

The Zwicky Transient Facility (ZTF) is a new robotic time-domain survey currently in progress using the Palomar 48-inch Schmidt Telescope. ZTF uses a 47 square degree field with a 600 megapixel camera to scan the entire northern visible sky at rates of ∼3760 square degrees/hour to median depths of g ~ 20.8 and r ~ 20.6 mag (AB, 5σ in 30 sec). We describe the Science Data System that is housed at IPAC, Caltech. This comprises the data-processing pipelines, alert production system, data archive, and user interfaces for accessing and analyzing the products. The real-time pipeline employs a novel image-differencing algorithm, optimized for the detection of point-source transient events. These events are vetted for reliability using a machine-learned classifier and combined with contextual information to generate data-rich alert packets. The packets become available for distribution typically within 13 minutes (95th percentile) of observation. Detected events are also linked to generate candidate moving-object tracks using a novel algorithm. Objects that move fast enough to streak in the individual exposures are also extracted and vetted. We present some preliminary results of the calibration performance delivered by the real-time pipeline. The reconstructed astrometric accuracy per science image with respect to Gaia DR1 is typically 45 to 85 milliarcsec. This is the RMS per-axis on the sky for sources extracted with photometric S/N ≥10 and hence corresponds to the typical astrometric uncertainty down to this limit. The derived photometric precision (repeatability) at bright unsaturated fluxes varies between 8 and 25 millimag. The high end of these ranges corresponds to an airmass approaching ∼2—the limit of the public survey. Photometric calibration accuracy with respect to Pan-STARRS1 is generally better than 2%. The products support a broad range of scientific applications: fast and young supernovae; rare flux transients; variable stars; eclipsing binaries; variability from active galactic nuclei; counterparts to gravitational wave sources; a more complete census of Type Ia supernovae; and solar-system objects.

Frank J. Masci↗

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

In this work, we perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass x-ray binaries (QLMXBs), photospheric radius expansion x-ray bursting sources, and x-ray timing observations of J0030+0451. With this dataset, the form of the prior distribution still has an impact on the posterior mass-radius curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an “intrinsic scattering” contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

79 ASTRONOMY AND ASTROPHYSICS↗

View Angle Effects on MODIS Snow Mapping in Forests

Binary snow maps and fractional snow cover data are provided routinely from MODIS (Moderate Resolution Imaging Spectroradiometer). This paper investigates how the wide observation angles of MODIS influence the current snow mapping algorithm in forested areas. Theoretical modeling results indicate that large view zenith angles (VZA) can lead to underestimation of fractional snow cover (FSC) by reducing the amount of the ground surface that is viewable through forest canopies, and by increasing uncertainties during the gridding of MODIS data. At the end of the MODIS scan line, the total modeled error can be as much as 50% for FSC. Empirical analysis of MODIS/Terra snow products in four forest sites shows high fluctuation in FSC estimates on consecutive days. In addition, the normalized difference snow index (NDSI) values, which are the primary input to the MODIS snow mapping algorithms, decrease as VZA increases at the site level. At the pixel level, NDSI values have higher variances, and are correlated with the normalized difference vegetation index (NDVI) in snow covered forests. These findings are consistent with our modeled results, and imply that consideration of view angle effects could improve MODIS snow monitoring in forested areas.

MODIS↗

Constraining Bedrock Groundwater Residence Times in a Mountain System With Environmental Tracer Observations and Bayesian Uncertainty Quantification

Groundwater residence time distributions provide fundamental insights on the hydrological processes within watersheds. Yet, observations that can constrain groundwater residence times over broad timescales remain scarce in mountain catchment studies. We use environmental tracers (CFC-12, SF 6 , 3 H, and 4 He) to investigate groundwater residence times along a hillslope in the East River Watershed, Colorado, USA. We develop a Bayesian inference framework that applies a Markov-chain Monte Carlo (MCMC) approach to estimate noble gas recharge temperature, elevation, and excess-air parameters and the resulting environmental tracer concentrations. MCMC is then used to propagate the environmental tracer uncertainties to estimates of groundwater mean residence times inferred with lumped parameter models. All samples contain 3 H, CFC-12, and SF 6 in addition to terrigenic 4 He, suggesting a mixture of water characterized by modern and premodern residence times. 4He exponential mean residence times range from hundreds of years at the upslope well to thousands of years at the toe-slope well assuming average crustal production rates. We find that binary mixing residence time distributions with separate young and old mixing fractions are needed to predict the 4 He, CFC-12, SF 6 , and 3 H observations, supporting the importance of flow path mixing in this bedrock system. Our findings that the fractured bedrock hosts groundwater with a mixture of residence times ranging from decades to millennia suggest variable recharge dynamics and flow path mixing along the hillslope and highlight the importance of characterizing groundwater systems with observations that are sensitive to transport over a broad range of residence times.

54 ENVIRONMENTAL SCIENCES↗

Nested Pebble Bed Blanket (NesPeB)

Recent advances in magnetic confinement fusion technology have attracted billions of dollars of investments in startups from venture capitals and corporations, resulting in the development of devices aiming to demonstrate net energy gain in a self-heated burning plasma, such as SPARC (under construction) and others. However, future fusion power plants must operate in regimes that will require technologies far beyond current experience. According to a National Academies of Science, Engineering, and Medicine report, to have nuclear fusion power plants contributing in a timely manner to the planned reduction of atmospheric carbon dioxide, a pilot plant should be built by 2035, and it should demonstrate fusion power production and the performance of the tritium fuel system (requiring a high enough tritium breeding) by 2040. A recognized key technology gap by [26] is the fusion first wall and blanket since no current blanket concept is considered satisfactory or has been built and proven. The first wall and blanket in magnetic fusion reactors form a vital and complex system, as it must satisfy different functions such as power extraction, tritium breeding, plasma containment, radiation shielding, and safety. The list of design requirements is even longer: high enough tritium production for fusion self-sufficiency, low material activation, decay heat and shutdown dose rates, high thermal efficiency, high-capacity factor, high magnets-divertor-vacuum vessel-first wall life, low corrosion, low cost, and intrinsically safe (requiring minimal licensing). Despite fifty-plus years of research, the first wall and blanket concepts proposed suffer from fundamental technical problems and immaturity (TRL=2-3) that jeopardize the timely delivery of a commercial fusion power plant. A fusion first-wall blanket has never been built nor tested, and a "winning", practical functioning design requires enough engineering margins (high enough tritium breeding considering the uncertainty, etc.), manufacturing simplicity, ease of continuous operation, maintenance, and low cost. A new, groundbreaking blanket concept called "Nested Pebble Bed Blanket" (NesPeB) was developed at ORNL under the successful ARPA-E GAMOW FERMI project (patent application allowed by the USPTO). The NesPeB blanket concept addresses current blanket concepts' shortcomings and technical immaturity, paving the way for accelerated delivery of fusion power plants. NesPeB is based on nested pebbles, which are binary-sized lithium-ceramic pebbles enclosed in "Beryllide" perforated and coated spherical shells, which are also binary-sized, stacked on top of each other, forming a "bed" and cooled by Nitrogen gas also "sweeping" the Helium and Tritium generated by the neutron irradiation of Lithium; the vacuum vessel plasma facing material is Molybdenum-96 and -97 with the first wall cooled by Helium while the divertor armor is made of Tungsten. The simulations of the NesPeB blanket using Fusion Reactors Models Integrator (FERMI) are encouraging as they estimate a tritium breeding ratio (TBR) greater than 1.2 using natural Lithium, acceptable pressure drop, and excellent heat transfer properties. Furthermore, the NesPeB blanket is not limited by magneto-hydro-dynamics (MHD) effects, is designed for online refueling, relies on existing tritium extraction technologies, has a simple construction, and limits the corrosion and chemical reactivity problems. NesPeB has the potential to be transformational and disruptive since it can solve all the main, challenging technical problems of fusion device blankets and accelerate a pilot plant delivery for 10 or more years.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The orbital inclination of A0620 - 00 measured polarimetrically

The mass of the degenerate primary in A0620 - 00 is inferred from its spectroscopic mass function to be not less than 3.2 solar masses, making it an excellent candidate for a black hole. The exact value of the mass depends on the orbital inclination. The inclination of a binary system can be determined from the shape of its Stokes parameter light curves if the linear polarization of the system varies as a function of orbital phase. A0620 - 00 over one 8-hour binary period was observed with the 4.5-m equivalent MMT. Its polarization in the visible is variable with orbital phase. The standard theory of Brown et al. (1978) was used to derive an orbital inclination of i = 57 deg (+20 deg, -50 deg), where the error is the 90-percent confidence interval. An inclination of i = 57 deg corresponds to a mass of the compact primary of 6.6 solar masses, but the large uncertainty in the measured value of the inclination allows the derived mass of A0620 - 00 to be as low as 3.8 solar masses. If this is taken to be the maximum mass of any degenerate configuration consistent with general relativity except a black hole, then the mass of A0620 - 00 is still not well enough determined to conclude that it must be a black hole.

Dolan, Joseph F.↗

Impact of Mergers on USA Parameter Estimation for Nonspinning Black Hole Binaries

We investigate the precision with which the parameters describing the characteristics and location of nonspinning black hole binaries can be measured with the Laser Interferometer Space Antenna (LISA). By using complete waveforms including the inspiral, merger and ringdown portions of the signals, we find that LISA will have far greater precision than previous estimates for nonspinning mergers that ignored the merger and ringdown. Our analysis covers nonspinning waveforms with moderate mass ratios, q > or = 1/10, and total masses 10(exp 5) < M/M_{Sun} < 10(exp 7). We compare the parameter uncertainties using the Fisher matrix formalism, and establish the significance of mass asymmetry and higher-order content to the predicted parameter uncertainties resulting from inclusion of the merger. In real-time observations, the later parts of the signal lead to significant improvements in sky-position precision in the last hours and even the final minutes of observation. For comparable mass systems with total mass M/M_{Sun} = approx. 10(exp 6), we find that the increased precision resulting from including the merger is comparable to the increase in signal-to-noise ratio. For the most precise systems under investigation, half can be localized to within O(10 arcmin), and 18% can be localized to within O(1 arcmin).

McWilliams, Sean T.↗

Modeling Macro- and Micro-Scale Turbulent Mixing and Chemistry in Engine Exhaust Plumes

Simulation of turbulent mixing and chemical processes in the near-field plume and plume-vortex regimes has been successfully carried out recently using a reduced gas phase kinetics mechanism which substantially decreased the computational cost. A detailed mechanism including gas phase HOx, NOx, and SOx chemistry between the aircraft exhaust and the ambient air in near-field aircraft plumes is compiled. A reduced mechanism capturing the major chemical pathways is developed. Predictions by the reduced mechanism are found to be in good agreement with those by the detailed mechanism. With the reduced chemistry, the computer CPU time is saved by a factor of more than 3.5 for the near-field plume modeling. Distributions of major chemical species are obtained and analyzed. The computed sensitivities of major species with respect to reaction step are deduced for identification of the dominant gas phase kinetic reaction pathways in the jet plume. Both the near field plume and the plume-vortex regimes were investigated using advanced mixing models. In the near field, a stand-alone mixing model was used to investigate the impact of turbulent mixing on the micro- and macro-scale mixing processes using a reduced reaction kinetics model. The plume-vortex regime was simulated using a large-eddy simulation model. Vortex plume behind Boeing 737 and 747 aircraft was simulated along with relevant kinetics. Many features of the computed flow field show reasonable agreement with data. The entrainment of the engine plumes into the wing tip vortices and also the partial detrainment of the plume were numerically captured. The impact of fluid mechanics on the chemical processes was also studied. Results show that there are significant differences between spatial and temporal simulations especially in the predicted SO3 concentrations. This has important implications for the prediction of sulfuric acid aerosols in the wake and may partly explain the discrepancy between past numerical studies (that employed parabolic or temporal approximations) and the measured data. Finally to address the major uncertainty in the near-field plume modeling related to the plume processing of sulfur compounds and advanced model was developed to evaluate its impact on the chemical processes in the near wake. A comprehensive aerosol model is developed and it is coupled with chemical kinetics and the axisymmetric turbulent jet flow models. The integrated model is used to simulate microphysical processes in the near-field jet plume, including sulfuric acid and water binary homogeneous nucleation, coagulation, non-equilibrium heteromolecular condensation, and sulfur-induced soot activation. The formation and evolution of aerosols are computed and analyzed. The computed results show that a large number of ultra-fine (0.3--0.6 nm in radius) volatile HSO4 - HO embryos are generated in the near-field plume. These embryos further grow in size by self coagulation and condensation. Soot particles can be activated by both heterogeneous nucleation and scavenging of H2SO4-H2O aerosols. These activated soot particles can serve as water condensation nuclei for contrail formation. Conditions under which ice contrails can form behind aircrafts are studied. The sensitivities of the threshold temperature for contrail formation with respect to aircraft propulsion efficiency, relative humidity, and ambient pressure are evaluated. The computed aerosol properties for different extent of fuel sulfur conversion to S(VI) (SO3 and H2SO4) in engine are examined and the results are found to be sensitive to this conversion fraction.

Menon, Suresh↗

Colloidal superionic conductors

Nanoparticles with highly asymmetric sizes and charges that self-assemble into crystals via electrostatics may exhibit behaviors reminiscent of those of metals or superionic materials. Here, we use coarse-grained molecular simulations with underdamped Langevin dynamics to explore how a binary charged colloidal crystal reacts to an external electric field. As the field strength increases, we find transitions from insulator (ionic state), to superionic (conductive state), to laning, to complete melting (liquid state). In the superionic state, the resistivity decreases with increasing temperature, which is contrary to metals, yet the increment decreases as the electric field becomes stronger. Additionally, we verify that the dissipation of the system and the fluctuation of charge currents obey recently developed thermodynamic uncertainty relation. Our results describe charge transport mechanisms in colloidal superionic conductors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope

If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A new interpretation of luminous blue stars

A major revision of current theoretical ideas about the brightest blue stars must be made if Carson's (1976) radiative opacities are adopted in stellar models. Unlike earlier opacities, these exhibit a large 'bump' due to CNO ionization, which leads to very strong central condensation, convective instability, and pulsational instability in hot diffuse stellar envelopes. Despite a number of theoretical uncertainties, the new picture of the structure of very luminous stars is reasonably successful in accounting for a variety of previously unexplained observations. The stellar models for the phase of core hydrogen burning predict large radii and rather cool effective temperatures for O stars and a spreading out of the main-sequence band in the H-R diagram toward luminous cool supergiants for masses higher than about 20 solar masses. In massive X-ray binary systems, circular orbits and supergiant-like visual companions are expected to be quite common. Long-period variability is predicted to exist for massive blue supergiants of luminosity class Ia. The models for helium stars predict large radii and rather cool effective temperatures for Wolf-Rayet stars, as well as multimodal pulsational instability and, possibly, surface turbulence for these stars.

Stothers, R.↗

Radial velocity measurements of the chromospherically-active stars (2): HD 28591 = V492 Per

From two sets of the spectroscopic observations covering a ten year period we have obtained 59 radial velocities of the chromospherically-active star HD 28591 = V492 Per. It is a G9III single-lined spectroscopic binary with a period of 21.2910 days and a circular orbit. The upsilon sin i of 24.6 km/sec, results in a minimum radius 10.3 solar radii. We estimate a distance of 165 +/- 40 pc and an orbital inclination of 65 +/- 25 degrees. The secondary is probably a mid to late-type K dwarf. The star is brighter than the limiting magnitude of the Bright Star Catalogue. The mean photometric and the orbital periods are identical within their uncertainties. Since the star fills a significant fraction of its Roche lobe, about 62%, the photometric light curve may be the result of starspots and a modest ellipticity effect.

Dadonas, V.↗

Search for 22 Na in novae supported by a novel method for measuring femtosecond nuclear lifetimes

Classical novae are thermonuclear explosions in stellar binary systems, and important sources of 26 Al and 22 Na. While γ rays from the decay of the former radioisotope have been observed throughout the Galaxy, 22 Na remains untraceable. Its half-life (2.6 yr) would allow the observation of its 1.275 MeV γ-ray line from a cosmic source. However, the prediction of such an observation requires good knowledge of its nucleosynthesis. The 22 Na(p, γ) 23 Mg reaction remains the only source of large uncertainty about the amount of 22 Na ejected. Its rate is dominated by a single resonance on the short-lived state at 7785.0(7) keV in 23 Mg. Here, we propose a combined analysis of particle-particle correlations and velocity-difference profiles to measure femtosecond nuclear lifetimes. The application of this method to the study of the 23 Mg states, places strong limits on the amount of 22 Na produced in novae and constrains its detectability with future space-borne observatories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The mass transfer rate in X1916-053 - It is driven by gravitational radiation?

A 50-minute period for a binary system harboring an X-ray burster would allow several alternatives for the mass-giving secondary, including an H-shell burning-plus-He degenerate core composite model. The burst properties of X1916-053 are presently used to argue against the He degenerate as well as the He main sequence solutions and to estimate whether, for any of the other solutions, the mass transfer rate could be consistent with that expected from gravitational radiation (GR). Within an uncertainty of a factor of 2, the transfer rate for the composite model solution is consistent with gravitational radiation, but enhancement by other mechanisms should be investigated.

Swank, J. H.↗

Improving Rain/No-Rain Detection Skill by Merging Precipitation Estimates from Different Sources

Rain/no-rain detection error is a key source of uncertainty in regional and global precipitation products that propagates into offline hydrological and land surface modeling simulations. Such detection error is difficult to evaluate and/or filter without access to high-quality reference precipitation datasets. For cases where such access is not available, this study proposes a novel approach for improved rain/no-rain detection. Based on categorical triple collocation (CTC) and a probabilistic framework, a weighted merging algorithm (CTC-M) is developed to combine noisy, but independent, precipitation products into an optimal binary rain/no-rain time series. Compared with commonly used approaches that directly apply the best parent product for rain/no-rain detection, the superiority of CTC-M is demonstrated analytically and numerically using spatially dense precipitation measurements over Europe. Our analysis also suggests that CTC-M is tolerant to a range of cross-correlated rain/no-rain detection errors and detection biases of the parent products. As a result, CTC-M will benefit global precipitation estimation by improving the representation of precipitation occurrence in gauge-based and multisource merged precipitation products.

Jianzhi Dong↗