# The Gaia-ESO Survey: Calibrating the lithium–age relation with open clusters and associations ★★★

## II. Expanded cluster sample and final membership selection.

M.L. Gutiérrez Albarrán¹, D. Montes¹, H.M. Tabernero¹,², J.I. González Hernández³,⁴, E. Marfil¹,⁵, A. Frasca⁶, A. C. Lanzafame⁷,⁶, A. Klutsch⁸,⁶, E. Franciosini⁹, S. Randich⁹, R. Smiljanic¹⁰, A. J. Korn¹¹, G. Gilmore¹², E. J. Alfaro¹³, T. Bensby¹⁶, K. Biazzo¹⁷, A. Casey²⁰, G. Carraro²³, F. Damiani¹⁵, S. Feltzing¹⁶, P. François²², F. Jiménez Esteban¹⁴, L. Magrini⁹, L. Morbidelli⁹, L. Prisinzano¹⁵, T. Prusti²¹, C. C. Worley¹⁹, S. Zaggia¹⁸ and GES builders

*(Affiliations can be found after the references)*

Received October 31, 2023; accepted January 27, 2024.

### ABSTRACT

**Context.** The Li abundance observed in pre-main sequence and main sequence late-type stars is strongly age-dependent, but also shows a complex pattern depending on several parameters, such as rotation, chromospheric activity, and metallicity. The best way to calibrate these effects, and with the aim of studying Li as an age indicator for FGK stars, is to calibrate coeval groups of stars, such as open clusters (OCs) and associations.

**Aims.** We present a considerable target sample of 42 OCs and associations —with an age range from 1 Myr to 5 Gyr— observed within the *Gaia*-ESO survey (GES), and using the latest data provided by GES iDR6 and the most recent release of *Gaia* that was then available, EDR3. As part of this study, we update and improve the membership analysis for all 20 OCs presented in our previous article.

**Methods.** We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES, complemented with detailed bibliographical searches, and based on numerous criteria: from radial velocity distributions, to the astrometry (proper motions and parallaxes) and photometry provided by *Gaia*, to gravity indicators (log g and the γ index), [Fe/H] metallicity, and Li content in diagrams of EW(Li) (Li equivalent widths) versus T<sub>eff</sub>.

**Results.** We obtain updated lists of cluster members for the whole target sample, as well as a selection of Li-rich giant contaminants obtained as an additional result of the membership process. Each selection of cluster candidates was thoroughly contrasted with numerous existing membership studies using data from *Gaia* to ensure the most robust results.

**Conclusions.** These final cluster selections will be used in the third and last paper of this series, which reports the results of a comparative study characterising the observable Li dispersion in each cluster and analysing its dependence on several parameters, allowing us to calibrate a Li–age relation and obtain a series of empirical Li envelopes for key ages in our sample.

**Key words.** Galaxy: open clusters and associations: general – Stars: late-type – Stars: abundances – Techniques: spectroscopic

## 1. Introduction

Lithium in the form of its main isotope ⁷Li (henceforth Li) is mainly observable through the spectral resonance lines of neutral Li, forming a doublet at 6707.76 Å and 6707.89 Å (e.g., Pallavicini et al. 1992; Soderblom 2010; Lyubimkov 2016; Randich & Magrini 2021). Li is a very fragile element that is easily destroyed in stellar interiors, burning at temperatures above ~ 2.5–3 x 10⁶ K, which corresponds to the temperature at the base of the convective zone of a solar-mass star on the zero-age main sequence (ZAMS; e.g., Pallavicini et al. 1990; Siess et al. 2000; Magrini et al. 2021). As a result, Li is slowly being depleted and its surface abundance decreases over time in solar-type and lower-mass FGK stars. The depletion of Li in these

stars is strongly mass (colour) and age dependent: Li depletion increases with decreasing mass at a given age, and the timescales for significant Li depletion typically range from 10–20 Myr in M-type stars to ~ 100 Myr in K-type stars, and ~ 1 Gyr in G-type stars (e.g., Jeffries et al. 2014; Soderblom et al. 2014; Lyubimkov 2016; Randich & Magrini 2021). For their part, Li-rich giants show an additional contribution of surface Li that can only be accounted for with non-standard mixing mechanisms (e.g., Casey et al. 2016; Lyubimkov 2016; Magrini et al. 2021). As it only survives in the outer layers of a star, Li is a very sensitive indicator of youth in late-type stars, a useful tracer of stellar evolution for pre-main sequence (PMS), main sequence (MS), and post-MS stars, and is also particularly relevant for the determination of the age of stellar open clusters (OCs) and associations (e.g., Sestito & Randich 2005; Bouvier et al. 2018; Dumont et al. 2021a,b; Randich & Magrini 2021). For a more comprehensive overview of Li as an age indicator and the mechanics of Li depletion, we refer the reader to the introductory chapter of the published thesis of Gutiérrez Albarrán (2022) and all references therein.

<footnote>* Based on observations collected with ESO telescopes at the La Silla Paranal Observatory in Chile, for the Gaia-ESO Large Public Spectroscopic Survey (188.B-3002, 193.B-0936).</footnote>

<footnote>** All tables in Appendix C are only available in electronic form at the CDS via anonymous ftp to http://cdsarc.u-strasbg.fr/ (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/.</footnote>
The present work is the second of a series of three papers where we study Li as an age indicator for PMS and MS FGKM late-type stars, with the ultimate aim being to calibrate an empirical Li–age relation. In this paper, we describe an expanded sample of 42 OCs and associations covering an age range from 1 Myr to 5 Gyr, using data from the Gaia-ESO Survey (GES) and from the Gaia mission (both briefly discussed below). For all target clusters, we here constrain the cluster membership to obtain the final lists of candidate members for all target clusters, updating and expanding the work already published in Gutiérrez Albarrán et al. (2020) (hereafter, Paper I). In Gutiérrez Albarrán et al. (2024, in prep.) (hereafter, Paper III), presented as a companion paper to the present work, we finally reach our objective of constraining an empirical Li–age relation by first presenting a comparative study to quantify the observable Li dispersion for all target clusters and to analyse its dependence on several GES parameters, finally obtaining a set of empirical Li envelopes for several key ages in our sample. The analyses and results of this work are additionally presented in full in the published thesis of Gutiérrez Albarrán (2022).

The Gaia-ESO Survey (GES – Gilmore et al. 2012; Randich et al. 2013; Gilmore et al. 2022; Randich et al. 2022; Hourihane et al. 2023)<fnref n="1" /> is a large public spectroscopic survey that systematically covers all major components of the Milky Way and provides a homogeneous and ambitious overview of the distribution of the kinematics, dynamical structure, and chemical compositions in the Galaxy. The GES uses the multi-object spectrograph FLAMES (Fibre Large Array Multi Element Spectrograph) at the Very Large Telescope (VLT at ESO, Chile) to obtain high-resolution spectra (R=50,000) with UVES (Ultraviolet and Visual Echelle Spectrograph), as well as medium-resolution spectra (R=5,000–20,000) with GIRAFFE (e.g., Smiljanic et al. 2014; Frasca et al. 2015; Lanzafame et al. 2015; Gilmore et al. 2022). The observations, which began in December 2011 and were completed in January 2018, provide high-quality, uniformly calibrated spectroscopy for about 10⁵ stars and a sample of 65 (plus 18–20 archive) OCs and star forming regions of all ages, metallicities, and stellar masses (Bragaglia et al. 2022; Gilmore et al. 2022; Randich et al. 2022). In addition to measured spectra, this exploration uses well-defined samples from photometric surveys, such as VISTA (Sutherland et al. 2015), 2MASS (Skrutskie et al. 2006) and a variety of photometric surveys of OCs that cover all the major components of the Galaxy (Gilmore et al. 2012; Bragaglia et al. 2022).

Gaia (Gaia Collaboration et al. 2016b)<fnref n="2" /> is an ongoing, ambitious global space-astrometry ESA mission that launched in 2013 and started scientific operations in mid-2014. Initially a five-year mission, in 2018 the Gaia mission was extended to 2020, and currently Gaia has been firmly extended until the end of 2022, with indicative approval until the end of 2025 (Gaia Collaboration et al. 2021). Its primary science goal is to chart a three-dimensional map of the Galaxy, providing photometry and astrometry of unprecedented precision for most stars brighter than G=20 mag, and to obtain low- and medium-resolution spectroscopy data for most stars brighter than G=17 mag (Gilmore et al. 2012; Gaia Collaboration et al. 2016b). Gaia provides a large stellar census that allows an extensive overview of the origin, composition, formation, and evolution of the Galaxy, producing a stereoscopic dataset that includes positions, proper motions, parallaxes, radial velocities, brightness, and astrophysical

parameters of about one billion stars and other astronomical objects in the Milky Way and throughout the Local Group.

Gaia, with its revolutionary high-precision astrometry and photometry, can be optimally complemented with GES and the superior spectroscopic capabilities of its large ground-based telescope to obtain high-quality data, providing a rich dataset yielding 3D spatial distributions, 3D kinematics, individual chemical abundances, and improved astrophysical parameters for all target objects (e.g., Gilmore et al. 2012; Beccari et al. 2018; Cantat-Gaudin et al. 2018; Randich et al. 2018; Soubiran et al. 2018; Cánovas et al. 2019; Bossini et al. 2019). We finish this introductory section by further acknowledging that, in our analysis, we also made use of several sets of rotational periods ($P_{rot}$) provided by missions such as CoRoT, Kepler, K2, and TESS.

This paper is organised as follows. In Sect. 2 we describe the expanded cluster sample upon which we base our analysis. In Sect. 3 we discuss the expanded and improved membership cluster analysis, describing our selection criteria, which take into account radial velocities (RVs), Gaia astrometry and photometry, gravity indicators, metallicity, and Li content, as well as various comparisons with previous studies from the literature using Gaia data (Subsect. 3.3). In Subsect. 3.4 we additionally describe the selection process for Li-rich contaminants obtained as an additional result of the membership analyses. We present the final lists of candidate members for the results of this project in Sect. 5.

## 2. Data

### 2.1. GES data

The GES consortium is structured in 19 working groups (WGs), WG0 to WG18, organised in a workflow and dedicated to different survey tasks, from data flow implementation and target selection, to the homogenisation of recommended parameters and the determination of abundances, preparation and documentation, and, finally, the delivery of the external data products to both ESO and the public archive (e.g., Gilmore et al. 2012; Lanzafame et al. 2015; Sacco et al. 2015; Gilmore et al. 2022). Among these WGs, WG10 and WG11 are focused on the spectroscopic analysis of the GIRAFFE and UVES FGKM late-type stars, respectively (e.g., Gilmore et al. 2012; Sacco et al. 2015), while WG12 is dedicated to the analysis of PMS stars in the fields of young clusters using both UVES and GIRAFFE data. Each WG is divided into several nodes, providing internally consistent results with various approaches, models and methodology<fnref n="3" /> (Smiljanic et al. 2014; Lanzafame et al. 2015; Gilmore et al. 2022). For more specific and extensive details on the methodology employed by each node, see Smiljanic et al. (2014) (WG11), Lanzafame et al. (2015) (WG12), and Gilmore et al. (2022) (all WGs). In addition, for an in-depth overview of the GES consortium and the data reduction process for GES spectra, we further refer the reader to Paper I, and also to chapter 2 of the thesis and all references therein.

Six analysis cycles and internal data releases have been carried out by GES (from iDR1 to iDR6), including homogenised recommended astrophysical parameters and elemental abundances derived from all the observations collected until the com-

<footnote n="3">³ Our UCM node is part of both WG11 and WG12, and we determined atmospheric parameters by employing the code STEPAR (Tabernero et al. 2019), based on equivalent widths (EWs) automatizing the use of MOOG code (Sneden et al. 2012) and measured with the automatic tool TAME (Tool for Automatic Measurement of Equivalent Widths – Kang & Lee 2012; Tabernero et al. 2019). Also see Paper I.</footnote>

<footnote n="1">¹ https://www.gaia-eso.eu/</footnote>

<footnote n="2">² https://sci.esa.int/web/gaia</footnote>
pletion of the survey in January 2018 (Randich et al. 2022). Sub-
sequent public releases to the ESO archive include DR2 (pub-
lished in July 2015), DR3 (December 2016 and May 2017), DR4
(December 2020), and, finally, DR6 was released to the public in
July 2023⁴. iDR4, which was made available to the consortium
in February 2016, was the data release we used in paper I for our
preliminary cluster membership analysis. For all the following
analysis presented in this paper we used the data provided by
the sixth and last internal data release of GES (iDR6)<fnref n="5" />. iDR6 is
a full release, including all observations from the beginning of
the Survey until its completion in January 2018 (Gilmore et al.
2022; Franciosini et al. 2022a; Randich et al. 2022). Similarly
to iDR5, the same line list and grid of synthetic spectra used for
iDR4 were employed in this final release. This final GES cata-
logue provides results for 114,325 targets and includes data on
velocities, stellar parameters, and abundances for up to 32 ele-
ments (Randich et al. 2022).

The *Gaia* mission has also provided three data releases so
far, from *Gaia* DR1 (released in September 2016) and *Gaia* DR2
(April 2018) to *Gaia* DR3, which was released in June 2022,
following the *Gaia* Early Data Release 3 (EDR3)<fnref n="6" /> (*Gaia Col-
laboration et al. 2016a, 2018b, 2021, 2023*). Pending further
approvals of mission extensions, additional releases will also
take place, with *Gaia* DR5 being currently anticipated to con-
tain all collected data. In combination with GES iDR6, for the
present project we made use of the proper motions, parallaxes
and photometry provided by the latest data release then avail-
able, namely, the first installment of the third intermediate data
release, *Gaia* EDR3, which is based on the data of the first 34
months of the mission. *Gaia* EDR3 provides, with increased
precision and accuracy, an updated source list, astrometry, and
broad photometry in the *G*, *G*<sub>BP</sub> and *G*<sub>RP</sub> bands (*Gaia Collabo-
ration et al. 2021*).

## 2.2. Cluster sample

Our present sample from GES iDR6 includes 114,325 UVES
and GIRAFFE spectra of 42 OCs ranging in ages from 1 Myr
to 4.5 Gyr. Here below we list the changes and improvements
of the updated and improved iDR6 sample used in this work, in
contrast to our former iDR4 sample of only 20 OCs and associ-
ations, as used in Paper I:

- The number of stars in the current iDR6 sample is consid-
erably larger than the 12,493 UVES and GIRAFFE spectra
we had at our disposal for the preliminary cluster calibration
results in Paper I. This larger number of stars in the fields of
many of our sample clusters improved our membership anal-
ysis and increased the number of final candidates in more
than one case.

- Several issues were also solved thanks to this larger sample, such as the fact that with our former iDR4 file only a few UVES Li values were listed for many of the intermediate-age and old clusters, especially in the 3000–4000 K temperature range (M-type stars, more difficult to measure in those age

ranges)<fnref n="7" />. Thus, we had to make additional use of the indi-
vidual Li measurements derived by the OACT (*Osservatorio
Astrofisico di Catania*) node, adding to our sample a number
of GIRAFFE stars that had no recommended *EW*(*Li*) values
in the iDR4 file<fnref n="8" />. We note that, even considering the list of
spectra that were deprecated in the working version of iDR6
released to the consortium in December 2020, we observed
no such lack with our present sample, and our former mem-
ber lists are as a result either typically enlarged, or at least
stay with approximately the same number of candidate stars.

- In addition to the larger number of spectra overall, we now count with 86 OCs in the current iDR6, 48 more than in iDR4. Our former iDR4 sample offered a total of 38 clusters, 26 of them being OCs, and the remaining 12 of them being globular, which we discarded. We do not consider any iDR6 globular cluster in this updated analysis either, as Li cannot be used as a youth indicator in those cases.

- In our current iDR6 sample, the data for the OCs suffering from the contamination of nebular lines, which can affect the RV distributions and therefore the final membership analysis (e.g., Bonito et al. 2013, 2020), were corrected by recalculating the RVs and applying an alternative different sky background subtraction (Gilmore et al. 2022)<fnref n="9" />. Thus, we were able to add to our sample a number of young and intermediate-age clusters we had to formerly discard because they exhibited high differential nebulosity (such as NGC 2264, NGC 2451 A and B, NGC 3532, NGC 6530, and Trumpler 14¹⁰).

- Finally, in this study we focus on FGK stars, and so we discarded all stars with $T_{\text{eff}} > 7500$ K from the cluster sample in all cases. Given the nature of our membership criteria (see Subsect. 3), for each cluster we also deprecated from our analysis all stars with no measurements of RVs (VRAD in the iDR6 file), either $T_{\text{eff}}$ (TEFF) or the newly measured infrared photometric temperatures (TEFF IRFM), and/or EW(Li) (be it the corrected EWs, EWC_LI, or the improved EWs with an additional veiling correction, taking into account the estimate of the continuum emission due to accretion, EW_LI_UNVEIL (e.g., Lanzafame et al. 2015)).

Out of the 86 OCs measured in the iDR6 data, we discarded
all clusters with less than 100 stars, in order to have a sufficient
minimum number of stars in each cluster to ensure a statisti-
cally significant membership analysis. On the other hand, we did
not take into account any old clusters in the iDR6 file with ages

<footnote n="7">⁷ This was due to the homogenisation process in iDR4, as only WG11 and WG12 provided homogenised *EWs*, while WG10 provided only abundances. As a result, *EWs* for GIRAFFE stars in intermediate-age and old clusters were not available in this release, and in some cases, *EWs* for young clusters were also missing if the node values proved inconsistent, which happened especially for M-type stars in WG12.</footnote>

<footnote n="8">⁸ The iDR4 issues have since been overcome thanks to the new homogeneous analysis done in iDR6 (Gilmore et al. 2022; Randich et al. 2022).</footnote>

<footnote n="9">⁹ In general terms (this is by no means a straightforward method), an initial estimate would firstly be formed from the background emission affecting the spectrum of a target object, taking into account both astronomical and atmospheric sources. This estimated background spectrum would then be subtracted from the input science spectrum on a pixel by pixel basis, which gives rise an output spectrum consisting of emission from the science object alone.</footnote>

<footnote n="10">¹⁰ The reason why this improved sky background subtraction could not be done for earlier releases to the level iDR6 has finally achieved is that, the Survey being fiber-fed, subtraction of the nebular sky background is not a straightforward procedure (Bonito et al. 2020)</footnote>

<footnote n="4">⁴ DR6 data products available from the ESO Science Archive Por-
tal at https://archive.eso.org/scienceportal/home?data_
collection=GAIAESO&publ_date=2023-07-02</footnote>

<footnote n="5">⁵ iDR6 data products released to the GES consortium at http://ges.ast.cam.ac.uk/GESwiki/GeSDR6/DataProducts</footnote>

<footnote n="6">⁶ See https://www.cosmos.esa.int/web/gaia/earlydr3 and https://www.cosmos.esa.int/web/gaia/dr3 to access *Gaia* EDR3 and DR3, respectively</footnote>
<table><caption>**Table 1.** Age estimates, reddening, distance to the Sun, and GES and Gaia membership studies from the literature for the 16 star forming regions and young clusters in our sample.</caption><thead><tr><th>Cluster name</th><th>Age (Myr)</th><th>E(B-V)<sup>a</sup> (dex)</th><th>Distance (kpc)</th><th>References Ages</th><th>References Distance</th><th>GES and Gaia membership studies</th></tr></thead><tbody><tr><td>NGC 6530</td><td>1–2</td><td>0.44 ± 0.10</td><td>1.33</td><td>3, 4, 18, 19, 20, 33</td><td>18, 19, 20, 21</td><td>2, 3, 16, 18, 19, 22</td></tr><tr><td>ρ Oph</td><td>1–3</td><td>0.76 ± 0.13</td><td>0.13 ± 0.01</td><td>1, 2, 3, 4, 5, 33, 61</td><td>1, 2, 6</td><td>1, 2, 3, 6</td></tr><tr><td>Trumpler 14</td><td>1–3</td><td>0.61 ± 0.10</td><td>2.90</td><td>3, 4, 13, 33</td><td>13, 14</td><td>2, 3, 15, 17</td></tr><tr><td>Cha I</td><td>2</td><td>0.18 ± 0.08</td><td>0.16 ± 0.02</td><td>3, 4, 7, 8, 9, 33</td><td>7, 8, 10, 11</td><td>2, 3, 7, 8, 10, 11<sup>b</sup>, 12</td></tr><tr><td>NGC 2244</td><td>4</td><td>0.49 ± 0.09</td><td>1.59</td><td>3, 23, 24, 33</td><td>21, 23, 24, 25</td><td>3, 15, 24, 26</td></tr><tr><td>NGC 2264</td><td>4</td><td>0.05 ± 0.05</td><td>0.76</td><td>2, 3, 4, 27, 28, 29, 30, 33</td><td>2, 15, 25, 28, 29</td><td>2, 3, 15, 31, 32</td></tr><tr><td>λ Ori</td><td>6</td><td>0.09 ± 0.04</td><td>0.41</td><td>3, 43</td><td>34</td><td>3, 15</td></tr><tr><td>Col 197</td><td>13</td><td>0.64 ± 0.07</td><td>0.80–0.90</td><td>3, 5, 13, 35, 36, 37</td><td>13, 25, 35, 36, 37</td><td>3, 15</td></tr><tr><td>γ Vel</td><td>10–20</td><td>0.04 ± 0.03</td><td>0.35–0.40</td><td>2, 3, 8, 10, 38, 39, 40, 41, 42</td><td>8, 10, 38, 40, 41, 42</td><td>2, 3, 8, 10, 38, 39, 40, 41, 42, 44, 45, 46</td></tr><tr><td>NGC 2232</td><td>18–32</td><td>0.04 ± 0.03</td><td>0.32</td><td>3, 5, 43, 47, 48, 64</td><td>13, 25, 48</td><td>3, 15, 43, 48</td></tr><tr><td>NGC 2547</td><td>20–45</td><td>0.06 ± 0.03</td><td>0.36 ± 0.02</td><td>2, 3, 5, 8, 32, 41, 42, 48, 49, 50, 51, 63, 64</td><td>2, 8, 41, 42, 48</td><td>2, 3, 8, 15, 32, 42, 48, 49, 52</td></tr><tr><td>IC 2391</td><td>36 ± 2<sup>c</sup></td><td>0.03 ± 0.01</td><td>0.16 ± 0.01</td><td>3, 48, 51, 53, 64</td><td>2, 48, 54, 55, 56</td><td>2, 3, 15, 49, 52, 62</td></tr><tr><td>IC 2602</td><td>35 ± 1<sup>c</sup></td><td>0.04 ± 0.02</td><td>0.15 ± 0.01</td><td>3, 5, 48, 51, 64</td><td>2, 48, 56</td><td>2, 3, 15, 49, 52, 62</td></tr><tr><td>IC 4665</td><td>38 ± 3<sup>c</sup></td><td>0.15 ± 0.02</td><td>0.36 ± 0.01</td><td>3, 5, 48, 51, 64</td><td>2, 48, 57, 58</td><td>2, 3, 15, 48, 49, 52, 62</td></tr><tr><td>NGC 2451 B</td><td>39 ± 1<sup>c</sup></td><td>0.10 ± 0.03</td><td>0.36</td><td>3, 42, 48, 49, 51, 59, 60</td><td>42, 48, 59</td><td>2, 3, 15, 42, 48, 49</td></tr><tr><td>NGC 2451 A</td><td>44 ± 2<sup>c</sup></td><td>0.02 ± 0.02</td><td>0.19</td><td>3, 42, 48, 49, 51, 59, 60</td><td>42, 48, 59</td><td>2, 3, 15, 42, 48, 49</td></tr></tbody></table>

**Notes.** (a) References for the E(B-V) reddening adopted from Jackson et al. (2021) for all clusters. (b) GES studies that reference the clusters and/or study them without taking membership analysis primarily into account. (c) Updated cluster ages using Gaia data, as listed by Bossini et al. (2019).

**References.** For the cluster ages, reddening and distances shown here we chose the most recent or most robust estimates for each cluster, but several other studies are further cited here, and age estimates are additionally briefly discussed in the individual notes of Appendix A: (1) Rigliaco et al. (2016); (2) Spina et al. (2017); (3) Jackson et al. (2021); (4) Randich et al. (2020); (5) Romano et al. (2021); (6) Cánovas et al. (2019); (7) Spina et al. (2014a); (8) Sacco et al. (2015); (9) López Martí et al. (2013); (10) Frasca et al. (2015); (11) Roccatagliata et al. (2018); (12) Galli et al. (2021); (13) Sampedro et al. (2017); (14) Mel'nik & Dambis (2017); (15) Cantat-Gaudin et al. (2018); (16) Castro-Ginard et al. (2020); (17) Damiani et al. (2017); (18) Wright et al. (2019); (19) Damiani et al. (2019); (20) Prisinzano et al. (2005); (21) Kuhn et al. (2019); (22) Prisinzano et al. (2019); (23) Mužić et al. (2019); (24) Michalska (2019); (25) Kharchenko et al. (2005); (26) Carrera et al. (2019); (27) Arancibia-Silva et al. (2020); (28) Bonito et al. (2020); (29) Gillen et al. (2020); (30) Venuti et al. (2019); (31) Maíz Apellániz (2019); (32) Jackson et al. (2016); (33) Randich et al. (2022); (34) Dib et al. (2018); (35) Dias et al. (2019); (36) Bonatto & Bica (2010); (37) Vande Putte et al. (2010); (38) Jeffries et al. (2014); (39) Spina et al. (2014b); (40) Franciosini et al. (2018); (41) Beccari et al. (2018); (42) Franciosini et al. (2022b); (43) Binks et al. (2022); (44) Damiani et al. (2014); (45) Prisinzano et al. (2016); (46) Cantat-Gaudin et al. (2019); (47) Liu & Pang (2019); (48) Pang et al. (2021); (49) Randich et al. (2018); (50) Sestito & Randich (2005); (51) Bossini et al. (2019); (52) Bravi et al. (2018); (53) Dumont et al. (2021b); (54) Platais et al. (2007); (55) De Silva et al. (2013); (56) Smiljanic et al. (2011); (57) Martín & Montes (1997); (58) Jeffries et al. (2009); (59) Silaj & Landstreet (2014); (60) Netopil & Paunzen (2013); (61) Kiman et al. (2021); (62) Gómez Garrido (2015); (63) Oliveira et al. (2003), (64) Jeffries et al. (2023).

