This White Paper presents the scientific rationale and instrument concept for HRMOS (High-Resolution Multi-Object Spectrograph), a next-generation instrument proposed for the ESO Very Large Telescope within the VLT 2030 roadmap. Current and planned facilities offer either multi-object spectroscopy or ultra-high spectral resolution, but not both. HRMOS fills this gap by combining very high spectral resolution, multi-object capability, and radial-velocity stability, enabling transformative studies in Galactic and extragalactic astrophysics. The baseline design provides a resolving power of R = 80000, radial-velocity precision of 10 m s-1 (goal: 5 m s-1), simultaneous observations of 50-60 targets, and broad optical coverage down to 385 nm. These capabilities enable precise measurements of elemental abundances, isotopic ratios, line profiles, and radial velocities for large stellar samples, including crowded fields, star clusters, the Galactic bulge, and nearby dwarf galaxies. HRMOS will address key questions on the age of the oldest stellar populations through nucleocosmochronology, the formation and survival of planetary systems, the assembly history of the Milky Way and satellites, the origin of the heaviest elements, stellar evolution, and the chemical and dynamical properties of the interstellar and circumgalactic medium. It will bridge large spectroscopic surveys and the next generation of extremely large telescopes, with strong synergies with 4MOST, Gaia, TESS, PLATO, the proposed Haydn mission, and future ELT instruments. Building on VLT/FLAMES heritage, HRMOS represents a strategic investment for European astronomy in the 2030s.

Context. The typically large distances, extinction, and crowding of Galactic supermassive star clusters (stellar clusters more massive than 10(4) M-circle dot) have so far hampered the identification of their very low mass members, required to extend our understanding of star and planet formation, and early stellar evolution, to the extremely energetic star-forming environment typical of starbursts. This situation has now evolved thanks to the James Webb Space Telescope (JWST), and its unmatched resolution and sensitivity in the infrared. Aims. In this paper, the third of the series of the Extended Westerlund 1 and 2 Open Clusters Survey (EWOCS), we present JWST/NIRCam and JWST/MIRI observations of the supermassive star cluster Westerlund 1. These observations are specifically designed to unveil the cluster members down to the brown dwarf mass regime, and to allow us to select and study the protoplane-tary disks in the cluster and to study the mutual feedback between the cluster members and the surrounding environment. Methods. Westerlund 1 was observed as part of JWST GO-1905 for 23.6 hours. The data have been reduced using the JWST calibration pipeline, together with specific tools necessary to remove artifacts, such as the 1/f random noise in NIRCam images. Source identification and photometry were performed with DOLPHOT. Results. The MIRI images show a plethora of different features. Diffuse nebular emission is observed around the cluster, which is typically composed of myriads of droplet-like features pointing toward the cluster center or the group of massive stars surrounding the Wolf-Rayet star W72/A. A long pillar is also observed in the northwest. The MIRI images also show resolved shells and outflows surrounding the M-type supergiants W20, W26, W75, and W237, the sgB[e] star W9 and the yellow hypergiant W4. Some of these shells have been observed before at other wavelengths, but never with the level of detail provided by JWST. The color-magnitude diagrams built using the NIRCam photometry show a clear cluster sequence, which is marked in its upper part by the 1828 NIRCam stars with X-ray counterparts. NIRCam observations using the F115W filter have reached the 23.8 mag limit with 50% completeness (roughly corresponding to a 0.06 M-circle dot brown dwarf).

