National Space Science Center

Abstract In December 2019, the International Association of Geomagnetism and Aeronomy (IAGA) Division V Working Group (V-MOD) adopted the thirteenth generation of the International Geomagnetic Reference Field (IGRF). This IGRF updates the previous generation with a definitive main field model for epoch 2015.0, a main field model for epoch 2020.0, and a predictive linear secular variation for 2020.0 to 2025.0. This letter provides the equations defining the IGRF, the spherical harmonic coefficients for this thirteenth generation model, maps of magnetic declination, inclination and total field intensity for the epoch 2020.0, and maps of their predicted rate of change for the 2020.0 to 2025.0 time period.

We report on the population properties of compact binary mergers inferred from gravitational-wave observations of these systems during the first three LIGO-Virgo observing runs. The Gravitational-Wave Transient Catalog 3 (GWTC-3) contains signals consistent with three classes of binary mergers: binary black hole, binary neutron star, and neutron star–black hole mergers. We infer the binary neutron star merger rate to be between 10 and 1700 Gpc − 3 yr − 1 and the neutron star–black hole merger rate to be between 7.8 and 140 Gpc − 3 yr − 1 , assuming a constant rate density in the comoving frame and taking the union of 90% credible intervals for methods used in this work. We infer the binary black hole merger rate, allowing for evolution with redshift, to be between 17.9 and 44 Gpc − 3 yr − 1 at a fiducial redshift ( z = 0.2 ). The rate of binary black hole mergers is observed to increase with redshift at a rate proportional to ( 1 + z ) κ with κ = 2. 9 − 1.8 + 1.7 for z ≲ 1 . Using both binary neutron star and neutron star–black hole binaries, we obtain a broad, relatively flat neutron star mass distribution extending from 1.2 − 0.2 + 0.1 to 2.0 − 0.3 + 0.3 M ⊙ . We confidently determine that the merger rate as a function of mass sharply declines after the expected maximum neutron star mass, but cannot yet confirm or rule out the existence of a lower mass gap between neutron stars and black holes. We also find the binary black hole mass distribution has localized over- and underdensities relative to a power-law distribution, with peaks emerging at chirp masses of 8.3 − 0.5 + 0.3 and 27.9 − 1.8 + 1.9 M ⊙ . While we continue to find that the mass distribution of a binary’s more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above approximately 60 M ⊙ , which would indicate the presence of a upper mass gap. Observed black hole spins are small, with half of spin magnitudes below χ i ≈ 0.25 . While the majority of spins are preferentially aligned with the orbital angular momentum, we infer evidenc

We perform decoy-state quantum key distribution between a low-Earth-orbit satellite and multiple ground stations located in Xinglong, Nanshan, and Graz, which establish satellite-to-ground secure keys with ∼kHz rate per passage of the satellite Micius over a ground station. The satellite thus establishes a secure key between itself and, say, Xinglong, and another key between itself and, say, Graz. Then, upon request from the ground command, Micius acts as a trusted relay. It performs bitwise exclusive or operations between the two keys and relays the result to one of the ground stations. That way, a secret key is created between China and Europe at locations separated by 7600 km on Earth. These keys are then used for intercontinental quantum-secured communication. This was, on the one hand, the transmission of images in a one-time pad configuration from China to Austria as well as from Austria to China. Also, a video conference was performed between the Austrian Academy of Sciences and the Chinese Academy of Sciences, which also included a 280 km optical ground connection between Xinglong and Beijing. Our work clearly confirms the Micius satellite as a robust platform for quantum key distribution with different ground stations on Earth, and points towards an efficient solution for an ultralong-distance global quantum network.