older than 2 Gyr, other than the ones we already analysed in Paper I, due to the fact that the Li content in GKM-type stars is strongly depleted in old clusters<fnref n="11" />. We also only analysed four additional clusters in the 0.8–2 Gyr age range, as we already counted with nine old clusters in our first sample selection with ages 0.8–4.5 Gyr, and preferred rather to add more new young and intermediate-age clusters to our updated study, adding to the age ranges where Li is increasingly relevant for the analysis at hand. This criteria leaves us with a total of 42 clusters constituting our current sample, including seven star forming regions (1–6 Myr) and nine young clusters (10–50 Myr), along with 13 intermediate-age clusters (90–575 Myr), and 13 old clusters (0.7–4.5 Gyr).

A number of membership studies have already been conducted, with potential members identified for all 42 clusters selected in the present paper. These studies, particularly those who made use of GES data and specifically iDR6 data, were of great help to evaluate the goodness of our membership analysis by comparing our final candidates with previous membership lists. We show all these studies in Tables 1 and 2. We divided the sample clusters into groups according to age, and show them in two tables for convenience due to their length: Table 1 lists the young (1–50 Myr) clusters, while Table 2 lists the intermediate-age (50–700 Myr) and old clusters (> 700 Myr). In the individual notes of Appendix A, where we detail all cluster membership results, we reference these studies in greater detail for each of the clusters. These tables also list the age estimates, distances and reddening values from the literature for all clusters, and we additionally refer to Table 3 for RVs from the literature, to Table 4 for proper motions (pmra and pmdec) and parallaxes (π) from Gaia, and to Table 5 for [Fe/H] metallicities from the literature.

### 3. Updated membership analysis

To obtain final lists of candidate members for the 42 clusters in our sample, we conducted a homogeneous and coherent analysis of their membership according to several criteria, which we now detail below:

- RV analysis (full description in Paper I and the published thesis): We studied the radial velocity distribution of each cluster with the aid of the RStudio environment for R by applying a two-sigma clipping procedure on the median, and then adopting a 2σ limit about the cluster mean yielded by the Gaussian fit to identify the most likely kinematic members. Depending on the specific cluster, this 2σ limit was sometimes enlarged to 3σ to further consider marginal members, as detailed in Appendix A. In Table 3 we show the fit parameters (mean values, associated dispersions and 2σ intervals) and the number of kinematic members resulting from the updated study of the distributions of RV for the cluster sample. We additionally report updated reference values estimated from the literature, and also present in the last two

<footnote n="11">¹¹ Furthermore, these clusters being so far, the already depleted GKM stars are also too faint to be observed by GES</footnote>
<table><caption>**Table 2.** Age estimates, reddening, distance to the Sun, and GES and Gaia membership studies from the literature for the 26 intermediate-age and old clusters in our sample.</caption><thead><tr><th>Cluster name</th><th>Age (Myr)</th><th>E(B-V)<sup>a</sup> (dex)</th><th>Distance (kpc)</th><th>References Ages</th><th>References Distance</th><th>GES and Gaia membership studies</th></tr></thead><tbody><tr><td>NGC 6405</td><td>94</td><td>0.14 ± 0.04</td><td>0.46</td><td>1, 2, 3</td><td>2, 3, 4, 5</td><td>1, 3, 6</td></tr><tr><td>Blanco 1</td><td>94 ± 5<sup>c</sup></td><td>-0.01 ± 0.03</td><td>0.23–0.24</td><td>1, 10, 15, 16, 17, 18, 53</td><td>17, 18</td><td>1, 18</td></tr><tr><td>NGC 6067</td><td>120</td><td>0.34 ± 0.04</td><td>1.4–1.7</td><td>1, 7, 8, 9, 10, 11</td><td>4, 5</td><td>1, 6, 11</td></tr><tr><td>NGC 6649</td><td>120</td><td>1.43 ± 0.05</td><td>1.8 ± 0.1</td><td>12, 13, 14</td><td>12, 14</td><td>1, 6</td></tr><tr><td>NGC 2516</td><td>125–138</td><td>0.11 ± 0.03</td><td>0.41</td><td>1, 9, 19, 21, 22</td><td>22, 23, 24</td><td>1, 6, 20, 22, 25, 26, 27, 28</td></tr><tr><td>NGC 6709</td><td>173 ± 34<sup>c</sup></td><td>0.27 ± 0.02</td><td>1.1</td><td>1, 10, 15, 29</td><td>4, 29</td><td>1, 6</td></tr><tr><td>NGC 6259</td><td>210</td><td>0.63 ± 0.09</td><td>1.9</td><td>1, 9, 30, 31</td><td>4, 29</td><td>1, 6</td></tr><tr><td>NGC 6705</td><td>280</td><td>0.40 ± 0.06</td><td>1.88</td><td>1, 9, 10, 26, 27</td><td>4, 24</td><td>1, 6, 26, 27, 32, 33, 34<sup>b</sup></td></tr><tr><td>Berkeley 30</td><td>300</td><td>0.51 ± 0.04</td><td>4.7–4.9</td><td>1, 4, 10, 30</td><td>4, 30</td><td>1, 6</td></tr><tr><td>NGC 6281</td><td>314</td><td>0.18 ± 0.02</td><td>0.47–0.51</td><td>1, 24, 30, 36</td><td>4, 30, 36</td><td>1, 6</td></tr><tr><td>NGC 3532</td><td>399 ± 5<sup>c</sup></td><td>0.05 ± 0.02</td><td>0.48–0.49</td><td>1, 10, 15, 28, 35</td><td>4, 28</td><td>1, 6, 28, 36, 37</td></tr><tr><td>NGC 4815</td><td>560</td><td>0.70 ± 0.07</td><td>2.40–2.90</td><td>1, 26, 27, 31, 37</td><td>24, 37</td><td>1, 6, 26, 27, 32, 33, 34, 37</td></tr><tr><td>NGC 6633</td><td>575</td><td>0.15 ± 0.04</td><td>0.39</td><td>1, 26, 27, 31, 38</td><td>4, 24</td><td>1, 6, 20, 26, 27</td></tr><tr><td>NGC 2477</td><td>700</td><td>0.31</td><td>1.4</td><td>39</td><td>3, 4, 6, 12, 39</td><td>6, 40</td></tr><tr><td>Trumpler 23</td><td>800</td><td>0.68 ± 0.04</td><td>2.20</td><td>1, 26, 27, 31, 43</td><td>24, 43</td><td>1, 26, 27, 42</td></tr><tr><td>Berkeley 81</td><td>860 ± 100</td><td>0.85 ± 0.04</td><td>3.00</td><td>1, 10, 26, 27, 31, 34, 45</td><td>24, 45</td><td>1, 6, 26, 27</td></tr><tr><td>NGC 2355</td><td>900</td><td>0.13 ± 0.03</td><td>1.80 ± 0.07</td><td>1, 4, 10, 30</td><td>42</td><td>1, 6</td></tr><tr><td>NGC 6802</td><td>900</td><td>0.79 ± 0.06</td><td>1.80</td><td>1, 10, 26, 27, 31, 44</td><td>24</td><td>1, 6, 26, 27, 44</td></tr><tr><td>NGC 6005</td><td>973 ± 4<sup>c</sup></td><td>0.49 ± 0.06</td><td>2.70</td><td>1, 15, 26, 27, 31</td><td>4, 24</td><td>1, 6, 26, 27</td></tr><tr><td>Pismis 18</td><td>1200 ± 400</td><td>+0.22 ± 0.04</td><td>2.20</td><td>27, 28, 31</td><td>24</td><td>1, 27, 28</td></tr><tr><td>Melotte 71</td><td>1294 ± 89<sup>c</sup></td><td>0.11</td><td>2.2–3.2</td><td>8, 15</td><td>4, 42</td><td>6</td></tr><tr><td>Pismis 15</td><td>1300</td><td>0.56 ± 0.05</td><td>2.6–2.9</td><td>1, 4, 30</td><td>4, 30</td><td>1, 6</td></tr><tr><td>Trumpler 20</td><td>1400</td><td>0.37 ± 0.03</td><td>3.00</td><td>1, 10, 26, 27, 41</td><td>24, 41</td><td>1, 6, 26, 27, 33, 41, 46</td></tr><tr><td>Berkeley 44</td><td>1600</td><td>0.90 ± 0.07</td><td>1.80–3.10</td><td>1, 10, 26, 27, 31, 47</td><td>24, 47</td><td>1, 6, 26, 27</td></tr><tr><td>NGC 2243</td><td>4000 ± 120</td><td>0.04 ± 0.04</td><td>4.50</td><td>1, 10, 27, 31, 48, 49, 50</td><td>24, 32, 49</td><td>1, 6, 27</td></tr><tr><td>M67</td><td>4000–4500</td><td>0.059</td><td>0.90</td><td>21, 38, 51, 52</td><td>4, 12</td><td>40</td></tr></tbody></table>

**Notes.** (a) References for the E(B-V) reddening adopted from Jackson et al. (2021) for all clusters except for NGC 2477 Rain et al. (2021). (b) GES studies that reference the clusters and/or study them without taking membership analysis primarily into account. (c) Updated cluster ages using Gaia data, as listed by Bossini et al. (2019).

**References.** For the cluster ages, reddening and distances shown here we chose the most recent or most robust estimates for each cluster. Several other studies are further cited here, and age estimates are additionally briefly discussed in the individual notes of Appendix A. (1) Jackson et al. (2021); (2) Kılıçoğlu et al. (2016); (3) Gao (2018); (4) Kharchenko et al. (2005); (5) Mel'nik & Dambis (2017); (6) Cantat-Gaudin et al. (2018); (7) Frinchaboy & Majewski (2008); (8) Netopil et al. (2016); (9) Randich et al. (2020); (10) Romano et al. (2021); (11) Rangwal et al. (2019); (12) Dib et al. (2018); (13) Liu & Pang (2019); (14) Alonso-Santiago et al. (2020); (15) Bossini et al. (2019); (16) Pang et al. (2021); (17) Gillen et al. (2020); (18) Zhang et al. (2020); (19) Binks et al. (2022); (20) Randich et al. (2018); (21) Dumont et al. (2021b); (22) Franciosini et al. (2022b); (23) Jeffries et al. (2001); (24) Dias et al. (2002); (25) Jackson et al. (2016); (26) Jacobson et al. (2016); (27) Magrini et al. (2017); (28) Fritzewski et al. (2019); (29) Vande Putte et al. (2010); (30) Sampedro et al. (2017); (31) Magrini et al. (2018); (32) Magrini et al. (2014); (33) Tautvaišienė et al. (2015); (34) Magrini et al. (2015); (35) Dobbie et al. (2012); (36) Fritzewski et al. (2021); (35) Hetem & Gregorio-Hetem (2019); (36) Joshi et al. (2016); (37) Friel et al. (2014); (38) Sestito & Randich (2005); (39) Rain et al. (2021); (40) Jadhav et al. (2021); (41) Donati et al. (2014b); (42) Buckner & Froebrich (2014); (43) Overbeek et al. (2017); (44) Tang et al. (2017); (45) Donati et al. (2014a); (46) Smiljanic et al. (2016); (47) Hayes & Friel (2014); (48) Heiter et al. (2014); (49) Jacobson et al. (2011); (50) Friel & Janes (1993); (51) Pallavicini et al. (2005); (52) Richer et al. (1998), (53) Jeffries et al. (2023).

columns of the table the mean values and associated disper-
sions of all final members of each cluster, after applying all
membership criteria and concluding our membership analy-
sis.

- **Gaia astrometry** (Subsect. 3.1): As the next step we reinforced the kinematic selection and obtain lists of probable astrometric members with the aid of the proper motions and parallaxes provided by Gaia EDR3, analysing the locus of probable members in pmra-versus-pmdec diagrams, and studying the parallax distributions with a Gaussian fit, similarly to the case of RVs. See Table 4 for the results of the study of parallax distributions for the cluster sample, including all fit parameters as well as the mean values and associated dispersions of all final candidates for each cluster,

after applying all criteria.

- **Gaia photometry** (Subsect. 3.2): We also made use of Gaia EDR3 photometry to analyse the goodness of the astrometric candidates and discard field contaminants in G-versus-G<sub>BP-RP</sub> colour-magnitude diagrams (CMDs).