OB stars are crucial for our understanding of Galactic structure, star formation, stellar feedback and multiplicity. In this paper we have compiled a census of all OB stars within 1 kpc of the Sun. We performed evolutionary and atmospheric model fits to observed SEDs compiled from astro-photometric survey data. We have characterized and mapped 24,706 O- and B-type stars (T-eff > 10,000 K) within 1 kpc of the Sun, whose overdensities correspond to well-studied OB associations and massive star-forming regions such as Sco-Cen, Orion OB1, Vela OB2, Cepheus and Circinus. We have assessed the quality of our catalogue by comparing it with spectroscopic samples and similar catalogues of OB(A) stars, as well as catalogues of OB associations, star-forming regions and young open clusters. Finally, we have also exploited our list of OB stars to estimate their scale height (76 +/- 1 pc), a local star formation rate of 2896(-1)(+417) M circle dot Myr(-1) and a local core-collapse supernova rate of similar to 15-30 per Myr. We extrapolate these rates to the entire Milky Way to derive a Galactic SFR of 0.67(-0.01)(+0.09 )M(circle dot) yr(-1) and a core-collapse supernova rate of 0.4-0.5 per century. These are slightly lower than previous estimates, which we attribute to improvements in our census of OB stars and changes to evolutionary models. We calculate a near-Earth core collapse supernova rate of similar to 2.5 per Gyr that supports the view that a nearby supernova explosion could have caused one of the recorded mass extinction events on Earth.

We present results from a large spectroscopic survey of the nearest young association to the Sun, Upper Scorpius, conducted using 2dF/HERMES on the Anglo-Australian Telescope. We use spectroscopic youth criteria such as Li-equivalent widths to identify >1000 pre-main sequence (PMS) members across the region and measure radial velocities, combining these with Gaia EDR3 5-parameter astrometry to obtain 6D kinematic information. We separate confirmed PMS association members into distinct kinematic groups and measure expansion and rotation trends in each. We find strong evidence for asymmetric expansion in several groups and derive expansion timescales from the greatest rates of expansion in each group. We also trace the past motion of these groups using an epicycle approximation and estimate the time since their most compact configuration. These kinematic properties are compared to literature ages and the star formation history of Upper Scorpius is discussed. We find evidence that a scenario in which star formation in the subgroups of Upper Scorpius proceeded independently, either by self-instability or external feedback from Upper Centaurus-Lupus, is more likely than a recently proposed "cluster chain" scenario in which these subgroups have triggered each other.

NGC 6302 is a spectacular bipolar planetary nebula (PN) whose spectrum exhibits fast outflows and highly ionized emission lines, indicating the presence of a very hot central star (similar to 220 000 K). Its infrared spectrum reveals a mixed oxygen and carbon dust chemistry, displaying both silicate and polycyclic aromatic hydrocarbon (PAH) features. Using the James Webb Space Telescope Mid-Infrared Instrument and Medium Resolution Spectrometer, a mosaic map was obtained over the core of NGC 6302, covering the wavelength range of 5-28 mu m and spanning an area of similar to 18.5 arcsec x 15arcsec. The spatially resolved spectrum reveals similar to 200 molecular and ionized lines from species requiring ionization potentials of up to 205 eV. The spatial distributions highlight a complex structure at the nebula's centre. Highly ionized species such as [Mg VII ] and [Si VII ] show compact structures, while lower ionization species such as H+ extend much farther outwards, forming filament-defined rims that delineate a bubble. Within the bubble, the H+ and H-2 emission coincide, while the PAH emission appears farther out, indicating an ionization structure distinct from typical photodissociation regions, such as the Orion Bar. This may be the first identification of a PAH formation site in a PN. This PN appears to be shaped not by a steady, continuous outflow, but by a series of dynamic, impulsive bubble ejections, creating local conditions conducive to PAH formation. A dusty torus surrounds the core, primarily composed of large (mu m-sized) silicate grains with crystalline components. The long-lived torus contains a substantial mass of material, which could support an equilibrium chemistry and a slow dust-formation process.

We present the first large-scale 3D kinematic study of ~2000 spectroscopically-confirmed young stars (<20 Myr) in 18 star clusters and OB associations (hereafter groups) from the combination of Gaia astrometry and Gaia-ESO Survey spectroscopy. We measure 3D velocity dispersions for all groups, which range from 0.61 to 7.4 km/s (1D velocity dispersions of 0.35 to 4.3 km/s). We find the majority of groups have anisotropic velocity dispersions, suggesting they are not dynamically relaxed. From the 3D velocity dispersions, measured radii and estimates of total mass we estimate the virial state and find that all systems are super-virial when only the stellar mass is considered, but that some systems are sub-virial when the mass of the molecular cloud is taken into account. We observe an approximately linear correlation between the 3D velocity dispersion and the group mass, which would imply that the virial state of groups scales as the square root of the group mass. However, we do not observe a strong correlation between virial state and group mass. In agreement with their virial state we find that nearly all of the groups studied are in the process of expanding and that the expansion is anisotropic, implying that groups were not spherical prior to expansion. One group, Rho Oph, is found to be contracting and in a sub-virial state (when the mass of the surrounding molecular cloud is considered). This work provides a glimpse of the potential of the combination of Gaia and data from the next generation of spectroscopic surveys.