We report results of a search for an isotropic gravitational-wave background (GWB) using data from Advanced LIGO's and Advanced Virgo's third observing run (O3) combined with upper limits from the earlier O1 and O2 runs. Unlike in previous observing runs in the advanced detector era, we include Virgo in the search for the GWB. The results of the search are consistent with uncorrelated noise, and therefore we place upper limits on the strength of the GWB. We find that the dimensionless energy density ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{GW}}\ensuremath{\le}5.8\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}9}$ at the 95% credible level for a flat (frequency-independent) GWB, using a prior which is uniform in the log of the strength of the GWB, with 99% of the sensitivity coming from the band 20--76.6 Hz; ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{GW}}(f)\ensuremath{\le}3.4\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}9}$ at 25 Hz for a power-law GWB with a spectral index of $2/3$ (consistent with expectations for compact binary coalescences), in the band 20--90.6 Hz; and ${\mathrm{\ensuremath{\Omega}}}_{\mathrm{GW}}(f)\ensuremath{\le}3.9\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}10}$ at 25 Hz for a spectral index of 3, in the band 20--291.6 Hz. These upper limits improve over our previous results by a factor of 6.0 for a flat GWB, 8.8 for a spectral index of $2/3$, and 13.1 for a spectral index of 3. We also search for a GWB arising from scalar and vector modes, which are predicted by alternative theories of gravity; we do not find evidence of these, and place upper limits on the strength of GWBs with these polarizations. We demonstrate that there is no evidence of correlated noise of magnetic origin by performing a Bayesian analysis that allows for the presence of both a GWB and an effective magnetic background arising from geophysical Schumann resonances. We compare our upper limits to a fiducial model for the GWB from the merger of compact binaries, updating the model to use the most recent data-driven population inference from the systems detected during O3a. Finally, we combine our results with observations of individual mergers and show that, at design sensitivity, this joint approach may yield stronger constraints on the merger rate of binary black holes at $z\ensuremath{\gtrsim}2$ than can be achieved with individually resolved mergers alone.

There is a growing appreciation that the environmental conditions that we call space weather impact the technological infrastructure that powers the coupled economies around the world. With that comes the need to better shield society against space weather by improving forecasts, environmental specifications, and infrastructure design. We recognize that much progress has been made and continues to be made with a powerful suite of research observatories on the ground and in space, forming the basis of a Sun–Earth system observatory. But the domain of space weather is vast – extending from deep within the Sun to far outside the planetary orbits – and the physics complex – including couplings between various types of physical processes that link scales and domains from the microscopic to large parts of the solar system. Consequently, advanced understanding of space weather requires a coordinated international approach to effectively provide awareness of the processes within the Sun–Earth system through observation-driven models. This roadmap prioritizes the scientific focus areas and research infrastructure that are needed to significantly advance our understanding of space weather of all intensities and of its implications for society. Advancement of the existing system observatory through the addition of small to moderate state-of-the-art capabilities designed to fill observational gaps will enable significant advances. Such a strategy requires urgent action: key instrumentation needs to be sustained, and action needs to be taken before core capabilities are lost in the aging ensemble. We recommend advances through priority focus (1) on observation-based modeling throughout the Sun–Earth system, (2) on forecasts more than 12 h ahead of the magnetic structure of incoming coronal mass ejections, (3) on understanding the geospace response to variable solar-wind stresses that lead to intense geomagnetically-induced currents and ionospheric and radiation storms, and (4) on developing a comprehensive specification of space climate, including the characterization of extreme space storms to guide resilient and robust engineering of technological infrastructures. The roadmap clusters its implementation recommendations by formulating three action pathways, and outlines needed instrumentation and research programs and infrastructure for each of these. An executive summary provides an overview of all recommendations.

Abstract Detections of gravitational waves are now starting to probe the mass distribution of stellar mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernovae. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45 is remarkably robust against variations in the metallicity (≈3 ), the treatment of internal mixing (≈1 ), and stellar wind mass loss (≈4 ), making it the most robust predictor for the final stages of the evolution of massive stars. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying the reaction rate within its 1 σ uncertainties shifts the location of the gap between 40 and 56 . We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a “standard siren” across the universe.

This paper explores the evolution of geoscientific inquiry, tracing the progression from traditional physics-based models to modern data-driven approaches facilitated by significant advancements in artificial intelligence (AI) and data collection techniques. Traditional models, which are grounded in physical and numerical frameworks, provide robust explanations by explicitly reconstructing underlying physical processes. However, their limitations in comprehensively capturing Earth's complexities and uncertainties pose challenges in optimization and real-world applicability. In contrast, contemporary data-driven models, particularly those utilizing machine learning (ML) and deep learning (DL), leverage extensive geoscience data to glean insights without requiring exhaustive theoretical knowledge. ML techniques have shown promise in addressing Earth science-related questions. Nevertheless, challenges such as data scarcity, computational demands, data privacy concerns, and the "black-box" nature of AI models hinder their seamless integration into geoscience. The integration of physics-based and data-driven methodologies into hybrid models presents an alternative paradigm. These models, which incorporate domain knowledge to guide AI methodologies, demonstrate enhanced efficiency and performance with reduced training data requirements. This review provides a comprehensive overview of geoscientific research paradigms, emphasizing untapped opportunities at the intersection of advanced AI techniques and geoscience. It examines major methodologies, showcases advances in large-scale models, and discusses the challenges and prospects that will shape the future landscape of AI in geoscience. The paper outlines a dynamic field ripe with possibilities, poised to unlock new understandings of Earth's complexities and further advance geoscience exploration.