- **Gravity indicators** (full description in Paper I and the thesis): We used the Kiel (log g-versus-T<sub>eff</sub>) diagram to identify and discard evolved outliers, such as Li-rich giant stars (see Subsect. 3.4) and other field contaminants. For all clusters we made use of the PARSEC isochrones (Bressan et al. 2012), with Z=0.019, and ages ranging from 50 Myr to 5 Gyr. In the case of young clusters, we mainly used the γ index defined by Damiani et al. (2014) to effectively discard
<table><caption>**Table 3.** Fit parameters and RV members for the target star forming regions and OCs.</caption><thead><tr><th rowspan="2">Cluster<sup>a</sup></th><th>RV<sup>b</sup></th><th colspan="2">2σ clipping</th><th rowspan="2">2σ membership<br>intervals</th><th rowspan="2">No. RV<br>members</th><th colspan="2">Final fit of member list</th></tr><tr><th>(km s<sup>-1</sup>)</th><th>⟨RV⟩ (km s<sup>-1</sup>)</th><th>σ (km s<sup>-1</sup>)</th><th>⟨RV⟩ (km s<sup>-1</sup>)</th><th>σ (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td>NGC 6530</td><td>-1.9 ± 7.2<sup>1</sup></td><td>0.3</td><td>2.1</td><td>[-3.8, 4.6]</td><td>477</td><td>0.0</td><td>2.5</td></tr><tr><td>ρ Oph</td><td>-7.0 ± 0.2<sup>2</sup></td><td>-6.4</td><td>2.1</td><td>[-10.6, -2.2]</td><td>48</td><td>-6.6</td><td>1.0</td></tr><tr><td>Trumpler 14</td><td>-10.5 ± 6.9<sup>3</sup></td><td>-7.9</td><td>2.9</td><td>[-13.7, -2.1]</td><td>231</td><td>-8.4</td><td>3.0</td></tr><tr><td>Cha I</td><td>14.6 ± 1.2<sup>4</sup></td><td>15.7</td><td>1.1</td><td>[13.5, 17.9]</td><td>102</td><td>15.6</td><td>1.1</td></tr><tr><td>NGC 2244</td><td>33.6 ± 0.7<sup>5</sup></td><td>30.7</td><td>2.7</td><td>[25.3, 36.1]</td><td>143</td><td>30.5</td><td>2.7</td></tr><tr><td>NGC 2264</td><td>20.2 ± 6.4<sup>3</sup></td><td>20.2</td><td>2.5</td><td>[15.2, 25.2]</td><td>618</td><td>19.9</td><td>2.2</td></tr><tr><td>λ Ori</td><td>27.5 ± 0.4<sup>5</sup></td><td>27.1</td><td>1.0</td><td>[25.1, 29.1]</td><td>208</td><td>27.3</td><td>1.0</td></tr><tr><td>Col 197</td><td>35.8 ± 2.3<sup>3</sup></td><td>20.8</td><td>1.1</td><td>[18.6, 23.0]</td><td>124</td><td>20.7</td><td>1.1</td></tr><tr><td>NGC 2232</td><td>25.4 ± 0.9<sup>3</sup></td><td>29.7</td><td>14.1</td><td>[-0.2, 56.2]</td><td>750</td><td>25.3</td><td>0.6</td></tr><tr><td>IC 2391</td><td>15.3 ± 0.2<sup>3</sup></td><td>15.4</td><td>2.3</td><td>[10.5, 19.7]</td><td>55</td><td>14.6</td><td>0.7</td></tr><tr><td>IC 2602</td><td>17.4 ± 0.2<sup>3</sup></td><td>13.2</td><td>10.2</td><td>[-8.0, 34.0]</td><td>309</td><td>17.5</td><td>0.6</td></tr><tr><td>IC 4665</td><td>-14.4 ± 0.8<sup>3</sup></td><td>-13.4</td><td>14.8</td><td>[-43.0, 16.2]</td><td>233</td><td>-14.0</td><td>0.7</td></tr><tr><td>NGC 6405</td><td>-9.2 ± 0.8<sup>6</sup></td><td>-9.0</td><td>6.5</td><td>[-22.0, 4.0]</td><td>251</td><td>-8.8</td><td>0.8</td></tr><tr><td>Blanco 1</td><td>5.8 ± 0.1<sup>6</sup></td><td>5.8</td><td>0.7</td><td>[4.4, 7.2]</td><td>141</td><td>5.8</td><td>0.6</td></tr><tr><td>NGC 6067</td><td>-39.9 ± 0.2<sup>1</sup></td><td>-38.1</td><td>3.3</td><td>[-44.7, -31.5]</td><td>209</td><td>-37.8</td><td>2.0</td></tr><tr><td>NGC 6649</td><td>-8.9 ± 0.5<sup>3</sup></td><td>-10.3</td><td>3.8</td><td>[-18.7, -1.9]</td><td>32</td><td>...</td><td>...</td></tr><tr><td>NGC 2516</td><td>23.8 ± 0.2<sup>3</sup></td><td>23.6</td><td>0.7</td><td>[22.3, 25.1]</td><td>460</td><td>23.8</td><td>0.8</td></tr><tr><td>NGC 6709</td><td>-7.8 ± 5<sup>1</sup></td><td>-7.7</td><td>10.3</td><td>[-28.3, 12.9]</td><td>322</td><td>-9.4</td><td>0.7</td></tr><tr><td>NGC 6259</td><td>-32.8 ± 0.4<sup>3</sup></td><td>-32.2</td><td>3.3</td><td>[-38.8, -25.6]</td><td>125</td><td>-32.5</td><td>2.5</td></tr><tr><td>NGC 6705</td><td>36.0 ± 0.2<sup>3</sup></td><td>35.6</td><td>1.9</td><td>[31.8, 39.4]</td><td>391</td><td>35.4</td><td>1.8</td></tr><tr><td>Berkeley 30</td><td>...</td><td>47.7</td><td>3.6</td><td>[40.5, 54.9]</td><td>78</td><td>48.1</td><td>3.0</td></tr><tr><td>NGC 6281</td><td>-5.0 ± 0.1<sup>3</sup></td><td>-4.8</td><td>1.9</td><td>[-8.6, -1.0]</td><td>82</td><td>-4.5</td><td>0.7</td></tr><tr><td>NGC 3532</td><td>5.4 ± 0.2<sup>3</sup></td><td>5.3</td><td>0.9</td><td>[3.5, 7.1]</td><td>518</td><td>5.4</td><td>0.9</td></tr><tr><td>NGC 4815</td><td>-29.8 ± 0.3<sup>3</sup></td><td>-29.0</td><td>3.4</td><td>[-36.4, -22.8]</td><td>68</td><td>-27.6</td><td>4.0</td></tr><tr><td>NGC 6633</td><td>-28.6 ± 0.1<sup>3</sup></td><td>-22.4</td><td>13.2</td><td>[-49.0, 0.0]<sup>e</sup></td><td>617</td><td>-28.3</td><td>0.9</td></tr><tr><td>NGC 2477</td><td>7.9 ± 0.1<sup>3</sup></td><td>8.1</td><td>1.0</td><td>[6.1, 10.1]</td><td>86</td><td>7.2</td><td>2.1</td></tr><tr><td>Trumpler 23</td><td>-61.4 ± 0.5<sup>3</sup></td><td>-61.3</td><td>1.8</td><td>[-57.7, -64.9]</td><td>51</td><td>-62.0</td><td>1.5</td></tr><tr><td>Berkeley 81</td><td>50.0 ± 0.7<sup>3</sup></td><td>47.7</td><td>2.2</td><td>[43.3, 52.1]</td><td>69</td><td>48.2</td><td>0.3</td></tr><tr><td>NGC 2355</td><td>36.9 ± 0.7<sup>3</sup></td><td>36.2</td><td>0.9</td><td>[34.4, 37.0]</td><td>119</td><td>36.3</td><td>0.9</td></tr><tr><td>NGC 6802</td><td>11.8 ± 0.4<sup>3</sup></td><td>13.0</td><td>1.8</td><td>[9.6, 16.6]</td><td>77</td><td>12.0</td><td>0.7</td></tr><tr><td>NGC 6005</td><td>-25.6 ± 0.5<sup>3</sup></td><td>-25.5</td><td>3.2</td><td>[-31.9, -19.1]</td><td>174</td><td>-24.6</td><td>1.1</td></tr><tr><td>Pismis 18</td><td>-28.5 ± 0.6<sup>3</sup></td><td>-29.7</td><td>2.8</td><td>[-35.3, -24.1]</td><td>41</td><td>-28.2</td><td>0.8</td></tr><tr><td>Melotte 71</td><td>51.3 ± 0.4<sup>3</sup></td><td>51.6</td><td>0.7</td><td>[50.2, 53.0]</td><td>71</td><td>51.1</td><td>0.4</td></tr><tr><td>Pismis 15</td><td>36.3 ± 0.7<sup>3</sup></td><td>35.0</td><td>0.1</td><td>[32.6, 37.4]</td><td>93</td><td>35.2</td><td>0.8</td></tr><tr><td>Trumpler 20</td><td>-39.8 ± 0.2<sup>3</sup></td><td>-39.8</td><td>1.4</td><td>[-42.6, -37]</td><td>451</td><td>-40.0</td><td>1.2</td></tr><tr><td>Berkeley 44</td><td>-7.6 ± 0.4<sup>3</sup></td><td>-8.8</td><td>0.4</td><td>[-10.4, -7.2]</td><td>39</td><td>-8.8</td><td>0.7</td></tr><tr><td>NGC 2243</td><td>59.6 ± 0.5<sup>3</sup></td><td>59.5</td><td>0.6</td><td>[58.3, 60.7]</td><td>469</td><td>59.6</td><td>0.6</td></tr><tr><td>M67</td><td>34.1 ± 0.1<sup>3</sup></td><td>34.0</td><td>0.8</td><td>[32.4, 35.6]</td><td>110</td><td>33.9</td><td>0.8</td></tr></tbody></table>

Notes. (a) As stated earlier in the text, regarding the clusters γ Vel, NGC 2547 and NGC 2451 A and B, we directly used the selections obtained by several studies listed in Table 1. (b) References for the mean cluster RVs: (1) Conrad et al. (2017); (2) Rigliaco et al. (2016); (3) Soubiran et al. (2018); (4) López Martí et al. (2013); (5) Carrera et al. (2019); (6) Gaia Collaboration et al. (2018b). (c) Final selection only consists of two members, final fit is not possible. (d) We could not find any prior RV measurements from the literature for this cluster. (e) The 2σ membership interval for this cluster has been obtained after excluding a large contaminant population in the middle of the distribution, as detailed in Appendix A.

giant contaminants as a first step of the membership process. Regarding the order of criteria, for the young clusters we discarded all giant contaminants using the γ index before performing the RV and astrometric analyses, due to their appreciable field contamination.

- Metallicity (full description in Paper I and the thesis): Using the spectroscopic index [Fe/H] as a proxy of metallicity,

an analysis of the metallicity distributions for each cluster (including all stars in the field before applying any other membership criteria) provided further confirmation of the membership of the candidate selections, and it was also taken into account to identify occasional rogue outliers that had been missed with all prior criteria. In Table 5 we present the updated results of the study of metallicity for the cluster sample, including all fit parameters, as well
<table><caption>**Table 4.** Fit parameters and parallax membership for the target star forming regions and OCs.</caption><thead><tr><th rowspan="2">Cluster<sup>a</sup></th><th>pmra<sup>b</sup></th><th>pmdec<sup>b</sup></th><th>Parallax (π)<sup>b</sup></th><th colspan="2">2σ clipping</th><th rowspan="2">2σ membership<br>intervals</th><th colspan="2">Final fit of member list</th></tr><tr><th>(mas yr<sup>-1</sup>)</th><th>(mas yr<sup>-1</sup>)</th><th>(mas)</th><th>⟨π⟩ (mas)</th><th>σ (mas)</th><th>⟨π⟩ (mas)</th><th>σ (mas)</th></tr></thead><tbody><tr><td>NGC 6530</td><td>1.32 ± 0.08<sup>1</sup></td><td>-2.07 ± 0.08<sup>1</sup></td><td>0.750<sup>2</sup></td><td>0.770</td><td>0.060</td><td>[0.65, 0.89]</td><td>0.769</td><td>0.062</td></tr><tr><td>ρ Oph</td><td>-1.70 ± 0.07<sup>1</sup></td><td>0.20 ± 0.07<sup>1</sup></td><td>7.692<sup>1</sup></td><td>7.237</td><td>0.120</td><td>[7.001, 7.481]</td><td>7.250</td><td>0.116</td></tr><tr><td>Trumpler 14</td><td>-6.54 ± 0.07<sup>1</sup></td><td>2.06 ± 0.07<sup>1</sup></td><td>0.340<sup>1</sup></td><td>0.390</td><td>0.050</td><td>[0.290, 0.490]</td><td>0.395</td><td>0.066</td></tr><tr><td>Cha I</td><td>-19.0 ± 5.00<sup>3</sup></td><td>2.00 ± 4.00<sup>3</sup></td><td>6.250<sup>1</sup></td><td>5.240</td><td>0.080</td><td>[5.080, 5.400]</td><td>5.247</td><td>0.087</td></tr><tr><td>NGC 2244</td><td>-1.70 ± 0.07<sup>1</sup></td><td>0.20 ± 0.07<sup>1</sup></td><td>0.590-0.710<sup>1</sup></td><td>0.670</td><td>0.060</td><td>[0.550, 0.790]</td><td>0.665</td><td>0.071</td></tr><tr><td>NGC 2264</td><td>-1.76 ± 0.08<sup>1</sup></td><td>-3.72 ± 0.07<sup>1</sup></td><td>1.250<sup>1</sup></td><td>1.380</td><td>0.040</td><td>[1.200, 1.600]</td><td>1.388</td><td>0.063</td></tr><tr><td>λ Ori</td><td>1.19 ± 0.51<sup>4</sup></td><td>-2.12 ± 0.39<sup>4</sup></td><td>2.462 ± 0.124<sup>4</sup></td><td>2.500</td><td>0.048</td><td>[2.404, 2.596]</td><td>2.492</td><td>0.060</td></tr><tr><td>Col 197</td><td>-5.81 ± 0.31<sup>4</sup></td><td>3.93 ± 0.39<sup>4</sup></td><td>1.135<sup>1</sup></td><td>1.020</td><td>0.050</td><td>[0.930, 1.130]</td><td>1.023</td><td>0.059</td></tr><tr><td>NGC 2232</td><td>-4.70<sup>4</sup></td><td>-1.80<sup>4</sup></td><td>3.067 ± 0.099<sup>4</sup></td><td>3.100</td><td>0.060</td><td>[2.980, 3.220]</td><td>3.111</td><td>0.058</td></tr><tr><td>IC 2391</td><td>-24.64 ± 0.88<sup>4</sup></td><td>23.32 ± 0.73<sup>4</sup></td><td>6.585<sup>5</sup></td><td>6.620</td><td>0.070</td><td>[6.490, 6.790]</td><td>6.613</td><td>0.086</td></tr><tr><td>IC 2602</td><td>-17.58 ± 0.83<sup>4</sup></td><td>10.70 ± 0.90<sup>4</sup></td><td>6.758<sup>5</sup></td><td>6.610</td><td>0.080</td><td>[6.460, 6.780]</td><td>6.629</td><td>0.084</td></tr><tr><td>IC 4665</td><td>-0.91 ± 0.28<sup>4</sup></td><td>-8.52 ± 0.28<sup>4</sup></td><td>2.872<sup>5</sup></td><td>2.900</td><td>0.060</td><td>[2.780, 3.020]</td><td>2.887</td><td>0.077</td></tr><tr><td>NGC 6405</td><td>1.31 ± 0.34<sup>4</sup></td><td>-5.85 ± 0.34<sup>4</sup></td><td>2.172<sup>4</sup></td><td>2.190</td><td>0.020</td><td>[2.150, 2.230]</td><td>2.185</td><td>0.022</td></tr><tr><td>Blanco I</td><td>18.74 ± 0.43<sup>4</sup></td><td>2.60 ± 0.44<sup>4</sup></td><td>4.210 ± 0.120<sup>4</sup></td><td>4.200</td><td>0.070</td><td>[4.060, 4.340]</td><td>4.225</td><td>0.062</td></tr><tr><td>NGC 6067</td><td>-1.91 ± 0.12<sup>4</sup></td><td>-2.59 ± 0.12<sup>4</sup></td><td>0.443 ± 0.065<sup>4</sup></td><td>0.470</td><td>0.028</td><td>[0.416, 0.528]</td><td>0.472</td><td>0.024</td></tr><tr><td>NGC 6649</td><td>-0.01 ± 0.18<sup>4</sup></td><td>-0.06 ± 0.18<sup>4</sup></td><td>0.467 ± 0.087<sup>4</sup></td><td>0.491</td><td>0.060</td><td>[0.371, 0.611]</td><td>...</td><td>...</td></tr><tr><td>NGC 2516</td><td>-4.75 ± 0.44<sup>4</sup></td><td>11.22 ± 0.35<sup>4</sup></td><td>2.417 ± 0.045<sup>4</sup></td><td>2.429</td><td>0.031</td><td>[2.367, 2.491]</td><td>2.428</td><td>0.033</td></tr><tr><td>NGC 6709</td><td>-1.76 ± 0.08<sup>1</sup></td><td>-3.72 ± 0.07<sup>1</sup></td><td>1.250<sup>1</sup></td><td>1.380</td><td>0.040</td><td>[1.200, 1.600]</td><td>0.916</td><td>0.033</td></tr><tr><td>NGC 6259</td><td>-1.02 ± 0.13<sup>4</sup></td><td>-2.89 ± 0.12<sup>4</sup></td><td>0.408 ± 0.057<sup>4</sup></td><td>0.432</td><td>0.061</td><td>[0.310, 0.550]</td><td>0.463</td><td>0.034</td></tr><tr><td>NGC 6705</td><td>-1.57 ± 0.16<sup>4</sup></td><td>-4.14 ± 0.17<sup>4</sup></td><td>0.427 ± 0.083<sup>4</sup></td><td>0.411</td><td>0.039</td><td>[0.333, 0.489]</td><td>0.404</td><td>0.054</td></tr><tr><td>Berkeley 30</td><td>-0.25 ± 0.17<sup>4</sup></td><td>-0.33 ± 0.13<sup>4</sup></td><td>0.164<sup>4</sup></td><td>0.209</td><td>0.138</td><td>[-0.074, 0.478]</td><td>0.201</td><td>0.144</td></tr><tr><td>NGC 6281</td><td>-1.86 ± 0.30<sup>4</sup></td><td>-4.02 ± 0.24<sup>4</sup></td><td>1.873 ± 0.080<sup>4</sup></td><td>1.880</td><td>0.030</td><td>[1.820, 1.940]</td><td>1.889</td><td>0.056</td></tr><tr><td>NGC 3532</td><td>-10.39 ± 0.40<sup>4</sup></td><td>5.18 ± 0.40<sup>4</sup></td><td>2.066 ± 0.062<sup>4</sup></td><td>2.086</td><td>0.025</td><td>[2.035, 2.136]</td><td>2.085</td><td>0.024</td></tr><tr><td>NGC 4815</td><td>-5.76 ± 0.11<sup>4</sup></td><td>-0.96 ± 0.10<sup>4</sup></td><td>0.261 ± 0.062<sup>4</sup></td><td>0.294</td><td>0.076</td><td>[0.142, 0.446]</td><td>0.294</td><td>0.084</td></tr><tr><td>NGC 6633</td><td>1.20 ± 0.33<sup>4</sup></td><td>-1.81 ± 0.30<sup>4</sup></td><td>2.525 ± 0.073<sup>4</sup></td><td>2.533</td><td>0.045</td><td>[-0.448, 3.432]<sup>d</sup></td><td>2.545</td><td>0.051</td></tr><tr><td>NGC 2477</td><td>-2.45 ± 0.01<sup>4</sup></td><td>0.88 ± 0.01<sup>4</sup></td><td>0.665 ± 0.037<sup>4</sup></td><td>0.688</td><td>0.030</td><td>[0.628, 0.784]</td><td>0.694</td><td>0.029</td></tr><tr><td>Trumpler 23</td><td>-4.18 ± 0.12<sup>4</sup></td><td>-4.75 ± 0.11<sup>4</sup></td><td>0.352 ± 0.059<sup>4</sup></td><td>0.353</td><td>0.040</td><td>[0.273, 0.433]</td><td>0.347</td><td>0.043</td></tr><tr><td>Berkeley 81</td><td>-1.20 ± 0.16<sup>4</sup></td><td>-1.83 ± 0.16<sup>4</sup></td><td>0.254 ± 0.090<sup>4</sup></td><td>0.261</td><td>0.054</td><td>[0.153, 0.369]</td><td>0.263</td><td>0.048</td></tr><tr><td>NGC 2355</td><td>-3.80 ± 0.14<sup>4</sup></td><td>-1.09 ± 0.13<sup>4</sup></td><td>0.497 ± 0.056<sup>4</sup></td><td>0.530</td><td>0.022</td><td>[0.486, 0.574]</td><td>0.529</td><td>0.015</td></tr><tr><td>NGC 6802</td><td>-2.81 ± 0.11<sup>4</sup></td><td>-6.44 ± 0.13<sup>4</sup></td><td>0.309 ± 0.067<sup>4</sup></td><td>0.335</td><td>0.028</td><td>[0.279, 0.391]</td><td>0.349</td><td>0.053</td></tr><tr><td>NGC 6005</td><td>-4.01 ± 0.12<sup>4</sup></td><td>-3.81 ± 0.12<sup>4</sup></td><td>0.362 ± 0.060<sup>4</sup></td><td>0.390</td><td>0.051</td><td>[0.288, 0.492]</td><td>0.327</td><td>0.135</td></tr><tr><td>Pismis 18</td><td>5.66 ± 0.11<sup>4</sup></td><td>-2.29 ± 0.11<sup>4</sup></td><td>0.332 ± 0.052<sup>4</sup></td><td>0.349</td><td>0.022</td><td>[0.305, 0.393]</td><td>0.342</td><td>0.100</td></tr><tr><td>Melotte 71</td><td>-2.45 ± 0.12<sup>4</sup></td><td>4.21 ± 0.11<sup>4</sup></td><td>0.434 ± 0.056<sup>4</sup></td><td>0.443</td><td>0.026</td><td>[0.391, 0.495]</td><td>0.348</td><td>0.100</td></tr><tr><td>Pismis 15</td><td>-5.30 ± 0.10<sup>4</sup></td><td>3.33 ± 0.12<sup>4</sup></td><td>0.354 ± 0.049<sup>4</sup></td><td>0.387</td><td>0.061</td><td>[0.265, 0.509]</td><td>0.401</td><td>0.058</td></tr><tr><td>Trumpler 20</td><td>-7.09 ± 0.09<sup>4</sup></td><td>0.18 ± 0.09<sup>4</sup></td><td>0.251 ± 0.045<sup>4</sup></td><td>0.272</td><td>0.030</td><td>[0.212, 0.332]</td><td>0.285</td><td>0.048</td></tr><tr><td>Berkeley 44</td><td>0.01 ± 0.13<sup>4</sup></td><td>-2.83 ± 0.14<sup>4</sup></td><td>0.303 ± 0.060<sup>4</sup></td><td>0.322</td><td>0.042</td><td>[0.238, 0.406]</td><td>0.341</td><td>0.062</td></tr><tr><td>NGC 2243</td><td>-1.28 ± 0.13<sup>4</sup></td><td>5.49 ± 0.13<sup>4</sup></td><td>0.211 ± 0.060<sup>4</sup></td><td>0.225</td><td>0.052</td><td>[0.121, 0.329]</td><td>0.220</td><td>0.048</td></tr><tr><td>M67</td><td>-10.97 ± 0.24<sup>4</sup></td><td>-2.95 ± 0.24<sup>4</sup></td><td>1.135 ± 0.051<sup>4</sup></td><td>1.161</td><td>0.026</td><td>[1.110, 1.214]</td><td>1.162</td><td>0.026</td></tr></tbody></table>

Notes. (a) As stated earlier in the text, regarding the clusters γ Vel, NGC 2547 and NGC 2451 A and B, we directly used the selections obtained by several studies listed in Table 1. (b) References for the mean cluster proper motions (pmra and pmdec) and parallaxes: (1) Kuhn et al. (2019); (2) Wright et al. (2019); (3) López Martí et al. (2013); (4) Cantat-Gaudin et al. (2018); (5) Kounkel & Covey (2019). (c) Final selection only consists of two members, final fit is not possible. (d) The 2σ membership interval for this cluster has been obtained after excluding a large contaminant population in the middle of the distribution, as detailed in Appendix A.

as the mean values and associated dispersions of all final candidates after concluding the membership analysis. We note that, in contrast to the study of RVs and parallaxes, we decided to adopt 3σ membership intervals, instead of a 2σ limit, as the starting point to identify the most likely metallicity members. We also note that we also accepted as final cluster candidates a number of marginal metallicity members deviating moderately from the starting 3σ limits, as they fully fulfilled the rest of criteria in our analysis. We refer the reader to the individual notes of Appendix A for more details on these marginal stars for each cluster, as well as listing, in the case of some of the clusters, a

smaller number of stars that were similarly accepted as robust cluster members, despite having [Fe/H] values that deviated more appreciably from the 3σ limit. It is thus apt to note that, in contrast to the most robust criteria in our analysis (particularly *Gaia* astrometry and CMDs), we consider metallicity to be among our less restrictive criteria and selected the final cluster members accordingly. This takes into account the inherent higher uncertainties related to [Fe/H] values derived from GIRAFFE spectra, even after the improvements achieved for the iDR6 data (e.g., Spina et al. 2014b; Gilmore et al. 2022), as well as the fact that the variations in stellar abundances also depend on several
<table><caption>**Table 5.** Fit parameters and metallicity membership for the target star forming regions and OCs.</caption><thead><tr><th rowspan="2">Cluster<sup>a</sup></th><th rowspan="2">[Fe/H]<br>(dex)</th><th rowspan="2">References<br>[Fe/H]</th><th colspan="2">3σ clipping</th><th rowspan="2">3σ membership<br>intervals<sup>d</sup></th><th colspan="2">Final fit of member list</th></tr><tr><th>⟨[Fe/H]⟩ (dex)</th><th>σ (dex)</th><th>⟨[Fe/H]⟩ (dex)</th><th>σ (dex)</th></tr></thead><tbody><tr><td>NGC 6530</td><td>-0.04 ± 0.01</td><td>1</td><td>0.00</td><td>0.08</td><td>[-0.24, 0.24]</td><td>-0.01</td><td>0.09</td></tr><tr><td>ρ Oph</td><td>-0.08 ± 0.02</td><td>2</td><td>0.02</td><td>0.09</td><td>[-0.25, 0.29]</td><td>-0.03</td><td>0.20</td></tr><tr><td>Trumpler 14</td><td>-0.03 ± 0.02</td><td>1</td><td>0.00</td><td>0.06</td><td>[-0.18, 0.18]</td><td>-0.02</td><td>0.07</td></tr><tr><td>Cha I</td><td>-0.07 ± 0.04</td><td>2, 3</td><td>-0.02</td><td>0.04</td><td>[-0.15, 0.10]</td><td>-0.04</td><td>0.05</td></tr><tr><td>NGC 2244</td><td>-0.23</td><td>4</td><td>-0.04</td><td>0.05</td><td>[-0.20, 0.11]</td><td>-0.05<sup>b</sup></td><td>0.05</td></tr><tr><td>NGC 2264</td><td>-0.06 ± 0.01</td><td>1, 2, 5</td><td>-0.03</td><td>0.06</td><td>[-0.21, 0.15]</td><td>-0.03</td><td>0.06</td></tr><tr><td>λ Ori</td><td>-0.01</td><td>4</td><td>-0.04</td><td>0.07</td><td>[-0.25, 0.17]</td><td>-0.03</td><td>0.06</td></tr><tr><td>Col 197</td><td>...</td><td>...</td><td>0.02</td><td>0.06</td><td>[-0.16, 0.20]</td><td>0.01</td><td>0.06</td></tr><tr><td>NGC 2232</td><td>0.04 ± 0.01</td><td>4</td><td>0.02</td><td>0.04</td><td>[-0.15, 0.15]</td><td>0.00</td><td>0.05</td></tr><tr><td>IC 2391</td><td>-0.03 ± 0.02</td><td>1, 2, 6, 7, 8</td><td>-0.03</td><td>0.03</td><td>[-0.10, 0.10]</td><td>-0.03</td><td>0.09</td></tr><tr><td>IC 2602</td><td>-0.02 ± 0.02</td><td>2, 7, 8</td><td>0.00</td><td>0.07</td><td>[-0.21, 0.21]</td><td>0.00</td><td>0.07</td></tr><tr><td>IC 4665</td><td>0.00 ± 0.02</td><td>2</td><td>-0.01</td><td>0.09</td><td>[-0.19, 0.17]</td><td>0.00</td><td>0.01</td></tr><tr><td>NGC 6405</td><td>0.07 ± 0.03</td><td>9, 10</td><td>0.03</td><td>0.12</td><td>[-0.33, 0.39]</td><td>0.00</td><td>0.06</td></tr><tr><td>Blanco 1</td><td>0.00</td><td>9, 12</td><td>0.01</td><td>0.02</td><td>[-0.05, 0.07]</td><td>0.01</td><td>0.02</td></tr><tr><td>NGC 6067</td><td>0.20 ± 0.08</td><td>9, 11</td><td>0.11</td><td>0.09</td><td>[-0.16, 0.38]</td><td>0.07</td><td>0.10</td></tr><tr><td>NGC 6649</td><td>...</td><td>...</td><td>0.03</td><td>0.19</td><td>[-0.54, 0.60]</td><td>...</td><td>...</td></tr><tr><td>NGC 2516</td><td>-0.02 ± 0.01</td><td>5, 13</td><td>0.02</td><td>0.05</td><td>[-0.13, 0.17]</td><td>0.00</td><td>0.06</td></tr><tr><td>NGC 6709</td><td>...</td><td>...</td><td>0.00</td><td>0.17</td><td>[-0.51, 0.51]</td><td>-0.08</td><td>0.07</td></tr><tr><td>NGC 6259</td><td>0.21 ± 0.04</td><td>11, 14</td><td>0.18</td><td>0.09</td><td>[-0.09, 0.45]</td><td>0.17</td><td>0.10</td></tr><tr><td>NGC 6705</td><td>0.16 ± 0.04</td><td>11, 13</td><td>0.06</td><td>0.08</td><td>[-0.18, 0.30]</td><td>0.06</td><td>0.08</td></tr><tr><td>Berkeley 30</td><td>0.10</td><td>15</td><td>-0.19</td><td>0.18</td><td>[-0.73, 0.35]</td><td>-0.20</td><td>0.19</td></tr><tr><td>NGC 6281</td><td>0.06</td><td>9, 17</td><td>0.07</td><td>0.20</td><td>[-0.53, 0.67]</td><td>0.08</td><td>0.15</td></tr><tr><td>NGC 3532</td><td>-0.07 ± 0.10</td><td>16</td><td>0.02</td><td>0.06</td><td>[-0.16, 0.20]</td><td>0.02</td><td>0.07</td></tr><tr><td>NGC 4815</td><td>0.11 ± 0.01</td><td>11, 13, 18</td><td>0.07</td><td>0.20</td><td>[-0.53, 0.67]</td><td>0.01</td><td>0.08</td></tr><tr><td>NGC 6633</td><td>-0.01 ± 0.11</td><td>11, 13, 19</td><td>-0.01</td><td>0.19</td><td>[-0.58, 0.56]</td><td>0.03</td><td>0.15</td></tr><tr><td>NGC 2477</td><td>0.07 ± 0.03</td><td>9, 17, 20</td><td>0.10</td><td>0.04</td><td>[-0.02, 0.22]</td><td>0.14</td><td>0.04</td></tr><tr><td>Trumpler 23</td><td>0.21 ± 0.04</td><td>11, 13, 21</td><td>0.16</td><td>0.10</td><td>[-0.14, 0.46]</td><td>0.19</td><td>0.07</td></tr><tr><td>Berkeley 81</td><td>0.22 ± 0.07</td><td>11, 13</td><td>0.18</td><td>0.11</td><td>[-0.15, 0.51]</td><td>0.18</td><td>0.10</td></tr><tr><td>NGC 2355</td><td>-0.11</td><td>4</td><td>-0.12</td><td>0.05</td><td>[-0.27, 0.03]</td><td>-0.13</td><td>0.05</td></tr><tr><td>NGC 6802</td><td>0.10 ± 0.02</td><td>11, 13, 22</td><td>0.03</td><td>0.15</td><td>[-0.42, 0.48]</td><td>0.05</td><td>0.13</td></tr><tr><td>NGC 6005</td><td>0.19 ± 0.02</td><td>11, 13</td><td>0.14</td><td>0.10</td><td>[-0.16, 0.44]</td><td>0.19</td><td>0.08</td></tr><tr><td>Pismis 18</td><td>0.22 ± 0.04</td><td>11, 13</td><td>0.06</td><td>0.11</td><td>[-0.27, 0.39]</td><td>0.13<sup>b</sup></td><td>0.02</td></tr><tr><td>Melotte 71</td><td>-0.27</td><td>9, 12</td><td>-0.33</td><td>0.24</td><td>[-1.08, 0.42]</td><td>-0.09</td><td>0.01</td></tr><tr><td>Pismis 15</td><td>0.01 ± 0.01</td><td>23</td><td>-0.10</td><td>0.22</td><td>[-0.76, 0.56]</td><td>-0.07</td><td>0.11</td></tr><tr><td>Trumpler 20</td><td>0.10 ± 0.05</td><td>11, 13</td><td>0.06</td><td>0.10</td><td>[-0.24, 0.36]</td><td>0.08</td><td>0.10</td></tr><tr><td>Berkeley 44</td><td>0.27 ± 0.06</td><td>11, 13</td><td>0.12</td><td>0.09</td><td>[-0.15, 0.39]</td><td>0.12</td><td>0.09</td></tr><tr><td>NGC 2243</td><td>-0.38 ± 0.04</td><td>11, 13, 17</td><td>-0.57</td><td>0.13</td><td>[-0.96, -0.18]</td><td>-0.61</td><td>0.13</td></tr><tr><td>M67</td><td>-0.01 ± 0.04</td><td>11, 17, 19, 24</td><td>0.00</td><td>0.05</td><td>[-0.15, 0.15]</td><td>0.00</td><td>0.05</td></tr></tbody></table>

Notes. (a) For γ Vel, NGC 2547 and NGC 2451 A and B, we used the selections listed in Table 1. (b) For these clusters, the final mean [Fe/H] values deviate more appreciably from the literature values. (c) We could not find any prior measurements from the literature for these clusters. (d) In the individual notes of Appendix A we also list the number of stars that we accepted as robust members despite deviating from these 3σ intervals.

(e) Final selection only consists of two members, final fit is not possible.

[Fe/H] references: Shown here are the most recent or robust estimates for each cluster, several studies are further cited in Appendix A: (1) Randich et al. (2018); (2) Spina et al. (2017); (3) Sacco et al. (2015); (4) Carrera et al. (2019); (5) Binks et al. (2022); (6) De Silva et al. (2013); (7) Dumont et al. (2021b); (8) Smiljanic et al. (2011); (9) Netopil & Paunzen (2013); (10) Kılıçoğlu et al. (2016); (11) Magrini et al. (2018); (12) Bossini et al. (2019); (13) Jacobson et al. (2016); (14) Casali et al. (2019); (15) Paunzen et al. (2010); (16) Fritzewski et al. (2019); (17) Heiter et al. (2014); (18) Friel et al. (2014); (19) Sestito & Randich (2005); (20) Rain et al. (2021); (21) Overbeek et al. (2017); (22) Tang et al. (2017); (23) Carraro et al. (2005); (24) Liu & Pang (2019).

astrophysical factors and processes, such as atomic diffusion or the influence of planets (e.g., Oh et al. 2018; Liu & Pang 2019; Yong et al. 2023). Finally, stellar Li itself also depends on metallicity (a dependence that will be further analysed in Paper III) (e.g., Jeffries 2014; Randich & Magrini

2021; Martos et al. 2023), which makes it important not to discard those cluster members that are iron-rich or iron-poor.

– Li content (full description in Paper I and the thesis): Probable candidates fulfilling the rest of criteria were consid-
ered Li members according to their locus in the EW(Li)-
versus-T<sub>eff</sub> diagrams, using as guides the upper Li enve-
lope of IC 2602 (35 Myr) (Montes et al. 2001; López-
Santiago et al. 2003), the upper (Neuhaeuser et al. 1997) and
lower (Soderblom et al. 1993) envelopes of the Pleiades (78–
125 Myr), and the upper envelope of the Hyades (750 Myr)
(Soderblom et al. 1993). As this work revolves around the
calibration of Li and its observable dispersion, this is one of
the final criteria so as to already count with a robust list of
probable cluster members and add the least bias to our study.
We refer to the individual notes of Appendix A to illustrate
some examples of the bias that we could add by filtering can-
didates on the basis of Li before other criteria. In addition,
for the star forming regions in our sample we additionally
considered as cluster candidates all kinematic and astromet-
ric members that are strong accretors with Hα10% > 270–
300 km s<sup>-1</sup>, disregardless of their Li content (Frasca et al.
2015; Sacco et al. 2017; Bonito et al. 2020)<sup>12</sup>. This is due to
the enhancement and/or underestimation that can be caused
by strong accretion (for more details, see Paper III).

- **Other Gaia membership studies** (Subsect. 3.3): Finally, we made use of additional studies conducted from Gaia DR1, DR2 and EDR3 data (e.g., Cantat-Gaudin et al. 2018; Randich et al. 2018; Soubiran et al. 2018; Bossini et al. 2019; Jackson et al. 2021) to better confirm the robustness of our final candidate selections. These studies were also of great help to decide on the membership of marginal members according to one of more criteria.

In this work, we improve and update all criteria presented
and described in Paper I (namely, the analysis of RV distribu-
tions, surface gravity indicators, Li content, metallicity and com-
parison with several Gaia studies from the literature), as well as
using the astrometry and photometry from Gaia EDR3 as new
criteria, upon which we further elaborate in Subsects. 3.1 and 3.2
in this section below. We additionally note that in our final analy-
sis we modify the order of the criteria applied in the membership
analysis in some cases, as briefly mentioned above, and we also
improve on several of our methods in comparison to our analy-
sis in Paper I. A complete in-depth description of all updated and
improved criteria in our final membership analysis can be found
in Chapter 2 of the thesis, and we also refer the reader to the indi-
vidual notes of Appendix A, where we offer a detailed discussion
of the membership analysis for all clusters in the sample. Finally,
we also refer to Appendix B, where we show the individual fig-
ures for each cluster corresponding to most of the membership
criteria.

We also note that for all clusters, we identified and discarded
a series of binary stars, both SB1 (single line spectroscopic bi-
naries) and SB2/3/4 (double and multiple line spectroscopic bi-
naries), which can add significant contamination to our analysis.
SB1s were excluded from our kinematic analysis, as they can
strongly affect the observed RV distributions, but we decided to
include those SB1 stars that had not already been discarded after
applying all kinematic and astrometric filters for the rest of our
membership criteria, seeing as Li measurements are not affected
by SB1 binaries, and they could still be of interest for our anal-
ysis. On the other hand, SB2/3/4s were fully discarded from our
data sample for all clusters and membership cirteria. All binary
stars were identified using the iDR6 data release metadata (via

the column *PECULI*) (e.g., Gilmore et al. 2022), as well as ex-
isting studies (Merle et al. 2017, 2020). SB1 and SB2/3/4 stars
are listed in the long tables of Appendix C.

3.1. Proper motions and parallaxes

<figure><img src="image_1.png" /><figcaption>**Fig. 1.** pmra-versus-pmdec proper motion diagram showing the final candidate selection (red squares) for NGC 2516, a 125–138 Myr intermediate-age cluster.</figcaption></figure>

In our earlier cluster calibration work using iDR4 data as
described in Paper I, we made use of several GES membership
studies using Gaia DR1 and DR2 to reinforce our final candidate
selections, thus making use in an indirect way of the precision
and robustness of the data provided by Gaia. One of the main
improvements of the present analysis is to use the proper mo-
tions and parallaxes provided by Gaia EDR3 to identify robust
astrometric candidates among the kinematic selection.

In order to be able to use Gaia data for our cluster sample,
we crossmatched the GES iDR6 files with Gaia EDR3 by mak-
ing use of the CDS Upload X-Match functionality in TOPCAT<sup>13</sup>,
which allows users to join a local table with any table provided
by VizieR<sup>14</sup>, or SIMBAD<sup>15</sup> (Wenger et al. 2000; Taylor 2005).
We chose a search radius of 5 arcsec for all clusters in our sam-
ple, and did not encounter any problems with duplicates. We
then constrained all the Gaia measurements according to a series
of quality indicators in order to ensure that the proper motions,
parallaxes and photometry we used throughout the membership
analysis are of sufficient quality regarding their precision, relia-
bility and consistency (Lindegren 2018; Lindegren et al. 2018).
We here used the renormalized unit weight error (RUWE – Lin-
degren et al. 2018; Gaia Collaboration et al. 2021), only con-
sidering for the astrometric study of each cluster the stars with
RUWE< 1.4 (Lindegren 2018; de la Torre Rojo 2020)<sup>16</sup> (also

<footnote n="12">¹² A tracer of accretion and youth indicator in young PMS stars, Hα10% refers to the width of the Hα emission line at 10% peak intensity.</footnote>

<footnote n="13">¹³ http://www.star.bris.ac.uk/~mbt/topcat/</footnote>

<footnote n="14">¹⁴ https://vizier.cds.unistra.fr/viz-bin/VizieR</footnote>

<footnote n="15">¹⁵ http://simbad.u-strasbg.fr/simbad/</footnote>

<footnote n="16">¹⁶ We refer to the following studies for more details on Gaia DR2 and EDR3 data processing and quality indicators (Cantat-Gaudin et al.</footnote>
see Subsect. 3.2 below for the quality indicators applied to Gaia photometry).

We studied the astrometric goodness of the kinematic members whose associated Gaia data fulfilled all these quality indicators by first plotting them in pmra-versus-pmdec diagrams and analysing the locus of probable members. As reference in order to identify this location for each cluster, we used the estimated proper motions from the literature listed in Table 4. As we can see in the example of Fig. 1 for intermediate-age cluster NGC 2516 (125–138 Myr), all astrometric members can be found clustered together in the expected location in the pmra-versus-pmdec diagram as estimated by the literature. This is one of our most restrictive criteria, and so we discarded from our analysis all the kinematic members that strayed significantly from the expected locus. As a limit, we applied the criterion of studies such as Cantat-Gaudin et al. (2018) to select sources with proper motions within a maximum of 0.5 mas yr⁻¹ of the centroid.

This criterion, combining precise measurements and accurate literature estimates, allowed us to optimally reinforce our list of kinematic members. We were also able to ascertain the astrometric membership of any marginal kinematic stars by studying their position in the proper motions diagrams, fully accepting them as potential candidates if they proved to be robust astrometric members. On the other hand, there were also a number of stars in each cluster to which we were not able to apply astrometric criteria, be it because there were no Gaia data available for those GES stars, or because the crossmatched Gaia data was filtered out with quality indicators. In those cases, we analysed the stars with the same criteria we used in Paper I (kinematics, gravity indicators, metallicity, and Li). We refer the reader to the individual notes of Appendix A for more details on the individual decisions for each cluster in these cases.

The next step in the astrometric analysis was to study the distributions of parallaxes (π) for each cluster. Similarly to the analysis of RV and metallicity distributions (as fully discussed in Paper I and the thesis), we fitted the initial parallax distributions for all clusters in the sample, making use of Gaussian curves, applying a 2σ clipping procedure on the median, and adopting a 2σ limit about the cluster mean yielded by the Gaussian fit to obtain the most likely parallax members. As was also the case with our RV and metallicity selections, we also considered a series of marginal parallax members that fully fulfilled the kinematic and proper motions criteria, and for these borderline members we similarly chose to enlarge the obtained 2σ intervals up to a certain threshold to assess them, typically resulting in a slightly larger 3σ interval (the exact value depends on each cluster, as is detailed in Appendix A).

In Table 4, we list all the mean parallax values, their associated dispersions and the resulting 2σ intervals for all 38 clusters analysed, alongside the estimated parallax values from the literature. The updated results for study of the distributions of RV and metallicity are also summarised in Tables 3 and 5, respectively (see Paper I and Chapter 2 of the thesis). As already discussed in Paper I for the kinematic analysis, the candidate selections for clusters γ Vel, NGC 2547 and NGC 2451 A and B were taken from several prior membership studies from the literature. However, in this case we also improved the final selections for these clusters by discarding several spurious stars that deviated appreciably in regards to our astrometry criteria. As in the case of RV

<figure><img src="image_2.png" /><figcaption>**Fig. 2.** Distributions of parallaxes and Gaussian fits for the intermediate-age cluster NGC 2516 (125–138 Myr). We display the histogram both for sources resulting from the 2σ clipping procedure on all the GES sources in this field (top panel), and for likely cluster members after applying all of our membership criteria (bottom panel).</figcaption></figure>

and metallicity distributions (similarly further discussed in both Paper I and the published thesis), we also fitted the parallax distributions of our final lists of candidates for each cluster, and compared the obtained final central mean parallaxes and their associated dispersions with those available in the literature for each cluster, finding all of our final estimates to be in agreement. Figure 2 additionally shows an example of the parallax distribution analysis for intermediate-age cluster NGC 2516, comparing the initial fit following the 2σ clipping procedure, from which we obtain a preliminary 2σ membership interval, with the final distribution of the parallaxes for the final candidates for the cluster.

We consider astrometry to be one of the most restrictive criteria in our analysis due to its reliability and precision, and thus we fully discarded from our analysis all stars that proved to be non-members according to proper motions, and all stars that fulfilled the proper motions criterion but proved to be non-members according to the study of parallaxes were similarly fully discarded. In contrast to this, other criteria such as gravity indicators and metallicity being less robust, we chose not to make them as restrictive in order to obtain the most probable lists of candidate members (as discussed in more detail in both Paper I and Chapter 2 of the thesis).

<footnote>2018; Gaia Collaboration et al. 2018a,b; Lindegren 2018; Lindegren et al. 2018; Randich et al. 2018; Gaia Collaboration et al. 2021; Riello et al. 2021; Torra et al. 2021)</footnote>
## 3.2. Colour–magnitude diagrams

<figure><img src="image_3.png" /><figcaption>**Fig. 3.** CMD showing the final candidate selection (red squares) of IC 2602, a 35 Myr-old cluster. We overplot the PARSEC isochrones with Z=0.019 for 30 Myr (green curve) and 40 Myr (grey curve).</figcaption></figure>

We made use of the photometry in the *G*, *G*<sub>BP</sub>, and *G*<sub>RP</sub> bands provided by Gaia EDR3 to reinforce the selections of kinematic and astrometric candidates and discard field contaminants by means of *M*<sub>G</sub>-versus-*G*<sub>BP</sub>-*G*<sub>RP</sub> colour–magnitude diagrams (CMDs; e.g., Riello et al. 2021). The absolute magnitudes *M*<sub>G</sub> were calculated using the apparent magnitudes *G* measured by Gaia, as well as the parallaxes (π) measured in milliarcseconds (mas). In addition to the criterion of RUWE< 1.4, already applied to astrometric criteria in Subsect. 3.1, to ensure the most precise dataset possible, we used a series of additional quality indicators to the Gaia photometry, applying the following filters on the relative magnitude errors on the *G*, *G*<sub>BP</sub> and *G*<sub>RP</sub> photometric bands: σ<sub>G</sub> < 0.022 mag, σ<sub>RP</sub> < 0.054 mag, and σ<sub>BP</sub> < 0.054 mag (Gaia Collaboration et al. 2018a; Lindegren et al. 2018; Duque Arribas 2020; de la Torre Rojo 2020; Gaia Collaboration et al. 2021).

The use of Gaia photometry in CMDs offers a more precise and robust way to assess the membership of our astrometric candidates and discard field contaminants than criteria based on spectroscopic gravity indicators, such as the analysis of the potential candidate list in Kiel diagrams (the *T*<sub>eff</sub>-versus-log g plane), or using the γ index for the young clusters in the sample (both criteria are described in detail in Paper I and the thesis). The reason why Kiel diagrams prove to be less reliable is that log g values are generally less precise, and also scarcer in our sample, and so, we consider the analysis of CMDs to be a far more restrictive and reliable criterion than gravity indicators<sup>17</sup> (as further discussed in Paper I and Chapter 2 of the thesis). We encountered several cases, to give an example, where a target star fully fulfilled the CMD criterion but appeared to deviate more appreciably in the Kiel diagram, and in such cases we relied more heavily on our CMD analysis due to its superior precision and reliability. Similarly to the astrometric criteria described above, we fully discarded any astrometric candidates

from the analysis when they deviated significantly from the expected trend in the CMD in a way that could not be explained by the existing inherent dispersion among the cluster members. We note, however, that most of our astrometric selections already being quite robust, these spurious astrometric candidates were not particularly common in our analysis.

For all clusters we made use of the PARSEC CMD isochrones (Bressan et al. 2012; de la Torre Rojo 2020)<sup>18</sup>, choosing the Gaia EDR3 photometric system, for ages ranging from 1 Myr to 5 Gyr, and with a metal fraction of Z=0.019 (except for the very low metallicity old cluster NGC 2243, where we used isochrones with Z=0.006). We also took the interstellar reddening and extinction into account when obtaining the isochrones for each cluster age by applying the corresponding *A*<sub>V</sub> (mag) extinction values. For each cluster, we calculated the *A*<sub>V</sub> extinction in the V band by means of the expression *A*<sub>V</sub> = 3.2*E*(B − V), using the *E*(B − V) reddening given by (Jackson et al. 2021) and listed in Tables 1 and 2. As an example, in Fig. 3 we present the CMD for the cluster IC 2602, a 35 Myr old young cluster, while the CMDs for all clusters in the sample are shown in Appendix B.

## 3.3. Comparison with Gaia studies

As part of the preliminary work on the membership analysis of the 42 OCs in the sample, we did an extensive research on each cluster, and listed all available membership studies from the literature that included lists of candidates. In Appendix A we delve in detail on the comparisons between all these membership studies and our own final cluster selections. Among these, several membership studies were particularly relevant for our work, conducted from Gaia DR1 (Randich et al. 2018), DR2 (Cantat-Gaudin et al. 2018; Soubiran et al. 2018; Bossini et al. 2019) and EDR3 data (Jackson et al. 2021).

We adopted the ages revised by Bossini et al. (2019) for nine of our sample clusters, and the RVs from Soubiran et al. (2018) for 28 of them, as reported in Tables 1 and 2 for the ages, and Table 3 for the RVs, respectively. We also observed that, judging by the measured Gaia ages by Bossini et al. (2019), as well as the empirical Li envelopes constructed by using our cluster candidates (see Paper III) it is possible that some of the former age estimates for the preselected intermediate-age and old clusters could be overestimates - NGC 6005, for example, had a former age estimate of 1.2 Gyr, while Bossini et al. (2019) gives a lower age of 973 ± 5 Myr, which is more in accordance with the Li envelope of our candidate selection. However, we note that we also decided not to use the age estimates by Bossini et al. (2019) for two clusters (NGC 2516 and NGC 6633). The reason for this is that we believe these ages to be overestimates as well, judging by both more recent age estimations, and once again our own candidate selections and obtained empirical Li envelopes (see Appendix A for more details on the sample cluster ages).

We used the other three studies cited above (Cantat-Gaudin et al. 2018; Randich et al. 2018; Jackson et al. 2021) as an additional tool to assess our own selections after concluding all membership analyses and applying all the criteria discussed in this section. All these works were of great aid to confirm and reinforce the robustness of our final candidate selections, and were also markedly useful to aid in the confirmation of marginal members in those cases when the membership criteria were not sufficient to fully confirm their membership to the clusters. For the preliminary work published in Paper I, we mainly made use

<footnote n="17"><sup>17</sup> Similarly to surface gravity indicators, [Fe/H] metallicity is not as robust a criterion as kinematics, astrometry and CMDs, as the [Fe/H] values derived from GIRAFFE spectra are widely dispersed and subject to larger uncertainties</footnote>

<footnote n="18"><sup>18</sup> http://stev.oapd.inaf.it/cgi-bin/cmd</footnote>
of the first two studies, conducted from Gaia DR1 and DR2 (Cantat-Gaudin et al. 2018; Randich et al. 2018), and for the final version of this work extensive use was made of the data in Jackson et al. (2021) in order to assess our final member selections. We also note the interest of later membership studies that were published after the completion of the present work Franciosini et al. (2022a); Prisinzano et al. (2022).

Jackson et al. (2021) combined GES iDR6 spectroscopic data with the astrometry provided by Gaia EDR3 to assign membership probabilities for a target sample of 63 OCs (including 39 out of 42 of the clusters in our sample (all except for NGC 2477, Melotte 71 and M67), as well as 7 globular clusters. Out of 43211 targets, Jackson et al. (2021) listed 13985 as highly probable cluster members, with $P > 0.9$, and an average membership probability of 0.993. Similarly to Cantat-Gaudin et al. (2018), the membership selection is purely kinematic and independent of photometry and chemistry and the final selection catalogues from both studies can be successfully combined with other photometric and spectroscopic criteria from GES.

For each of the clusters considered in these studies, see the individual notes in Appendix A for in-depth details regarding the comparison between the candidates listed in these Gaia studies and our final member selections. In the online tables described in Appendix C we also include for reference these Gaia membership selections alongside the columns listing the results of our membership analysis criteria.

### 3.4. Identification of Li-rich giant contaminants

<figure><img src="image_4.png" /></figure>

**Fig. 4.** EW(Li)-versus-$T_{eff}$ diagram for the Li-rich giant outliers in the field of the young (red squares), intermediate-age (green squares), and old clusters (blue squares). Also shown are the upper envelope of EW(Li) for IC 2602 (35 Myr, red), the upper and lower envelopes of the Pleiades (78–125 Myr, grey), and the upper envelope of the Hyades (750 Myr, turquoise).

As discussed in more detail in Paper I, gravity indicators help identify giant contaminants in the field of the clusters by plotting the sample stars in the Kiel diagram and the $(\gamma, T_{eff})$ plane. Given their interest<fnref n="19" /> (e.g., Smiljanic et al. 2018; Magrini et al. 2021), as a parallel result of the membership analysis we also listed

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Fig. 5.** $\gamma$ index and log g as a function of $T_{eff}$ for the Li-rich giant outliers (red squares) in the field of the young (top), and intermediate-age (green squares), and old (blue squares) clusters (bottom).

some of these outliers for future study, specifically potential Li-rich giants with $A(Li) > 1.5$. We consider as likely giants any source with log g < 3.5 (Spina et al. 2014a, 2017) and/or with $\gamma > 1.01$ (Damiani et al. 2014; Sacco et al. 2015; Casey et al. 2016; Spina et al. 2017). We also consider Li-rich giants to have $T_{eff} < 5200$ K (Casey et al. 2016; Spina et al. 2017) and, in the case of stars in the field of young clusters, a lack of $H\alpha$ emission, given that this is a youth indicator for PMS stars (Casey et al. 2016). In Figs. 4 and 5 we show all Li-rich giant outliers obtained in the field of each of the 42 clusters in our sample as a result of the membership analysis in diagrams of EW(Li), log g, and $\gamma$ as a function of $T_{eff}$, for both the young clusters and the intermediate-age and old clusters in the sample. All Li-rich giants are also listed in Tables 6 and 7 in Sect. 4, as well as in the long tables of Appendix B.

We note that the classification of Li-rich giant stars in this work is only preliminary. We find a large number of potential Li-rich giants in the field of some clusters (e.g., IC 2602) and, while these stars fulfil the adopted criteria ($T_{eff} < 5200$ K and $A(Li) > 1.5$), given the rare nature of these objects, further confirmation would be required to accept them as bona fide Li-rich giants. It is also worth noting that we selected all Li-rich giant out-

<footnote n="19">¹⁹ Comprising 1–2% of FGK giants and supergiants, most are still not well understood, requiring non-standard evolution models to account for the fresh Li detected on their surface (see the thesis and references therein).</footnote>
liers according to the filters on gravity criteria (γ index for young clusters, and log g for intermediate-age and old clusters) detailed above. In the case of several clusters, we do find some inconsistencies when plotting Li-rich giant outliers in CMDs, with the preselected Li-rich giants sometimes appearing to be non-giants according to photometric data. We detail these instances in the individual notes of Appendix A. For the moment, we decided to limit our classification of Li-rich giant contaminants to gravity indicators, as was also done in all cited instances in the literature, but this issue does reinforce the necessity to further confirm the goodness of all selected Li-rich giant outliers in this study.

## 4. Cluster member selections

We finish this section with the final results from the membership analysis of all 42 sample clusters, as summarised in Table 6. For each cluster, we report i) the number of stars from the iDR6 sample observed with both UVES and GIRAFFE; ii) those with measured values of EW(Li); iii) the number of stars selected as candidate members; and iv) the number of Li-rich giant outliers obtained as a parallel result during the membership analysis. Readers are directed to the individual notes of Appendix A, where we offer a detailed and in-depth discussion of the membership analysis for each cluster, as well as commenting on features of interest regarding individual stars in the selection, and comparing our candidate lists with former membership studies. The full tables resulting from our membership analysis, divided into young, intermediate-age and old age ranges, are provided in online form, and described in Appendix C. In these tables we list all membership criteria in our analysis and the final selections of candidate members for each of the 42 OCs analysed.

