Nature Astronomy (@natastron.nature.com)

Bowshocks driven by the pole-on molecular jet of outbursting protostar SVS 13 Abstract Outflows play a key role in the star and planet formation processes. Some outflows show discrete clumps of cold molecular gas moving at extremely high velocities (EHVs) of ~100 km s−1, known as ‘molecular bullets’, that are likely closely associated with their primary driving agent. Here we present ALMA CO (J = 3–2) observations of a bright EHV molecular bullet that reveal its morphology in detail down to scales of 30 au and its kinematic structure across the entire intermediate velocity range (~30–100 km s−1). These provide important insights into how outflows transfer mass and momentum to the surrounding medium. The observed channel maps display several sequences of ring-like features whose velocity increases and size decreases with projected distance from the driving source, each sequence tracing a thin, bow-shaped shell culminating on axis in a bright EHV head. The shape, kinematics and mass of each shell all agree remarkably well with the simplest textbook models of momentum-conserving bowshocks produced by a time-variable EHV jet. The dynamical timescale between consecutive shells is of a few decades, with the latest ejection event coinciding with the protostar optical/infrared outburst observed in ~1990. The very strong evidence for bowshock-driven entrainment induced by jet variability revealed by this work suggests that accretion bursts, and therefore variations in the disk snowlines, should occur on decade timescales, which could substantially impact grain growth and planet formation. Access options Get Nature+, our best-value online-access subscription Prices may be subject to local taxes which are calculated during checkout The ALMA raw data are available from the ALMA Science Archive (https://almascience.eso.org/aq/) using the project identifiers 2015.1.01229.S (for Cycle 3 data) and 2016.1.101305.S (for Cycle 4 data). The reduced data is available via Zenodo at https://doi.org/10.5281/zenodo.17249853 (ref. 96). This work has made extensive use of the open-source Python libraries NumPy97, SciPy86, Matplotlib98, Photutils99, Astropy100,101,102 and LMFIT103. The code used to compute the channel maps of the bowshock model is available at GitHub104. References Frank, A. et al. in Protostars and Planets VI (eds Beuther, H. et al.) Ch. 20 (Univ. of Arizona Press, 2014). Krumholz, M. R. & Federrath, C. The role of magnetic fields in setting the star formation rate and the initial mass function. Front. Astron. Space Sci. 6, 7 (2019). Manara, C. F. et al. in Protostars and Planets VII (eds Inutsuka, S. et al.) Ch. 15 (ASP, 2023). Pascucci, I. et al. in Protostars and Planets VII (eds Inutsuka, S. et al.) Ch. 16 (ASP, 2023). Bachiller, R. Bipolar molecular outflows from young stars and protostars. Annu. Rev. Astron. Astrophys. 34, 111–154 (1996). Bally, J. Protostellar outflows. Annu. Rev. Astron. Astrophys. 54, 491–528 (2016). Anglada, G., Rodríguez, L. F. & Carrasco-González, C. Radio jets from young stellar objects. Astron. Astrophys. Rev. 26, 3 (2018). Lee, C.-F. Molecular jets from low-mass young protostellar objects. Astron. Astrophys. Rev. 28, 1 (2020). Ray, T. P. & Ferreira, J. Jets from young stars. New. Astron. Rev. 93, 101615 (2021). Ray, T. P. et al. Outflows from the youngest stars are mostly molecular. Nature 622, 48–52 (2023). Federman, S. A. et al. Investigating protostellar accretion-driven outflows across the mass spectrum: JWST NIRSpec integral field unit 3–5 μm spectral mapping of five young protostars. Astrophys. J. 966, 41 (2024). Reipurth, B. & Bally, J. Herbig-Haro flows: probes of early stellar evolution. Annu. Rev. Astron. Astrophys. 39, 403–455 (2001). Anglada, G. et al. Proper motions of the jets in the region of HH 30 and HL/XZ Tau: evidence for a binary exciting source of the HH 30 jet. Astron. J. 133, 2799–2814 (2007). Reipurth, B. et al. The HH 24 complex: jets, multiple star formation, and orphaned protostars. Astron. J. 165, 209 (2023). Bachiller, R., Cernicharo, J., Martin-Pintado, J., Tafalla, M. & Lazareff, B. High-velocity molecular bullets in a fast bipolar outflow near L 1448/IRS 3. Astron. Astrophys. 231, 174–186 (1990). Bachiller, R., Martin-Pintado, J. & Planesas, P. High-velocity molecular jets and bullets from IRAS 03282+3035. Astron. Astrophys. 251, 639–648 (1991). Hatchell, J., Fuller, G. A. & Ladd, E. F. Hot molecular bullets in HH 111 and CEP E. Astron. Astrophys. 346, 278–284 (1999). Zapata, L. A. et al. A highly collimated, young, and fast CO outflow in OMC-1 South. Astrophys. J. Lett. 630, L85–L88 (2005). Lefloch, B., Cernicharo, J., Reipurth, B., Pardo, J. R. & Neri, R. Anatomy of HH 111 from CO observations: a bow-shock-driven molecular outflow. Astrophys. J. 658, 498–508 (2007). Lee, C.-F. et al. Rotation and outflow motions in the very low-mass class 0 protostellar system HH 211 at subarcsecond resolution. Astrophys. J. 699, 1584–1594 (2009). Santiago-García, J., Tafalla, M., Johnstone, D. & Bachiller, R. Shells, jets, and internal working surfaces in the molecular outflow from IRAS 04166+2706. Astron. Astrophys. 495, 169–181 (2009). Hirano, N. et al. Extreme active molecular jets in L1448C. Astrophys. J. 717, 58–73 (2010). Podio, L. et al. First image of the L1157 molecular jet by the CALYPSO IRAM-PdBI survey. Astron. Astrophys. 593, L4 (2016). Tychoniec, Ł. et al. Chemical and kinematic structure of extremely high-velocity molecular jets in the Serpens Main star-forming region. Astron. Astrophys. 632, A101 (2019). de A. Schutzer, A. et al. SOLIS. XVI. Mass ejection and time variability in protostellar outflows: Cep E. Astron. Astrophys. 662, A104 (2022). Dutta, S. et al. ALMA survey of Orion Planck galactic cold clumps (ALMASOP): molecular jets and episodic accretion in protostars. Astron. J. 167, 72 (2024). Moraghan, A., Lee, C.-F., Huang, P.-S. & Vaidya, B. A study of the wiggle morphology of HH 211 through numerical simulations. Mon. Not. R. Astron. Soc. 460, 1829–1838 (2016). Podio, L. et al. The CALYPSO IRAM-PdBI survey of jets from Class 0 protostars. Exploring whether jets are ubiquitous in young stars. Astron. Astrophys. 648, A45 (2021). Yoshida, T., Hsieh, T.-H., Hirano, N. & Aso, Y. Multi-epoch submillimeter array observations of the L1448C(N) protostellar SiO jet. Astrophys. J. 906, 112 (2021). Rodriguez, L. F., Carral, P., Ho, P. T. P. & Moran, J. M. Anisotropic mass outflow in regions of star formation. Astrophys. J. 260, 635–646 (1982). Lada, C. J. Cold outflows, energetic winds, and enigmatic jets around young stellar objects. Annu. Rev. Astron. Astrophys. 23, 267–317 (1985). Konigl, A. & Pudritz, R. E. in Protostars and Planets IV (eds Mannings, V. et al.) Ch. 28 (Univ. of Arizona Press, 2000). Levreault, R. M. A search for molecular outflows toward pre–main-sequence objects. Astrophys. J. Suppl. 67, 283–371 (1988). Raga, A. & Cabrit, S. Molecular outflows entrained by jet bowshocks. Astron. Astrophys. 278, 267–278 (1993). Tafalla, M. et al. Anatomy of the internal bow shocks in the IRAS 04166+2706 protostellar jet. Astron. Astrophys. 597, A119 (2017). Rabenanahary, M., Cabrit, S., Meliani, Z. & Pineau des Forêts, G. Wide-angle protostellar outflows driven by narrow jets in stratified cores. Astron. Astrophys. 664, A118 (2022). Li, Z.-Y. & Shu, F. H. Interaction of wide-angle MHD winds with flared disks. Astrophys. J. 468, 261–268 (1996). Wang, L.-Y., Shang, H. & Chiang, T.-Y. Ejection history of the IRAS 04166+2706 molecular jet. Astrophys. J. 874, 31 (2019). Shang, H., Liu, C.-F., Krasnopolsky, R. & Wang, L.-Y. A unified model for bipolar outflows from young stars: kinematic signatures of jets, winds, and their magnetic interplay with the ambient toroids. Astrophys. J. 944, 230 (2023). Pudritz, R. E. & Norman, C. A. Centrifugally driven winds from contracting molecular disks. Astrophys. J. 274, 677–697 (1983). Yvart, W., Cabrit, S., Pineau des Forêts, G. & Ferreira, J. Molecule survival in magnetized protostellar disk winds. II. Predicted H2O line profiles versus Herschel/HIFI observations. Astron. Astrophys. 585, A74 (2016). de Valon, A. et al. Modeling the CO outflow in DG Tauri B: swept-up shells versus perturbed MHD disk wind. Astron. Astrophys. 668, A78 (2022). Fernández-López, M. et al. A ringed pole-on outflow from DO Tauri revealed by ALMA. Astron. J. 159, 171 (2020). Zhang, Y. et al. An episodic wide-angle outflow in HH 46/47. Astrophys. J. 883, 1 (2019). Liu, C.-F. et al. ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): nested morphological and kinematic structures of outflows revealed in SiO and CO emission. Astrophys. J. 979, 17 (2025). Lee, C.-F., Li, Z.-Y., Shang, H. & Hirano, N. Magnetocentrifugal origin for protostellar jets validated through detection of radial flow at the jet base. Astrophys. J. Lett. 927, L27 (2022). Louvet, F. et al. The HH30 edge-on T Tauri star. A rotating and precessing monopolar outflow scrutinized by ALMA. Astron. Astrophys. 618, A120 (2018). Strom, S. E., Vrba, F. J. & Strom, K. M. Infrared surveys of dark-cloud complexes. II. The NGC 1333 region. Astron. J. 81, 314–316 (1976). Bachiller, R. & Cernicharo, J. Extremely high-velocity emission from molecular jets in NGC 6334I and NGC 1333 (HH 7-11). Astron. Astrophys. 239, 276–286 (1990). Masson, C. R., Mundy, L. G. & Keene, J. The extremely high velocity CO flow in HH 7–11. Astrophys. J. Lett. 357, L25–L28 (1990). Bachiller, R., Gueth, F., Guilloteau, S., Tafalla, M. & Dutrey, A. The origin of the HH 7-11 outflow. Astron. Astrophys. 362, L33–L36 (2000). Chen, X., Arce, H. G., Zhang, Q., Launhardt, R. & Henning, T. Rotating bullets from a variable protostar. Astrophys. J. 824, 72 (2016). Lefèvre, C. et al. CALYPSO view of SVS 13A with PdBI: multiple jet sources. Astron. Astrophys. 604, L1 (2017). Ortiz-León, G. N. et al. The Gould’s Belt Distances Survey (GOBELINS). V. Distances and kinematics of the Perseus molecular cloud. Astrophys. J. 865, 73 (2018). Gaia Collaborationet al. Gaia Data Release 3. Summary of the content and survey properties. Astron. Astrophys. 674, A1 (2023). Anglada, G., Rodríguez, L. F. & Torrelles, J. M. Discovery of a subarcsecond radio binary associated with the SVS 13 star in the HH 7-11 region. Astrophys. J. Lett. 542, L123–L126 (2000). Anglada, G. et al. A single circumstellar disk in the SVS 13 close binary system. Astrophys. J. Lett. 605, L137–L140 (2004). Diaz-Rodriguez, A. K. et al. The physical properties of the SVS 13 protobinary system: two circumstellar disks and a spiraling circumbinary disk in the making. Astrophys. J. 930, 91 (2022). Eisloeffel, J. et al. An optical/infrared outburst in the exciting source of HH 7-11. Astrophys. J. Lett. 383, L19–L22 (1991). Hodapp, K. W. & Chini, R. The launch region of the SVS 13 outflow and jet. Astrophys. J. 794, 169 (2014). Hartigan, P., Holcomb, R. & Frank, A. Proper motions and shock wave dynamics in the HH 7-11 stellar jet. Astrophys. J. 876, 147 (2019). Liseau, R., Sandell, G. & Knee, L. B. G. The structure of the molecular outflow near SSV 13 and HH 7-11 in the NGC 1333 region. Astron. Astrophys. 192, 153–164 (1988). Plunkett, A. L. et al. CARMA observations of protostellar outflows in NGC 1333. Astrophys. J. 774, 22 (2013). Gardner, C. L., Jones, J. R. & Hodapp, K. W. Numerical simulation of the SVS 13 microjet and bow shock bubble. Astrophys. J. 830, 113 (2016). Cieza, L. A. et al. Imaging the water snow-line during a protostellar outburst. Nature 535, 258–261 (2016). Houge, A. & Krijt, S. Collisional evolution of dust and water ice in protoplanetary discs during and after an accretion outburst. Mon. Not. R. Astron. Soc. 521, 5826–5845 (2023). Simon, M. N. et al. Tracing slow winds from T Tauri stars via low-velocity forbidden line emission. Astrophys. J. 831, 169 (2016). Hartigan, P., Morse, J. A., Heathcote, S. & Cecil, G. Observations of entrainment and time variability in the HH 47 jet. Astrophys. J. Lett. 414, L121–L124 (1993). Birney, M. et al. A kinematical study of the launching region of the blueshifted HH 46/47 outflow with SINFONI K-band observations. Astron. Astrophys. 692, A143 (2024). Jhan, K.-S. et al. ALMA Survey of Orion Planck Galactic Cold Clumps (ALMASOP): deriving inclination angle and velocity of the protostellar jets from their SiO knots. Astrophys. J. Lett. 931, L5 (2022). Fischer, W. J. et al. in Protostars and Planets VII (eds Inutsuka, S. et al.) Ch. 10 (ASP, 2023). Masson, C. R. & Chernin, L. M. Properties of jet-driven molecular outflows. Astrophys. J. 414, 230–241 (1993). Ostriker, E. C., Lee, C.-F., Stone, J. M. & Mundy, L. G. A ballistic bow shock model for jet-driven protostellar outflow shells. Astrophys. J. 557, 443–450 (2001). Tabone, B., Raga, A., Cabrit, S. & Pineau des Forêts, G. Interaction between a pulsating jet and a surrounding disk wind. A hydrodynamical perspective. Astron. Astrophys. 614, A119 (2018). Sperling, T. et al. Probing the hidden atomic gas in Class I jets with SOFIA. Astron. Astrophys. 642, A216 (2020). Cabrit, S. in Star-Disk Interaction in Young Stars, Vol. 243, (eds Bouvier, J. & Appenzeller, I.) Ch. 21 (Cambridge Univ. Press, 2007). Königl, A., Romanova, M. M. & Lovelace, R. V. E. Are the outflows in FU Orionis systems driven by the stellar magnetic field? Mon. Not. R. Astron. Soc. 416, 757–766 (2011). Ferreira, J., Pelletier, G. & Appl, S. Reconnection X-winds: spin-down of low-mass protostars. Mon. Not. R. Astron. Soc. 312, 387–397 (2000). Zanni, C. & Ferreira, J. MHD simulations of accretion onto a dipolar magnetosphere. II. Magnetospheric ejections and stellar spin-down. Astron. Astrophys. 550, A99 (2013). de Valon, A. et al. ALMA reveals a large structured disk and nested rotating outflows in DG Tauri B. Astron. Astrophys. 634, L12 (2020). Matzner, C. D. & McKee, C. F. Bipolar molecular outflows driven by hydromagnetic protostellar winds. Astrophys. J. Lett. 526, L109–L112 (1999). Federrath, C., Schrön, M., Banerjee, R. & Klessen, R. S. Modeling jet and outflow feedback during star cluster formation. Astrophys. J. 790, 128 (2014). Rohde, P. F. et al. Evolution of Hubble wedges in episodic protostellar outflows. Mon. Not. R. Astron. Soc. 483, 2563–2580 (2019). Offner, S. S. R. & Arce, H. G. Investigations of protostellar outflow launching and gas entrainment: hydrodynamic simulations and molecular emission. Astrophys. J. 784, 61 (2014). Raga, A. C., Canto, J., Binette, L. & Calvet, N. Stellar jets with intrinsically variable sources. Astrophys. J. 364, 601–610 (1990). Virtanen, P. et al. SciPy 1.0: fundamental algorithms for scientific computing in Python. Nat. Methods 17, 261–272 (2020). Romero-Muñoz, L., Garcia-Villanueva, M. & Gómez-Suárez, C. A tutorial on the total least squares method for fitting a straight line and a plane. Rev. Cienc. Ingen. Inst. Technol. Super. Coatzacoalcos 1, 167–173 (2014). Reynolds, N. K. et al. The disk orientations of Perseus protostellar multiples at ~8 au resolution. Astrophys. J. 963, 164 (2024). Estalella, R. & Anglada, G. Introducción a la Física del Medio Interestelar (Publicacions i Edicions de la Universitat de Barcelona, 2008). Frerking, M. A., Langer, W. D. & Wilson, R. W. The relationship between carbon monoxide abundance and visual extinction in interstellar clouds. Astrophys. J. 262, 590–605 (1982). Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, N. J. & Lee, C. W. MAMBO mapping of Spitzer c2d small clouds and cores. Astron. Astrophys. 487, 993–1017 (2008). Palla, F. & Stahler, S. W. The pre-main-sequence evolution of intermediate-mass stars. Astrophys. J. 418, 414–425 (1993). Enoch, M. L., Evans, N. J., Sargent, A. I. & Glenn, J. Properties of the youngest protostars in Perseus, Serpens, and Ophiuchus. Astrophys. J. 692, 973–997 (2009). Osorio, M. et al. Imaging the inner and outer gaps of the pre-transitional disk of HD 169142 at 7 mm. Astrophys. J. Lett. 791, L36 (2014). Blazquez-Calero, G. et al. CO(J=3–2) spectral cubes of the reduced ALMA observations of SVS13 outflow. Zenodo https://doi.org/10.5281/zenodo.17249853 (2025). Harris, C. R. et al. Array programming with NumPy. Nature 585, 357–362 (2020). Hunter, J. D. Matplotlib: a 2D graphics environment. Comput. Sci. Eng. 9, 90–95 (2007). Astropy Collaboration. et al. Astropy: a community Python package for astronomy. Astron. Astrophys. 558, A33 (2013). Astropy Collaboration. et al. The Astropy Project: building an open-science project and status of the v2.0 core package. Astron. J. 156, 123 (2018). Astropy Collaboration. et al. The Astropy Project: sustaining and growing a community-oriented open-source project and the latest major release (v5.0) of the core package. Astrophys. J. 935, 167 (2022). Newville, M., Stensitzki, T., Allen, D. B. & Ingargiola, A. Lmfit: non-linear least-square minimization and curve-fitting for Python. Zenodo https://doi.org/10.5281/zenodo.11813 (2014). Blazquez-Calero, G. et al. BowshockPy: a Python package for the generation of synthetic spectral channel maps of a jet-driven bowshock model. GitHub https://github.com/gblazquez/bowshockpy (2025). Stephens, I. W. et al. Alignment between protostellar outflows and filamentary structure. Astrophys. J. 846, 16 (2017). Rodríguez, L. F., Anglada, G. & Curiel, S. The nature of the radio continuum sources embedded in the HH 7-11 region and its surroundings. Astrophys. J. Suppl. 125, 427–438 (1999). Khanzadyan, T. et al. A multi-epoch near-infrared study of the HH 7-11 protostellar outflow. Mon. Not. R. Astron. Soc. 338, 57–66 (2003). We dedicate this work to the memory of our dear colleagues and friends, Alejandro (Alex) Raga and Robert Estalella. Their passion and brilliance were central to the theoretical foundations of this research and played a crucial role in shaping the Article as it stands today. This manuscript owes much of its depth and rigour to their insight, dedication and collaboration. Although they did not live to see its publication, their spirit and intellectual legacy live on in every page. We remain deeply grateful for the time we shared and the science we built together. We thank B. Tabone for his useful comments and suggestions on this work. G.A., G.B.-C., I.d.G.-M., A.K.D.-R., G.A.F., J.F.G. and M.O. acknowledge financial support from grant nos. PID2020-114461GB-I00, PID2023-146295NB-I00 and CEX2021-001131-S, funded by MCIN/AEI/10.13039/501100011033. G.B.-C., G.A.F. and M.O. acknowledge financial support from Junta de Andalucia (Spain) grant no. P20-00880 (FEDER, EU). G.B-C acknowledges support from grant no. PRE2018-086111, funded by MCIN/AEI/ 10.13039/501100011033 and by ‘ESF Investing in your future’ and thanks ESO Science Support Discretionary Fund for their financial support under the 2024 SSDF 06 project. S.C. acknowledges support from Conseil Scientifique of Observatoire de Paris and from the Programme National de Physique et Chimie du Milieu Interstellaire (PCMI) of CNRS/INSU co-funded by CEA and CNES. A.-K.D.R acknowledges support from STFC grant no. ST/T001488/1. G.A.F. also acknowledges support from the Collaborative Research Centre 956, funded by the Deutsche Forschungsgemeinschaft (DFG) project ID 184018867, the DFG for funding through SFB 1601 ‘Habitats of massive stars across cosmic time’ (sub-project B1) and the University of Cologne and its Global Faculty programme. R.E. acknowledges partial financial support from grant nos. PID2020-117710GB-I00 and CEX2019-000918-M funded by MCIN/ AEI /10.13039/501100011033. J.M.T. acknowledges support from grant no. PID2023-146675NB funded by MCIN/AEI/10.13039/501100011033 and by the programme Unidad de Excelencia María de Maeztu CEX2020-001058-M. L.F.R. acknowledges support from grant no. CBF-2025-I-2471 of SECIHTI, Mexico. L.A.Z. acknowledges financial support from grant nos. CONACyT-280775, UNAM-PAPIIT IN110618 and IN112323, México. This work makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.01229.S, ADS/JAO.ALMA#2016.1.01305.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan) and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This publication use images based on observations made with the NASA/ESA Hubble Space Telescope and obtained from the Hubble Legacy Archive, which is a collaboration between the Space Telescope Science Institute (STScI/NASA), the Space Telescope European Coordinating Facility (ST-ECF/ESA) and the Canadian Astronomy Data Centre (CADC/NRC/CSA). G.B.-C. led the data reduction, analysis, modelling, interpretation and writing of the paper. G.A., S.C., M.O., A.C.R., G.A.F. and R.E. contributed important inputs to the methodology, modelling and interpretation of the results and to the writing of paper sections. G.A., J.F.G. and A.K.D.-R. contributed to the data reduction. G.A. led the ALMA observation proposal, conceived and prepared with contributions from M.O., J.F.G., A.K.D.-R., J.M.T., L.F.R., E.M., I.d.G.-M. and P.T.P.H. All authors participated in discussions of the results and contributed to the paper preparation and revision. Peer review information Nature Astronomy thanks Somnath Dutta, Klaus Hodapp and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Extended data Observed outflow features surrounding SVS 13 at different scales (see their observational properties in Supplementary Table 1). Large boxes indicate zoomed-in regions. The synthesized beams are shown in the bottom right corner of the panels. Velocities are relative to VLA 4B (vLSR = + 9.3 km/s)58. Top left: CARMA CO(J=1-0) map of the large-scale SHV molecular outflow63,105. The emission has been integrated from − 11.6 to − 3.1 km/s (blue lobe) and from 0.7 to 8 km/s (red lobe). The location of the EHV molecular bullets51,52 (triangles), Herbig-Haro objects (squares), and radio sources proposed as YSOs58,106 (plus signs, where SVS 13 is the largest) is indicated. Top right: HST optical image (F606W filter) of the region, showing the outflow cavity, with the SMA map of the CO(2-1) emission of the EHV bullets52 (with velocities from − 161 to +129 km/s) overlapped in blue contours. Bottom left: Our ALMA CO(3-2) map (not corrected for the primary beam response) of the central region of the molecular outflow (beam=0.527″ × 0.333″, PA=2.7 deg) is shown in contours. The emission has been integrated from − 126.4 to − 6.5 km/s (blue lobe) and from 6.5 to 60.8 km/s (red lobe). Contours are 3, 6, 10, 15, 21, 28, and 37 times 0.41 Jy/beam km/s (blue lobe) and 0.32 Jy/beam km/s (red lobe).The H2 arcuate features60 are plotted as arcs. Bottom right: Close-up toward Bullet 1. Our ALMA CO(3-2) map of the blueshifted emission (beam=0.537″ × 0.333″, PA=2.7 deg), integrated from − 9.3 to − 126.4 km/s (thus including the SHV, IHV, and EHV components), is shown in contours. Contours are 3, 6, 10, 16, and 25 times 0.61 Jy/beam km/s. The colored filled areas represent the heads of the detected families of rings (see Fig. 3). The green squares represent the location of two CO clumps (where c2 coincides in position and velocity with the SiO knot reported by ref. 53). Image of the velocity-integrated intensity of the SO(88 − 77) line observed with the ALMA 12-m array using a low spectral resolution spectral window dedicated to continuum observation, with a channel spacing of 13.6 km/s. The emission has been integrated in the LOS velocity range from -9.8 to -118.6 km/s relative to the velocity of VLA 4B (+9.3 km/s)58. The positions of the two protostars of the SVS 13 binary56 are indicated by plus signs. The synthesized beam, shown in the bottom left corner, is 0.54″ × 0.36″ (PA=8.37 deg). Contours are -3, 3, 5, 8, 13, 20, 30, 50, 80, 140, 260 times 0.04 Jy/beam km/s. The image has not been corrected by the primary beam response. A sample of spectral channel images of the CO(J=3-2) emission observed by ALMA with a synthesized beam of 0.554″ × 0.352″ (PA=3.1 deg), where natural weighting has been used. The obtained data cover a range of velocities wider than the high angular resolution data shown in Fig. 2. The positions of the two protostars of the SVS 13 binary56 are indicated by plus signs. The LOS velocity, relative to the velocity of VLA 4B (vLSR = + 9.3 km/s)58, is shown in the top left corner of each image. The width of each of the spectral channels shown in the figure is 2.12 k/ms, which corresponds to the average of 10 native channels. The rms of the images is 8 mJy/beam. Images have not been corrected by the primary beam response. The synthesized beam is plotted as an ellipse in the bottom left corner of the first image. Position angle as a function of the dynamical time for the heads of the families of rings (plotted as triangles and labeled F I to F VI; see Methods), as well as for other features associated with the SVS 13 outflow (we plot atomic features as circles, and molecular as rectangles), such as HH objects61,107 (plotted as squares and labeled HH7, HH11), [FeII] jet (triangle labeled as Microjet), and H2 arcs60 (squares labeled as HC1, HC2, HC3). Colors indicate features with position angles around 160 deg (in blue) and in the range of 120-140 deg (in gray). The data for the families of rings (F I-VI) are presented as mean values, and error bars represent the standard deviation (see ‘Dynamical times of the heads of the families of rings’ in Methods). For the rest of the objects we represent the values and uncertainties as reported in the literature52,60,61,107. We assumed inclination angles from 22 to 25 deg. See Methods and Supplementary Table 1 for details on the calculations. Epoch is given in Julian years, with epoch 2000 corresponding to the standard definition of Julian epoch J2000.0. Dynamical times correspond to the date of our high-angular resolution ALMA observations (epoch 2016.69). The bowshock is seen in a reference system moving at the velocity vjet of the working surface. The cylindrical jet beam has a diameter 2rjet, and the impinging ambient gas moves to the left at a velocity vjet − vamb in this reference system. We show a cylindrical coordinate system (x*,r), where r is the cylindrical radius and x* the distance measured from the head of the working surface towards the outflow source. The working surface ejects material sideways at a velocity v0 (which is approximately equal to the post-cooling region sound speed of ~ 10 km/s). This sideways ejection interacts with the impinging ambient gas, forming a thin shell bowshock that has a well defined locus rb(x*), and locally has a slope \(\tan (\alpha )={\rm{dr}}_{\rm{b}}/{{\rm{dx}}}^{*}\). Top left: Geometry of the bowshock shell, where black vectors show velocities of the shell, blue vectors their LOS projections, and the shaded band outlines the spatial extent of the emitting region in the velocity range of a given channel map. Top center: Sketch of the corresponding channel image. Due to projection effects, the observed ring emission looks wider on the side closest to the star, while the far side appears narrower. Top right: Observed channel image (rotated so that the top faces SVS 13) illustrating the asymmetry of the rings, with the side closest in projection to the star appearing brighter and wider than the opposite side. Middle row: At low optical depths, in the side closest to the star, the emission appears spread over a wider range of radii and with lower intensity than on the far side, where it appears narrower but with higher intensity (left). When convolved with an elongated beam (center), the intensity becomes almost uniform throughout the ring, except near the positions where the beam axis is tangent to the ring (a well-known beam-filling-factor effect95). When noise is added (right), no substantial asymmetries are detected. Bottom row: When the optical depth is high enough (left), the side closer to the star becomes almost as bright, but more extended, than the far side, where any increase in intensity is hindered by opacity saturation. When convolved with the beam (center), the intensity drops in the far (narrower) side of the ring, because of the smaller beam-filling factor, producing an asymmetry similar to that observed. A model image with noise (right panel) shows that a qualitative agreement with observations (top panel) can be achieved. In the modeling, the assumed velocity dispersion is 2 km/s, the inclination angle is 20 deg, and the channel width is 0.53 km/s. Decomposition of the observed CO (J = 3–2) channel map emission of Bullet 1 in SVS 13 into elliptical rings. The LOS velocity, relative to VLA 4B (vLSR = +9.3 km s−1 (ref. 58)), ranges from −0.9 km s−1 to −102.5 km s−1. Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law. https://bit.ly/3XU2i7k