A measurement of the cosmic ray positron fraction e+/(e++e−) in the energy range of 1–30 GeV is presented. The measurement is based on data taken by the AMS-01 experiment during its 10 day Space Shuttle flight in June 1998. A proton background suppression on the order of 106 is reached by identifying converted bremsstrahlung photons emitted from positrons.

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including DNS formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova natal kick distributions, remnant mass prescriptions, and common-envelope energetics. We compute the likelihood of observing the orbital period–eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) on to a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian |$\sigma =30\, \rm km\, s^{-1}$|⁠) and high (⁠|$\sigma =265\, \rm km\, s^{-1}$|⁠) components is preferred over a single high-kick component. We conclude that the observed DNS mass distribution can place strong constraints on model assumptions.

years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time that an experiment directly measures the cosmic ray protons up to ~100 TeV with high statistics. The measured spectrum confirms the spectral hardening at ~300 GeV found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee and sheds new light on the origin of Galactic cosmic rays.

As observed in Thomson-scattered white light, coronal mass ejections (CMEs) are manifest as large-scale expulsions of plasma magnetically driven from the corona in the most energetic eruptions from the Sun. It remains a tantalizing mystery as to how these erupting magnetic fields evolve to form the complex structures we observe in the solar wind at Earth. Here, we strive to provide a fresh perspective on the post-eruption and interplanetary evolution of CMEs, focusing on the physical processes that define the many complex interactions of the ejected plasma with its surroundings as it departs the corona and propagates through the heliosphere. We summarize the ways CMEs and their interplanetary CMEs (ICMEs) are rotated, reconfigured, deformed, deflected, decelerated and disguised during their journey through the solar wind. This study then leads to consideration of how structures originating in coronal eruptions can be connected to their far removed interplanetary counterparts. Given that ICMEs are the drivers of most geomagnetic storms (and the sole driver of extreme storms), this work provides a guide to the processes that must be considered in making space weather forecasts from remote observations of the corona.

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/.

We perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (i) identify which predictions are robust against model uncertainties and assess their implications, (ii) investigate which physical processes leave a clear imprint and may therefore be constrained observationally, and (iii) provide a suite of publicly available model predictions to allow for the use of kinematic constraints from the Gaia mission. We find that 22+26−8% of all massive binary systems merge prior to the first core-collapse in the system. Of the remainder, 86+11−9% become unbound because of the core-collapse. Remarkably, this rarely produces runaway stars (observationally defined as stars with velocities above 30 km s−1). These are outnumbered by more than an order of magnitude by slower unbound companions, or “walkaway stars”. This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. For stars more massive than 15 M⊙, we estimate that 10+5−8% are walkaways and only 0.5+1.0−0.4% are runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies; however, the low runaway fraction we find is in tension with observed fractions of about 10%. Thus, astrometric data on presently single massive stars can potentially constrain the physics of massive binary evolution. Finally, we show that the high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks, and we propose this as a potentially stringent test for the explosion mechanism. We also discuss companions remaining bound that can evolve into X-ray and gravitational wave sources.

Since the beginning of the 21st century, the pace of lunar exploration has accelerated, with more than a dozen probes having undertaken scientific exploration of the Moon. Prominent among these have been the robotic "Chang'E" (CE) missions of the China Lunar Exploration Program (CLEP). We discuss technological and scientific goals and achievements for the four completed, and four planned, CE missions, and longer-term goals and plans of the CLEP beyond the CE missions. The exploration plan is flexible and iterative, with an emphasis on international cooperation.

A new three‐dimensional asymmetric magnetopause model has been developed for corrected GSM coordinates and parameterized by the solar wind dynamic and magnetic pressures ( P d + P m ), the interplanetary magnetic field (IMF) B z , and the dipole tilt angle. On the basis of the magnetopause crossings from Geotail, IMP 8, Interball, TC1, Time History of Events and Macroscale Interactions during Substorms (THEMIS), Wind, Cluster, Polar, Los Alamos National Laboratory (LANL), GOES, and Hawkeye, and the corresponding upstream solar wind parameters from ACE, Wind, or OMNI, this model is constructed by the Levenberg‐Marquardt method for nonlinear multiparameter fitting step‐by‐step over the divided regions. The asymmetries of the magnetopause and the indentations near the cusps are appropriately described in this new model. In addition, the saturation effect of IMF B z on the subsolar distance and the extrapolation for the distant tail magnetopause are also considered. On the basis of this model, the power law index for the subsolar distance versus P d + P m is a bit less than −1/6, the northward IMF B z almost does not influence the magnetopause, and the dipole tilt angle is very important to the north–south asymmetry and the location of indentations. In comparison with the previous empirical magnetopause models based on our database, the new model improves prediction capability to describe the three‐dimensional structure of the magnetopause. It is shown that this new model can be used to quantitatively study how P d + P m compresses the magnetopause, how the southward IMF B z erodes the magnetopause, and how the dipole tilt angle influences the north–south asymmetry and the indentations.