We also show our final selections in the following figures: Figure 6 shows the EW(Li)-versus-T<sub>eff</sub> diagrams subdivided into young, intermediate-age, and old clusters. These figures also include for reference the upper envelope of EW(Li) for IC 2602 (35 Myr, in red), the upper and lower envelopes of the Pleiades (78–125 Myr, in grey), and the upper envelope of the Hyades (750 Myr, turquoise). In Figure 6 we also show the representative average errors in T<sub>eff</sub> and EW(Li) for all members of the young, intermediate-age and old clusters. These average errors amount to 68 K and 12 mÅ for the young clusters; 69 K and 7 mÅ for the intermediate-age clusters; and 65 K and 8 mÅ for the old clusters. On the other hand, Fig. 7 shows the γ-versus-T<sub>eff</sub> diagram for the young clusters in our sample, as well as the Kiel diagram for all intermediate-age and old clusters in the sample. Additionally, Appendix B shows all the individual figures for the cluster sample, including both candidate members as well as Li-rich contaminants of interest.

Finally, in Table 7 we show some further results of our membership analysis for the 42 sample clusters. As in Table 6, we show the number of stars in the field of each cluster from the initial iDR6 sample and the number of candidate stars for all clusters, as well as the Li-rich giant contaminants. With these results we derived percentages of the candidate members and outlier contaminants, which can be used to rank the target clusters and different age ranges in terms of the percentage of candidate members and contaminants identified in each case. Regarding the cluster members, these percentages are considered firstly with respect to all stars in the field of each cluster, and also with respect to all stars that present Li in the initial sample. However, we only present percentages for the Li-rich giant outliers with respect to all stars in the field, given that these were selected taking A(Li) and not EW(Li) into account.

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** EW(Li)-versus-T<sub>eff</sub> diagrams for the candidate members of the young (1–50 Myr; top panel), intermediate-age (50–700 Myr; middle panel), and old clusters (> 700 Myr; bottom panel). Open squares indicate improbable EW(Li) values for some members. Average errors in T<sub>eff</sub> (K) and EW(Li) (mÅ) are also shown. Reference envelopes for IC 2602, Pleiades and Hyades shown as described in Fig. 4.</figcaption></figure>
<table><caption>**Table 6.** Main results for the 42 OCs analysed, indicating, for each cluster, the number of stars from the sample detected in UVES and GIRAFFE; the number of stars with EW(Li) values; the number of stars selected as candidate members (including RV, astrometric and Li members, as well as the final members); and the number of Li-rich giant contaminants.</caption><thead><tr><th rowspan="2">Cluster<sup>a</sup></th><th colspan="2">UVES</th><th colspan="2">GIRAFFE</th><th rowspan="2">RV</th><th colspan="3">Membership</th><th rowspan="2">Li-rich G outliers</th></tr><tr><th>All stars</th><th>With Li</th><th>All stars</th><th>With Li</th><th>Gaia</th><th>Li</th><th>Final</th></tr></thead><tbody><tr><td>NGC 6530</td><td>52</td><td>5</td><td>1931</td><td>1325</td><td>470</td><td>200</td><td>359</td><td>343</td><td>11</td></tr><tr><td>ρ Oph</td><td>23</td><td>20</td><td>288</td><td>277</td><td>48</td><td>29</td><td>37</td><td>44</td><td>2</td></tr><tr><td>Trumpler 14</td><td>43</td><td>11</td><td>1859</td><td>1045</td><td>228</td><td>55</td><td>165</td><td>159</td><td>11</td></tr><tr><td>Cha I</td><td>47</td><td>36</td><td>660</td><td>623</td><td>102</td><td>54</td><td>88</td><td>87</td><td>6</td></tr><tr><td>NGC 2244</td><td>8</td><td>6</td><td>444</td><td>385</td><td>143</td><td>79</td><td>123</td><td>116</td><td>1</td></tr><tr><td>NGC 2264</td><td>113</td><td>70</td><td>1740</td><td>1539</td><td>621</td><td>423</td><td>507</td><td>503</td><td>10</td></tr><tr><td>λ Ori</td><td>116</td><td>103</td><td>720</td><td>675</td><td>207</td><td>142</td><td>163</td><td>161</td><td>3</td></tr><tr><td>Col 197</td><td>8</td><td>3</td><td>401</td><td>363</td><td>123</td><td>86</td><td>104</td><td>92</td><td>7</td></tr><tr><td>γ Vel</td><td>79</td><td>50</td><td>1183</td><td>1140</td><td>...</td><td>...</td><td>...</td><td>234</td><td>1</td></tr><tr><td>NGC 2232</td><td>47</td><td>28</td><td>1822</td><td>1722</td><td>750</td><td>56</td><td>84</td><td>68</td><td>9</td></tr><tr><td>NGC 2547</td><td>54</td><td>33</td><td>423</td><td>372</td><td>...</td><td>...</td><td>...</td><td>151</td><td>0</td></tr><tr><td>IC 2391</td><td>48</td><td>26</td><td>386</td><td>374</td><td>56</td><td>27</td><td>37</td><td>35</td><td>6</td></tr><tr><td>IC 2602</td><td>131</td><td>89</td><td>1721</td><td>1651</td><td>309</td><td>43</td><td>59</td><td>55</td><td>24</td></tr><tr><td>IC 4665</td><td>34</td><td>30</td><td>533</td><td>514</td><td>233</td><td>29</td><td>51</td><td>33</td><td>0</td></tr><tr><td>NGC 2451 A/B</td><td>90</td><td>70</td><td>1566</td><td>1537</td><td>...</td><td>...</td><td>...</td><td>106</td><td>10</td></tr><tr><td>NGC 6405</td><td>21</td><td>12</td><td>680</td><td>486</td><td>251</td><td>80</td><td>53</td><td>51</td><td>1</td></tr><tr><td>Blanco 1</td><td>37</td><td>31</td><td>426</td><td>373</td><td>142</td><td>119</td><td>101</td><td>98</td><td>0</td></tr><tr><td>NGC 6067</td><td>27</td><td>16</td><td>753</td><td>327</td><td>209</td><td>126</td><td>60</td><td>56</td><td>0</td></tr><tr><td>NGC 6649</td><td>6</td><td>0</td><td>277</td><td>62</td><td>42</td><td>21</td><td>4</td><td>2</td><td>0</td></tr><tr><td>NGC 2516</td><td>51</td><td>32</td><td>708</td><td>645</td><td>460</td><td>378</td><td>379</td><td>376</td><td>0</td></tr><tr><td>NGC 6709</td><td>10</td><td>10</td><td>720</td><td>590</td><td>322</td><td>71</td><td>53</td><td>49</td><td>10</td></tr><tr><td>NGC 6259</td><td>16</td><td>14</td><td>478</td><td>264</td><td>125</td><td>71</td><td>39</td><td>35</td><td>0</td></tr><tr><td>NGC 6705</td><td>49</td><td>31</td><td>1017</td><td>579</td><td>391</td><td>313</td><td>142</td><td>139</td><td>0</td></tr><tr><td>Berkeley 30</td><td>14</td><td>13</td><td>318</td><td>144</td><td>78</td><td>44</td><td>24</td><td>24</td><td>0</td></tr><tr><td>NGC 6281</td><td>16</td><td>7</td><td>304</td><td>214</td><td>82</td><td>38</td><td>23</td><td>23</td><td>1</td></tr><tr><td>NGC 3532</td><td>67</td><td>51</td><td>1094</td><td>809</td><td>518</td><td>411</td><td>323</td><td>384</td><td>11</td></tr><tr><td>NGC 4815</td><td>14</td><td>12</td><td>204</td><td>92</td><td>68</td><td>50</td><td>30</td><td>29</td><td>0</td></tr><tr><td>NGC 6633</td><td>57</td><td>38</td><td>1605</td><td>1463</td><td>617</td><td>35</td><td>62</td><td>67</td><td>9</td></tr><tr><td>NGC 2477</td><td>11</td><td>10</td><td>114</td><td>0</td><td>86</td><td>71</td><td>9</td><td>9</td><td>0</td></tr><tr><td>Trumpler 23</td><td>16</td><td>15</td><td>165</td><td>70</td><td>51</td><td>41</td><td>25</td><td>23</td><td>0</td></tr><tr><td>Berkeley 81</td><td>14</td><td>14</td><td>265</td><td>159</td><td>69</td><td>42</td><td>25</td><td>24</td><td>0</td></tr><tr><td>NGC 2355</td><td>11</td><td>11</td><td>197</td><td>149</td><td>119</td><td>129</td><td>87</td><td>86</td><td>0</td></tr><tr><td>NGC 6802</td><td>13</td><td>13</td><td>184</td><td>69</td><td>77</td><td>51</td><td>36</td><td>32</td><td>2</td></tr><tr><td>NGC 6005</td><td>19</td><td>19</td><td>541</td><td>298</td><td>174</td><td>112</td><td>55</td><td>49</td><td>3</td></tr><tr><td>Pismis 18</td><td>10</td><td>10</td><td>134</td><td>72</td><td>41</td><td>24</td><td>12</td><td>10</td><td>0</td></tr><tr><td>Melotte 71</td><td>9</td><td>9</td><td>111</td><td>0</td><td>71</td><td>64</td><td>4</td><td>4</td><td>0</td></tr><tr><td>Pismis 15</td><td>11</td><td>11</td><td>322</td><td>201</td><td>91</td><td>66</td><td>33</td><td>31</td><td>1</td></tr><tr><td>Trumpler 20</td><td>29</td><td>26</td><td>1184</td><td>404</td><td>451</td><td>367</td><td>116</td><td>104</td><td>0</td></tr><tr><td>Berkeley 44</td><td>7</td><td>7</td><td>86</td><td>73</td><td>39</td><td>33</td><td>31</td><td>30</td><td>0</td></tr><tr><td>NGC 2243</td><td>27</td><td>27</td><td>634</td><td>576</td><td>469</td><td>446</td><td>289</td><td>289</td><td>8</td></tr><tr><td>M67</td><td>42</td><td>36</td><td>95</td><td>85</td><td>109</td><td>89</td><td>96</td><td>96</td><td>0</td></tr></tbody></table>

