A global view on star formation: The GLOSTAR Galactic plane survey. IX. Radio Source Catalog III: 2 ◦ < ℓ < 28 ◦ , 36 ◦ < ℓ < 40 ◦ , 56 ◦ < ℓ < 60 ◦ and | b | < 1 ◦ , VLA B-configuration

,


Introduction
The GLObal view of STAR formation in the Milky Way (GLOSTAR) 1 ) survey covers 145 deg 2 of the northern Galactic plane, namely the region −2 o < ℓ < 60 o & |b| < 1 o and the details of the GLOSTAR survey have been described in an overview paper (Brunthaler et al. 2021) that, in particular discusses the results for a "pilot region"(28 • < ℓ < 36 • ).Radio emission of young stellar objects in the Galactic centre region (−2 o < ℓ < 2 o & |b| < 1.0 • ) have been discussed in Nguyen et al. (2021).Chakraborty et al. (2020) have studied the differential sources count for the pilot region.The detections and properties of the 6.7 GHz CH 3 OH maser sources found in the Cygnus X region have been reported by Ortiz-León et al. (2021), and by Nguyen et al. (2022) for the region −2 o < ℓ < 60 o & |b| < 1.0 • .Dokara et al. (2021Dokara et al. ( , 2023) ) have presented the population and properties of Galactic supernova remnants (SNRs) detected in the GLOSTAR survey.The 4.8 GHz H 2 CO absorption in the Cygnus X region using both the Effelsberg 100-m telescope and VLA are presented by Gong et al. (2023).
Since the GLOSTAR survey utilizes data from the VLA in two different configurations with quite different angular resolutions, two sub-catalogs have been created -a B-configuration catalog at high angular resolution (∼ 1 ′′ ) and a D-configuration catalog at low angular resolution (∼ 18 ′′ ).The sub-catalogs for the Galactic center (−2 o < ℓ < 2 o ) and the Cygnus X region in the B-configuration and the D-configuration catalog from 2 o < ℓ < 28 o , 36 o < ℓ < 60 o , and |b| < 1.0 • are in preparation (for example, Medina et al.; Ortiz-Leon et al., in prep.), while the catalogs of the pilot region (28 o < ℓ < 36 o ) have already been published for both the D-and B-configuration data (Medina et al. 2019;Dzib et al. 2023).Table 1 summarizes the observations, sky coverage, and context of the papers of the GLOSTAR survey.The region 40 o < ℓ < 56 o has not been observed in the Bconfiguration due to limited observing time allotted for the survey.In this paper, we construct and discuss the catalog extracted from the GLOSTAR B-configuration images for the remaining region covering 2 o < ℓ < 28 o , 36 o < ℓ < 40 o , 56 o < ℓ < 60 o , and |b| < 1.0 • .This paper is organized as follows: Section 2 describes the details of the observations, data reduction, and the noise map of the survey.Section 3 presents and describes the extraction of the sources and their properties, reliability, completeness level, the in-band spectral index determination, source classification, and clustered/fragmented sources.In Section 4, we present the catalog and discuss the properties and Galactic distribution of the sources in our high-reliability catalog (these are sources above a 7σ detection threshold).In Section 5 we compare the properties of our GLOSTAR catalog with other radio continuum surveys such as CORNISH (Co-Ordinated Radio 'N' Infrared Survey for High-mass star formation, Hoare et al. 2012;Purcell et al. 2013), MAGPIS (The Multi-Array Galactic Plane Imaging Survey, White et al. 2005;Helfand et al. 2006), and THOR(The H i OH, Recombination line survey of the Milky Way, Bihr et al. 2015Bihr et al. , 2016;;Beuther et al. 2016;Wang et al. 2018).We also discuss the properties of H ii region candidates, planetary nebula candidates, extragalactic source candidates, and variable sources detected in the GLOSTAR survey.We present a summary of this work and highlight our conclusions in Section 6.

Observations
As a part of the GLOSTAR survey, the full Stokes continuum observations were conducted with the VLA in B-configuration and covered two portions the C-band (4-8 GHz), namely 4.  GHz to avoid strong radio frequency interference (RFI) around 4.1 GHz and 6.3 GHz, with 16 spectral windows and 64 channels, each channel having a bandwidth of 2 MHz.The synthesized beam in B configuration at C-band is ∼ 1 ′′ .0 and the FWHM primary beam size is ∼ 6 ′ .5 (Brunthaler et al. 2021).With a typical integration time of ∼15 seconds per pointing, the total observation time per a 2 • × 1 • region is ∼ 5 h and the root mean square (rms) noise in the images is expected to be ∼ 0.08 mJy beam −1 .The phase calibrators for this work, listed in Table 2, were observed every 5-10 minutes to correct the amplitude and phase of the interferometer data for atmospheric and instrumental effects.The absolute flux density scale was calibrated by comparing the observations of the standard flux calibrators J1331+305 (3C286) and J0137+3309 (3C48) with their models provided by the NRAO (Perley & Butler 2017).The observation and instrument parameters of the B-configuration data presented in this work are summarized in Table 2.

Data reduction, calibration and imaging
The data calibration and imaging pipelines were performed in a semi-automatic manner using the OBIT 2 package (Cotton 2008) with scripts written in python, which made use of the ObitTalk interface to access tasks from the Astronomical Image Processing Software package (AIPS) 3 tasks (Greisen 2003).The full details of the data reduction pipelines for the whole survey are described in the GLOSTAR overview paper Brunthaler et al. (2021).The continuum data calibration and imaging are presented in Medina et al. (2019) for D-configuration data and in Dzib et al. (2023) for B-configuration data.As a summary with additional remarks for the above references, we illustrate the pipeline logic of the calibration and imaging processes in the Appendix Fig. A.1,which includes the OBIT tasks used for each step.For completeness, the process is briefly described below.
Raw data collected with the VLA were calibrated and edited following 12 major steps as shown in the left-panel of Fig A.1.The downloaded data in archival science data model (ASDM) format were converted to AIPS format using OBIT task BDFIn.The initial flagging table consists of online flags (i.e., bad data collected during the observing time) and new flags from shadowed antennas, outliers in the time and frequency domains, as well as the time domain RMS filtering of calibrator data, using OBIT tasks UVFlag, MednFlag, and AutoFlag.The following tasks were used to correct the parallactic angle (OBIT/CLCor), amplitude variations in system temperature, antenna delays relative to the reference antenna (OBIT/Calib & CLCal), amplitude and phase in channels (OBIT/BPass & BPCal), amplitude and phase in time (OBIT/Calib, SNSmo & GetJy).The various calibration steps were each followed by an editing step looking for deviant solutions and flagging the corresponding data.After the first calibration pass, the corrected flags were kept and the calibration redone applying the previous flags.Diagnostic plots of various stages were made for each IF (spectral window as defined by Cotton 2008), polarization, antenna, and baseline, including amplitude/phase/delay vs time for all sources and amplitude/phase vs frequency for calibrators.Each plot was visually inspected to mark unflagged bad data such as large phase scatters, errant amplitudes, system-temperature spikes, and large delays.These flags were manually added to the special editing lists that were automatically applied in the pipeline.Final flagging tables and calibration solutions were produced and applied to all the data.The calibrated uv datasets in AIPS format were of the GLOSTAR survey.Details of the B-configuration observations of the pilot region (28 • < ℓ < 36 • ) are presented in Dzib et al. (2023).Full observation details of the GLOSTAR survey are described in Brunthaler et al. (2021).
converted to uvfits format before proceeding to the next step of imaging.
As shown in the right panel of Fig.
A.1, OBIT/MFImage was used on the calibrated uvfits data for wide-band and wide-field imaging.The imaging procedure first did a shallow CLEAN to get a crude sky model to use for flagging RFI, and then did a deeper, multifrequency CLEAN with potential self calibration.The wideband imaging (OBIT/MFImage) in Stokes I for each pointing was made by forming a weighted combined image of the sub-band images, with the field of view of the primary beam at a given frequency.The bandpass was divided into 9 sub-bands, and a narrow sub-band image was made for each sub-band.A pixel-by-pixel spectrum was fitted to the sub-band images to determine the spectral index.A weighted fit to the spectrum of each pixel was used to estimate the flux density at the central frequency of 5.8 GHz.To improve the dynamic range, selfcalibration was performed for fields with peak brightnesses that exceed a given threshold.During the imaging/ deconvolution process of the B-configuration data, the baseline range was restricted to be larger than 50 kλ (corresponding to an angular size < 4 ′′ ), in order to reject emission from poorly mapped extended structures, which are numerous in the Galactic plane.This has a minor impact on the overall sensitivity, as discussed in Dzib et al. (2023).Each field of view consists of multiple facets and all facets are cleaned in parallel.The CLEAN facet images are combined into a single plane for each sub-band.A final 1 • × 2 • image was made by the mosaic process (Brunthaler et al. 2021).Here, a circular restoring Gaussian beam with an FWHM of 1 ′′ .0 and a pixel size of 0 ′′ .25 was used.Here the first few rows of pointings from the two neighboring images were also used.