Context. With a mass exceeding several 10^4 solar masses and a rich and dense population of massive stars, supermassive young star clusters represent the most massive star-forming environment that is dominated by the feedback from massive stars and gravitational interactions among stars. Aims. In this paper we present the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, which aims to investigate the influence of the starburst environment on the formation of stars and planets, and on the evolution of both low and high mass stars. The primary targets of this project are Westerlund 1 and 2, the closest supermassive star clusters to the Sun. Methods. The project is based primarily on recent observations conducted with the Chandra and JWST observatories. Specifically, the Chandra survey of Westerlund 1 consists of 36 new ACIS-I observations, nearly co-pointed, for a total exposure time of 1 Msec. Additionally, we included 8 archival Chandra/ACIS-S observations. This paper presents the resulting catalog of X-ray sources within and around Westerlund 1. Sources were detected by combining various existing methods, and photon extraction and source validation were carried out using the ACIS-Extract software. Results. The EWOCS X-ray catalog comprises 5963 validated sources out of the 9420 initially provided to ACIS-Extract, reaching a photon flux threshold of approximately 2x10^-8 photons/cm^2/s. The X-ray sources exhibit a highly concentrated spatial distribution, with 1075 sources located within the central 1 arcminute. We have successfully detected X-ray emissions from 126 out of the 166 known massive stars of the cluster, and we have collected over 71000 photons from the magnetar CXO J164710.20-455217

Agricultural ponds are a significant source of greenhouse gases, contributing to the ongoing challenge of anthropogenic climate change. Nations are encouraged to account for these emissions in their national greenhouse gas inventory reports. We present a remote sensing approach using open-access satellite imagery to estimate total methane emissions from agricultural ponds that account for (1) monthly fluctuations in the surface area of individual ponds, (2) rates of historical accumulation of agricultural ponds, and (3) the temperature dependence of methane emissions. As a case study, we used this method to inform the 2024 National Greenhouse Gas Inventory reports submitted by the Australian government, in compliance with the Paris Agreement. Total annual methane emissions increased by 58% from 1990 (26 kilotons CH4 year-1) to 2022 (41 kilotons CH4 year-1). This increase is linked to the water surface of agricultural ponds growing by 51% between 1990 (115 kilo hectares; 1,150 km2) and 2022 (173 kilo hectares; 1,730 km2). In Australia, 16,000 new agricultural ponds are built annually, expanding methane-emitting water surfaces by 1,230 ha yearly (12.3 km2 year-1). On average, the methane flux of agricultural ponds in Australia is 0.238 t CH4 ha-1 year-1. These results offer policymakers insights into developing targeted mitigation strategies to curb these specific forms of anthropogenic emissions. For instance, financial incentives, such as carbon or biodiversity credits, can mobilize widespread investments toward reducing greenhouse gas emissions and enhancing the ecological and environmental values of agricultural ponds. Our data and modeling tools are available on a free cloud-based platform for other countries to adopt this approach.

The spatial-kinematic structure of 48 young star clusters and associations is investigated. Moran's I statistic is used to quantify the degree of kinematic substructure in each region, and the results are compared to those expected assuming the hierarchical or monolithic models of star cluster formation. Of the observed regions, 39 are found to have significant kinematic substructure, such that they are compatible with the hierarchical model and incompatible with the monolithic model. This includes multiple regions whose Q parameter shows the region to be centrally concentrated and clustered. The remaining nine are compatible with both models. From this, it is concluded that the kinematic substructure of the observed star clusters represents strong evidence in favour of the hierarchical model of star cluster formation over the monolithic model.