Notes. (a) Regarding the clusters γ Vel, NGC 2547 and NGC 2451 A and B, we directly used the selections obtained by several studies listed in Table 1.
<table><caption>**Table 7.** Results for the 42 star forming regions and OCs in the sample, indicating, for each cluster: all stars (both UVES and GIRAFFE) from the GES sample and the number of stars detected with measured EW(Li) values, the number of stars selected as candidate members, and Li-rich giant field contaminants. Regarding member stars, we provide their percentages with respect to all GES stars and to those with a EW(Li) measurement in the field of each cluster.</caption><thead><tr><th rowspan="2">Cluster</th><th colspan="3">iDR6 stars</th><th colspan="2">Members</th><th colspan="2">Li-rich giants</th></tr><tr><th>All</th><th>With Li</th><th>#</th><th>%(All)</th><th>%(with Li)</th><th>#</th><th>%(All)</th></tr></thead><tbody><tr><td>NGC 6530</td><td>1983</td><td>1330</td><td>343</td><td>17.3</td><td>25.8</td><td>11</td><td>0.6</td></tr><tr><td>ρ Oph</td><td>311</td><td>297</td><td>44</td><td>13.8</td><td>14.5</td><td>2</td><td>0.6</td></tr><tr><td>Trumpler 14</td><td>1902</td><td>1056</td><td>159</td><td>8.4</td><td>15.1</td><td>11</td><td>0.6</td></tr><tr><td>Cha I</td><td>707</td><td>659</td><td>87</td><td>12.3</td><td>13.2</td><td>6</td><td>0.8</td></tr><tr><td>NGC 2244</td><td>452</td><td>391</td><td>116</td><td>25.7</td><td>29.7</td><td>1</td><td>0.2</td></tr><tr><td>NGC 2264</td><td>1853</td><td>1609</td><td>503</td><td>27.1</td><td>31.3</td><td>10</td><td>0.5</td></tr><tr><td>λ Ori</td><td>836</td><td>778</td><td>161</td><td>19.3</td><td>20.7</td><td>3</td><td>0.4</td></tr><tr><td>Col 197</td><td>409</td><td>366</td><td>92</td><td>22.5</td><td>25.1</td><td>7</td><td>1.7</td></tr><tr><td>γ Vel</td><td>1262</td><td>1190</td><td>234</td><td>18.5</td><td>19.7</td><td>1</td><td>0.1</td></tr><tr><td>NGC 2232</td><td>1869</td><td>1750</td><td>68</td><td>3.6</td><td>3.9</td><td>9</td><td>0.5</td></tr><tr><td>NGC 2547</td><td>477</td><td>405</td><td>151</td><td>31.7</td><td>37.3</td><td>0</td><td>0.0</td></tr><tr><td>IC 2391</td><td>434</td><td>400</td><td>35</td><td>8.1</td><td>8.8</td><td>6</td><td>1.4</td></tr><tr><td>IC 2602</td><td>1852</td><td>1740</td><td>55</td><td>3.0</td><td>3.2</td><td>24</td><td>1.3</td></tr><tr><td>IC 4665</td><td>567</td><td>544</td><td>33</td><td>5.8</td><td>6.1</td><td>0</td><td>0.0</td></tr><tr><td>NGC 2451 A/B</td><td>1656</td><td>1607</td><td>106</td><td>6.4</td><td>6.6</td><td>10</td><td>0.6</td></tr><tr><td>NGC 6405</td><td>701</td><td>498</td><td>51</td><td>7.3</td><td>10.2</td><td>1</td><td>0.1</td></tr><tr><td>Blanco 1</td><td>463</td><td>404</td><td>96</td><td>20.7</td><td>23.8</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6067</td><td>780</td><td>343</td><td>56</td><td>7.2</td><td>16.3</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6649</td><td>283</td><td>62</td><td>2</td><td>0.7</td><td>3.2</td><td>0</td><td>0.0</td></tr><tr><td>NGC 2516</td><td>759</td><td>677</td><td>376</td><td>49.5</td><td>55.5</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6709</td><td>730</td><td>600</td><td>49</td><td>6.7</td><td>8.2</td><td>10</td><td>1.4</td></tr><tr><td>NGC 6259</td><td>494</td><td>278</td><td>35</td><td>7.1</td><td>12.6</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6705</td><td>1066</td><td>610</td><td>139</td><td>13.0</td><td>22.8</td><td>0</td><td>0.0</td></tr><tr><td>Berkeley 30</td><td>332</td><td>157</td><td>24</td><td>7.2</td><td>15.3</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6281</td><td>320</td><td>221</td><td>23</td><td>7.2</td><td>10.4</td><td>1</td><td>0.3</td></tr><tr><td>NGC 3532</td><td>1145</td><td>860</td><td>384</td><td>33.5</td><td>44.7</td><td>11</td><td>1.0</td></tr><tr><td>NGC 4815</td><td>218</td><td>104</td><td>29</td><td>13.3</td><td>27.9</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6633</td><td>1662</td><td>1501</td><td>67</td><td>4.0</td><td>4.5</td><td>9</td><td>0.5</td></tr><tr><td>NGC 2477</td><td>125</td><td>10</td><td>9</td><td>7.2</td><td>90.0</td><td>0</td><td>0.0</td></tr><tr><td>Trumpler 23</td><td>165</td><td>85</td><td>23</td><td>13.9</td><td>27.1</td><td>0</td><td>0.0</td></tr><tr><td>Berkeley 81</td><td>279</td><td>173</td><td>24</td><td>8.6</td><td>13.9</td><td>0</td><td>0.0</td></tr><tr><td>NGC 2355</td><td>208</td><td>160</td><td>86</td><td>41.4</td><td>53.8</td><td>0</td><td>0.0</td></tr><tr><td>NGC 6802</td><td>197</td><td>82</td><td>32</td><td>16.2</td><td>39.0</td><td>2</td><td>1.0</td></tr><tr><td>NGC 6005</td><td>560</td><td>317</td><td>49</td><td>8.8</td><td>15.5</td><td>3</td><td>0.5</td></tr><tr><td>Pismis 18</td><td>142</td><td>82</td><td>10</td><td>7.0</td><td>12.2</td><td>0</td><td>0.0</td></tr><tr><td>Melotte 71</td><td>120</td><td>9</td><td>4</td><td>3.4</td><td>44.5</td><td>0</td><td>0.0</td></tr><tr><td>Pismis 15</td><td>333</td><td>211</td><td>31</td><td>9.3</td><td>14.7</td><td>1</td><td>0.3</td></tr><tr><td>Trumpler 20</td><td>1213</td><td>430</td><td>104</td><td>8.6</td><td>24.2</td><td>0</td><td>0.0</td></tr><tr><td>Berkeley 44</td><td>93</td><td>80</td><td>30</td><td>32.3</td><td>37.5</td><td>0</td><td>0.0</td></tr><tr><td>NGC 2243</td><td>661</td><td>603</td><td>289</td><td>43.7</td><td>48.0</td><td>8</td><td>1.2</td></tr><tr><td>M67</td><td>131</td><td>121</td><td>96</td><td>73.3</td><td>79.3</td><td>0</td><td>0.0</td></tr></tbody></table>

## 5. Summary

This work is a large-scale project, which we started in Paper I (Gutiérrez Albarrán et al. 2020) and concluded in Paper III (Gutiérrez Albarrán et al. 2024, in prep.), with the main objective of studying Li as an age indicator and calibrating a Li–age relation for PMS, ZAMS, and early MS FGK late-type stars. With a considerably expanded sample of 42 OCs spanning a wide age range from 1–3 Myr to 4–4.5 Gyr, in the present work we made

use of the most recent available GES-derived data provided by iDR6, as well as of the high-precision data provided by Gaia EDR3, in order to conduct a thorough membership analysis and obtain updated and expanded selections of candidates for all target clusters. We summarise our analysis and the main results of this work as follows:

- With the present work, we considerably enlarged and improved our target sample, with 42 young, intermediate-age and old clusters provided by GES iDR6 at our disposal, as
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Gravity index $\gamma$ as a function of $T_{\text{eff}}$ for the members for all young clusters of the sample (top panel), and Kiel diagrams for all intermediate-age and old clusters for the candidate members of the intermediate-age (50–700 Myr; middle panel) and old clusters (> 700 Myr; bottom panel). We overplot the PARSEC isochrones in a similar age range, with a metallicity of $Z=0.019$.</figcaption></figure>

compared to our former sample of 20 clusters using data from GES iDR4 from Paper I. The new homogeneous analysis done in iDR6 also offers marked improvements to both our cluster sample and our consequent analysis, including a significantly larger number of stars measured in the field of many of the clusters, the recalculation of parameters for higher accuracy, and accounting for the contamination of background nebular lines (as detailed in Sect. 2).

- We carried out extensive preliminary bibliographical research on each of the 42 sample clusters, compiling a thorough selection of literature data on ages, distances, and reddening values (see Tables 1 and 2 in Sect. 2), as well as RVs (Table 3 in Sect. 3), proper motions, parallaxes (Table 4 in Sect. 3), [Fe/H] metallicity (Table 5 in Sect. 3), previous Li measurements, and numerous existing membership studies (as also presented in Tables 1 and 2, as well as in Appendix A).

- We performed exhaustive and detailed membership analyses (see Sect. 3), firstly studying the RV distributions to obtain likely kinematic candidates, and complementing this with a thorough study of the proper motions and parallaxes provided by Gaia, a marked improvement in regards to our former analyses in Paper I. Gravity indicators such as log g and the $\gamma$ index were helpful in order to discard field giant contaminants and confirm the membership of the astrometric selections, also significantly improved by Gaia photometry in CMDs. We complemented all these criteria with a study of the distribution of [Fe/H] metallicity, which helped reinforce our cluster selections and discard further rogue contaminants. Finally, we used Li as a final criterion by plotting the candidates in EW(Li) versus $T_{\text{eff}}$ diagrams.

- We obtained selections of robust candidates for all 42 sample clusters, some of which had not been previously studied in detail by the GES consortium (see Table 6 in Sect. 4), and we also discuss which individual clusters and age ranges present the highest percentages of members in our target sample in Table 7. All individual figures are displayed in Appendix B, and this paper also includes descriptions for the associated online long tables of results, described in Appendix C, in which we list the final candidates for all clusters and specify which membership criteria are fulfilled by each of the target stars for each cluster.

- We made use of a number of studies from the literature conducted from Gaia DR1, DR2, and EDR3 data to do an in-depth assessment of our candidate selections after concluding our membership analysis, typically obtaining coherent results, as we found our lists of cluster members to be in general agreement with these previous GES studies (see Appendix A). In addition, we also fitted the distributions of RV, parallaxes, and [Fe/H] metallicity values of all final cluster selections, as displayed in Tables 3 (for RVs), 4 (for parallaxes), and 5 (for metallicity), obtaining mean values that are consistent with the literature for the majority of clusters, with very few exceptions.

- Given their importance for our understanding of stellar Li, as an additional result of our membership analysis we also selected a series of preliminary Li-rich giant outliers for 23 out of the 42 target clusters (see Sect. 3.4, as well as Table 6 in Sect. 4 and Appendix C). Given the scarcity of these objects, we require further confirmation to accept all stars listed as Li-rich giant contaminants in this study.

- The selections of cluster candidates obtained by means of recent GES and Gaia data can be used by the scientific community in multiple applications, from the calibration of a Li-age relation —the aim of this particular work—to further
analysis of (the same or other) cluster members, or the char-
acterisation of a variety of stellar parameters and other areas
of interest regarding OCs and the formation and evolution of
stars in the Galaxy.

- In the companion paper to this work, Paper III, we conclude our calibration of an empirical Li–age relation and obtain a set of empirical Li envelopes for several key ages in our sample. This final paper will also describe our future work.