The noise level
The root mean square (RMS) noise map determined from the B-configuration image of the region covered by the GLOSTAR .25.High noise levels are found in the fields associated with bright emissions (the star-forming regions, H ii regions, and other bright radio-emitting sources), observed at low declinations (close to the Galactic center), observed in bad weather conditions, and located at the edge area of the survey.High noise striping in Galactic Latitude is due to changes in observing conditions between scans.For B-configuration data of the "pilot region" of the GLOSTAR survey (28  (Fig. 1).The bin width is 0.002 mJy beam −1 .The red, blue and black histograms represent values determined for the region discussed here (2 Dzib et al. (2023), and the combined area 2 • < ℓ < 40 • , 56 • < ℓ < 60 • area.The median RMS noise level is ∼ 0.08 mJy beam −1 and 95% of the fields have a noise level < 0.13 mJy beam −1 .survey, (i.e., the total area minus the Galactic center region and the Cygnus X region) is presented in Fig. 1.The data points represent RMS noise values determined over a quadratic areas with a side length of 3 ′ .25 (i.e., half of the average primary beam FWHM, see Brunthaler et al. 2021).As expected, the noise is found to be high in fields associated with bright and extended emission (e.g., star-forming complexes, H ii regions, and other bright radio-emitting sources), close to the edge of the survey regions, and those observed at low elevations or in bad weather (e.g., Helfand et al. 2006;Hoare et al. 2012;Purcell et al. 2013;Bihr et al. 2015).For instance, three of the fields show noise levels of 1 σ > 0.9 mJy, all of which are observed at relatively low declinations (δ < −17 • ) and located in massive star-forming complexes such as W28 A (G005.89−0.39),W31 C (G010.6−0.4), and W33 (G012.8−0.2), with associated masers, compact and UC H ii regions and outflows (e.g., Liu et al. 2010;Beuther et al. 2011;Qiu et al. 2012;Immer et al. 2014;Wyrowski et al. 2016;Urquhart et al. 2018;Yang et al. 2022a).
Fig. 2 presents the distributions of the RMS noise level for the B-configuration images (2 Dzib et al. (2023).The 1 σ noise level varies spatially over the fields, from ∼ 0.05 mJy beam −1 to ∼ 1 mJy beam −1 , with a median RMS noise of ∼ 0.08 mJy beam −1 .About 95% of the fields show RMS values < 0.13 mJy beam −1 , and only ∼1% of the fields have noise levels exceeding the median value by a factor of 3 (i.e.> 0.24 mJy beam −1 ).
As presented in Fig. 2, the region covered in this work shows a statistically higher noise level than the already published pilot region, and the median RMS values are ∼ 0.079 mJy beam −1 and ∼ 0.065 mJy beam −1 , respectively.This could be due to the fact that this work covers a larger sky area (68 deg 2 or 79% of the total fields) with more high-noise fields than the pilot region (16 deg 2 ).For instance, this paper covers the inner part of the GLOSTAR area that has lower declinations (δ < −4 • ) compared to the pilot region (δ > −4 • ), and has been observed at lower elevations, resulting in higher levels of noise.This is also seen in the CORNISH survey that shows a significantly higher noise for lower declination fields (Purcell et al. 2013).From Fig. 2, we can see the noise distribution of this work is bimodal with a second peak of smaller amplitude at the high noise side, which is mainly attributed to the high noise level in the inner regions with lower declinations.
In summary, the B-configuration continuum images of the GLOSTAR survey show a spatially varying noise level of ∼ 0.05 − 0.13 mJy beam −1 for 95% of the covered area.The typical RMS noise ∼ 0.08 mJy beam −1 is consistent with the theoretical prediction.

Source catalog construction
The source catalog presented in this work is constructed following the same strategy that was used by Dzib et al. (2023) and Medina et al. (2019) for source extraction and estimation of physical parameters.

Source extraction, fluxes, and sizes
The radio continuum sources are initially extracted using a 5 σ local noise threshold determined with the software BLOB-CAT (Hales et al. 2012).BLOBCAT produces a catalog of BLOBS that contains the properties for every detected source, including peak pixel coordinates, peak and integrated flux, RMS noise levels, number of pixels comprising each source.The effective sizes can be determined from the total number of pixels comprising each source N pix and the pixel size 0 ′′ .25 as R e f f = √ A/π, where the source area in arc seconds, A, is N pix × 0 ′′ .25 × 0 ′′ .25.

Astrometry
In order to estimate the astrometric accuracy of our data, we determined the position offsets of VLBA calibrators between GLOSTAR and Petrov (2021) which have (sub)milliarcsecond astrometic accuracy.We also compared the GLOSTAR with the CORNISH positions (Hoare et al. 2012) which has an angular resolution of 1 ′′ .5 and a position accuracy of 0 ′′ .1.As seen in the upper-panel of Fig. 3, we find no significant position offset between the GLOSTAR B-configuration and VLBA calibrator positions for the 44 common sources with mean offsets ± standard deviations of -0 ′′ .09 ± 0 ′′ .13 in Galactic longitude and -0 ′′ .01 ± 0 ′′ .09 in Galactic latitude, respectively.Similarly, as shown in the bottom panel of Fig. 3, we find no significant offsets between GLOSTAR and CORNISH for 699 compact common sources, with 0 ′′ .09 ± 0 ′′ .17 and 0 ′′ .05 ± 0 ′′ .14 in the Galactic longitude and latitude directions, respectively.Based on the above analysis, we conclude that the systematic positional uncertainty of our GLOSTAR B-configuration positions is ≲0 ′′ .1, which is consistent with that determined for the pilot region in Dzib et al. (2023).

Spurious source estimation
A total of 13689 sources are detected above 5 σ in the Bconfiguration images of the GLOSTAR region discussed in this paper (2 It is to be noted that the number of false detections is strongly related to the signal-to-noise ratio (S/N) of sources and a significant fraction of 5 σ detections are likely to be spurious or caused by residual sidelobes of nearby strong sources (Helfand et al. 2006;Purcell et al. 2013).To estimate the number of the spurious sources expected in the GLOSTAR survey, we follow the strategy used by CORNISH (Purcell et al. 2013) and run a source finder on the inverted data to extract negative detections (i.e., spurious sources), which can be done by using the tool Aegean (Hancock et al. 2012).As the cumulative distribution of spurious sources shown in Fig. 4, the false detections decrease rapidly for increasing signal-to-noise ratios.For sources below 5 σ, more than twothirds are likely spurious sources, whereas the fraction of false detections falls to about half for sources between 5 σ and 6 σ.The false detections decrease to 186 at 6.1σ.The red dashed line in Fig. 4 shows that the population of false detections is close to what is expected from Gaussian statistics, by fitting the negative detections with f (σ) = 1 − er f (σ/ √ 2) (where er f (σ/ √ 2) is the Gaussian error function as outlined in Purcell et al. 2013).Based on this, for sources above 7 σ, less than 5 spurious sources are expected to be detected.Therefore, we use 7 σ as a threshold to split the catalog into two: a high-reliability catalog with 5497 sources above 7 σ and a low-reliability catalog with 7917 sources between 5 σ and 7 σ.The low-reliability 5−7σ threshold

Galactic latitude offsets ('')
GLOSTAR-CORNISH Fraction of detections G06.5 G07.5 G15.5 G16.5 G22.5 Total The mean RMS values for the high, low, and typical noise regions are 0.12 mJy beam −1 , 0.06 mJy beam −1 , and 0.08 mJy beam −1 , respectively.95% completeness limit is reached at a flux density of ∼0.35 mJy for the best case G22.5 (red line) and ∼0.8 mJy for the worst case G06.5.Typically, the GLOSTAR survey in B-configuration is 95% complete to point sources at the ∼0.6 mJy level (i.e., the chosen 7 σ threshold for the catalog in Sect.3.3), as shown in black line for all the selected fields.See the completeness in Sect.3.4 for details.

Completeness
The completeness limit of the catalog can be investigated by running a source finder for added artificial point sources in the empty sky, that is, fields with few detections above 5 σ.The noise distributions vary spatially as outlined in Sect.2.3, and hence to investigate the completeness limit of the whole survey, we chose several regions with a size of 0.5 • × 0.5 • that sampled all noise regimes such as high noise (∼ 0.12 mJy beam −1 ), low noise (∼ 0.06 mJy beam −1 ) and typ-ical noise (∼ 0.08 mJy beam −1 ).Then, one thousand artificial point sources are injected into each selected region.The position and peak flux densities of the injected sources are randomly produced from a uniform distribution, with a flux density range of 0.05 mJy − 1.0 mJy.Following the same source extraction strategy described in Sect.3.1, BLOBCAT was used to extract these added artificial sources.After 10 iterations of injection and extraction processes, the extracted results are compared with the injected source parameters to estimate the completeness.Fig. 5 shows the fraction of the recovered artificial sources as a function of the peak intensity of the added artificial sources for fields with different noise levels.The completeness fraction varies among these fields, and fields with lower noise levels tend to reach the same completeness limit at lower flux density values.We find that a 95% completeness limit corresponds to ∼0.35 mJy (red line in Fig. 5) for the best case of G22.5 and ∼0.8 mJy (lime line in Fig. 5) for the worst case of G06.5.Typically, above ∼0.6 mJy (i.e., the 7 σ-threshold detection limit for the high-reliability catalog in Sect.3.3), the survey is 95% complete for point sources, as shown in the black line in Fig. 5 for all the selected fields.

Spectral index determination
As mentioned in Sect.2, the GLOSTAR survey covers a wide bandwidth from 4-8 GHz, which was split into 9 sub-bands and each sub-band was imaged separately.We are thus able to determine the in-band spectral indices α by ( 1) extracting the peak flux density of each source within the sub-bands and ( 2) fitting the peaks to the formula S ν ∝ ν α based on the scipy function curve_ f it.The uncertainties in the peak flux densities are taken into account in the fitting process to estimate the uncertainty for the spectral index.When we are not able to extract the peak flux densities in all 9 sub-bands for some sources due to high noise or RFI in some sub-bands, the spectral indices are determined by fitting the remaining data points.More than 98% of the 7 σ sources have reliable data suitable for spectral fitting in at least five sub-bands.Fig. 6 gives an example of the spectral index (α) fitting process for two bright sources with all 9 sub-bands and two faint sources with 5 and 6 sub-bands.With this strategy, we measured in-band spectral indices for 5430 sources in the 7σthreshold catalog.The in-band spectral index of the 5430 sources ranges from -2.85 to 2.72, as listed in Table 5.The uncertainty in spectral index σ α ranges from 0.01 to 2.63 with mean and median values of ∼0.2 for both, and is strongly correlated with the signal-to-noise ratio of the sources.Given that 97% of the 7-σ sources are compact with Y factor < 2.0 (defined as the ratio between integrated flux density and peak flux density), the measured spectral index is not expected to be affected significantly by spatial filtering of the interferometer (due to missing short spacings).

Source classification
As outlined in the surveys of CORNISH (Purcell et al. 2013), THOR (Wang et al. 2018), and the GLOSTAR pilot (Dzib et al. 2023), the classification processes of the GLOSTAR sources were conducted based on multiwavelength counterparts and/or emission properties from the Galactic plane surveys such as the near-infrared (NIR) UKIDSS survey at J H K bands (Lucas et al. 2008), the GLIMPSE survey at 3−8 µm (Churchwell et al. 2009) far-infrared (FIR) Hi-GAL at 70 − 500 µm (Molinari et al. 2010), and the submillimeter (Submm) ATLASGAL survey at 870 µm (Schuller et al. 2009;Contreras et al. 2013;Csengeri et al. 2014;Urquhart et al. 2014Urquhart et al. , 2018Urquhart et al. , 2022)).To search for counterparts, we adopt the beam size of these surveys as the search radius, such as 0.8 ′′ for UKIDSS, 2 ′′ for GLIMPSE, 6 ′′ for WISE, and 6 ′′ for MIPSGAL at 24 µm.We also carried out a visual inspection of the emission in different bands to look for association with extended structure.This was done by constructing threecolor images from UKIDSS (red K-band 2.2 µm, green H-band 1.65 µm and blue J-band 1.2 µm), GLIMPSE (red 8.0 µm, green 4.5 µm, and blue 3.6 µm), and WISE (red 22.0 µm, green 12.0 µm and blue 4.6 µm) surveys, as well as images from the MIPSGAL, Hi-GAL, and ATLASGAL surveys.
The GLOSTAR sources are then classified as five types of candidates: H ii regions, radio stars, planetary nebulae, extragalactic sources, and others, using the criteria below: H ii regions: As ionized gas regions around massive stars are located in dense molecular clouds, H ii regions are bright in the Submm, FIR, and MIR (Churchwell 2002;Anderson et al. 2019).H ii regions are also bright in the GLIMPSE three-color images due to emission from the associated polycyclic aromatic hydrocarbons (PAHs) (Purcell et al. 2013;Tsai et al. 2009).Since they are deeply embedded in molecular clouds, young H ii regions may still be dark or weak in NIR and even in some MIR bands (Hoare et al. 2012;Murphy et al. 2010;Yang et al. 2021).Therefore, radio sources associated with Submm and/or FIR emission are classified as H ii regions.
Planetary nebulae (PNe): As ionized gas regions around young white dwarf stars (Bobrowsky et al. 1998), PNe and H ii regions have similar emission properties at infrared and radio wavelengths (Anderson et al. 2012).Compared to H ii regions, the spectral energy distributions (SEDs) of PNe tends to peak at shorter MIR wavelengths and fall off steeply at FIR (Anderson et al. 2012;Purcell et al. 2013), which often makes PNe undetectable in the submm range (Urquhart et al. 2013(Urquhart et al. , 2018)).Although some nearby PNe are detectable in the FIR (Hi-GAL) and the submm (ATLASGAL), their emission typically tends to be fainter at longer wavelengths (Anderson et al. 2012;Purcell et al. 2013).Since the SEDs typically peak in the MIR range, PNe tend to appear red in WISE and GLIMPS three-color images.Moreover, PNe are likely to be isolated point-like sources in GLIMPSE and UKIDSS images due to the absence of molecular clouds (Zhang & Kwok 2009;Hoare et al. 2012).
Radio stars: Both thermal and non-thermal radio emission from radio stars have been observed, that could arise from evolved OB stars, active stars, or active binaries, as discussec by Hoare et al. (2012).In the submm and FIR ranges, emssion from radio stars tends to be weak or absent.In the WISE, GLIMPSE, and UKIDSS three-color images of , radio stars tend to appear as blue and point-like sources (Hoare et al. 2012;Urquhart et al. 2013).Also, a GLOSTAR source is classified as a radio star if it has an NIR UKIDSS counterpart (offset < 0 ′′ .8) that is suggested to be a star in the UKIDSS catalog (Lucas et al. 2008).
Extragalactic sources: extragalactic sources are usually not seen in the Submm, FIR, MIR, and NIR Galactic plane surveys (Hoare et al. 2012;Urquhart et al. 2013).Some of them are associated with very faint diffuse and/or point-like counterparts (offset < 0 ′′ .8) in the NIR UKIDSS, and the counterparts are classified as galaxies by Lucas et al. (2008).
Others: Radio sources that cannot be categorized as one of the above four types are classified as type "others".One type in this class is the photodissociation region (PDR), namely, the interface between the ionized region and the molecular cloud, which are normally extended and bright at GLIMPSE, usually with weak or no emission in the FIR and Submm (Hoare et al. 2012).
An overall flow chart of the classification process is shown in Fig. 7 and the example sources are shown in Fig. 8.The presence or absence of emission at other wavelengths, along with the source classification is indicated in Table 3.In total, among the 5437 7 σ-threshold sources, we identify candidates of 251 H ii regions, 784 radio stars, 282 PNe, 4080 extragalactic sources, and 29 Others.Among sources classified as the type "Others", 11 are likely to be PDR region candidates (i.e., normally associated with extended emission at MIR, and only weak or no emission at FIR and SMM wavelengths; Hoare et al. 2012) while the rest are unidentified.As expected, most sources (∼75%) are classified as extragalactic candidates.
To quantify the quality of the classification, we use a matching radius of 2" for determining SIMBAD counterparts (Wenger et al. 2000), as outlined in the pilot paper of GLOSTAR (Dzib et al. 2023).Among the 5437 classified sources, 1010 are found to have SIMBAD counterparts, and more than half (542/1010) are classified as the SIMBAD type 'Radio'.For the 251 H ii regions candidates, 182 show SIMBAD counterparts, 95% (173/182) of which have SIMBAD types that are consistent with the expected properties of H ii regions at different wavelengths, such as SIMBAD types of NIR/MIR/IR sources, molecular clouds, YSOs, star formation regions, dense cores, millimetric/submillimetric sources, and H ii regions.Similarly, 166 PNe, 119 radio stars and 416 extragalactic sources have SIM-BAD counterparts, with our classification being consistent with the SIMBAD type for 95% (158/166), 92% (111/119) and 97% (406/416) of PNe, radio stars and extragalactic sources respectively.In addition, the consistency of classification between this work and the CORNISH survey are 100% for H ii regions and PNe, and 98% for extragalactic sources.This supports the validity of our classification criteria.
We noted that our classification criteria may misclassify a radio source if it displays the same multiband emission properties as the above source types but does not belong to any of these categories.For example, three extragalactic source candidates are found to be associated with three pulsars within a 1.1 ′′ radius from the SIMBAD database, such as G016.8052-01.0011(e.g., PSR J1825-1446, Hobbs et al. 2004;Wang et al. 2020), G023.2721+00.2979(e.g., PSR J1832-0827, Wang et al. 2001;Yao et al. 2017), and G023.3856+00.0631(e.g., PSR J1833-0827, Tian et al. 2007;Jankowski et al. 2019).Therefore, further investigation is required to understand the nature of these candidates.

Clustered sources
Small clusters of radio sources are identified in the catalog using a friends-of-friends method (e.g., Purcell et al. 2013), that is, a source is associated with a cluster if it is located within 12 ′′ of any other member in the catalog.Among the total of 5497 sources, 60 sources are likely to be artifacts or semi-ring-like fragmented components that have been removed from the catalog.570 sources are found to be associated with 258 clusters.Among the 258 clusters, the majority (∼84%; 216/258) harbors two radio sources.About 16% (42/258) clusters harbor more than two radio sources, with ∼13% (33/258) harboring three radio sources and ∼3% (9/258) having more than three sources.All cluster members of the B-configuration are usually detected as one source in the D-configuration (beam=18 ′′ ) of the GLOSTAR with the D-configuration names from Medina et al. (2019) and Medina et al., in prep., being shown in Column 12 of Table 3.As mentioned in Sect.2.2, the imaging process was restricted to have baselines of more than 50 kλ, due to which the survey is optimized to detect emission with angular size scales up to 4 ′′ .Hence, the extended sources with complex structures tend to be resolved and decomposed into multiple compact components, presented as group sources in clusters.To distinguish between adjacent but unrelated sources (top panel of Fig. 9) and resolved sources, we manually inspected each of the clustered sources.About 155 clustered members are associated with 56 clusters that are likely to trace resolved/fragmented sources, as shown in the middle-panel of Fig. 9.All the clustered sources that trace resolved sources are listed in Table 4.Some clustered sources show one or two semi-ring-like structures around the central compact bright component, suggesting that a single source that is fragmented as seen in the bottom panel of Fig 9 and Figure 3 of Dzib et al. (2023).The semi-ring-like fragments and the associated central compact component are regarded as one single radio source.It should be noted that the flux density is trustworthy only for the detected compact sources, and the flux of the extended sources, especially the fragmented cluster sources are underestimated due to the lack of short baselines.In the future, the combined B+D images from the GLOSTAR survey will be used for the reliable measurement of both compact and extended emissions of these resolved or fragmented sources.The imaging of the combined D+B configuration data of the entire survey is ongoing and will be published in a subsequent paper.

Results
We found 5497 sources above a 7 σ threshold in the region (2 We visually inspected the 7 σ detections and exclude 18 artifacts and 42 semi-ring-like fragmented components (see Sect. 3.7), giving a final catalog of 5437 sources.Among these 5437 sources, 251 are likely to be H ii regions, 784 are candidates of radio stars, 282 are PNe candidates, and 4080 are likely to have extragalactic origins.The remaining 40 that cannot be classified into any of the above four types are regarded as "Other", including 11 PDRs and 29 "unclear", as mentioned in Sect.3.6 and listed in Table 3.

Catalog description
The catalog contains 17 columns for each source, as presented in Table 3. Columns 1-8 are determined or derived by the source finder tool BLOBCAT as described in Sect.3.1, and correspond to the Galactic name of the GLOSTAR source, Galactic longitude ℓ and latitude b, the signal-to-noise ratio (S/N), peak flux density S peak and its uncertainty ∆S peak , integrated flux density S int and its uncertainty ∆S int .The Y factor (defined as the ratio between the integrated flux density and the peak flux density, Y factor =S int /S peak ) is listed in Column 9.The source effective radius R eff , determined by the pixel size and the total number of pixels composing the source as outlined in Sect.3.1, and the spectral index α, obtained by fitting the peak flux densities of the sub-band images (see Set. 3.5), are presented in Columns 10-11.Column 12 gives the corresponding D-configuration name from the GLOSTAR survey (Medina et al. 2019, Medina et al., in prep.), which includes cluster sources as discussed in Sect.3.7.The information about the presence or absence of counterparts in the NIR, MIR, FIR and Submm (see Sect. 3.6) are displayed in Columns 13-16.Column 17 gives the source classification based on the multiwavelength properties as outlined in Sect.3.6.All the columns are displayed in Table 3 for a small portion of the catalog, with full catalog being available at CDS and the GLOSTAR website 1 .

Galactic distribution
To show the Galactic distributions for the whole B-configuration catalog, we combined our work with the catalog published in Dzib et al. (2023).Fig. 10 presents the distributions of GLOSTAR sources as a function of Galactic Longitude (ℓ) and Latitude (b).The full sample, the extragalactic, and Galactic sources are shown in light purple, purple, and dark purple, respectively.
As shown in the upper-panel of Fig. 10, the number of sources per 0.1 • latitude bin increase gradually toward zero latitude for the total sample (light purple), which is due to the expected peak at b ∼ 0.0 • of the Galactic sources (dark purple).The extragalactic sources (purple) show a relatively flat distribution with Galactic latitude, which is also seen in Hoare et al. (2012) for the CORNISH survey.We note that the distribution of Galactic sources in latitude is slightly asymmetric and skewed toward b > 0 • , which is also seen for the Galactic SNR distribution in GLOSTAR (Dokara et al. 2021) and in THOR (Anderson et al. 2017).This asymmetry is supported by a Shapiro-Wilk test of the Galactic source distribution in latitude that gives a statistical confidence of p-value≪0.001.We found that the source counts of H ii regions peaks at latitude b ∼ 0.0 • , as noted by Urquhart et al. (2013), while the distributions of radio stars and planetary nebulae are relatively flat along latitude from −1.0 • to 1.0 • .
The lower-panel of Fig. 10 shows the source counts in 2 • bins of Galactic longitudes.The source counts are seen to have peaks at certain longitudes such as ℓ ≈ 10 • and ℓ ≈ 26 • .These correspond to locations of star formation regions such as W31, SFC1-4 (Thompson et al. 2006;Murray 2011) and W42 (Gao et al. 2019;Liu et al. 2019;Feng et al. 2021).The Kolmogorov−Smirnov (K-S) tests between the Galactic sources and the extragalactic sources give a statistical confidence p-value ≪ 0.001, which confirms that their distributions are significantly different in both Galactic longitude and latitude.The source count of extragalactic sources is seen to drop in regions with high noise such as |b| > 0.9 o at the survey edges and 2 • < ℓ <8 • at low declination, as shown in Fig. 1.Considering the high noise at 2 • < ℓ <8 • , the number of extragalactic sources per 2 • longtitude bin still decrease gradually toward longitude zero, which is also seen in the distribution of non-classified sources with extragalactic origin in THOR survey (Figure 13 of Wang et al. 2018).The number of Galactic sources per 2 • longtitude bin increase toward longitude zero, which is also seen for the distributions of the resolved sources with galactic origin in the CORNISH survey (Figure 18 of Purcell et al. (2013)) and the Galactic sources (PNe and H ii regions) in the THOR survey (Figure 13 of Wang et al. 2018).

Source Properties
In this section, we present the source properties for the high reliability (7 σ−threshold) catalog.
Infrared counterpart Submm Class Gname summary of the source properties for the total catalog in Table 3, the extragalactic sources, and the Galactic sources.The typical values of flux density, effective radius, and spectral index of the sources in the catalog are S int ∼ 1.0 mJy, R eff ∼ 0 ′′ .8, and α ∼ −0.52, indicating that the detected sources are typically compact and weak, and are dominated by extragalactic sources with non-thermal emission.As shown in Table 5, the effective radius of the sample ranges from 0.45 ′′ to 2.69 ′′ , with a mean and median value of 0.86 ′′ and 0.80 ′′ respectively.As expected, the Galactic sources have larger sizes with a higher mean R eff compared to the extragalactic sources, and the two distributions are significantly different as suggested by the K-S test (p-value≪0.001).The effective source sizes are calculated from the total number of pixels comprising each source obtained from the BLOB-CAT software (see Sect. 3.1).As discussed for the THOR survey (Bihr et al. 2016;Wang et al. 2018), the effective radius is not a good parameter to distinguish between resolved and unresolved sources.Hence, we use the parameter Y factor (i.e., S int /S peak in Table 3), namely, the ratio between the integrated (in units of mJy) and peak flux density (in units of mJy/beam), to divide sources into three subsamples: extended sources (Y factor > 2.0), compact sources (1.1 < Y factor ≤ 2.0) and unresolved/point-like sources (Y factor ≤ 1.1).This is identical to what was used in the pilot region of the GLOSTAR B-configuration (Dzib et al. 2023) and the D-configuration (Medina et al. 2019), as well as in previous work (e.g., Bihr et al. 2016;Yang et al. 2019).As expected, the majority of the detected sources (97%=5255/5437 with Y factor ≤ 2) are compact (2053 with 1.1 < Y factor ≤ 2.0) and unresolved/point-like sources (3202 with Y factor ≤ 1.1), due to the fact that the images are optimized for detection of compact emission (see Sect.2.2).In the above three categories of unresolved, compact, and extended sources, 79%, 73%, and 32% of them are extragalactic sources candidates, respectively.3) and the published pilot catalog in Dzib et al. (2023) .
About 98.6% of the extragalactic sources and 91% of the Galactic sources are compact and unresolved with Y factor ≤ 2.0.The remaining 1.4% of extragalactic sources are extended with Y factor > 2.0, which are expected to be associated with radio galaxy lobes, while the extended Galactic sources (9% of the Galactic sources) are found to be preferentially associated with H ii regions and PNe.

Peak and integrated flux density
The peak and integrated flux densities are obtained as outputs from the source extraction tool, BLOBCAT (see Sect. 3.1).Fig. 12 shows the distributions of peak and integrated flux density for the high reliability catalog (light purple), the extragalactic sources (purple), and the Galactic sources (dark purple).The decline in the source count for flux densities below 0.65 mJy is due to the non-uniform noise distribution and the resulting variation in the 7 σ detection limit over the survey region.It is to be noted that there is a significant population of unresolved sources with Y factor < 1.0, which means that their integrated flux densities (S int ) are lower than their peak intensities (S peak ).While for most of the sources, the integrated flux densities and the peak intensities are still consistent within 3σ, there are a few sources where the discrepancy is significant.This could happen when the unresolved sources are ( 1) not fitted with enough pixels by BLOBCAT (i.e., fitted area less than beam); (2) located in negative side lobes from nearby bright sources; and (3) not cleaned properly.These are also seen in other survey catalogs generated using BLOBCAT, such as that of THOR (Bihr et al. 2016;Wang et al. 2018) and the GLOSTAR pilot region (Dzib et al. 2023;Medina et al. 2019).Therefore, for these unresolved sources with a Y factor < 1.0, we used their peak flux densities for further analysis discussed in the following sections.3 and the published pilot catalog in Dzib et al. (2023) .

Spectral index
between the Galactic sources and the extragalactic sources (p-value≪0.001)suggests that they are significantly different.
The spectral index is a common and useful tool for distinguishing between thermal and non-thermal radio emission, broadly corresponding to positive and negative spectral indices, respectively.The catalog is dominated by non-thermal radiation, as it includes 74% of sources with α < −0.1.The spectral index of extragalactic sources (purple) and the total sample (light purple) peaks at α ∼ −0.6, as seen in the pilot paper (Dzib et al. 2023).Among the 4080 extragalactic sources, about 77% have negative spectral indices with α < −0.1, while 23% show spectral indices from thermal emission with α > −0.1.The spectral index of Galactic sources (dark purple in Fig. 13) peaks at ∼ −0.2, with 67% having α < −0.1 (mainly radio stars) and the remaining 33% having α > −0.1 (mainly H ii regions and planetary nebulae).This shows that the radio emission from both Galactic and extragalactic sources is dominated by non-thermal radiation.It is to be noted that the uncertainty in the spectral index is strongly correlated with the logarithm of single-to-noise ratio (S/N), with a correlation coefficient of ρ =-0.89 and p-value ≪ 0.001.Hence, the in-band spectral indices of weak sources just above the detection limit of 7σ are less reliable, with there being some cases of two adjacent pixels having very different spectral indices.A similar observation was made by Rosero et al.  3 and the published pilot catalog in Dzib et al. (2023) .( 2016) for weak and compact radio emission in high-mass starforming regions.

Comparison with other radio surveys
The comparison between GLOSTAR and other radio surveys allows us to discuss the consistency of peak and integrated flux densities, and positions of the detected sources.Due to the differences in spatial filtering, angular resolution, and uv coverage, the comparisons of radio properties are valid for compact sources detected in these surveys.In this section, we compare the properties of compact sources in our catalog with other surveys such as CORNISH (e.g., Hoare et al. 2012;Purcell et al. 2013), MAG-PIS (e.g., White et al. 2005;Helfand et al. 2006), and THOR (e.g., Bihr et al. 2016;Beuther et al. 2016;Wang et al. 2018).

The CORNISH survey
The CORNISH survey has the most similar observation setup and sky coverage to the B-configuration of the GLOSTAR survey, making it an excellent resource for a comprehensive inspection of the source properties in this study.CORNISH (Hoare et al. 2012;Purcell et al. 2013) used the VLA in B and BnA configurations at 5 GHz to conduct a Galactic plane survey from 10 • < ℓ < 65 • and |b| < 1 • , with a resolution of 1.5 ′′ and a median RMS noise of ∼ 0.4 mJy beam −1 .A total of 2638 highreliability CORNISH sources are detected above the 7σ limit (∼ 2.5 mJy beam −1 ), containing extended emission on scales no greater than 14 ′′ .Compared to CORNISH, GLOSTAR in Bconfiguration has a similar angular resolution of 1.0 ′′ and better sensitivity of ∼ 0.08 mJy beam −1 (i.e., the 7σ detection limit ∼ 0.56 mJy beam −1 ), but a poorer sampling of extended emission with the images being mostly insensitive to emission on scales >4 ′′ (see Sect. 2.2).Despite their differences in the uv coverage, the properties of compact and unresolved sources are expected to be common in both catalogs.The comparison between GLOSTAR and CORNISH for the pilot region (28 Dzib et al. (2023).
Within the overlap region of this work (10 1210 high-reliability sources above 7σ, including 742 compact sources (defined as sources with angular sizes≤1.8′′ in Purcell et al. 2013).In the overlap region, GLOSTAR detected 4381 sources above 7σ, highlighting the improvement in sensitivity.Among these 742 compact CORNISH sources, we find a match of 669 GLOSTAR sources, using a circular matching threshold of 1.8 ′′ , giving a match rate of ∼ 90%.The match rate is similar (90%=259/290) even if we only consider sources that are unresolved by CORNISH (i.e., sizes ≤ 1 ′′ .5).In spite of the improved sensitivity of the GLOSTAR survey, there are 73 sources that are detected by CORNISH but not by GLOSTAR.An examination of the location of these sources reveals that about 54% (40/73) are located at the survey edges with |b| > 1.0 that are not covered by GLOSTAR.The remaining sources are likely to be variable radio sources as discussed later in Sect.5.6.The catalog discrepancies between CORNISH and GLOSTAR have also been discussed in the GLOSTAR pilot papers by Dzib et al. (2023) and Medina et al. (2019).Fig. 14 shows the comparison of flux densities of CORNISH and GLOSTAR for the 669 compact sources that have been detected in both surveys.We measured a mean± standard deviation of the flux ratio of 1.14 ± 0.3, demonstrating a good agreement between the two surveys.The consistency of flux measurements between CORNISH and GLOSTAR was also seen in the Dconfiguration catalog of the pilot region (see Fig. 11 in Medina et al. 2019).Because CORNISH (beam=1 ′′ .5) and GLSOTAR (beam=1 ′′ .0) have similar angular resolutions, we have also used these compact matching sources to infer the astrometric accuracy of the GLOSTAR survey, indicating a position uncertainty of ≲ 0 ′′ . 1 for this catalog (see Fig. 3 and Sect.3.2 for details).
Among the compact sources detected in both surveys, COR-NISH has classified 17 sources as H ii regions, 32 as planetary nebulae, and 577 as extragalactic sources.Among these COR-NISH sources with types, GLOSTAR found all 17 H ii regions and 32 planetary nebulae, as well as 569 extragalactic sources.This gives a classification match rate of 100% for H ii regions and PNe, and 98% for extragalactic sources.This demonstrates a high level of agreement in classification between the two surveys.
For the 5 GHz catalog in White et al. (2005), within the overlapping region of the GLOSTAR B-configuration (2 • < ℓ < 40 • and |b| < 0.4 • ), MAGPIS detected 2345 sources above 5.5 σ with angular size measurements, including 935 compact sources (defined as sources with angular sizes < 6 ′′ in MAGPIS White et al. 2005).In the same region, GLOSTAR detected 6216 sources above 7σ.Among the 935 compact MAGPIS sources, we find a match of 684 GLOSTAR sources above 5σ using a matching radius of 5 ′′ .The remaining sources that are detected by MAGPIS but not by GLOSTAR are extended with angular sizes greater than 4 ′′ (i.e., the largest angular scale structure detected by GLOSTAR in B-configuration).Due to the differences in the adopted detection thresholds and the uv coverage between MAGPIS and GLOSTAR, among the 684 5σ matches, 663 are above 7σ in the high-reliability catalog of GLOSTAR.From the top panel of Fig. 15, we can see that the measured flux densities of the unresolved sources (i.e., Y factor < 1.1) in GLOSTAR and MAGPIS are in good agreement, as was also seen in (Medina et al. 2019) for the D-configuration catalog of the pilot region.The outliers that have Y factor < 2 could be from the variable radio source sample such as G031.0777+00.1703 in Dzib et al. (2023) which is the outlier point located at the bottom-right of Fig. 15.The extended sources with Y factor > 2 are responsible for the outliers that show higher flux densities in MAGPIS compared to GLOSTAR, which is mainly attributed to differences in uv coverage between the two surveys.For the 1.4 GHz catalog in Helfand et al. (2006), within the overlapping region (5 • < ℓ < 32 • and |b| < 0.8 • ), MAGPIS detects 3149 sources above 5 σ, 1153 of which are compact with angular sizes < 6 ′′ sources.Using a matching radius of 5 ′′ , GLOSTAR detects 860 sources above 7 σ at 5 GHz.To compare the flux densities at the same observing frequency with MAGPIS 1.4 GHz, we extrapolated the GLOSTAR 5 GHz flux densities to the 1.4 GHz flux densities according to the spectral indices of GLOSTAR catalog.To make the 1.4 GHz extrapolated flux densities reliable, we select the 484 compact sources that have low uncertainties in their spectral indices (i.e., σ α < 0.2, where 0.2 is the typical value of σ α as outlined in Sec.3.5).The bottom panel of Fig. 15 shows the comparison of 1.4 GHz flux densities of MAGPIS and GLOSTAR for the 484 compact sources detected by both surveys.This suggests that the extrapolated 1.4 GHz flux densities from GLOSTAR agree with the 1.4 GHz MAGPIS fluxes.Considering the differences in uv coverage and the observing frequency, the flux measurements of the two surveys are consistent.

The THOR survey
The THOR survey (Bihr et al. 2015(Bihr et al. , 2016;;Beuther et al. 2016;Wang et al. 2018)) observed HI, OH, recombination lines, and continuum of the Galactic plane between 14.5 • < ℓ < 67.4 • and |b| < 1.25 • , using the VLA in C-configuration at 1-2 GHz.The continuum images of THOR have a 1σ noise level of 0.3 − 1.0 mJy beam −1 and a typical angular resolution of ∼ 25 ′′ .The detection threshold of THOR is set as 5σ, and a total catalog of 10387 sources is detected above the threshold.Given the large differences in the observation setup and uv coverage between THOR and GLOSTAR in B-configuration, we can only roughly discuss the similarities and differences between the two surveys.
The comparison between GLOSTAR and THOR for the pilot region (28 ) can be found in Dzib et al. (2023).
Within the region of overlap between THOR and the GLOSTAR area presented in this work (14.5 • < ℓ < 28 • , 36 • < ℓ < 40 • , 56 • < ℓ < 60 • and |b| < 1.0 • ), THOR detected 4083 sources above 5σ as reported in Wang et al. (2018).In the overlapped region, GLOSTAR detected 3650 sources above 7σ detection level.Using a matching radius of 5 ′′ , we found 2363 common sources above 5σ and 2001 matches above 7σ in the GLOSTAR survey.Among these 7σ matches, 1764 are regarded as point sources in the THOR survey, 1390 of which have had spectral indices measured by both surveys.Fig. 16 shows the distribution of the difference in the measured spectral index (α Bconf − α THOR ) for these sources.The mean difference in spectral index between the two surveys is −0.04 ± 0.03 (the error is estimated from a Gaussian fit to the distribution in Fig. 16), with a standard deviation of 0.78.The mean value reduces to −0.02 if the matched sources are also compact (Y factor < 2) in GLOSTAR, and down to zero if the sample is restricted to be sources that are classified as extragalactic in Sect.3.6.The spectral indices of compact sources measured in the two surveys are thus consistent.A similar result was found in the catalog of the pilot region (Dzib et al. 2023).Fig. 16 also shows the presence of a number of sources that have significant differences in spectral index measurement between THOR and GLOSTAR, such as abs(α Bconf − α THOR ) > 2. This could be due to: (1) the big differences in beam sizes of THOR (∼ 25 ′′ ) and GLOSTAR in B-configuration (∼ 1 ′′ ); ( 2) the measurement uncertainties of spectral index measured in THOR (Bihr et al. 2016) and in this work; (3) a turnover in the radio spectra between THOR and GLOSTAR.
Among the 1764 7σ compact matches, only 130 sources have been classified by THOR, including 45 H ii regions and 57 PNe candidates.The classification in THOR is consistent with our classification for 38 out of the 45 H ii regions and 53 out of the 57 PNe.This gives a classification match rate of 84% (38/45) for H ii regions and 92% (53/57) for PNe.Differences in classification can be caused by the differences in the matching radius used for comparison in IR surveys between the two surveys, as discussed in the GLOSTAR catalog paper of the pilot region (Dzib et al. 2023).

H ii region candidates
As discussed in Sect.3.6, radio sources with submm and FIR emission are classified as H ii region candidates.These H ii region candidates trace radio emission in star formation regions (SFR).In this paper, we identified 251 H ii region candidates 5 Among these H ii regions, 138 are identified/detected by previous work using the CORNISH, THOR, and the SIMBAD database.Therefore, 113 H ii regions are newly identified in this work.
The H ii region candidates of this work are compact and show a mean effective angular size of 1.2 ′′ , ranging from 0.55 ′′ to 2.69 ′′ .This indicates that the majority of them belong to the category of the most compact H ii regions (e.g, Hoare et al. 2007), such as hyper-compact H ii (HC H ii) regions and ultra-compact H ii (UC H ii) regions (e.g., Kurtz 2005;Yang et al. 2021;Liu et al. 2021;Patel et al. 2023).Combined with the compact H ii region candidates in Table 1 of Dzib et al. (2023), we obtain a sample of 390 H ii region candidates in GLOSTAR B-configuration.
The distribution of the spectral index α for the 390 H ii region candidates is shown in the top panel of Fig. 17, and α ranges from −3.70 and 2.96, with a mean value of −0.53.Considering the uncertainties of the spectral indices σ α as outlined in Sec.3.5, we choose to discuss the H ii region candidates that have a reliable spectral index, that is, σ α < 0.2 (where 0.2 refers to the typical value of σ α in the catalog).This gives a sample of 183 H ii region candidates, with -2.33<α<1.58,as the hatched histogram shown in the top panel of Fig. 17  have reported positive and negative spectral indices for H ii region candidates in the CORNISH survey (Kalcheva et al. 2018) and H ii region candidates in the GLOSTAR pilot region (Dzib et al. 2023).
Theoretically, the radio continuum emission of an H ii region is thermal and has a spectral index α (S ν ∝ ν α ) varying from +2 (optically thick) at low frequency to −0.1 (optically thin) at high frequency.We use the spectral index α to divide H ii region candidates into two subsamples: the α < −0.1 group (66% = 120/183) and the α > −0.1 group (34% = 63/183 ), with mean values of α ∼ −0.7 and α ∼ 0.6, respectively.Considering 2 times of the uncertainties in spectral index (2 σ α ) for the two samples, the fraction of the α > −0.1 sample increase to about 50%.The mean spectral index of 0.6 for the α > −0.1 group refers to the thermal emission from H ii regions, which is similar to the mean of α ∼ 0.6 observed at 5 GHz for the young H ii regions (HC H ii and UC H ii) sample in Yang et al. (2019)    suggest the existence of H ii regions with a mix of optically thin and thick components along the line of sight.Intriguingly, the majority (∼ 66%) of the sample belongs to the α < −0.1 group with a mean spectral index of −0.7, which indicates that a substantial portion of radio emission in these H ii region candidates is non-thermal.A number of observational studies have reported the existence of H ii regions with a mixture of thermal and non-thermal radiation (e.g., Wang et al. 2018;Meng et al. 2019;Padovani et al. 2019) and dominated nonthermal emission (e.g., Wilner et al. 1999;van der Tak & Menten 2005;Rosero et al. 2019).The H ii region with radio continuum α < −0.5 are considered to be dominated by non-thermal emission (Kobulnicky & Johnson 1999).Considering that 74 out of the 120 H ii region candidates have α < −0.5, we suggest these are dominated by non-thermal emission, while the remaining 46 candidates with −0.5 < α < −0.1 are likely to be associated with a mixture of thermal and non-thermal radiation.
Given that the H ii region candidates are identified by radio emission in star formation regions (see Sect. 3.6), it is possible that the non-thermal radio emission originates from the processes such as radio jets, shocks and outflows from high-mass (van der Tak & Menten 2005) and low-mass stars (Gómez et al. 2005).We find that the α < −0.1 group is more likely to be located in clusters (see Sect. 3.7) and be associated with molecular outflows in Yang et al. 2018Yang et al. , 2022a compared compared to the α > −0.1 sample, implying that the non-thermal emission arises from lo-calized spots that are seen only when the large scale emission is filtered out by the interferometer.From Fig. 17, we can see that the α < −0.1 sample shows significantly lower values of peak flux density compared to the α > −0.1 sample, indicating that the α < −0.1 sample are relatively compact and weak.Thus, these clustered non-thermal emission spots are likely to be radio jets and outflows located in the vicinity of H ii regions and in star formation regions (e.g., Wang et al. 2012;Purser et al. 2016;Rosero et al. 2016;Liu et al. 2017;Qiu et al. 2019).This is consistent with the findings of Wang et al. (2022) who suggested that most of the radio sources in the Cygnus region are radio jets and winds originating from massive young stellar objects.In Fig. 18, we displayed an example of the two H ii region candidates (G024.7897and G024.7888) that shows positive and negative in-band spectral index derived by fitting the 8 sub-bands radio images of the GLOSTAR.The H ii region candidate G024.7897 (Y factor = 1.02) with α = 1.43 ± 0.02 is supposed to be a "real" H ii region, which shows extended green emission (as defined by Cyganowski et al. 2008) in Fig. C, and is associated with a maser-emitting UC H ii region in Hu et al. (2016).In contrast, the nearby H ii region candidate G024.7888 (Y factor = 1.97) with α = −0.55 ± 0.04 is likely to be the nonthermal emission from radio jets or outflows in massive starforming regions.Some of the sources in the catalog were confirmed as non-thermal sources through VLBI observations by Dzib et al. (2016).We note that there are many H ii regions like G024.7897 that are surrounded with at least one non-thermal source like G024.7888, as shown in Fig. C.1, which is consistent with the findings in Gómez et al. (2005) who detected a cluster of non-thermal sources around a young and compact UC H ii region with VLA observations and suggested these non-thermal clusters are originated from low-mass, pre-main-sequence stars.
In summary, from the α > −0.1 sample, we confirm the existence of the H ii regions with a mixture of optically thin and thick thermal emission components.From the α < −0.1 sample, we find that a large fraction of compact H ii region candidates are associated with non-thermal emission, suggesting that these candidates can be radio jets, winds and outflows from high-mass and low-mass young stellar objects.This further indicates that there is a significant amount of relativistic electrons that exist in star-forming regions.Further investigation is required to confirm the nature of non-thermal emissions in massive star-forming regions.

Planetary nebula candidates
We identified 282 planetary nebulae (PNe) based on the classification process described in Sect.3.6.Among these, 155 are identified/detected by previous work and 127 PNe are new detections.
These PNe candidates are compact and their effective size range from 0.53 ′′ to 2.46 ′′ with a mean and standard deviation of 1.0 ′′ and 0.3 ′′ respectively.About 54% of the sources have effective sizes less than 1.0 ′′ .The in-band spectral indices α of the PNe candidates range from −2.75 to 1.9, with a mean value of −0.52.Combining with the 68 PNe candidates identified in the Table 1 of Dzib et al. (2023) for the pilot region, we obtained a sample of 350 PNe candidates for the GLOSTAR Bconfiguration catalog.
As in H ii regions, from the 162 PNe candidates with σ α < 0.2, two subsamples are obtained: (1) PNe candidates with α < −0.1 (62%=100/162) and ( 2) PNe candidates with α > −0.1 (38%=62/162).This is similar to observations from previous surveys such as CORNISH (Irabor et al. 2018) and the pilot region of GLOSTAR (Dzib et al. 2023), who found that PNe were found with radio emission showing both positive and negative spectral indices.This suggests that both thermal and non-thermal emission components are associated with PNe candidates.Theoretically, PNe are thus expected to have radio continuum from thermal free-free emission with −0.1 < α < 2 (e.g., Gómez et al. 2005;Tafoya et al. 2009;Qiao et al. 2016).Thus, similar to the H ii region candidates, the PNe candidates with spectral index greater than and less than -0.1 are expected to be associated with thermal and non-thermal emission, respectively.
Observationally, a low spectral index threshold of α < −0.5 between 1 GHz and 20 GHz was suggested to identify nonthermal emission from planetary nebulae detected in the AT20G survey (i.e., the Australia Telescope 20 GHz survey, Chhetri et al. 2015).Planetary nebulae with both non-free-free emission components (e.g., Casassus et al. 2007;Irabor et al. 2018) and even non-thermal dominated emission (e.g., Suárez et al. 2015;Cerrigone et al. 2017;Anglada et al. 2018) have been reported before.Non-thermal emission is found in young PNe that are associated with masers (Cohen et al. 2006) or that are in the formation stage (Cerrigone et al. 2017).In our sample of sources with α < −0.1, 59 PNe candidates have α < −0.5 and thus are likely to be dominated by non-thermal emission, while the remaining 41 PNe candidates with −0.5 < α < −0.1 are thought to have both thermal and non-thermal components.The α < −0.5 candidates show significantly lower Y factor and flux densities com-pared to sources of the α > −0.5 sample.This suggests that the α < −0.5 sources are denser and weaker and thus likely associated with radio jets.In short, we found a large fraction of PNe are associated with compact non-thermal emissions likely from radio jets, suggesting that these PNe candidates are very young or in the formation stage.

Extragalactic source candidates
In this paper, we identified 4080 extragalactic sources based on the classification process outlined in Sect.3.6.This is consistent with the expected number of ∼ 4137 estimated from the FIRST survey (White et al. 1997) and ∼ 4100 calculated from Equation (A2) of Anglada et al. (1998) using the 5 GHz source counts of Condon (1984) and the typical detection level (∼ 0.6 mJy beam −1 ) of this work.Among these extragalactic sources, 1905 have been identified/detected previously in the CORNISH and THOR surveys and the SIMBAD database.Therefore, 2175 extragalactic sources are newly identified by this work.
Combining with the 1157 extragalactic source candidates in the Table 1 of Dzib et al. (2023) for the pilot region, we have a sample of 5237 extragalactic source candidates for the GLOSTAR B-configuration catalog.These extragalactic source candidates have effective sizes ranging from 0.45 ′′ to 3.24 ′′ , with a mean value of 0.86 ′′ and a standard deviation of 0.25 ′′ .As expected, the effective size of the extragalactic source candidates is systematically smaller than that of the Galactic sources (i.e., H ii region and PNe candidates).Some extragalactic source candidates show extended radio emissions that are likely to be associated with radio galaxy lobes as suggested in the CORNISH survey (Purcell et al. 2013).The in-band spectral index α of the 4836 extragalactic source candidates are measured, ranging from −3.24 to 2.72, with a mean value of −0.58.As expected, the majority of extragalactic source candidates (77%=3725/4836) have α < −0.1, indicative of non-thermal synchrotron emission.The remaining 23% extragalactic candidates show α > −0.1, which could be due to the spectral turnover due to synchrotron selfabsorption and free-free absorption (e.g., Bicknell et al. 2018;Shao et al. 2022) of gigahertz peaked spectrum sources (GPS, e.g., Bicknell et al. 1997) or high-frequency peakers (HFP, e.g., Dallacasa et al. 2000), as discussed in the pilot field GLOSTAR article (Dzib et al. 2023).Given that extragalactic radio sources with α < −1.3 are expected to be high-redshift galaxies with ultra steep spectra (Wang et al. 2018), 13% (698/5237) of extragalactic candidates with α < −1.3 in this work are thus likely to be in this category.We note that there are more than 200 extragalactic candidates associated with NIR ∼ 2µm UKIDSS counterparts that are suggested to be galaxies by Lucas et al. (2008), which are likely to be radio galaxies in the Zone of Avoidance (e.g., Marleau et al. 2008) or quasars behind the Galactic plane (e.g., Fu et al. 2022).

Radio star candidates
In this paper, we have identified 784 radio star candidates based on the classification process in Sect.3.6, which refers to a group of radio sources that are point-like and blue in the three-color images at near infrared UKIDSS, GLIMPSE, and WISE, but with weak or no emission at submm and FIR wavelength (due to the dispersion of natal molecular clouds), as shown in panel (c) of  3 and the published pilot catalog (28 • < ℓ < 36 • and |b| < 1 • ) in Table 1 of Dzib et al. (2023).The bin sizes are 0.1 • and 2.0 • for the upper and lower panel, respectively.The blank region in Galactic longitude 40 • < ℓ < 56 • is not covered by the GLOSTAR survey in B-configuration.
the SIMBAD database.Therefore, 581 radio stars are newly detected and identified by this work.
Combining the 81 radio star candidates in Table 1 of Dzib et al. (2023), a sample of 865 radio star candidates are obtained for the GLOSTAR B-configuration.These radio star candidates have effective sizes ranging from 0.51 ′′ to 3.04 ′′ , with a mean value of 0.83 ′′ .The in-band spectral index α of the radio star candidates ranges from -2.38 to 1.94, with a mean value of -0.41.As in H ii regions, from the 278 radio stars with σ α < 0.2, there are 61% (169/278) showing non-thermal emissions with α < −0.1 and 39% (109/278) showing thermal emissions with α > −0.1.This is expected as both the thermal and non-thermal radio emission from radio stars have been observed (Hoare et al. 2012).
Fig. 19 shows the Galactic distributions for radio star candidates in the GLOSTAR B-configuration catalog.The number of radio star candidates per 0.1 latitude bin is relatively flat and the source counts decrease near the edges of b range possibly due to the higher noise level as discussed in Sec.3.4.The source counts of radio star candidates per 2 • longitude bin are found to increase toward low longitudes, which is consistent with the increasing of the Galactic sources toward low longitudes as discussed in Sec.4.2.Given that there might be radio-emitting sources showing the same multiband emission properties, the nature of these radio star candidates are needed to be explored further.

Variable sources
Variable radio sources are defined based on the ratios of peak flux densities between the GLOSTAR and CORNISH surveys for compact sources (i.e., Y factor < 2.0), as outlined in Dzib et al. (2023).A source is considered to be a variable source if ( 1) it shows a flux density ratio between GLOSTAR and CORNISH larger than 2.0; (2) it is detected in CORNISH at 7σ but is not detected by GLOSTAR at 7σ; (3) it is a GLOSTAR source with a peak flux density higher than the CORNISH detection limit (2.7 mJy beam −1 ) but is not detected by CORNISH.
Using the criteria listed above, we identified 245 variable radio sources in this work.Together with the 49 variable sources identified in Dzib et al. (2023) for the pilot region, there are a total of 294 variable sources in the GLOSTAR B-configuration catalog.Table 6 lists the source names, peak flux densities, and the suggested source types of all these variable radio sources.The majority of variable sources in the current work (76%=186/245) are found to be extragalactic in origin or infrared quiet based on the classification of this work and CORNISH.The catalog of variable radio sources in this paper includes 12 H ii regions, 5 planetary nebulae, and 4 PDRs.The variability of H ii regions (e.g., Yang et al. 2022b;Dzib et al. 2023) and planetary nebulae (e.g., Cerrigone et al. 2011;Suárez et al. 2015;Yang et al. 2022b;Dzib et al. 2023) have been reported earlier, which are interesting targets for further exploration.

Conclusions
The GLOSTAR survey covers 145 square degrees of the Galactic plane and observes both spectral lines and continuum in the 4-8 GHz using the VLA in B-and D-configurations and the Effelsberg 100-m telescope.In this paper, we present a catalog of continuum sources in the 68 square degrees of the Galactic plane (2 o < ℓ < 28 o , 36 o < ℓ < 40 o , 56 o < ℓ < 60 o & |b| < 1 o ) observed with the VLA in B-configuration.The images have an angular resolution of 1 ′′ , a typical 1σ noise level of ∼ 0.08 mJy beam −1 , and are restricted to detect extended emission on angular scales smaller than 4 ′′ .The high and low reliability catalogs from this region are presented in Tables 3 and B.1, respectively.The main results are summarised below.dominated by non-thermal radio emission.The spectral indices of compact sources between GLOSTAR and THOR are consistent.3. We classified all sources above 7 σ-threshold based on the presence or absence of counterparts/emissions in Galactic plane surveys at infrared and submillimeter wavelengths, as well as the SIMBAD database.We identified candidates of 251 H ii regions, 784 radio stars, 282 PNe, 4402 extragalactic sources, and 29 others (11 PDRs and 18 unclear).The consistency in classification between GLOSTAR and COR-NISH survey is 100% for H ii regions and PNe, and 98% for extragalactic sources.4. A significant fraction of candidates of the H ii region and PNe show a spectral index α < −0.1 (or even α < −0.5), suggesting that there are significant numbers of non-thermal emission sources corresponding to radio jets, winds, and outflows in the vicinity of young H ii regions and PNe.As expected, the majority of extragalactic candidates (77%) show α < −0.1, indicating non-thermal emission.The remaining 23% showing α > −0.1 are likely to be the gigahertz peak spectrum sources or high-frequency peakers.For the above three source types, the α < −0.1 group is more likely to be associated with clustered sources compared to the α > −0.1 group.5. We compared the measured flux densities between GLOSTAR and other radio surveys (CORNISH, MAGPIS, and THOR), and find a high level of agreement in the flux density and spectral index for compact sources detected in these surveys.We identified 245 variable radio sources listed in Table 6 by comparing the fluxes between the GLOSTAR and CORNISH survey, and most of these variable sources are found to be infrared and submillimeter quiet.
To date, the GLOSTAR survey using VLA in its Bconfiguration has the highest resolution (∼ 1.0 ′′ ) and best sen-sitivity (1σ ∼ 0.08 mJy beam −1 ) for the Galactic plane in the C band (4-8 GHz).In this work, we present a continuum catalog of 13354 sources ≥ 5σ and 5437 sources ≥ 7σ detection level.From the high-reliability 7σ−threshold catalog, we identified the largest sample of candidates of H ii regions, PNe, extragalactic sources, and variable sources in the Galactic plane.It is worth noting that a significant fraction of radio emission associated with massive star formation regions are non-thermal, indicating that relativistic electrons commonly exist.All the catalogs and data are available online at the GLOSTAR website https://glostar.mpifr-bonn.mpg.de.

Fig. 1 .
Fig. 1.RMS noise map of the 2 • < ℓ < 40 • , 56 • < ℓ < 60 • , & |b| < 1 • region observed in B-configuration of the GLOSTAR survey.Each field is made by sampling on a size of 3 ′.25.High noise levels are found in the fields associated with bright emissions (the star-forming regions, H ii regions, and other bright radio-emitting sources), observed at low declinations (close to the Galactic center), observed in bad weather conditions, and located at the edge area of the survey.High noise striping in Galactic Latitude is due to changes in observing conditions between scans.For B-configuration data of the "pilot region" of the GLOSTAR survey (28 • < ℓ < 36 • , |b| < 1 • ), seeDzib et al. (2023).

Fig. 5 .
Fig. 5. Completeness fraction of recovered artificial sources as a function of the peak flux density of the total added artificial sources.The completeness of the selected fields with different typical noises are shown in different colors.The mean RMS values for the high, low, and typical noise regions are 0.12 mJy beam −1 , 0.06 mJy beam −1 , and 0.08 mJy beam −1 , respectively.95% completeness limit is reached at a flux density of ∼0.35 mJy for the best case G22.5 (red line) and ∼0.8 mJy for the worst case G06.5.Typically, the GLOSTAR survey in B-configuration is 95% complete to point sources at the ∼0.6 mJy level (i.e., the chosen 7 σ threshold for the catalog in Sect.3.3), as shown in black line for all the selected fields.See the completeness in Sect.3.4 for details.

Fig. 6 .
Fig. 6.Example of the peak flux as a function of the in-band frequency for bright sources (top panel) and faint sources (bottom panel).Each data point (circle or triangle) refers to the peak intensity at each subband, with the error of the peak measurement.The solid and dashed lines in the top panel show the best fit to the peak flux densities from 9 sub-bands for two bright sources: G020.0808-00.1354(α = 1.1 ± 0.07) and G021.7700+00.8423(α = −0.92±0.04),respectively.The solid and dashed lines in the bottom panel present the spectral fitting of 5 and 6 sub-bands for two faint sources: G011.9357+00.2871(α = −1.7 ± 0.36) and G036.8880+01.0288(α = 2.5 ± 0.6), respectively.

Fig. 9 .
Fig. 9. Example of cluster sources in three cases.Top panel: the cluster sources in normal case.Middle-panel: example of clusters consists of over-resolved sources are listed in Table 4. Bottom panel: example of clusters with semi-ring-like fragments.In this case, the sources with semi-ring-like structures have been removed from the catalog.The positions of 7σ cluster sources are noted by blue and white pluses.The Galactic name of each figure refers to the source in the center position.The cyan contours are emission from B-configuration images (this work), starting at 5σ with 5σ increment.The background of each image presents the B+D configuration of GLOSTAR (Brunthaler et al. 2021).The FWHM beams of B (1 ′′ .0) and B+D configuration (4 ′′ .0) are presented by the white circles in the lower-right and lower-left corners of each image.

Fig. 11
Fig. 11  shows the source counts as a function of source effective sizes (see Sect 3.1) for the high reliability catalog of GLOSTAR B-configuration (light purple), with the extragalactic and Galactic sources being shown in purple and dark purple respectively.As shown in Table5, the effective radius of the sample ranges from 0.45 ′′ to 2.69 ′′ , with a mean and median value of 0.86 ′′ and 0.80 ′′ respectively.As expected, the Galactic sources have larger sizes with a higher mean R eff compared to the extragalactic sources, and the two distributions are significantly different as suggested by the K-S test (p-value≪0.001).

Fig. 11 .
Fig. 11.Distribution of source effective radius for the 6894 sources in the 7 σ−threshold catalog of this work (Table3) and the published pilot catalog inDzib et al. (2023)

Fig. 13 Fig. 12 .
Fig. 13  shows the distribution of spectral index for the full high reliability sample of this work (light purple), extragalactic (purple), and Galactic sources (dark purple).The measured in-band spectral index of this work is consistent with that measured in THOR for compact sources that are detected in both surveys, as discussed later in Sect.5.1.3.The K-S test of the spectral index

Fig. 14 .
Fig. 14.Comparison of flux densities between GLOSTAR and COR-NISH.Top panel: The comparison of flux densities for 669 compact sources detected by both CORNISH and GLOSTAR catalogs.The error bar of each point shows the uncertainty of flux.Bottom panel: The histogram of the flux density ratios of compact sources between COR-NISH and GLOSTAR.The black line displays the Gaussian fit to the histogram.The dashed line presents the mean value of the distribution, with mean and standard deviation of the flux ratio from the Gaussian fit are 1.14 and 0.3, respectively.

Fig. 15 .
Fig. 15.Comparison of flux densities between GLOSTAR and MAG-PIS.Top panel: The comparison of measured flux densities for 663 MAGPIS compact sources at 5 GHz (White et al. 2005) common to GLOSTAR.These compact sources are defined as sources with angular sizes < 6 ′′ in MAGPIS.The dashed line means the flux densities are the same in MAGPIS and GLOSTAR.At the top, we show the color bar for the Y factor (defined as S int /S peak ) of the GLOSTAR detections, indicating the emission of sources in the GLOSTAR image are unresolved (defined as Y factor < 1.1 in Sect.4.3.1),compact (1.1 < Y factor < 2.0) or extended (Y factor > 2.0).Bottom panel: The comparison of flux densities of 484 compact MAGPIS at 1.4 GHz that are also detected by GLOSTAR at 5 GHz.The 1.4 flux densities from GLOSTAR is extrapolated from the 5 GHz flux densities and spectral indices of the GLOSTAR catalog.

Fig. 16 .
Fig. 16.The distribution for the difference (α Bconf − α THOR ) in the measured spectral indices for 1390 compact THOR sources common to GLOSTAR in B-configuration.The black solid line means the Gaussian fit for the distribution, giving a mean value of −0.04 ± 0.03 (the error of the mean is from the Gaussian fit) and a standard deviation of 0.78.
. Previous studies5 In Appendix C, we present the MIR emission around the H ii region candidates.

Fig. 17 .
Fig. 17.Distributions of spectral index and peak flux densities of H ii region candidates in GLOSTAR B-configuration.Top panel: The histogram of the measured spectral index for the total 390 H ii region candidates of GLOSTAR B-configuration.The hatched histogram represents the 183 candidates with σ α < 0.2.The bin size is 0.25.The 183 candidates are divided into two subsamples: (1) 120 H ii regions with α and (2) 63 H ii regions with α > −0.1, indicating non-thermal and thermal emissions, respectively.Bottom panel: The cumulative distribution of the peak flux densities for the α < −0.1 sample (orange line) and the α > −0.1 sample (black line).The orange and black vertical dashed lines show the mean values of 9.5 mJy beam −1 and 31.6 mJy beam −1 for the two subsamples, respectively. who

Fig. 18 .
Fig. 18.Example of two H ii region candidates with positive and negative in-band spectral index, as discussed in Sect.5.2.The top two rows show the 8 sub-bands of the GLOSTAR image.The bottom-left panel shows the averaged image at 5.8 GHz used to extract the source, and the bottom-right panel shows the in-band spectral index fitting for the two compact H ii regions candidates: the solid line for G024.7897(Y factor = 1.02) with α = 1.43 ± 0.02 and the dashed line for G024.7888(Y factor = 1.97) with α = −0.55 ± 0.04.The FWHM beam of GLOSTAR in B-configuration (1.0 ′′ ) is indicated by the white circles in the lower-left corner of each image.There is a clear trend that the fluxes increase and decrease as the increasing frequencies in the sub-bands for G024.7897 and G024.7888 respectively.

Table 1 .
Summary of the GLOSTAR papers

Table 2 .
Summary of VLA continuum observations in B-configuration for the region ( Dzib et al. (2023).
Table 5 displays a statistical

Table 3 .
GLOSTAR B-configuration catalog for 2 The description of each column is displayed in Sect.4.1.Please noted that the flux measurements might be unreliable for sources which have S int < S peak , that is, Y factor = S int /S peak < 1.0, and please see Sect.4.3.2 for details.Column (12) lists the corresponding D-configuration counterparts(Medina et al., In.Prep.), which indicates the cluster sources as discussed in Sect.3.7.Symbol ✓ in Column 13-16 indicates that the sources show counterparts and/or emission at NIR (i.e., the UKIDSS and/or GLIMPSE surveys), MIR (i.e., the WISE and/or MIPSGAL surveys), FIR (i.e., the Hi-GAL survey), Submm (i.e., the ATLASGAL survey) as discussed in Sect.3.6.Symbol − refers to no association of counterpart and/or emission.Column (17) for source classification as described in Sect.3.6:EgC= Extragalactic source candidate, HII = H ii region candidate, Radio star, PN = planetary nebula candidate, Other = source that cannot be classified as one of the above four types, such as PDR (Photodissociation region) or source with no clear classification.Only a small portion of the data is provided here.And the full version will be available in electronic form at the CDS via anonymous ftp to cdsarc.cds.unistra.fr(130.79.128.5)orvia https://cdsarc.cds.unistra.fr/cgi-bin/qcat?J/A+A/.

Table 5 .
(2)mary of source properties of the 7 σ−threshold catalog.In Columns(2)(3)(4)(5)we give the minimum (x min ), maximum (x max ), mean ± standard deviation (x mean ± x std ), and median values (x med ) of these parameters for the total sample, the extragalactic sample and the Galactic sample.

Table 6 .
Catalog of 245 variable radio sources identified by this work Source names in column (1) are taken from GLOSTAR or CORNISH if no GLOSTAR detections.Symbol − refers to no measurements.Only a small portion of the data is provided here.The full table will be available in electronic form at the CDS via anonymous ftp tocdsarc.cds.unistra.fr(130.79.128.5)orvia https://cdsarc.cds.unistra.fr/cgi-bin/qcat?J/A+A/.