Context. Wolf-Rayet (WR) stars are massive evolved stars that exhibit particularly fast and dense stellar winds. Although they constitute a very short phase near the end of a massive star's life, they play a crucial role in the evolution of massive stars and have a substantial impact on their surrounding environment. Aims. We present the most comprehensive and deepest X-ray study to date of the properties of the richest Wolf-Rayet population observed in a single stellar cluster, Westerlund 1 (Wd1). By examining the X-ray signatures of WR stars, we aim to shed light on the hottest plasma in their stellar winds and gain insights into whether they exist as single stars or within binary systems. Methods. This work is based on 36 Chandra observations obtained from the "Extended Westerlund 1 and 2 Open Clusters Survey" (EWOCS) project, plus 8 archival Chandra observations. The overall exposure depth Ms) and baseline of the EWOCS observations extending over more than one year enable us to perform a detailed photometric, colour, and spectral analysis, as well as to search for short- and long-term periodicity. Results. In X-rays, we detect 20 out of the 24 known Wolf-Rayet stars in Wd1 down to an observed luminosity of similar to 7 x 10(29) erg s(-1) (assuming a distance of 4.23 kpc to Wd1), with 8 WR stars being detected in X-rays for the first time. Nine stars show clear evidence of variability over the year-long baseline, with clear signs of periodicity. The X-ray colours and spectral analysis reveal that the vast majority of the WR stars are hard X-ray sources (kT >= 2.0 keV). The Fe XXV emission line at similar to 6.7 keV, which commonly originates from the wind-wind collision zone in binary systems, is detected for the first time in the spectra of 17 WR stars in Wd1. In addition the similar to 6.4 keV fluorescent line is observed in the spectra of three stars, which are among the very few massive stars exhibiting this line, indicating that dense cold material coexists with the hot gas in these systems. Overall, our X-ray results alone suggest a very high binary fraction (>= 80%) for the WR star population in Wd1. When combining our results with properties of the WR population from other wavelengths, we estimate a binary fraction of >= 92%, which could even reach unity. This suggests that either all the most massive stars are found in binary systems within Wd1, or that binarity is essential for the formation of such a rich population of WR stars.

Context. Star clusters, composed of stars born from the same molecular cloud, serve as invaluable natural laboratories for understanding the fundamental processes governing stellar formation and evolution. Aims. This study aims to investigate correlations between the Mean Interdistance (), Mean Closest Interdistance () and Median Weighted Central Interdistance () with the age of star clusters, examining their evolutionary trends and assessing the robustness of these quantities as possible age indicators. Methods. We selected a sample of open clusters in the solar region and with a representative number of members (e.g. well populated and without outliers). The interdistances are derived from the spatial distribution of member stars within a cluster. Their evolution over time allows us to use them as age indicators for star clusters. Results. Our investigation reveals a high-significant correlation between the interdistances and cluster age. Considering the full sample of clusters between 7 and 9 kpc, the relationship is very broad. This is due to uncertainties in parallax, which increase with increasing distance. In particular, we must limit the sample to a maximum distance from the Sun of about 200 pc to avoid artificial effects on cluster shape and on the spatial distribution of their stars along the line of sight. Conclusions. By conservatively restraining the distance to a maximum of similar to 200 pc, we have established a relationship between the interdistances and the age of the clusters. In our sample, the relationship is mainly driven by the internal expansion of the clusters and is marginally affected by external perturbative effects. Such relation might enhance our comprehension of cluster dynamics and might be used to derive cluster dynamical ages.

The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate

ABSTRACT WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, saw first light in late 2022. WEAVE comprises a new 2-deg field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable ‘mini’ integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366–959 nm at R ∼ 5000, or two shorter ranges at $R\sim 20\, 000$. After summarizing the design and implementation of WEAVE and its data systems, we present the organization, science drivers, and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy’s origins by completing Gaia’s phase-space information, providing metallicities to its limiting magnitude for ∼3 million stars and detailed abundances for ∼1.5 million brighter field and open-cluster stars; (ii) survey ∼0.4 million Galactic-plane OBA stars, young stellar objects, and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey ∼400 neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionized gas in z < 0.5 cluster galaxies; (vi) survey stellar populations and kinematics in ${\sim} 25\, 000$ field galaxies at 0.3 ≲ z ≲ 0.7; (vii) study the cosmic evolution of accretion and star formation using >1 million spectra of LOFAR-selected radio sources; and (viii) trace structures using intergalactic/circumgalactic gas at z > 2. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator.

The zero-age main sequence (ZAMS) is a critical phase for stellar angular momentum evolution, as stars transition from contraction-dominated spin-up to magnetic wind-dominated spin-down. We present the first robust observational constraints on rotation for FGK stars at ≈40 Myr. We have analyzed TESS light curves for 1410 members of five young open clusters with ages between 25 and 55 Myr: IC 2391, IC 2602, NGC 2451A, NGC 2547, and Collinder 135. In total, we measure 868 rotation periods, including 96 new, high-quality periods for stars around 1 M _⊙ . This is an increase of ten times the existing literature sample at the ZAMS. We then use the τ ^2 method to compare our data to models for stellar angular momentum evolution. Although the ages derived from these rotation models do not match isochronal ages, we show that these observations can clearly discriminate between different models for stellar wind torques. Finally, τ ^2 fits indicate that magnetic braking and/or internal angular momentum transport significantly impact rotational evolution even on the pre-main sequence.

Artificial waterbodies can help tackle the ongoing freshwater biodiversity crisis by providing new habitat for aquatic wildlife. Farm dams are among the most abundant artificial waterbodies in agricultural landscapes, yet general guidelines for maximising their ecological potential are yet to emerge. Here we used a continental-scale dataset of frog species near 8,800 Australian farm dams to ask: What characteristics promote higher frog species richness at farm dams? The highest values of species richness were at older (>20 years) farm dams of intermediate sizes (0.1 ha in surface area), with smaller rainfall catchments (<10 ha), and near other freshwater systems or conservation sites. By identifying quantifiable features improving the ecological value of farm dams, this work helps identify “win-win” outcomes for agricultural productivity and conservation. In the future, “biodiversity credit” policies could promote large-scale implementation by rewarding farmers who invest in improving the condition of their farm dams.

The Cygnus OB2 association is the largest concentration of young and massive stars within 2 kpc of the Sun, including an estimated ∼65 O-type stars and hundreds of OB stars. The Chandra Cygnus OB2 Legacy Survey is a large imaging program undertaken with the Advanced CCD Imaging Spectrometer on board the Chandra X-ray Observatory. The survey has imaged the central 0.5 deg ^2 of the Cyg OB2 association with an effective exposure of ∼120 ks and an outer 0.35 deg ^2 area with an exposure of ∼60 ks. Here we describe the survey design and observations, discuss the data reduction and source detection, and present a catalog of ∼8000 X-ray point sources. The survey design employs a grid of 36 heavily (∼50%) overlapping pointings, a method that overcomes Chandra's low off-axis sensitivity and produces a highly uniform exposure over the inner 0.5 deg ^2 . The full X-ray catalog is described here and is made available online.