Acknowledgements. We acknowledge financial support from the Universidad Complutense de Madrid (UCM) and the Agencia Estatal de Investigación (AEI/10.13039/501100011033) of the Ministerio de Ciencia e Innovación and the ERDF “A way of making Europe” through projects PID2019-109522GB-C5[4] and PID2022-137241NBC4[4]. We acknowledge the support from INAF and Ministero dell'Istruzione, dell'Università e della Ricerca (MIUR) in the form of the grant "Premiale VLT 2012". T.B. was supported by grant No. 2018-04857 from the Swedish Research Council. J.I.G.H. acknowledges financial support from the Spanish Ministry of Science, Innovation and Universities (MI- CIU) under the 2003 Ramón y Cajal program RYC-2013-14875, and also from the Spanish Ministry project MICIU AYA2017-86389-P. E.J.A. acknowledges financial support from the State Agency for Research of the Spanish MICIU through the "Center of Excellence Severo Ochoa" award to the Instituto de Astrofísica de Andalucá (CEX2021-001131-S). E. M. acknowledges financial support through a "Margarita Salas" postdoctoral fellowship from Universidad Complutense de Madrid (CT18/22), funded by the Spanish Ministerio de Univer- sidades with NextGeneration EU funds. F.J.E. acknowledges support from ESA through the Faculty of the European Space Astronomy Centre (ESAC) - Funding reference 4000139151/22/ES/CM. Based on data products from obser- vations made with ESO Telescopes at the La Silla Paranal Observatory under programme focusID 188.B-3002. These data products have been processed by the Cambridge Astronomy Survey Unit (CASU) at the Institute of Astronomy, University of Cambridge, and by the FLAMES/UVES reduction team at INAF- Osservatorio Astrofísico di Arcetri. These data have been obtained from the GES Data Archive, prepared and hosted by the Wide Field Astronomy Unit, Institute for Astronomy, University of Edinburgh, which is funded by the UK Science and Technology Facilities Council. The results presented here benefit from dis- cussions held during GES workshops and conferences supported by the ESF (European Science Foundation) through the GREAT Research Network Pro- gramme. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www. cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement. This publication makes use of the VizieR database (Ochsenbein et al. 2000) and the SIMBAD database (Wenger et al. 2000), both operated at CDS, Centre de Données astronomiques de Strasbourg, France. This research also made use of the WEBDA database, operated at the Department of Theoretical Physics and Astrophysics of the Masaryk University, and the interactive graphical viewer and editor for tabular data TOPCAT (Tay- lor 2005). For the analysis of the distributions of RV and metallicity we used RStudio Team (2015). Integrated Development for R. RStudio, Inc., Boston, MA (http://www.rstudio.com/). Finally, we would like to thank the anonymous referee for helpful comments and suggestions.

References

Alonso-Santiago, J., Negueruela, I., Marco, A., et al. 2017, MNRAS, 469, 1330
Alonso-Santiago, J., Negueruela, I., Marco, A., et al. 2020, A&A, 644, A136
Arancibia-Silva, J., Bouvier, J., Bayo, A., et al. 2020, A&A, 635, L13
Armstrong, J. J., Wright, N. J., Jeffries, R. D., et al. 2020, MNRAS, 494, 4794
Bailey, J. I., Mateo, M., White, R. J., et al. 2018, MNRAS, 475, 1609
Balachandran, S. 1995, ApJ, 446, 203
Balog, Z., Kiss, L. L., Vinkó, J., et al. 2009, ApJ, 698, 1989
Barrado, D., Stelzer, B., Morales-Calderón, M., et al. 2011, A&A, 526, A21
Barrado y Navascués, D. 2004, A&A, 422, 631
Barrado y Navascués, D., Stauffer, J. R., Briceño, C., et al. 2001, ApJS, 134, 103
Barsony, M., Haisch, K. E., Marsh, K. A., et al. 2012, ApJ, 751, 22
Bayo, A., Barrado, D., Huélamo, N., et al. 2012, A&A, 547, A80
Beccari, G., Boffin, H. M. J., Jerabkova, T., et al. 2018, MNRAS, 481, L11
Binks, A. S., Jeffries, R. D., Sacco, G. G., et al. 2022, MNRAS, 513, 5727
Bonatto, C. & Bica, E. 2010, A&A, 516, A81
Bonito, R., Prisinzano, L., Guarcello, M. G., et al. 2013, A&A, 556, A108
Bonito, R., Prisinzano, L., Venuti, L., et al. 2020, A&A, 642, A56
Bossini, D., Vallenari, A., Braglia, A., et al. 2019, A&A, 623, A108
Bouvier, J., Barrado, D., Moraux, E., et al. 2018, A&A, 613, A63

Braglia, A., Alfaro, E., Flaccomio, E., et al. 2022, A&A, 659, A200
Bravi, L., Zari, E., Sacco, G. G., et al. 2018, A&A, 615, A37
Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127
Brucalassi, A., Koppenhoefer, J., Saglia, R., et al. 2017, A&A, 603, A85
Buckner, A. S. M. & Froebrich, D. 2014, MNRAS, 444, 290
Cánovas, H., Cantero, C., Cieza, L., et al. 2019, A&A, 626, A80
Cantat-Gaudin, T., Jordi, C., Vallenari, A., et al. 2018, A&A, 618, A93
Cantat-Gaudin, T., Mapelli, M., Balaguer-Núñez, L., et al. 2019, A&A, 621, A115
Cargile, P. A., James, D. J., & Jeffries, R. D. 2010, ApJ, 725, L111
Carlberg, J. K. 2014, AJ, 147, 138
Carraro, G., Geisler, D., Baume, G., et al. 2005, MNRAS, 360, 655
Carrera, R., Braglia, A., Cantat-Gaudin, T., et al. 2019, A&A, 623, A80
Casali, G., Magrini, L., Tognelli, E., et al. 2019, A&A, 629, A62
Casamiquela, L., Carrera, R., Balaguer-Núñez, L., et al. 2018, A&A, 610, A66
Casey, A. R., Ruchti, G., Masseron, T., et al. 2016, MNRAS, 461, 3336
Castro-Ginard, A., Jordi, C., Luri, X., et al. 2020, A&A, 635, A45
Conrad, C., Scholz, R. D., Kharchenko, N. V., et al. 2017, A&A, 600, A106
Damiani, F., Klutsch, A., Jeffries, R. D., P., et al. 2017, MNRAS, 603, A81
Damiani, F., Prisinzano, L., Micela, G., et al. 2014, A&A, 566, A50
Damiani, F., Prisinzano, L., Micela, G., et al. 2019, A&A, 623, A25
de la Torre Rojo, A. 2020, Master's thesis in Astrophysics. UCM, Universidad Complutense de Madrid. 2019-2020. msc.ucm.delatorrojo.pdf
De Silva, G. M., D'Orazi, V., Melo, C., et al. 2013, MNRAS, 413, 1005
de Wit, W. J., Bouvier, J., Palla, F., et al. 2006, A&A, 448, 189
Dias, W. S., Alessi, B. S., Moitinho, A., et al. 2002, A&A, 389, 871
Dias, W. S., Monteiro, H., Lépine, J. R. D., et al. 2019, MNRAS, 486, 5726
Dib, S., Schmeja, S., & Parker, R. J. 2018, MNRAS, 473, 849
Dobbie, P. D., Day-Jones, A., Williams, K. A., et al. 2012, MNRAS, 423, 2815
Dodd, R. J. 2004, MNRAS, 355, 959
Donati, P., Beccari, G., Braglia, A., et al. 2014a, MNRAS, 437, 1241
Donati, P., Cantat Gaudin, T., Braglia, A., et al. 2014b, A&A, 561, A94
Ducourant, C., Teixeira, R., Krone-Martins, A., et al. 2017, A&A, 597, A90
Dumont, T., Charbonnel, C., Palacios, A., et al. 2021a, A&A, 654, A46
Dumont, T., Palacios, A., Charbonnel, C., et al. 2021b, A&A, 646, A48
Duque Arribas, C. 2020, Master's thesis in Astrophysics. UCM, Universidad Complutense de Madrid. 2019-2020. msc.ucm.duquearribas.pdf
Eigenbrod, A., Mermilliod, J. C., Clariá, J. J., et al. 2004, A&A, 423, 189
Elliott, P., Bayo, A., Melo, C. H. F., et al. 2016, A&A, 590, A13
Esplin, T. L., Luhman, K. L., Faherty, J. K., et al. 2017, AJ, 154, 46
Franciosini, E., Randich, S., de Laverny, P., et al. 2022a, A&A, 668, A49
Franciosini, E. & Sacco, G. G. 2011, A&A, 530, A150
Franciosini, E., Sacco, G. G., Jeffries, R. D., et al. 2018, A&A, 616, L12
Franciosini, E., Tognelli, E., Degl'Innocenti, S., et al. 2022b, A&A, 659, A85
Frasca, A., Biazzo, K., Lanzafame, A. C., et al. 2015, A&A, 575, A4
Friel, E. D., Donati, P., Braglia, A., et al. 2014, A&A, 563, A117
Friel, E. D. & Janes, K. A. 1993, A&A, 267, 75
Frinchaboy, P. M. & Majewski, S. R. 2008, AJ, 136, 118
Fritzewski, D. J., Barnes, S. A., James, D. J., et al. 2019, A&A, 622, A110
Fritzewski, D. J., Barnes, S. A., James, D. J., et al. 2020, A&A, 641, A51
Fritzewski, D. J., Barnes, S. A., James, D. J., et al. 2021, A&A, 652, A60
Gaia Collaboration, Busiaux, C., van Leeuwen, F., et al. 2018a, A&A, 616, A10
Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2016a, A&A, 595, A2
Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018b, A&A, 616, A1
Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2021, A&A, 649, A1
Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016b, A&A, 595, A1
Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1
Galli, P. A. B., Bouy, H., Olivares, J., et al. 2021, A&A, 646, A46
Gao, X. 2018, AJ, 156, 121
Geller, A. M., Latham, D. W., & Mathieu, R. D. 2015, AJ, 150, 97
Gillen, E., Hillenbrand, L. A., Stauffer, J., et al. 2020, MNRAS, 495, 1531
Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25
Gilmore, G., Randich, S., et al. 2022, A&A, 666, A120
Gómez Garrido, M. 2015, Master's thesis in Astrophysics. UCM, Universidad Complutense de Madrid. 2014-2015. eprint/33629/
Gómez Garrido, M., Montes, D., Gutiérrez Albarrán, M. L., et al. 2017, Highlights on Spanish Astrophysics IX, Proceedings of the XII Scientific Meeting of the Spanish Astronomical Society held on July 18-22, 2016, in Bilbao, Spain, ISBN 978-84-617-8931-3, p. 506
Gutiérrez Albarrán, M. L. 2022, PhD thesis in Astrophysics. UCM, Universidad Complutense de Madrid. 2014-2022. 2022PhDT........6G
Gutiérrez Albarrán, M. L., Montes, D., Gómez Garrido, M., et al. 2020, A&A, 643, A71
Gutiérrez Albarrán, M. L., Montes, D., Tabernero, H. M., et al. 2024, in prep., A&A
Hatzidimitriou, D., Held, E. V., Tognelli, E., et al. 2019, A&A, 626, A90
Hayes, C. R. & Friel, E. D. 2014, AJ, 147, 69
Heiter, U., Soubiran, C., Netopil, M., et al. 2014, A&A, 561, A93
Hernández, J., Morales-Calderón, M., Calvet, N., et al. 2010, ApJ, 722, 1226
Hetem, A. & Gregorio-Hetem, J. 2019, MNRAS, 490, 2521

Hobbs, L. M. & Pilachowski, C. 1986, ApJ, 311, L37

Hourihane, A., Francois, P., Worley, C. C., et al. 2023, A&A, 676, A129

Hünsch, M., Randich, S., Hempel, M., et al. 2004, A&A, 418, 539

Ilin, E., Schmidt, S. J., Poppenhäger, K., et al. 2021, A&A, 645, A42

Irwin, J., Hodgkin, S., Aigrain, S., et al. 2007, MNRAS, 377, 741

Jackson, R. J. & Jeffries, R. D. 2010, MNRAS, 407, 465

Jackson, R. J., Jeffries, R. D., Randich, S., et al. 2016, A&A, 586, A52

Jackson, R. J., Jeffries, R. D., Wright, N. J., et al. 2021, MNRAS, 509, 1664

Jacobson, H. R., Friel, E. D., Jílková, L., et al. 2016, A&A, 591, A37

Jacobson, H. R., Friel, E. D., & Pilachowski, C. A. 2011, AJ, 141, 58

Jadhav, V. V., Pennock, C. M., Subramaniam, A., et al. 2021, MNRAS, 503, 236

Jeffries, R. D. 1997, MNRAS, 292, 177

Jeffries, R. D. 2014, 65, 289

Jeffries, R. D., Jackson, R. J., & Binks, A. S. 2023, MNRAS, 526, 1260

Jeffries, R. D., Jackson, R. J., Cottaar, M., et al. 2014, A&A, 563, A94

Jeffries, R. D., Jackson, R. J., James, D. J., et al. 2009, MNRAS, 400, 317

Jeffries, R. D., James, D. J., & Thurston, M. R. 1998, MNRAS, 300, 550

Jeffries, R. D., Thurston, M. R., & Hambly, N. C. 2001, A&A, 375, 863

Jeffries, R. D., Totten, E. J., Harmer, S., et al. 2002, MNRAS, 336, 1109

Jones, B. F., Fischer, D., & Soderblom, D. R. 1999, AJ, 117, 330

Joshi, Y. C., Dambis, A. K., Pandey, A. K., et al. 2016, ApJ, 593, A116

Juarez, A. J., Cargile, P. A., James, D. J., et al. 2014, ApJ, 795, 143

Kang, W. & Lee, S.-G. 2012, MNRAS, 425, 3162

Kharchenko, N. V., Piskunov, A. E., Röser, S., et al. 2005, A&A, 438, 1163

Kılçoğlu, T., Monier, R., Richer, J., et al. 2016, AJ, 151, 49

Kiman, R., Faherty, J. K., Cruz, K. L., et al. 2021, AJ, 161, 277

Kounkel, M. & Covey, K. 2019, AJ, 158, 122

Kuhn, M. A., Hillenbrand, L. A., Sills, A., et al. 2019, ApJ, 870, 32

Lanzafame, A. C., Frasca, A., Damiani, F., et al. 2015, A&A, 576, A80

Lim, B., Sung, H., Kim, J. S., et al. 2016, ApJ, 831, 116

Lindegren, L. 2018, Proceedings from the third Gaia DPAC consortium meeting. GAIA-C3-TN-LU-LL-124-01

Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2

Liu, L. & Pang, X. 2019, ApJS, 245, 32

Lodieu, N., de Wit, W. J., Carraro, G., et al. 2011, A&A, 532, A103

López Martí, B., Jiménez-Esteban, F., Bayo, A., et al. 2013, A&A, 556, A144

López-Santiago, J., Montes, D., Fernández-Figueroa, M. J., et al. 2003, A&A, 411, 489

Luhman, K. L., Allen, L. E., Allen, P. R., et al. 2008, ApJ, 675, 1375

Luhman, K. L. & Muench, A. A. 2008, ApJ, 684, 654

Lyubimkov, L. S. 2016, Astrophysics, 59, 411

Magrini, L., Randich, S., Donati, P., et al. 2015, A&A, 580, A85

Magrini, L., Randich, S., Kordopatis, G., et al. 2017, A&A, 603, A2

Magrini, L., Randich, S., Romano, D., et al. 2014, A&A, 563, A44

Magrini, L., Smiljanic, R., Lagarde, N., et al. 2021, The Messenger, 185, 18

Magrini, L., Vincenzo, F., Randich, S., et al. 2018, A&A, 618, A102

Maíz Apellániz, J. 2019, A&A, 630, A119

Manzi, S., Randich, S., de Wit, W. J., et al. 2008, A&A, 479, 141

Martín, E. L. & Montes, D. 1997, A&A, 318, 805

Martos, G., Meléndez, J., Rathsam, A., & Carvalho Silva, G. 2023, MNRAS, 522, 3217

Mel'nik, A. M. & Dambis, A. K. 2017, MNRAS, 472, 3887

Merle, T., Van der Swaelmen, M., Van Eck, S., et al. 2020, A&A, 635, A155

Merle, T., Van Eck, S., Jorissen, A., et al. 2017, A&A, 608, A95

Messina, S., Desidera, S., Lanzafame, A. C., et al. 2011, A&A, 532, A10

Michalska, G. 2019, MNRAS, 487, 3505

Montes, D., López-Santiago, J., Fernández-Figueroa, M. J., et al. 2001, A&A, 379, 976

Mužić, K., Scholz, A., Peña Ramírez, K., et al. 2019, ApJ, 881, 79

Netopil, M. & Paunzen, E. 2013, A&A, 557, A10

Netopil, M., Paunzen, E., Heiter, U., et al. 2016, A&A, 585, A150

Neuhaeuser, R., Torres, G., Sterzik, M. F., et al. 1997, A&A, 325, 647

Ochsenbein, F., Bauer, P., & Marcout, J. 2000, A&AS, 143, 23

Oh, S., Price-Whelan, A. M., Brewer, J. M., et al. 2018, ApJ, 854, 138

Oliveira, J. M., Jeffries, R. D., Devey, C. R., et al. 2003, MNRAS, 342, 651

Overbeek, J. C., Friel, E. D., Donati, P., et al. 2017, A&A, 598, A68

Pace, G., Castro, M., Meléndez, J., et al. 2012, A&A, 541, A150

Pallavicini, R., Pasquini, L., & Randich, S. 1997, Mem. Soc. Astron. Italiana, 68, 933

Pallavicini, R., Randich, S., Giampapa, M., & Cutispoto, G. 1990, The Messenger, 62, 51

Pallavicini, R., Randich, S., & Giampapa, M. S. 1992, A&A, 253, 185

Pallavicini, R., Randich, S., & Sestito, P. 2005, 560, 867

Pang, X., Li, Y., Yu, Z., et al. 2021, ApJ, 912, 162

Pasquini, L., Brucalassi, A., Ruiz, M. T., et al. 2012, A&A, 545, A139

Pasquini, L., Randich, S., & Pallavicini, R. 1997, A&A, 325, 535

Paunzen, E., Heiter, U., Netopil, M., et al. 2010, A&A, 517, A32

Platais, I., Girard, T. M., Vieira, K., et al. 2011, MNRAS, 413, 1024

Platais, I., Melo, C., Mermilliod, J. C., et al. 2007, A&A, 461, 509

Postnikova, E. S., Chupina, N. V., & Vereshchagin, S. V. 2019, INASAN Science Reports, 3, 336

Postnikova, E. S., Elsanhoury, W. H., Sariya, D. P., et al. 2020, Research in Astronomy and Astrophysics, 20, 016

Prisinzano, L., Damiani, F., Kalari, V., et al. 2019, A&A, 623, A159

Prisinzano, L., Damiani, F., Micela, G., et al. 2005, A&A, 430, 941

Prisinzano, L., Damiani, F., Micela, G., et al. 2007, A&A, 462, 123

Prisinzano, L., Damiani, F., Micela, G., et al. 2016, A&A, 589, A70

Prisinzano, L., Damiani, F., Sciortino, S., et al. 2022, A&A, 664, A175

Prisinzano, L., Micela, G., Sciortino, S., et al. 2012, A&A, 546, A9

Rain, M. J., Ahumada, J. A., & Carraro, G. 2021, A&A, 650, A67

Randich, S., Aharpour, N., Pallavicini, R., et al. 1997, A&A, 323, 86

Randich, S., Gilmore, G., & Gaia-ESO Consortium. 2013, The Messenger, 154, 47

Randich, S., Gilmore, G., et al. 2022, A&A, 666, A121

Randich, S. & Magrini, L. 2021, Frontiers in Astronomy and Space Sciences, 8, 6

Randich, S., Pallavicini, R., Meola, G., et al. 2001, A&A, 372, 862

Randich, S., Pasquini, L., Franciosini, E., et al. 2020, A&A, 640, L1

Randich, S., Primas, F., Pasquini, L., et al. 2002, A&A, 387, 222

Randich, S., Tognelli, E., Jackson, R., et al. 2018, A&A, 612, A99

Rangwal, G., Yadav, R. K. S., Durgapal, A., et al. 2019, MNRAS, 490, 1383

Richer, H. B., Fahlman, G. G., Rosvick, J., et al. 1998, ApJ, 504, L91

Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3

Rigliaco, E., Wilking, B., Meyer, M. R., et al. 2016, A&A, 588, A123

Robrade, J. & Schmitt, J. H. M. M. 2007, A&A, 461, 669

Roccatagliata, V., Sacco, G. G., Franciosini, E., et al. 2018, A&A, 617, L4

Romano, D., Magrini, L., Randich, S., et al. 2021, A&A, 653, A72

Sacco, G. G., Franciosini, E., Randich, S., et al. 2008, A&A, 488, 167

Sacco, G. G., Jeffries, R. D., Randich, S., et al. 2015, A&A, 574, L7

Sacco, G. G., Spina, L., Randich, S., et al. 2017, A&A, 601, A97

Sampedro, L., Dias, W. S., Alfaro, E. J., et al. 2017, MNRAS, 470, 3937

Sestito, P. & Randich, S. 2005, A&A, 442, 615

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593

Silaj, J. & Landstreet, J. D. 2014, A&A, 566, A132

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Smiljanic, R., Franciosini, E., Bragliał, A., et al. 2018, A&A, 617, A4

Smiljanic, R., Franciosini, E., Randich, S., et al. 2016, A&A, 591, A62

Smiljanic, R., Korn, A. J., Bergemann, M., et al. 2014, A&A, 570, A122

Smiljanic, R., Randich, S., & Pasquini, L. 2011, A&A, 535, A75

Sneden, C., Bean, J., Ivans, I., Lucatello, S., & Sobeck, J. 2012, "MOOG: LTE line analysis and spectrum synthesis"

Soderblom, D. R. 2010, ARA&A, 48, 581

Soderblom, D. R., Hillenbrand, L. A., Jeffries, R. D., et al. 2014, "Ages of Young Stars", Protostars and Planets VI, H. Beuther, R. S. Klessen, C. P. Dullemond, and T. Henning (eds.), University of Arizona Press, Tucson, 219

Soderblom, D. R., Jones, B. F., Balachandran, S., et al. 1993, AJ, 106, 1059

Soubiran, C., Cantat-Gaudin, T., Romero-Gómez, M., et al. 2018, A&A, 619, A155

Spina, L., Randich, S., Magrini, L., et al. 2017, A&A, 601, A70

Spina, L., Randich, S., Palla, F., et al. 2014a, A&A, 568, A2

Spina, L., Randich, S., Palla, F., et al. 2014b, A&A, 567, A55

Stauffer, J. R., Hartmann, L. W., Prosser, C. F., et al. 1997, ApJ, 479, 776

Subramaniam, A. & Sagar, R. 1999, AJ, 117, 937

Sutherland, W., Emerson, J., Dalton, G., et al. 2015, A&A, 575, A25

Tabernero, H. M., Marfil, E., Montes, D., et al. 2019, A&A, 628, A131

Tang, B., Geisler, D., Friel, E., et al. 2017, A&A, 601, A56

Tautvaišienė, G., Drazdauskas, A., Mikolaitis, Š., et al. 2015, A&A, 573, A55

Taylor, M. B. 2005, 347, 29

Terndrup, D. M., Pinsonneault, M., Jeffries, R. D., et al. 2002, ApJ, 576, 950

Torra, F., Castañeda, J., Fabricius, C., et al. 2021, A&A, 649, A10

Umezu, M. & Saio, H. 2000, MNRAS, 316, 307

Vande Putte, D., Garnier, T. P., Ferreras, I., et al. 2010, MNRAS, 407, 2109

Venuti, L., Damiani, F., & Prisinzano, L. 2019, A&A, 621, A14

Vereshchagin, S. V., Tutukov, A. V., & Chupina, N. V. 2019, INASAN Science Reports, 3, 302

Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9

Wright, N. J., Drake, J. J., Mamajek, E. E., et al. 2011, ApJ, 743, 48

Wright, N. J., Jeffries, R. D., Jackson, R. J., et al. 2019, MNRAS, 486, 2477

Yong, D., Liu, F., Ting, Y.-S., et al. 2023, MNRAS, 526, 2181

Zhang, Y., Tang, S.-Y., Chen, W. P., et al. 2020, ApJ, 889, 99
<footnote n="1">¹ Departamento de Física de la Tierra y Astrofísica and IPARCOS-UCM (Instituto de Física de Partículas y del Cosmos de la UCM), Facultad de Ciencias Físicas, Universidad Complutense de Madrid, E-28040, Madrid, Spain e-mail: mlgutierrez@ucm.es</footnote>

<footnote n="2">² Centro de Astrobiología (CSIC/INTA), Instituto Nacional de Técnica Aeroespacial, Ctra de Torrejón a Ajalvir, km 4, 28850 Torrejón de Ardoz, Madrid, Spain</footnote>

<footnote n="3">³ Instituto de Astrofísica de Canarias (IAC), E-38205 La Laguna, Tenerife, Spain</footnote>

<footnote n="4">⁴ Universidad de la Laguna, Dept. Astrofísica, E-38206 La Laguna, Tenerife, Spain</footnote>

<footnote n="5">⁵ Hamburger Sternwarte, Gojenbergsweg 112, 21029 Hamburg, Germany</footnote>

<footnote n="6">⁶ INAF - Osservatorio Astrofisico di Catania, via S. Sofia, 78, 95123 Catania, Italy</footnote>

<footnote n="7">⁷ Università di Catania, Dipartimento di Fisica e Astronomia, Sezione Astrofisica, Via S. Sofia 78, I-95123 Catania, Italy</footnote>

<footnote n="8">⁸ Institut für Astronomie und Astrophysik, Eberhard Karls Universität, Sand 1, D-72076 Tübingen, Germany</footnote>

<footnote n="9">⁹ INAF - Osservatorio Astrofisico di Arcetri, Largo E. Fermi 5, 50125, Firenze, Italy</footnote>

<footnote n="10">¹⁰ Nicolaus Copernicus Astronomical Center, Polish Academy of Sciences, ul. Bartycka 18, 00-716, Warsaw, Poland</footnote>

<footnote n="11">¹¹ Observational Astrophysics, Division of Astronomy and Space Physics, Department of Physics and Astronomy, Uppsala University, Box 516, SE-751 20 Uppsala, Sweden</footnote>

<footnote n="12">¹² Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, United Kingdom</footnote>

<footnote n="13">¹³ Instituto de Astrofísica de Andalucía (CSIC), Glorieta de la Astronomía s/n, Granada 18008, Spain</footnote>

<footnote n="14">¹⁴ Centro de Astrobiología, CSIC-INTA, Camino bajo del castillo s/n, E-28692, Villanueva de la Canàda, Madrid, Spain</footnote>

<footnote n="15">¹⁵ INAF - Osservatorio Astronomico di Palermo, Piazza del Parla-mento, 1, 90134 Palermo, Italy</footnote>

<footnote n="16">¹⁶ Lund Observatory, Division of Astrophysics, Department of Physics, Lund University, Box 43, SE-22100 Lund, Sweden</footnote>

<footnote n="17">¹⁷ INAF - Rome Astronomical Observatory (OAR), Via di Frascati, 33, I-00044, Monte Porzio Catone, Italy</footnote>

<footnote n="18">¹⁸ INAF - Osservatorio Astronomico di Padova, Vicolo dell'Osservatorio 5, 35122 Padova, Italy</footnote>

<footnote n="19">¹⁹ School of Physical and Chemical Sciences — Te Kura Matū, University of Canterbury, Private Bag 4800, Christchurch 8140, New Zealand</footnote>

<footnote n="20">²⁰ 114, 10 College Walk, School of Physics and Astronomy, Monash University VIC 3800, Australia</footnote>

<footnote n="21">²¹ European Space Agency (ESA), European Space Research and Technology Centre (ESTEC), Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands</footnote>

<footnote n="22">²² GEPI, Observatoire de Paris, CNRS, Université Paris Diderot, 5 Place Jules Janssen, 92190 Meudon, France</footnote>

<footnote n="23">²³ Dipartimento di Fisica e Astronomia, Università di Padova, Vicolo dell'Osservatorio 3, 35122 Padova, Italy</footnote>
# Appendix A: Cluster selections: Individual notes

In this appendix we present an in depth account of each of the 42 clusters in the sample, including general information about the different age estimations of the cluster sample, as well as a detailed discussion on membership analysis and final selections from Sect. 3, especially featuring relevant particular cases, as well as a comparison of the final candidates with a series of studies from the literature. In Paper III, we further comment on particular cases in regards to the study of rotation, activity and metallicity. We additionally refer to the individual figures of Sect. 3 in Appendix B.

## Appendix A.1: Star forming regions (age ≤ 6 Myr) and young OCs (age ≤ 50 Myr)

### NGC 6530

NGC 6530 is a 1–2 Myr star forming region (Prisinzano et al. 2005; Damiani et al. 2019; Wright et al. 2019; Randich et al. 2020; Jackson et al. 2021) associated to the HII region M8 (Lagoon Nebula).

Of the final 359 Li candidates for NGC 6530 we discarded 16 stars according to the following criterion: We decided to discard all stars which we could not analyse using the astrometric and photometric Gaia data (be it because there were not Gaia data available for those stars, or because said data did not pass the quality indicators) if they were additionally listed as clear non-members by Jackson et al. (2021). All of these 16 stars are listed as clear non-members with a probability of 0.00–0.16, and one of them also has a metallicity value deviating appreciably from the mean of the cluster. As a result of the analysis of NGC 6530 we found 14 strong accreting stars, all of which are Li members and classified as high-probability members ($P > 0.90$) by Jackson et al. (2021). As for kinematics, 44 stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating up to $13.0–15.0 \text{ km s}^{-1}$ from the mean of the cluster. We accepted the stars which deviate from the limits established by the $2\sigma$ interval by ~10 km s$^{-1}$ given that they fulfil the rest of membership criteria (especially regarding available Gaia data, as well as gravity indicators and lithium) (see Paper I). All of these stars are additionally listed as high-probability members by several studies (e.g., Damiani et al. 2019; Prisinzano et al. 2019; Wright et al. 2019; Jackson et al. 2021). On the other hand, we also accepted as final candidates four stars with no measured RV values in iDR6, due to the fact that they fully fulfil the astrometry criteria (as well as the other criteria), and are also listed as members by these studies (e.g., Damiani et al. 2019; Prisinzano et al. 2019; Wright et al. 2019; Jackson et al. 2021). Finally, regarding metallicity, 26 stars in our final selection are marginal members which deviate moderately from the limits of the initial $3\sigma$ membership interval, with values up to 0.40–0.80 dex from the mean of the cluster. As discussed in Sect. 3 (also see Paper I), we accepted these marginal stars as final cluster members, seeing as they fulfil all prior criteria (including the more robust criteria of kinematics and astrometry), and all of these stars are also listed as high-probability members by various studies (e.g., Damiani et al. 2019; Prisinzano et al. 2019; Wright et al. 2019; Jackson et al. 2021). In addition, five of the stars in our final selection additionally present [Fe/H] values which deviate more appreciably from the $3\sigma$ limit, up to 0.93–1.68 dex from the cluster mean. As further discussed in both Paper I and the updated analysis presented in the thesis, however, we also listed these stars as final candidates, given that they fully fulfil the rest of criteria, including the more restric-

tive criteria, and they are all additionally listed as members by Jackson et al. (2021) with $P=0.99$.

In regards to previous selections, we found the following number of common stars with our selection: 52 stars in Prisinzano et al. (2007) and Prisinzano et al. (2012), 148 in Spina et al. (2017), 255 in Damiani et al. (2019), 320 in Prisinzano et al. (2019), 286 in Wright et al. (2019), and one in Castro-Ginard et al. (2020). Finally, Jackson et al. (2021) includes 331 out of our 343 candidates. Of the remaining 12 stars, 10 of them were not included in their analysis, and two of them are listed with probabilities of 0.26–0.33 (Jackson et al. (2021) considers high-probability members to have $P > 0.90$ and definite non-members to have $P< 0.10$). We accepted these two stars given that they fulfilled all of our available membership criteria, and were furthermore included in several other studies (e.g., Damiani et al. 2019; Prisinzano et al. 2019; Wright et al. 2019). As already discussed in Subsect. 3.4, we also note that we seem to find some inconsistencies when plotting the Li-rich giant outliers selected in the field of this cluster in the CMD diagram, with a couple of them appearing among the cluster candidates. However, when plotting them in the $\gamma$ index-versus-$T_{\text{eff}}$ diagram the expected clear distinction is found between the non-giant cluster candidates and the Li-rich giant outliers.

### Rho Oph

The $\rho$ Ophiuchi ($\rho$ Oph) molecular cloud complex, consisting of two major regions of dense gas and dust (L1688 and L1689), is an 1–3 Myr star forming region (Barsony et al. 2012; Rigliaco et al. 2016; Spina et al. 2017; Cánovas et al. 2019; Randich et al. 2020; Kiman et al. 2021; Jackson et al. 2021).

In addition to the final 37 Li members, for $\rho$ Oph we also considered for completeness seven strong accreting stars which presented no values of EW(Li) in the iDR6 file. We also note that one of the stars in our final selection (16270456-2442140) presents a [Fe/H] value which deviates appreciably from the $3\sigma$ limit (0.45 dex from the cluster mean). We listed it as a final candidate given that it fully fulfils all the other criteria, and is additionally listed as a member by Cánovas et al. (2019) and Jackson et al. (2021) (with $P=0.99$). Of the final members in $\rho$ Oph, 29 belong to the star forming region L1688 (Rigliaco et al. 2016), and 15 of them are also strong accretors with $H\alpha 10\% > 270–300 \text{ km s}^{-1}$. Of these 15 accreting stars, only one (16273311-2441152) does not pass our gravity criteria. This star, with $\gamma=1.022$ and $A(\text{Li})=3.2$, could be listed as a potential Li-rich giant contaminant, but due to the fact that it is also a strong accretor (with $H\alpha 10\% = 469.13 \text{ km s}^{-1}$, we counted it as a likely member of $\rho$ Oph, also in accordance with Rigliaco et al. (2016) and Jackson et al. (2021). Regarding other membership studies from the literature, Rigliaco et al. (2016) has 42 of our final candidates in their selection, while Ducourant et al. (2017) lists 10 common members, Spina et al. (2017) also lists 10 common stars for $\rho$ Oph, and Cánovas et al. (2019) obtained 35 stars in common with our selection. Finally, we have 42 stars in common with Jackson et al. (2021). The remaining two stars in our selection for $\rho$ Oph were not included in the analysis of Jackson et al. (2021), but are however listed as members by Rigliaco et al. (2016) and fulfil all our membership criteria.

### Trumpler 14

Trumpler 14 is a 1–3 Myr star forming region (Sampedro et al. 2017; Randich et al. 2020; Jackson et al. 2021) associated to the Carina Nebula, one of the most massive HII regions in the Galaxy. Of the final 165 Li candidates in Trumpler 14 we
discarded six stars we could not analyse using the astrometric and photometric Gaia data, and they were additionally listed as non-members by Jackson et al. (2021) (most of them are clear non-members with a probability of 0.00, one of them is listed with P=0.27, but it is similarly discarded). One of these six stars also has a metallicity value deviating appreciably from the mean of the cluster. We also found five strong accreting stars among our final members, all of which are Li members and classified as high-probability members by Jackson et al. (2021). We also note that we discarded an F-type star with an EW(Li) that seemed to be higher than we would expect for a star forming region (150 mÅ). This star had no recorded measurements of either accretion or Hα that could perhaps explain its EW(Li), and it was additionally listed as a non-member by Damiani et al. (2017). Regarding kinematics, five stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating up to $13.0–15.0 \text{ km s}^{-1}$ from the mean of the cluster. We accepted the stars which deviate from the limits established by the $2\sigma$ interval by ~10 km s$^{-1}$, given that they fulfil the rest of membership criteria (especially regarding available Gaia data, as well as gravity indicators and lithium. All of these stars are additionally listed as high-probability members by Jackson et al. (2021). In addition, we also accepted as final candidates 24 stars with no measured RV values in iDR6, due to the fact that they fully fulfil the astrometry criteria (as well as the other criteria), and are also listed as members by Jackson et al. (2021) with high probabilities in the range of 0.85–1.00. Finally, we also accepted as final cluster candidates 11 marginal members deviating from the $3\sigma$ metallicity limit up to 0.26 dex from the mean of the cluster. All of these stars fulfil all other criteria, and 10 of them are listed as high-probability members by Jackson et al. (2021) (with P=0.86–0.99).

Regarding other membership selections from the literature, we found the following number of common stars with our selection: 63 stars in Spina et al. (2014b), as well as 65 in Damiani et al. (2017), 51 stars in Cantat-Gaudin et al. (2018), and one star in Castro-Ginard et al. (2020). Finally, Jackson et al. (2021) includes 153 out of our 159 candidates (the remaining six were not included in their analysis). We also note that Jackson et al. (2021) includes several stars which we initially discarded for not fulfilling our gravity criteria of $\gamma < 1.01$ for non-giants (see Subsect. 3.4). The reason for this is that Jackson et al. (2021) considers a less restrictive criterion of $\gamma < 1.33$. We also note that we seem to find some inconsistencies when plotting the Li-rich giant outliers selected in the field of this cluster in the CMD diagram, with eight of them appearing among the cluster candidates. However, when plotting them in the $\gamma$ index-versus-$T_{\text{eff}}$ diagram the expected distinction is found between the non-giant cluster candidates and the Li-rich giant outliers.

## Chamaeleon I

Chamaeleon I (Cha I) is a ~2 Myr star forming region (López Martí et al. 2013; Spina et al. 2014a; Sacco et al. 2015; Randich et al. 2020; Jackson et al. 2021), composed of two subclusters (Luhman et al. 2008; Sacco et al. 2017; Roccatagliata et al. 2018), with a shift in velocity of ~1 km s$^{-1}$, and ages in the range of 5–6 Myr and 3–4 Myr, respectively.

Of the final 88 Li candidates for Cha I we discarded two stars: The first one (11080297-7738425) had no Gaia data available (thus, astrometry could not be analysed) and its [Fe/H] value deviated appreciably from the mean of the distribution, beyond the extended $3\sigma$ membership interval. The second star (11095340-7634255), also a non-member according to metal-

licity, we considered as a final non-member because we could not analyse this star using the astrometric and photometric Gaia data and it was additionally listed as a clear non-member with a probability of 0.00 by Jackson et al. (2021). As a result of the analysis of Cha I we found 35 strong accreting stars, of which only one (11104959-7717517) is not a Li member, probably due to possible veiling suppressing the absorption Li line. We classified this strong accretor as an additional likely member, and it is additionally considered as a member of Cha I by Jackson et al. (2021). We also decided to classify as final cluster members four very young active stars which exhibit little dispersion with respect to the rest of candidates, but also present values of EW(Li) which are slightly larger than in the case of most member stars (841–948 mÅ). Two of them, the ones with the largest EW(Li)s, are additionally strong accretors with $H\alpha$10% > 270–300 km s$^{-1}$, and so their larger EW(Li)s are a potential effect of either accretion-induced enhancement, or, more probably in this case, a result of a strong veiling correction, where derived EWs may be too high if the veiling factor is overestimated. For more details on chromospheric activity and accretion mechanisms in regards to Li depletion, we refer the reader to Paper III. We also note that three of the stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating up to 3.7–10.2 km s$^{-1}$ from the mean of the cluster. We accepted the stars which deviate from the limits established by the $2\sigma$ interval by ~10 km s$^{-1}$, given that they fulfil the rest of membership criteria (especially regarding gravity indicators and lithium, as only one of them has available Gaia data), and all three of them are additionally listed as candidates by several studies (e.g. Esplin et al. 2017; Sacco et al. 2017; Galli et al. 2021; Jackson et al. 2021). Regarding metallicity, we also accepted seven marginal members (three of which are also strong accreting stars), deviating up to 0.45 dex from the mean of the cluster; as well as one additional star (11091172-7729124) which deviated more appreciably from the $3\sigma$ limit (0.80 dex from the cluster mean). All of these stars fulfil all other criteria, and are listed as high-probability members with P=0.99 by Jackson et al. (2021).

Regarding previous selections, we found the following number of common stars with our selection: nine stars in Robrade & Schmitt (2007), 23 in Luhman & Muench (2008), and 47 in López Martí et al. (2013). We found 12 common stars in Spina et al. (2014a), all UVES members in their selection, except for 10555973-7724399 and 11092378-7623207, for which several parameters are not released in the iDR6 catalogue. We also found 88 of our candidate stars in the member list of Esplin et al. (2017), 75 in Sacco et al. (2017), and 66 in Galli et al. (2021). Finally, Jackson et al. (2021) includes 82 out of our 87 candidates (the remaining five were not included in their analysis). Regarding field contaminants, one Li-rich giant (11000515-7623259) is listed in Casey et al. (2016).

## NGC 2244

NGC 2244 is a 2–4 Myr star forming region (Michalska 2019; Mužić et al. 2019; Jackson et al. 2021), associated with the Rosette Nebula in the Perseus Arm of the Galaxy.

Of the final 123 Li candidates for NGC 2244, we discarded 11 stars which we could not analyse using the astrometric and photometric Gaia data, and which were additionally listed as non-members by Jackson et al. (2021) (all but one are definite non-members with a probability of 0.00, while the remaining one is listed with P=0.29, and is similarly discarded). We also found 10 strong accreting stars among our final members, and all but one are Li members and classified as high-probability
members by Jackson et al. (2021) (the remaining one has no measured EW(Li) in DR6 and is not included in Jackson et al. (2021)). Regarding kinematics, one star in our final selection was a RV non-member according to the initial $2\sigma$ interval, with a RV deviating 6.0 km s$^{-1}$ from the mean of the cluster. We accepted it given that it fulfilled the rest of membership criteria, and it was also listed as a high-probability member by Jackson et al. (2021). On the other hand, we also accepted as final candidates nine stars with no measured RV values in iDR6, due to the fact that they fully fulfil the astrometry criteria (as well as the other criteria), and are also listed as members by Jackson et al. (2021) with probabilities in the range of 0.53–0.99. As for metallicity, we accepted 22 marginal members which deviate up to 0.35 dex from the mean of the cluster, and also one additional star (06312952+0454342) deviating more appreciably from the $3\sigma$ limit (0.45 dex from the cluster mean). All of these stars fulfil the rest of criteria, and are listed as high-probability members with P=0.99 Jackson et al. (2021). Regarding other membership selections from the literature, we found the following number of common stars with our selection: 78 stars in Cantat-Gaudin et al. (2018), one star in Carrera et al. (2019), and 39 stars in Michalska (2019). Finally, Jackson et al. (2021) includes 114 out of our 116 candidates (the remaining two were not included in their analysis). We also note that we included several stars with a marginal metallicity, deviating slightly from the $2\sigma$ interval, as they fulfilled the rest of criteria and were furthermore reinforced by high probabilities listed by Jackson et al. (2021).

### NGC 2264

NGC 2264 is a 4 Myr star forming region (Jackson et al. 2021), the dominant component of the Mon OB1 association in the Monoceros constellation. Several studies list ages in the 1–5 Myr range (Spina et al. 2017; Cantat-Gaudin et al. 2018; Venuti et al. 2019; Arancibia-Silva et al. 2020; Bonito et al. 2020; Gillen et al. 2020; Randich et al. 2020).

Of the final 507(535) Li candidates for NGC 2264, we discarded 28 stars on the basis of not being able to analyse them using astrometric and photometric Gaia data, and all of these 28 stars were additionally listed as non-members by Jackson et al. (2021) (most are clear non-members with a probability of 0.00, while a few are listed with probabilities up to P=0.35, and are similarly discarded). In the analysis of NGC 2264 we found 144 strong accreting stars with $H\alpha$10% > 270–300 km s$^{-1}$ among our final members, and a high number of the remaining candidates show considerable values of accretion as well, albeit under that limit. Of these 144 strong accretors, all are Li members and classified as high-probability members by Jackson et al. (2021), except for 22 stars which show higher levels of EW(Li)s than expected (marked as ‘Y?’ in the table of Appendix C), and two stars which do not appear to be Li members, but are similarly included as final members due to their strong accretion. These two stars furthermore fulfil the rest of our criteria and are included as high probability members by Jackson et al. (2021). As was also discussed earlier with Cha I, we also considered as members all of the 28 strong accretors which present values of EW(Li) notably larger than the rest of candidates (800–1270 mÅ). Furthermore, 20 of them are listed as high probability members by Jackson et al. (2021) and several are also included in several others membership studies (e.g., Lim et al. 2016; Cantat-Gaudin et al. 2018; Maíz Apellániz 2019; Venuti et al. 2019). Also see Paper III for more details on the influence of accretion for this cluster specifically. As for kinematics, three star in our final selection were RV non-members according to the initial $2\sigma$ interval, with

RVs deviating up to 9.0 km s$^{-1}$ from the mean of the cluster. As in the case of other cluster analyses, we accepted them as final members, as they fulfilled the rest of membership criteria, and were also listed as high-probability members by Jackson et al. (2021). On the other hand, we also accepted as final candidates 22 stars with no measured RV values in iDR6, due to the fact that they fully fulfilled all astrometric criteria (as well as the other criteria), and are also listed as members by Jackson et al. (2021) with probabilities in the range of 0.67–0.99. Finally, we also accepted 12 marginal metallicity members (three of which are also strong accreting stars), deviating from the $3\sigma$ limit up to 0.20–0.40 dex from the mean of the cluster. All of these stars fulfil all other criteria, and are listed as high-probability members with P=0.96–0.99 by Jackson et al. (2021).

Regarding previous membership selections from the literature, we found the following number of common stars with our selection in several studies: 197 stars in Jackson et al. (2016), 69 stars in Lim et al. (2016), 130 stars in Spina et al. (2017), 116 stars in Cantat-Gaudin et al. (2018), 198 stars in Maíz Apellániz (2019), and 226 stars in Venuti et al. (2019). Finally, Jackson et al. (2021) includes 479 out of our 503 candidates. Of the remaining candidates in our selection, 21 were not included in their analysis, and the other four (which fully fulfil all our criteria) are listed with probabilities in the range of 0.29–0.38 by Jackson et al. (2021). We also note that we seem to find some inconsistencies when plotting the Li-rich giant outliers selected in the field of this cluster in the CMD diagram, with one of them appearing among the cluster candidates. However, when plotting them in the $\gamma$ index-versus-$T_{\text{eff}}$ diagram the expected distinction is found between the non-giant cluster candidates and the Li-rich giant outliers.

### $\lambda$ Ori

$\lambda$ Ori (Collinder 69) is a young cluster with an age range of 5–12 Myr (Jackson et al. 2021; Binks et al. 2022), the oldest of the associations included in the $\lambda$ Ori star forming region (Barrado et al. 2011).

Of the final 163 Li candidates we discarded two stars because we were not able to analyse them using astrometric and photometric Gaia data, and they were listed as non-members with P=0.00–0.01 by Jackson et al. (2021). In the analysis we also found 34 strong accreting stars among our final candidates. Of these 34 strong accretors, all of them are Li members, and 17 are also listed as high-probability members by Jackson et al. (2021) (the rest are not included in their study). In addition, we list four stars with EW(Li) values slightly larger than the rest of candidates (870–950 mÅ). Similarly to the cases of Cha I and NGC 2264, we consider them as members given that they exhibit little dispersion with respect to the rest of candidates, and we explain their larger apparent values of Li for being young active stars, with high values of both accretion, and in some cases also chromospheric $H\alpha$. Two of them, the ones with the largest EW(Li) values, are additionally strong accretors with $H\alpha$10% > 270–300 km s$^{-1}$, and so, once again, their larger EW(Li)s could also be an effect of either accretion-induced enhancement or the result of a strong veiling correction. As for kinematics, 105 stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating 7–10 km s$^{-1}$ from the mean of the cluster. As in the case of other cluster analyses, we accepted them as final members, as they fulfilled the rest of membership criteria, and all stars except for 14 (which are not included in their study) were also listed as high-probability members by Jackson et al. (2021). Regarding previous membership selections from the
literature, we found the following number of common stars with our selection in several studies: 10 stars in Sacco et al. (2008), 70 stars in Hernández et al. (2010), 30 stars in Barrado et al. (2011), 16 stars in Franciosini & Sacco (2011), 51 stars in Bayo et al. (2012), and 94 stars in Cantat-Gaudin et al. (2018). Finally, Jackson et al. (2021) includes 116 out of our 161 candidates (the 45 remaining stars were not included in their study).

* Col 197

Collinder 197 (Col 197 or Cr 197) is a low-population young cluster with an age of 12–14 Myr (Bonatto & Bica 2010; Vande Putte et al. 2010; Sampedro et al. 2017; Dias et al. 2019; Jackson et al. 2021; Romano et al. 2021), embedded in the HII region of Gum 15.

Similarly to other clusters, of the final 104 Li candidates for Col 197 we discarded 11 stars as we were not able to analyse them using astrometric and photometric Gaia data, and they were listed as non-members with P=0.00–0.01 by Jackson et al. (2021). We also decided to discard one star, in spite of fulfilling all our criteria, for being listed as a non-member by both Cantat-Gaudin et al. (2018) and Jackson et al. (2021). In addition, we also accepted as final candidates six stars with no measured RV values in iDR6, due to the fact that they fully fulfil the astrometry criteria (as well as the other criteria), and are also listed as members by Jackson et al. (2021) with probabilities in the range of 0.70–1.00. In the analysis we also found 11 strong accreting stars among our final candidates. Of these 11 strong accretors, all of them are Li members, and eight are also listed as members by Jackson et al. (2021). Finally, we also accepted 10 marginal metallicity members (one of which is also a strong accretor), violating from the $3\sigma$ limit up to 0.38–0.50 dex from the mean of the cluster. All of these stars fulfil all other criteria, and are listed as members by Cantat-Gaudin et al. (2018) and/or Jackson et al. (2021) (in the case of 5 of these stars, with P=0.99).

Regarding previous membership selections from the literature, we found 79 common stars in Cantat-Gaudin et al. (2018), and Jackson et al. (2021) includes 65 out of our 92 candidates. Of the remaining candidates in our selection, eight stars were not included in their study, and in an untypical way in all our membership comparisons with this work, we found the remaining 19 stars listed as non-members with P=0.00–0.02. All of these 19 stars seem to fulfil all of our criteria, including Gaia, so we did not discard them from our final selection, but the disagreement is worth noting, even more so considering that this is the only cluster in our sample of 42 clusters whose final RV mean estimation additionally does not agree with the values given by the literature.

* Vela OB2 association: γ Velorum and NGC 2547

The γ Velorum (γ Vel) and NGC 2547 clusters are two of the PMS clusters around the Vela OB2 association (γ² Vel), each of them composed of two kinematically distinct populations. The age of both γ Vel A and B is typically listed as 10–20 Myr, with γ Vel A being 1–2 Myr older than γ Vel B (e.g., Jeffries et al. 2014; Spina et al. 2014b; Frasca et al. 2015; Sacco et al. 2015, 2017; Spina et al. 2017; Beccari et al. 2018; Franciosini et al. 2018; Armstrong et al. 2020; Jackson et al. 2021). The age of NGC 2547 is listed as either 35–45 Myr (with Jeffries et al. (2023) recently obtaining an age of $39.2_{-1.3}^{+1.6}$ Myr), or a slightly younger range of 20–35 Myr. NGC 2547 B is younger than NGC 2547 A, and coeval with γ Vel, with an age of 10–20 Myr (e.g., Oliveira et al. 2003; Soderblom et al. 2014; Sacco et al. 2017; Spina et al. 2017; Beccari et al. 2018; Randich et al.

2018; Bossini et al. 2019; Jackson et al. 2021; Pang et al. 2021; Romano et al. 2021).

The cluster membership selections for both γ Vel and NGC 2547 consist of 234 stars in γ Vel (99 in γ Vel A and 84 in γ Vel B, as well as 51 additional candidate stars that are not associated to a specific population), and 151 stars in NGC 2547 (138 in NGC 2547 A and 13 in NGC 2547 B). As mentioned in Sect. 2, for this study we used the membership selections obtained by a series of former studies in the literature: Regarding γ Vel, we first used Jeffries et al. (2014), the study specifying the members of γ Vel A and B, as well as a series of other GES studies (Damiani et al. 2014; Spina et al. 2014b; Frasca et al. 2015; Prisinzano et al. 2016; Cantat-Gaudin et al. 2019; Jackson et al. 2021). For NGC 2547 we similarly used several membership studies from the literature, some of which list the membership probabilities of each candidate to Pop. A/B Sacco et al. (2015); Randich et al. (2018); Cantat-Gaudin et al. (2018), while others do not Spina et al. (2017); Jackson et al. (2021). Gaia studies Randich et al. (2018), Cantat-Gaudin et al. (2018) and Jackson et al. (2021) offer updated membership probabilities for γ Vel A Jackson et al. (2021) and NGC 2547 (all three of them). Randich et al. (2018) listed a series of new members with respect to Sacco et al. (2015). Although their membership probabilities for each population are generally consistent with each other, a small number of stars were associated with different populations in Sacco et al. (2015) and Randich et al. (2018). In this case, we adopted the membership from Randich et al. (2018), as this is the most recent study of the two. For NGC 2547 we also compared the member stars in Sacco et al. (2015), Cantat-Gaudin et al. (2018) and Randich et al. (2018) with the candidate members in Bravi et al. (2018). For all our adopted candidates, Bravi et al. (2018) found high probabilities, namely ranging from 60 to 100%, of being RV members of the cluster. Finally, Jackson et al. (2021) offer updated probabilities for both clusters, and we decided to also prioritize the membership probabilities of this work with respect to earlier studies, for those cases where they might appear to disagree.

While we mainly used the membership lists offered by the aforementioned works, we also note that we did make use of the Gaia EDR3 data to analyse the proper motions, parallaxes and the position in the CMD for all stars marked as members for γ Vel and NGC 2547, in order to improve the final selections. We thus discarded a series of stars which deviated appreciably from the locus of the rest of members in the pmra-versus-pmdec diagram, as well as those stars which proved to be parallax non-members, and those which also deviated appreciably from the rest in the CMDs. As to field contaminants, we also find two of the γ Vel Li-rich giants in our list (J08095783-4701385 and J08102116-4740125), as well as one of the Li-rich giants from NGC 2547 (J08110403-4852137), in Casey et al. (2016). Another one, (J08110403-4658057) in γ Vel, is listed in Smiljanic et al. (2018).

* NGC 2232

NGC 2232 is a young cluster in Monoceros with an estimated age range of 18–38 Myr, with the most recent studies veering towards younger ages in the 18–25 Myr range (Liu & Pang 2019; Jackson et al. 2021; Pang et al. 2021; Romano et al. 2021; Binks et al. 2022). Jeffries et al. (2023) recently obtained an age of $29.3_{-1.0}^{+1.2}$ Myr.

Of the final 68 Li candidates for NGC 2232 we discarded 16 stars, listed as non-members with P=0.00–0.01 by Jackson et al. (2021), and not including any astrometric or photometric Gaia data. One of them additionally has [Fe/H] value deviating
appreciably from the mean of the cluster, beyond the extended
3σ membership interval. In our final selection we note as a
particular case a star with a EW(Li) value which is higher than
the rest of candidates (841–948 mÅ), but without exhibiting a
large dispersion with respect to the rest of the candidates. We
consider it as a final member given that it fulfilled the other
criteria as well as being listed in several studies as a candidate
of the cluster (Cantat-Gaudin et al. 2018; Jackson et al. 2021;
Binks et al. 2022). As for metallicity, we accepted one marginal
member deviating 0.24 dex from the mean of the cluster, as
it fulfilled all other criteria, and is listed as a high-probability
member with P=0.98 by Jackson et al. (2021). Regarding
previous membership selections from the literature, we found
46 stars in common with Cantat-Gaudin et al. (2018), as well
as 49 stars in Pang et al. (2021), and 58 stars in Binks et al.
(2022). Finally, Jackson et al. (2021) includes 56 out of our 68
candidates (the remaining stars were not included in their study,
except for one which was listed with P=0.30).

* IC 2391

IC 2391 is an open cluster in the constellation Vela, with an
age in the 20–60 Myr range (e.g., Stauffer et al. 1997; Randich
et al. 2001; Barrado y Navascués 2004; Platais et al. 2007; Smil-
janic et al. 2011; De Silva et al. 2013). Bossini et al. (2019) gave
an age of 36±2 Myr using Gaia data, while LDB estimates result
in an age range of 50–58 Myr, and Jeffries et al. (2023) recently
proposed an age of 60.1<sup>+26.8</sup><sub>-9.2</sub> Myr.

Of the final 37 Li candidates for IC 2391, we discarded two
stars from the final member selection. These stars could not be
analysed with any Gaia data (gravity, RV and Li were the only
criteria available), and they were additionally classified as non-
members with a probability of 0.00 by Jackson et al. (2021).
We also note that we added one star as part of our final se-
lection which initially was a RV non-member according to the
initial 2σ interval, with a RV deviating 6.5 km s⁻¹ from the
mean of the cluster. We accepted it as it deviated by less than
10 km s⁻¹ from the mean RV, fulfilled the rest of criteria and is
listed as a member by Jackson et al. (2021). Finally, we accepted
a marginal metallicity member deviating 0.18 dex from the mean
of the cluster, as well as one additional star (08384138-5255039)
which deviated more appreciably from the 3σ limit (0.30 dex
from the cluster mean). Both of these stars fulfil all other crite-
ria, and are listed as members by several studies (e.g, Randich
et al. 2018; Bravi et al. 2018; Jackson et al. 2021).

Regarding previous selections, we found the following
number of common candidates in the membership lists of these
non-GES studies: 19 stars in Barrado y Navascués et al. (2001)
and Barrado y Navascués (2004); three stars in Randich et al.
(2001); four stars in Dodd (2004); six stars in Platais et al.
(2007), Messina et al. (2011) and De Silva et al. (2013); and
three stars in Elliott et al. (2016). We also refer to Gómez
Garrido (2015) and Gómez Garrido et al. (2017) for a mem-
bership study of this cluster, alongside IC 2602 and IC 4665.
Regarding GES studies, we firstly found four common stars
in Spina et al. (2017) and 14 stars in Bravi et al. (2018). The
latter study derived RV membership probabilities and lists of
candidate members for this cluster (alongside IC 2602, IC 4665
and NGC 2547) using iDR4 data. Comparing our final selection
of 35 members with the list of 53 candidate stars of Bravi
et al. (2018), we find 30 kinematic candidates and 14 final
members in common (all of the common members stars have
high RV membership probabilities of at least 0.95, except for
two stars in the 0.7–0.9 range). We note that, for many stars
Bravi et al. (2018) used values of EW(Li) and/or T<sub>eff</sub> which

were derived from one of the WG12 nodes and do not appear
in our sample. As a result, we excluded these stars from our
membership analysis and only consider those stars in Bravi
et al. (2018) with EW(Li) values in our iDR6 sample. Our
mean RV and σ for IC 2391 also agree with the estimates
in Bravi et al. (2018) (see Table 3). Finally, a series of more
recent studies using both GES and Gaia data also list updated
membership probabilities for IC 2391 (Cantat-Gaudin et al.
2018; Randich et al. 2018; Jackson et al. 2021; Pang et al.
2021). We have 28 common stars as members with Randich
et al. (2018), as well as 24 stars with Cantat-Gaudin et al.
(2018), 29 with Jackson et al. (2021) (the remaining six were
not included in their analysis), and 23 with Pang et al. (2021).
As we mention in Subsect. 3.3, in this study we relied more
heavily on Jackson et al. (2021), the most extensive and one of
the most recent studies using both GES and Gaia data. Thus, we
decided to discard as a final member a marginal star according
to our analysis given that Jackson et al. (2021) listed it as a
non-member, even though Randich et al. (2018) considered it
a candidate. To finish, for completeness we also mention here
a series of recent studies that used Gaia-DR2 data to study
the spacial-kinematic distribution and cluster membership of
IC 2391 (Postnikova et al. 2019, 2020; Vereshchagin et al. 2019).

* IC 2602

IC 2602 is an open cluster in the constellation Carina, with
an age of 35 ± 1 Myr according to Bossini et al. (2019), while
previous studies found ages ranging from 30–67 Myr (Randich
et al. 1997; Stauffer et al. 1997; Randich et al. 2001; Smiljanic
et al. 2011; Jackson et al. 2021; Pang et al. 2021; Romano et al.
2021). Ages in the range of 46–53 Myr were also derived from
the LDB method (Dobbie et al. 2012; Jeffries et al. 2023), and
Jeffries et al. (2023) proposed an slightly older estimation of
45.8<sup>+4.7</sup><sub>-3.3</sub> Myr for this cluster.

Of the final 59 Li candidates for IC 2602, we discarded four
stars from the final member selection. Two of them were also
discarded via the criteria of RV and metallicity. In addition, sim-
ilarly to other clusters, all four of these stars could not be anal-
ysed with any Gaia data (gravity, RV, Li, and sometimes metal-
licity were the only criteria available), and they were additionally
classified as non-members with a probability of 0.00 by Jackson
et al. (2021) (all four stars were also listed as non-members by
Randich et al. (2018)). Finally, we accepted two marginal metal-
licity members, deviating moderately from the 3σ limit (0.30–
0.40 dex from the mean of the cluster), both fulfilling all other
criteria. One of them (10452826-6413450) is also listed as a
high-probability member with P=0.99 by Jackson et al. (2021).

Comparing our selection with former studies, we found
three stars in common with non-GES Randich et al. (2001), and
we also reinforced our membership selection with 28 additional
members of IC 2602 with measured EW(Li)s listed by this
study. We also made use of 44 additional members with EW(Li)
measurements from Jeffries et al. (2009). All these member stars
which were not measured by GES were of particular interest to
construct the empirical lithium envelope for IC 2602 (with the
envelope of Montes et al. (2001) as a base) as well, especially in
the region of the LDB (Lithium Depletion Boundary, see Paper
III). We also refer to Gómez Garrido (2015) and Gómez Garrido
et al. (2017) for a membership study of this cluster, alongside
IC 2391 and IC 4665. Regarding GES studies, we firstly found
11 common stars in Spina et al. (2017), as well as 55 kinematic
candidates and 27 final members in common with the list of
101 candidates in Bravi et al. (2018). As in the case of IC 2391,
many of the stars in Bravi et al. (2018) have no EW(Li) values
in the iDR6 sample, and therefore we excluded them in our own
membership analysis. The mean RV and σ derived in Bravi et al.
(2018) for IC 2602 are also in agreement with the ones obtained
in this paper. As for GES studies that include Gaia data, we
have 36 common stars listed as members with Randich et al.
(2018), as well as 29 stars with Cantat-Gaudin et al. (2018),
and 49 with Jackson et al. (2021) (the remaining six were not
included in their analysis).

* IC 4665

IC 4665 is an open cluster in the constellation Ophiuchus.
Bossini et al. (2019) calculated an age of 38 ± 3 Myr for this
cluster, in agreement with the 35–43 Myr range reported in sev-
eral studies (Martín & Montes 1997; Jeffries et al. 2001; Jack-
son et al. 2021; Pang et al. 2021; Romano et al. 2021). How-
ever, while they also cite an LDB age of ~ 32 Myr, Jeffries et al.
(2023) recently obtained a slightly older age of 52.9<sup>+6.8</sup><sub>-5.2</sub> Myr for
this cluster, proposing that the published LDB age should be in-
terpreted as a lower limit.

Of the final 51 Li candidates for IC 4665, we discarded 18
stars from the final member selection. Similarly to the case of
other clusters, all of these stars could not be analysed with any
Gaia data, and they were all classified as non-members by Jack-
son et al. (2021) (15 stars with a probability of 0.00, the remain-
ing three in the potential but not definite non-member range of
0.17-0.49). Randich et al. (2018) also lists seven of these 18 stars
as non-members (the rest were not considered in their analysis).

Regarding previous selections, we found the following
number of common candidates in the membership lists of these
non-GES studies: six stars in de Wit et al. (2006) and Jeffries
et al. (2009); 11 stars in Manzi et al. (2008); and 17 stars
in Lodieu et al. (2011). We also reinforced our membership
selection with 23 additional members of IC 4665 with mea-
sured EW(Li)s listed by Manzi et al. (2008). We also refer to
Gómez Garrido (2015) and Gómez Garrido et al. (2017) for
a membership study of this cluster, alongside IC 2391 and
IC 2602. As for GES studies, we found 15 common stars in
Spina et al. (2017), as well as 30 kinematic candidates and 19
final members in common with the list of 122 candidates in
Bravi et al. (2018). Our mean RV and σ are also consistent
with the previous estimates in this work. Finally, regarding GES
studies that include Gaia data, we have 11 common stars listed
as members with Randich et al. (2018), as well as 16 stars with
Cantat-Gaudin et al. (2018), 23 stars with Pang et al. (2021),
and 29 candidate stars in common with Jackson et al. (2021). In
this latter comparison, three out of the remaining four stars were
not included in the analysis of citejackson2021, and the other
one was listed with a probability of 0.14, but seeing as it fully
fulfilled all our criteria, including Gaia astrometry, we accepted
it as a final member. In spite of using the same data sets and very
similar membership criteria, it is probable that the reason for
these individual deviations between our two analyses originated
from the existing differences regarding, for example, criteria
limits or the weight of the factors implicated in the decision to
accept or discard stars (especially if they seem to be marginal
stars according to the individual membership criteria) as final
candidates.

* NGC 2451 A and B

NGC 2451 A and NGC 2451 B are two open clusters pro-
jected along the same line of sight. The age of both NGC 2547 A
and B is generally listed in the 50–80 Myr range (e.g., Balog
et al. 2009; Netopil & Paunzen 2013; Silaj & Landstreet 2014;
Randich et al. 2018; Franciosini et al. 2022b; Jackson et al. 2021;

Pang et al. 2021), while individual ages are also listed, such as
the estimations of Bossini et al. (2019), who calculated an age of
44 ± 2 Myr for NGC 2547 A and 39 ± 1 for NGC 2547 B. Thus,
the nearer cluster NGC 2547 A seems to be slightly older than
the more distant cluster NGC 2547 B Hünsch et al. (2004).

The cluster membership selections for both NGC 2451 A
and NGC 2451 B consist of 42 stars in NGC 2451 A and 64
stars in NGC 2451 B. For this study we used the membership
selections obtained by a series of former studies in the litera-
ture: Silaj & Landstreet (2014) only gives a candidate list for
NGC 2451 A, while the rest include final member lists for both
NGC 2451 A and NGC 2451 B (indicated as 'A' and 'B' in
the table of Appendix C) (Spina et al. 2017; Randich et al.
2018; Cantat-Gaudin et al. 2018; Jackson et al. 2021). Gaia
studies Randich et al. (2018), Cantat-Gaudin et al. (2018) and
Jackson et al. (2021) offer updated membership probabilities for
NGC 2451 A and NGC 2451 B. Although their membership
probabilities for each population are generally consistent with
each other, a small number of stars were considered as members
by Randich et al. (2018) and Cantat-Gaudin et al. (2018) but
were listed as definite non-members in Jackson et al. (2021).
In these cases, we adopted the membership from Jackson et al.
(2021), as this is the most recent study. As seen with other
clusters, we decided to prioritize the membership probabilities
of Jackson et al. (2021) with respect to earlier studies, for those
cases were they might appear to disagree. While we mainly
used the membership lists offered by the aforementioned works,
we also note that we did make use of the Gaia EDR3 data to
analyse the proper motions, parallaxes and the position in the
CMD for all stars marked as members for NGC 2451 A/B, in
order to improve the final selections. We thus discarded a series
of stars which deviated appreciably from the locus of the rest
of members in the pmra-versus-pmdec diagram, as well as
those stars which proved to be parallax non-members, and those
which also deviated appreciably from the rest in the CMDs.
Most of these discarded stars were also listed as non-members
by Jackson et al. (2021), in contrast to them being considered
as members by earlier studies, which reinforced our decision to
discard them on the basis of our astrometric criteria.

Appendix A.2: Intermediate-age clusters (age= 50–700 Myr)

* NGC 6405

NGC 6405 (also known as M6) is a 80–100 Myr cluster
(Kılıçoğlu et al. 2016; Netopil et al. 2016; Gao 2018; Jackson
et al. 2021), located between the local arm and the Sagittarius
arm of the Galaxy.

Of the final 53 Li candidates for NGC 6405, we discarded
two stars: one of them, additionally listed as a clear non-member
by Jackson et al. (2021), we could not analyse using the
astrometric and photometric Gaia data. The remaining one was
similarly listed as a non-member by both Cantat-Gaudin et al.
(2018) and Jackson et al. (2021), and it was also a marginal
non-member according to our parallax criterion. We also
accepted four marginal metallicity members, deviating from the
3σ limit up to 0.35 dex from the mean of the cluster. All of these
stars fulfil all other criteria, and are listed as high-probability
members with P=0.89–0.99 by Jackson et al. (2021). Regarding
other membership selections from the literature, we found 46
stars in common with Cantat-Gaudin et al. (2018), and 45 star
in Gao (2018). Jackson et al. (2021) includes 50 out of our 51
candidates. The remaining one is listed with P=0.28. Seeing as
it fully fulfilled all our criteria, including Gaia astrometry, we
accepted it as a final member.

### * Blanco 1

Blanco 1 is a cluster located towards the South Galactic Pole, with an age in the range 90–150 Myr (Gillen et al. 2020; Zhang et al. 2020; Jackson et al. 2021; Pang et al. 2021; Romano et al. 2021). Bossini et al. (2019) calculated an age estimation of 94 ± 5 Myr, while studies based on the LDB method give slighter higher ages: 114 ± 10 Myr (Juarez et al. 2014) and 132 ± 24 Myr Cargile et al. (2010). On the other hand, Jeffries et al. (2023) recently proposed a younger age of $73.5_{-6.9}^{+12.9}$ Myr for this cluster.

Of the final 101 Li candidates for Blanco I, we discarded three stars because we were not able to analyse them using astrometric and photometric Gaia data, and they were listed as non-members with $P=0.00$ by Jackson et al. (2021). As for kinematics, nine stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating $4–8$ km s$^{-1}$ from the mean of the cluster. As in the case of other cluster analyses, we accepted them as final members, and all stars were also listed as high-probability members by Jackson et al. (2021). On the other hand, we also accepted as final candidates two stars with no measured RV values in iDR6, due to the fact that they fully fulfil the astrometry criteria (as well as the other criteria), and are also listed as members by several studies (Platais et al. 2011; Zhang et al. 2020; Jackson et al. 2021). Finally, we also accepted seven marginal metallicity members, deviating from the $3\sigma$ limit up to 0.16 dex from the mean of the cluster. All of these stars fulfil all other criteria, and are listed as high-probability members by Jackson et al. (2021), with $P=0.99–1.00$.

Regarding previous membership selections from the literature, we found 58 stars in common with Platais et al. (2011), one star in Juarez et al. (2014), and 80 stars in Zhang et al. (2020). Finally, Jackson et al. (2021) includes 92 out of our 98 candidates. Four of the remaining ones were not included in their study, one is listed with $P=0.37$, and the last one is listed as a non-member with $P=0.00$. The latter two stars fully fulfil all our membership criteria, so we accepted them as final members. As in the case of IC 4665, it is probable that the reason for these individual deviations between our two analyses originated from the existing differences regarding criteria limits or the weight of the factors implicated in the decision to accept or discard stars.

### * NGC 6067

NGC 6067 is an open cluster superimposed on the Norma star cloud, with an age in a 60–150 Myr range, with the most recent studies citing ages of 100–120 Myr (Frinchaboy & Majewski 2008; Netopil et al. 2016; Alonso-Santiago et al. 2017; Magrini et al. 2018; Rangwal et al. 2019; Randich et al. 2020; Jackson et al. 2021; Romano et al. 2021).

Of the final 60 Li candidates for NGC 6067, we discarded four stars with no astrometric and photometric Gaia data, further listed as non-members with $P=0.00$ (for two stars) and $P=0.19–0.36$ by Jackson et al. (2021). As for kinematics, two stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating $13–14$ km s$^{-1}$ from the cluster mean. Listed as high-probability members by Jackson et al. (2021) and fulfilling the rest of criteria, we accepted them as final members. Regarding metallicity, we also accepted six marginal members (deviating up to 0.20 dex from the mean of the cluster), as well as one additional star (16130269-5408556) which fulfilled all criteria except for metallicity, deviating more appreciably from the $3\sigma$ limit (0.36 dex from the cluster mean). All of these stars are further listed as high-probability members by Jackson et al. (2021), with $P=0.97–0.99$. As for previous

membership selections from the literature, we found three stars in common with Frinchaboy & Majewski (2008), eight stars in Alonso-Santiago et al. (2017), and 42 stars in Cantat-Gaudin et al. (2018). Jackson et al. (2021) includes all 56 of our final candidates.

### * NGC 6649

NGC 6649 is a heavily-reddened ($E(B-V)=1.43 \pm 0.05$ Myr) open cluster located in the first Galactic quadrant with an age in a 50–126 Myr range, close to the age of the Pleiades, with the most recent studies citing ages of 120 Myr (Kharchenko et al. 2005; Dib et al. 2018; Liu & Pang 2019; Alonso-Santiago et al. 2020; Jackson et al. 2021).

Of the final 3 Li candidates for NGC 6649, we discarded one star with no astrometric and photometric Gaia data, also listed as a non-member with $P=0.00$ by Jackson et al. (2021). As for kinematics, one star in our final selection was a RV non-member according to the initial $2\sigma$ interval, with a RV value deviating $10$ km s$^{-1}$ from the mean of the cluster. Similarly to other clusters in this case, we accepted it as a final member, and it is also listed as a high-probability member by Jackson et al. (2021). Regarding previous membership selections from the literature, we found no stars in common with Cantat-Gaudin et al. (2018), and one star out of our two final candidates in common with Jackson et al. (2021) (the other one was not included in their study).

### * NGC 2516

NGC 2516 is an open cluster in the southern constellation of Carina. The age of 251 ± 3 Myr in Bossini et al. (2019) is older than the previous age estimates ranging from 110 to 150 Myr (e.g., Jacobson et al. 2016; Magrini et al. 2017; Fritzewski et al. 2020; Randich et al. 2020; Dumont et al. 2021b; Franciosini et al. 2022b; Jackson et al. 2021; Binks et al. 2022).

Of the final 379(383) Li candidates for NGC 2516, we discarded seven stars which we were not able to analyse using astrometric and photometric Gaia data, and all of them were further listed as non-members with $P=0.00–0.10$ by Jackson et al. (2021). The four stars we classified as possible Li members (marked as ‘Y?’ in the table of Appendix C) present values of $EW(Li)$ slightly larger than the rest of candidates (300–325 mÅ), especially regarding what we could expect from the age of the cluster. However, they do not exhibit appreciable dispersion with respect to the other members, all four stars fulfil the rest of membership criteria, and all of them are also listed as high probability members by Jackson et al. (2016), and so we accepted them as final members of the cluster. As well as possible scatter due to errors in the measures, we could also explain the higher observable values of $EW(Li)$ of these stars considering the fact that they also exhibit either higher values of rotation and/or chromospheric activity, and thus, a potential slower rate of Li depletion (see more on rotation and activity in regards to Li depletion in Paper III, where we analyse in detail the dependence with several parameters). As for kinematics, 23 stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating $5–8$ km s$^{-1}$ from the mean of the cluster. Fulfilling the rest of cirteria and listed as high-probability members with $P=0.99–1.00$ by Jackson et al. (2021), we accepted all of these stars as members. In addition, we accepted as final members 63 marginal metallicity members, deviating 0.12–0.30 dex from the mean of the cluster, as well as eight additional stars with values which deviate more appreciably from the $3\sigma$ limit, up to 0.40 dex from the cluster mean. As with prior clusters, all of these stars fulfil the rest of criteria and are also listed as members by several
studies (Randich et al. 2018; Cantat-Gaudin et al. 2018). Finally, all but one of these stars are further listed as high-probability members by Jackson et al. (2021), with P=0.98–1.00.

Regarding previous membership selections from the literature, we found the following number of common stars with our selection in a series of both non-GES and GES studies: Firstly, we found 326 stars in Jeffries et al. (2001), 25 stars in Terndrup et al. (2002), 55 stars in Irwin et al. (2007), 44 stars in Jackson & Jeffries (2010), six stars in Wright et al. (2011), and two stars in Heiter et al. (2014). As for GES studies, we found 56 common stars in Jacobson et al. (2016), 15 stars in Jackson et al. (2016), 13 stars in Magrini et al. (2017), 49 stars in Bailey et al. (2018), 210 stars in Cantat-Gaudin et al. (2018), 332 stars in Randich et al. (2018), and 226 stars in Fritzewski et al. (2019). Finally, Jackson et al. (2021) includes 373 out of our 379 candidates (the remaining six stars were not included in their analysis). We also note that, in order to further reinforce the membership of the stars in the field of this cluster, we also made use of the additional non-GES members with EW(Li) values of Jeffries et al. (1998).

### NGC 6709

NGC 6709 is an open cluster situated towards the center of the Galaxy in the constellation Aquila, with an age in a 150–190 Myr range (Subramaniam & Sagar 1999; Vande Putte et al. 2010; Jackson et al. 2021; Romano et al. 2021). Bossini et al. (2019) calculated an age estimation of $173 \pm 34$ Myr.

Of the final 53 Li candidates for NGC 6709, we discarded four stars with no astrometric and photometric Gaia data, all four of them also listed as non-members with P=0.00 by Jackson et al. (2021). Regarding previous membership selections from the literature, we found 45 stars in common with Cantat-Gaudin et al. (2018), and Jackson et al. (2021) includes 46 out of our 49 of final candidates. Of the remaining three, two stars were not included in their study, and the last one is listed with P=0.31. Given that said star fully fulfils our membership criteria, including astrometry, we included it in our final selection. We also note that we seem to find some inconsistencies when plotting the Li-rich giant outliers selected in the field of this cluster in the CMD diagram, with two of them, for example, appearing very close among the non-giant cluster candidates. However, when plotting them in the Kiel diagram the expected distinction is found between the non-giant cluster candidates and the Li-rich giant outliers.

### NGC 6259

NGC 6259 is an open cluster with an age of 210 Myr (Sampedro et al. 2017; Magrini et al. 2018; Casali et al. 2019; Randich et al. 2020; Jackson et al. 2021).

Of the final 39 Li candidates for NGC 6259, we discarded four stars with no astrometric or photometric data, additionally listed as non-members with P=0.00 by Jackson et al. (2021). We also accepted two marginal metallicity members (deviating up to 0.28 dex from the mean of the cluster), which fulfilled all other criteria and were listed as members by Jackson et al. (2021). Regarding previous membership selections from the literature, we found 20 stars in common with Cantat-Gaudin et al. (2018), and 34 out of our 35 members in Jackson et al. (2021). The remaining star in our final selection is listed in the latter study with P=0.32. Seeing as it fully fulfils all our membership criteria, we accepted it as a final member.

### NGC 6705

NGC 6705 (also known as M11) is a rich open cluster in the constellation Scutum, with an estimated age of 280–300 Myr

(Jacobson et al. 2016; Magrini et al. 2017; Randich et al. 2020; Jackson et al. 2021; Romano et al. 2021).

Of the final 142 Li candidates for NGC 6705, we discarded three stars with no photometric or astrometric data, also listed as non-members by Jackson et al. (2021). As for kinematics, eight stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating up to $10–13$ km s$^{-1}$ from the mean of the cluster. Also listed as high-probability members by Jackson et al. (2021), as in the case of other clusters, we accepted all of these stars as final members. In addition, we also accepted 14 marginal metallicity members (deviating up to 0.30–0.40 dex from the mean of the cluster), which fulfilled the rest of criteria, as well as one additional star (18504954-0614585) which fulfilled all criteria except for log g and metallicity, deviating appreciably from the starting $3\sigma$ limit (0.95 dex from the cluster mean). All of these stars are further listed as high-probability members by Jackson et al. (2021), with P=0.99–1.00.

Comparing our final selection with existing membership studies from the literature for this cluster, we found the following number of common stars with a number of studies: Firstly, we found 21 common stars in Magrini et al. (2014), as well as 27 stars in Tautvaišienė et al. (2015) and Jacobson et al. (2016), 15 stars in Sampedro et al. (2017) and Magrini et al. (2017), and 44 stars in Cantat-Gaudin et al. (2018) and 11 in Casamiquela et al. (2018). Finally, Jackson et al. (2021) includes 136 out of our final 139 candidates (the remaining three stars in our selection were not included in their study).

### Berkeley 30

Berkeley 30 is an open cluster with an age of 300–313 Myr (Kharchenko et al. 2005; Sampedro et al. 2017; Jackson et al. 2021; Romano et al. 2021).

Among the final 24 Li candidates for this cluster, we included one marginal metallicity member (deviating 0.20 dex from the mean of the cluster), as well as two stars with [Fe/H] values deviating more appreciably from the $3\sigma$ limit (0.83 dex from the cluster mean). All three of them are additionally listed as members by Cantat-Gaudin et al. (2018), and by Jackson et al. (2021), with P=0.94–0.99. Regarding previous membership selections from the literature, we found five stars in common with Cantat-Gaudin et al. (2018), and Jackson et al. (2021) includes all 24 of our final candidates.

### NGC 6281

NGC 6281 is an open cluster with an age of 314–316 Myr (Dias et al. 2002; Joshi et al. 2016; Sampedro et al. 2017; Jackson et al. 2021).

Among the final 24 Li candidates for NGC 6281, we accepted as a final candidate one star, listed as a high-probability member by Jackson et al. (2021), in spite of being a RV non-member according to the initial $2\sigma$ interval, with a RV value deviating $5$ km s$^{-1}$ from the mean of the cluster. Regarding previous membership selections from the literature, we found the following number of stars in common with several studies: Firstly, we found two stars in common with Frinchaboy & Majewski (2008) and Heiter et al. (2014), as well as 19 common stars in Cantat-Gaudin et al. (2018). Finally, Jackson et al. (2021) includes all 23 of our final candidates.

### NGC 3532

NGC 3532 is a very rich southern open cluster in a crowded Galactic field in Carina. Intermediate in age between the Pleiades (78–125 Myr) and the Hyades (750 Myr), it has an age
estimate of 399 ± 5 Myr according to Bossini et al. (2019), while other studies give an age range of 300–399 Myr (Dobbie et al. 2012; Fritzewski et al. 2019; Jackson et al. 2021; Romano et al. 2021).

Of the final 323 Li candidates for NGC 3532, we discarded a small number of stars according to the following criteria: Firstly, we discarded seven stars which we were not able to analyse using astrometric and photometric Gaia data, and which were further listed as non-members by Jackson et al. (2021). We also discarded one SB2 star, as well as another star we also could not analyse with astrometric criteria and deviated appreciably in the Kiel diagram. Lastly, we discarded two final stars with no Gaia data analysis available which was marked as a non-member by Fritzewski et al. (2019). As for kinematics, we note that 42 stars in our final selection were RV non-members according to the initial $2\sigma$ interval, with RVs deviating up to 8–11 km s$^{-1}$ from the mean of the cluster. As in the case of other cluster analyses, we accepted all of them as final members, and all 42 stars were also listed as high-probability members by Jackson et al. (2021). Regarding metallicity, we accepted as final candidates 76 marginal members (deviating up to 0.48 dex from the mean of the cluster), as well as two additional stars which deviate more appreciably from the $3\sigma$ limit, up to 0.83 dex from the cluster mean. All of these stars fulfil all other criteria, and are further listed as members by several studies (e.g., Cantat-Gaudin et al. 2018; Fritzewski et al. 2019, 2021), with all but eight of these stars being also considered as high-probability members by Jackson et al. (2021), with $P=0.98–1.00$.

Regarding previous membership selections from the literature, we found the following number of stars in common with these studies: Firstly, we found one star in common with Heiter et al. (2014), as well as 310 stars in Cantat-Gaudin et al. (2018), 216 stars in Fritzewski et al. (2019), 27 stars in Hetem & Gregorio-Hetem (2019), and 132 stars in Fritzewski et al. (2021). Finally, Jackson et al. (2021) includes 382 out of our 384 final candidates (the remaining two stars in our final selection were not included in their study).

We will now discuss the 72 stars that we classified as possible Li members, marked as ‘Y?’ in the table of Appendix C and represented as open squares in all $EW(Li)$-versus-$T_{eff}$ diagrams to differentiate them from the rest of the candidates (see Appendix B). These 72 stars are late-K and M-type stars with $T_{eff} < 4100$ K, and all of them consistently present values of $EW(Li)$ which are considerably larger than we would expect for KM stars in NGC 3532, a cluster with an age of 300–399 Myr. However, similarly to other clusters such as NGC 6633 (see below), the apparently overestimated $EW(Li)$ values of these stars can be explained by taking into account the inherent difficulty of obtaining accurate Li measurements for KM type stars in this age range: At low temperatures ($T_{eff} \le 4250$), the $EW$s can be thought of as pseudo-$EW$s derived by integrating over an interval which depends on $v\sin i$ values, and so they include several additional components. As a result, Li measurements for this age range do not go to zero even when Li would already be depleted, and they furthermore actually increase with decreasing $T_{eff}$ as the additional components become stronger (Gilmore et al. 2022; Franciosini et al. 2022a; Randich et al. 2022). And, while we see such $EW$s for KM stars very clearly in intermediate-age clusters such as NGC 3532 and NGC 6633, all clusters reaching such low $T_{eff}$s can show this effect, starting with younger clusters such as NGC 2547 A and NGC 2451 A/B.

Given that the $EW(Li)$ values of these stars are for these reasons not truly representative of the age of this cluster, we decided not to take them into account when creating the empirical enve-

lope for this cluster (see Paper III). However, we also consider all these stars to be robust members of NGC 3532, given that i) They do not exhibit appreciable dispersion with respect to each other or the other members, making it highly improbable for all 93 stars are spurious outliers; ii) All of them consistently fulfil the rest of membership criteria (including kinematics, astrometry, gravity and metallicity criteria); and iii) All of them are listed as high probability members by Jackson et al. (2016) with $P=0.97–1.00$, and several appear as candidates in other studies as well (Cantat-Gaudin et al. 2018; Fritzewski et al. 2019, 2021).