This paper describes the scientific objectives and payloads of Tianwen-1, China’s first exploration mission to Mars. An orbiter, carrying a lander and a rover, lifted-off in July 2020 for a journey to Mars where it should arrive in February 2021. A suite of 13 scientific payloads, for in-situ and remote sensing, autonomously commanded by integrated payload controllers and mounted on the orbiter and the rover will study the magnetosphere and ionosphere of Mars and the relation with the solar wind, the atmosphere, surface and subsurface of the planet, looking at the topography, composition and structure and in particular for subsurface ice. The mission will also investigate Mars climate history. It is expected that Tianwen-1 will contribute significantly to advance our scientific knowledge of Mars.

Solar radiation, which varies over multiple temporal scales, modulates remarkably the evolution of the ionosphere. The solar activity dependence of the ionosphere is a key and fundamental issue in ionospheric physics, providing information essential to understanding the variations in the ionosphere and its processes. Selected recent studies on solar activity effects of the ionosphere are briefly reviewed in this report. This report focuses on (1) observations of solar irradiance at X-ray and extreme ultraviolet wavelengths and the outstanding problems of solar proxies, in the view of ionospheric studies, (2) new findings and improved representations of the features of the solar activity dependence of ionospheric key parameters and the corresponding physical processes, (3) possible phenomena in the ionosphere under extremely high and low solar activity conditions that are unique, as indicated by historical solar datasets and the deep solar minimum of solar cycle 23/24, and (4) statistical studies and model simulations of the ionosphere response to solar flares. The above-mentioned studies provide new clues for comprehensively explaining basic processes in the ionosphere and improving the prediction capability of ionospheric models and related applications.

Abstract We present the first catalog of very-high-energy and ultra-high-energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory. The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array from 2021 March to 2022 September and 933 days of data recorded by the Kilometer Squared Array from 2020 January to 2022 September. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering decl. from −20° to 80°. In total, the catalog contains 90 sources with an extended size smaller than 2° and a significance of detection at >5 σ . Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ( E > 100 TeV) emission at >4 σ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog.

A post sunset bubble manifested by total electron content depletion was observed at midlatitudes (∼30°–34°N, ∼130°–134°E) during the main phase of a storm on 12 February 2000. With loss of lock and the rate of the total electron content index maps, the bubble was seen to bifurcate at its early growth phase. The upward drift speed was observed ∼300 m/s at ∼2150 km, and decreasing with increasing altitude and time. The bubble had unusually large latitudinal extension reaching midlatitude of 36.5°N (31.5°N magnetic latitude), indicating an apex height of ∼2500 km. In process of the evolution, the bubble drifted eastward at a speed of ∼50 m/s. The F region peak height and density obtained by a meridional ionosonde chain suggested a prompt penetrating magnetospheric electric field helped to trigger the super bubble.

We present an updated all-particle energy spectrum of primary cosmic rays in a wide range from 10 14 to 10 17 eV using 5.5 × 10 7 events collected from 2000 November through 2004 October by the Tibet-III air-shower array located 4300 m above sea level (an atmospheric depth of 606 g cm −2 ). The size spectrum exhibits a sharp knee at a corresponding primary energy around 4 PeV. This work uses increased statistics and new simulation calculations for the analysis. We discuss our extensive Monte Carlo calculations and the model dependencies involved in the final result, assuming interaction models QGSJET01c and SIBYLL2.1, and heavy dominant (HD) and proton dominant (PD) primary composition models. Pure proton and pure iron primary models are also examined as extreme cases. A detector simulation was also performed to improve our accuracy in determining the size of the air showers and the energy of the primary particle. We confirmed that the all-particle energy spectra obtained under various plausible model parameters are not significantly different from each other, which was the expected result given the characteristics of the experiment at high altitude, where the air showers of the primary energy around the knee reach near-maximum development, with their features dominated by electromagnetic components, leading to a weak dependence on the interaction model or the primary mass. This is the highest statistical and the best systematics-controlled measurement covering the widest energy range around the knee energy region.