We investigate the uncertainties of fitted X-ray model parameters and fluxes for relatively faint Chandra ACIS-I source spectra. Monte-Carlo (MC) simulations are employed to construct a large set of 150,000 fake X-ray spectra in the low photon count statistics regime (from 20 to 350 net counts) using the XSPEC spectral model fitting package. The simulations employed both absorbed thermal (APEC) and non-thermal (power-law) models, in concert with the Chandra ACIS-I instrument response and interstellar absorption. Simulated X-ray spectra were fit assuming a wide set of different input parameters and C-statistic minimization criteria to avoid numerical artifacts in the accepted solutions. Results provide an error estimate for each parameter (absorption, NH, plasma temperature, kT, or power-law slope, Gamma, and flux, and for different background contamination levels. The distributions of these errors are studied as a function of the 1 sigma quantiles and we show how these correlate with different model parameters, net counts in the spectra and relative background level. Maps of uncertainty in terms of the 1 sigma quantiles for parameters and flux are computed as a function of spectrum net counts. We find very good agreement between our estimated X-ray spectral parameter and flux uncertainties and those recovered from spectral fitting for a subset of the X-ray sources detected in the Chandra Cygnus OB2 Legacy Survey diagnosed to be Association members and that have between 20 and 350 net counts. Our method can provide uncertainties for spectral parameters whenever formal X-ray spectral fits cannot be well-constrained, or are unavailable, and predictions useful for computing Chandra ACIS-I exposure times for observation planning.

Habitat loss is a key factor in the ongoing freshwater biodiversity crisis. A promising way to help tackle the rapid decline in freshwater biodiversity is to improve the potential for artificial wetlands to provide habitat for aquatic wildlife. Farm dams (also known as agricultural ponds) are among the most abundant waterbodies in agricultural landscapes and can act as "oases" against droughts for many species. Despite their prominent role in agriculture, predictive models to evaluate their ecological potential are yet to emerge. Here we use a continental-scale data set of 104,013 audio recordings from citizen scientists to identify and locate 107 species of frogs near 8800 Australian farm dams. Frog species are among the most threatened taxa on earth and we asked: What characteristics promote higher frog species richness at farm dams? We found that the highest values of frog species richness were at old (>20 years) farm dams of intermediate size (0.1 ha in surface area), with small or medium rainfall catchments (<10 ha), and situated near other freshwater systems or conservation sites. The relationships shown here are highly generalisable and applicable on a continental scale. By identifying quantifiable features improving the ecological value of farm dams, we help identify "win-win" outcomes for agricultural productivity and conservation. In the future, "biodiversity credit" policies could incentivise large-scale ecological restoration by rewarding individuals who invest in enhancing their farm dams to support and promote local biodiversity.

We analyze the X-ray spectra of the ∼8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (this focus issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter sample we fit the spectra with absorbed isothermal models. In order to limit uncertainties, for most of the fainter Cygnus OB2 members in this sample we constrain the spectral parameters to characteristic ranges defined from the brightest stars. For X-ray sources with <20 net photons we adopt a conversion factor from detected photon flux to intrinsic flux. This was defined, building on the results for the previous sample, as a function of the 20% quantile of the detected photon energy distributions, which we prove to also correlate well with extinction. We then use the X-ray extinction from the spectral fits to constrain the ratio between optical and X-ray extinction toward Cyg OB2, finding it consistent with standard “Galactic” values, when properly accounting for systematics. Finally, we exploit the large number of sources to constrain the average coronal abundances of several elements, through two different ensemble analyses of the X-ray spectra of low-mass Cyg OB2 members. We find the pattern of abundances to be largely consistent with that derived for the young stellar coronae in the Orion Nebula Cluster.

ABSTRACT OB associations are important probes of recent star formation and Galactic structure. In this study, we focus on the Auriga constellation, an important region of star formation due to its numerous young stars, star-forming regions, and open clusters. We show using Gaia data that its two previously documented OB associations, Aur OB1 and OB2, are too extended in proper motion and distance to be genuine associations, encouraging us to revisit the census of OB associations in Auriga with modern techniques. We identify 5617 candidate OB stars across the region using photometry, astrometry and our SED fitting code, grouping these into five high-confidence OB associations using HDBSCAN. Three of these are replacements to the historical pair of associations – Aur OB2 is divided between a foreground and a background association – while the other two associations are completely new. We connect these OB associations to the surrounding open clusters and star-forming regions, analyse them physically and kinematically, constraining their ages through a combination of 3D kinematic traceback, the position of their members in the HR diagram and their connection to clusters of known age. Four of these OB associations are expanding, with kinematic ages up to a few tens of Myr. Finally, we identify an age gradient in the region spanning several associations that coincides with the motion of the Perseus spiral arm over the last ∼20 Myr across the field of view.