Another factor with which we could potentially explain some of the higher $EW(Li)$ values for these stars involves the study of their levels of rotation and activity. It is possible to explain higher observable values of $EW(Li)$ considering the fact that they might exhibit either higher values of rotation and/or higher values of chromospheric activity, as we already saw in the case of clusters such as NGC 2516, where undepleted members due to magnetic activity were found both in the literature and possibly in our final selection as well (see above). Due to the Li-rotation and Li-activity anti-correlations, these stars would consequently deplete lithium in a slower rate than their less active and/or slower rotating counterparts, thus exhibiting higher values of $EW(Li)$ as a result. Studies such as Pallavicini et al. (1990) concluded that, while less common, this effect can also be observed in K-type stars in older clusters with intermediate ages between the Pleiades (78–125 Myr) and the Hyades (750 Myr), such as both NGC 3532 (399 Myr) and NGC 6633 (575 Myr).

This cluster does indeed include a series of K-type stars exhibiting high rotation and chromospheric activity, and thus, some of these faster rotators and the more active stars could have consequently depleted Li more slowly, resulting in higher $EW(Li)$s than would be expected for late K-type stars at this age. However, while this may be a factor for some of the 93 KM stars here discussed, not all of them exhibit sufficiently high rotations and activity to be able to explain their observable Li as a consequence of these effects, many of them actually being among the slowest rotating stars in the final selection for this cluster. And so, for clusters such as NGC 3532 and NGC 6633 (see below), we believe that the most probable explanation for the $EW$s of KM members points to the GES measurement process detailed above.

## NGC 4815

NGC 4815 is a heavily populated open cluster in the constellation of Musca, with an age of 500–570 Myr (Friel et al. 2014; Jacobson et al. 2016; Magrini et al. 2017, 2018; Jackson et al. 2021).

Similarly to our analysis of others clusters in the sample, of the final 30 Li candidates we discarded one star for not being able to analyse it using astrometric and photometric Gaia data, and because it was additionally listed as a non-member by Jackson et al. (2021). As for kinematics, we accepted five stars in our final selection (also listed as high-probability members by Jackson et al. (2021)) as final members in spite of being RV non-members, with RVs deviating up to 14 km s$^{-1}$ from the mean of the cluster. Regarding previous membership selections from the literature, we found the following number of stars in common with several studies: Firstly, we found five stars in common with Friel et al. (2014), Magrini et al. (2014), Tautvaišienė et al. (2015) and Jacobson et al. (2016), as well as three common stars in Magrini et al. (2017), and seven stars in Cantat-Gaudin et al. (2018). Finally, Jackson et al. (2021) includes all 29 of our final candidates. We also note that two of the final candidates in our sample (12572442-6455173 and