Free Access
Issue
A&A
Volume 653, September 2021
Article Number A35
Number of page(s) 14
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/202140500
Published online 13 September 2021

© ESO 2021

1. Introduction

Local galaxy observations reveal the large-scale structure (LSS) of the Universe to be a rich network of filaments, walls, nodes, and voids (de Lapparent et al. 1986; Bond et al. 1996; Alpaslan et al. 2014). According to the Λ cold dark matter (ΛCDM) cosmological model these structures are built of two main elements: baryonic and dark matter. The former exists in the form of stars, gas, and dust; these can be traced at different wavelengths using large sky surveys. Dark matter, however, cannot be observed directly, although it is gravitationally dominant. Therefore, we rely on visible baryonic matter observations to indirectly trace the underlying dark matter distribution.

One of the methods used to connect baryonic and dark matter involves measurements of the galaxy two-point correlation function (2pCF; Peebles 1980). This powerful statistical tool describes the spatial distribution of galaxies and has been extensively used in the past to quantify their clustering and its various dependences: on luminosity (Norberg et al. 2001; Pollo et al. 2006; de la Torre et al. 2007; Meneux et al. 2009; Zehavi et al. 2011; Marulli et al. 2013; Guo et al. 2015; Farrow et al. 2015), stellar mass (Meneux et al. 2008; Marulli et al. 2013; Beutler et al. 2013; Dolley et al. 2014; Skibba et al. 2015; Durkalec et al. 2018), star formation rate (SFR; Hartley et al. 2010; Mostek et al. 2013), colour (Zehavi et al. 2005; Coil et al. 2008; Skibba et al. 2014), and spectral type (Norberg et al. 2002; Meneux et al. 2006). The general conclusion from all these studies is that galaxy clustering strongly depends on galaxy properties. More luminous, massive, redder, and early-type galaxies exhibit stronger clustering, and tend to exist in denser regions of the universe, than their less massive, bluer, and later-type counterparts.

These observations can be explained in the framework of the ΛCDM cosmology and hierarchical model of structure formation. Small density fluctuations in the early universe evolved under gravity to form the present LSS (Springel et al. 2005). The initially stronger over-densities evolved faster, thus resulting in the formation of self-bound clumps of dark matter called dark matter haloes. Such haloes provided the gravitational potential to trap the baryonic matter and thereby form galaxies at their centres (Press & Schechter 1974; White & Rees 1978). Therefore, it is expected that the properties of the parent halo play a major role in defining galaxy properties such as luminosity, stellar mass, colour, and SFR. The halo mass is believed to greatly influence the halo clustering and hence the properties of hosted galaxies (e.g. Zheng et al. 2005; More et al. 2009; Gu et al. 2016); this is frequently referred to as ‘halo bias’. Large and massive haloes have potentials strong enough to form bigger, more massive, and more luminous galaxies. However, it has also been shown that the clustering of haloes has a dependence on properties other than halo mass–in large part the halo assembly history that is commonly referred to as ‘halo assembly bias’ (Zentner et al. 2014; Croton et al. 2007; Mao et al. 2018). This means that halo properties correlate with environment (Sheth & Tormen 2004). These two dependences, that is those between halo properties and galaxy properties, and between halo properties and environment, prompt a correlation between galaxy properties and environment (see Wechsler & Tinker 2018, for a review).

Studies of these dependences between galaxy properties and environment are crucial to understand structure formation in the universe, and over the past decade there has been remarkable progress in the development of galaxy formation models describing connections between dark matter haloes and their galaxies (see Somerville & Davé 2015, for a review). Methods used in these models include numerical hydrodynamic techniques (e.g. McCarthy et al. 2012; Vogelsberger et al. 2013; Kannan et al. 2014), semi-analytic models (SAM; e.g. Baugh 2006; Benson 2012; Linke et al. 2020), and even empirical methods in which physical constraints are taken entirely from observations (e.g. Yang et al. 2012; Moster et al. 2020; Grylls et al. 2020). We are now able to simulate the physics of galaxy formation and, to some extent, link galaxy properties to the host halo properties. For example, simulations have been used, such as UNIVERSEMACHINE (Behroozi et al. 2019), which reasonably parametrises the galaxy growth-halo assembly correlation, and SHARK (Lagos et al. 2018), which agrees fairly well with observations (Bravo et al. 2020). But there is still a need for improvement in understanding the process of galaxy formation (see Naab & Ostriker 2017, for a review). There has not yet been a method that would perfectly reconstruct the observed dependence of galaxy properties on halo properties and environment and the dependences of galaxy clustering on galaxy properties. Hence, a better understanding of how different galaxy properties trace the environment is needed to establish better constraints on galaxy formation and evolution models. It would be preferred if this understanding came from galaxy observations rather than simulations.

There is a problem that has to be faced here: galaxy clustering strength depends on the photometric passband in which the galaxies are selected for measurement (Milliard et al. 2007; de la Torre et al. 2007; Zehavi et al. 2011; Skibba et al. 2014). In other words, it is common that different works report various clustering strengths for galaxies selected using different methods. For example, in the optical range, Zehavi et al. (2005) measured correlation functions (CFs) in volume-limited samples of galaxies from the Sloan Digital Sky Survey (SDSS) selected in bins of the r-band absolute magnitude. These authors observed that galaxies that are brighter in the r band show stronger clustering than the fainter galaxies. Similar behaviour of clustering was observed for the B band from the VIMOS survey (Marulli et al. 2013), g band from PRIMUS (Skibba et al. 2014), and K band from HiZELS (Sobral et al. 2010). Measurements in the u band, however, show an opposite trend. Galaxies luminous in that band tend to exist in low-density regions of the universe, whereas their u-band fainter counterparts are preferentially found in high-density locations (Deng 2012). Additionally, measurements based on farther infrared (IR) indicate stronger clustering than that measured for galaxies observed at optical wavelengths (Oliver et al. 2004; Pollo et al. 2013a,b). Moreover, Heinis et al. (2004) and Milliard et al. (2007) reported weaker clustering of ultra-violet (UV) galaxies in the local universe compared to optical and IR galaxies. All these results imply that galaxies selected based on different properties trace the local environment differently.

The main aim of this paper is to demonstrate how different galaxy properties trace the small-scale galaxy clustering. In particular, we show which property can be a better tracer of galaxy environment. Environment has been defined in many different ways in the past (Muldrew et al. 2012); in this work we define it as the galaxy over-density around the object.

We also study the possible influence of selection methods on clustering results. In particular we show which photometric passband best serves as a proxy for stellar mass in the absence thereof, and how survey flux limitations can influence clustering results. Frequently, in the literature, luminosity and stellar mass are considered to be one-to-one correlated: more luminous galaxies are assumed to be more massive (Blanton & Moustakas 2009). In particular, near-infrared (NIR) luminosity is known to be a good proxy of the galaxy stellar mass (Kochanek et al. 2001). Mid-infrared (MIR) fluxes, particularly those with 3.4 μm and 4.6 μm wavelengths, are also reliable tracers of galaxy stellar mass (Jarrett et al. 2013; Cluver et al. 2014). For instance, K-band selected samples are used to construct stellar mass limited samples with high completeness (van Dokkum et al. 2006; Taylor et al. 2009). However, it is yet unclear if such a sample can be a perfect proxy for clustering measurements. Better understanding on this issue gives us a better idea on the cautions to be taken while these kinds of proxy samples are used for clustering studies.

In addition, galaxy surveys are inevitably flux-limited. Therefore, extra care has to be taken while working with stellar mass selected samples extracted from such surveys. They tend to miss galaxies that are massive enough to pass the mass selection, but not luminous enough to reach the flux limit of the survey (Meneux et al. 2008, 2009; Marulli et al. 2013). This effect makes such samples incomplete and hence can influence the galaxy clustering measurements. In this work we try to understand how the flux limit of a survey affects the measured clustering and what steps are needed to account for resulting inaccuracies.

In our work we use the Galaxy and Mass Assembly spectroscopic survey (GAMA; Driver et al. 2009). We choose GAMA over SDSS owing to its high completeness (> 98.5%) down to rpetro < 19.8 (2 mag fainter than SDSS). Moreover, GAMA does not suffer from fibre collisions that affect the close galaxy pairs in SDSS (Robotham et al. 2010). It also provides reliable measurements of absolute magnitudes in a wide wavelength range, stellar masses, and SFR. The GAMA survey has been used for various aspects of environmental effects, in particular the impact of group, cluster, local, and large-scale environment on galaxy properties (Wijesinghe et al. 2012; Burton et al. 2013; Brough et al. 2013; McNaught-Roberts et al. 2014; Alpaslan et al. 2015; Davies et al. 2016, 2019a,b; Schaefer et al. 2017, 2019; Grootes et al. 2017; Barsanti et al. 2018; Wang et al. 2018; Vázquez-Mata et al. 2020).

Our methods rely on marked statistics tools (Stoyan & Stoyan 1994), in particular on measurements of the galaxy marked correlation function (MCF), which has been proven to be very sensitive to the environment (Sheth & Tormen 2004; Skibba et al. 2013). In this method, each galaxy is assigned a ‘mark’ that is defined as any measurable property of the galaxy. The MCF accounts for the clustering of positions of the marks (i.e. galaxies of a given property). Hence, MCF measurements with different galaxy properties as marks help us to study how these different properties trace the galaxy clustering, particularly on small scales (Sheth 2005). Marked statistics has been widely used to show that closer pairs of galaxies are more luminous, redder, and older than pairs that have larger separations (Beisbart & Kerscher 2000; Skibba et al. 2006; Sheth et al. 2006). Additionally, Sheth et al. (2005) showed that these observations are in qualitative agreement with the semi-analytic galaxy formation models. Gunawardhana et al. (2018) used marked statistics on a set of luminosity- and stellar-mass-selected galaxy samples from GAMA with SFR, specific SFR (sSFR), and (g − r)rest colour as marks. They observed that sSFR is a better tracer of interactions between star-forming galaxies than colour. Riggs et al. (2021) uses MCF to explore the clustering of galaxy groups in GAMA.

In our study we compute the MCFs on stellar-mass-selected samples in the redshift range 0.1 < z < 0.16 using absolute magnitudes in the u, g, r, J, and K bands, stellar mass, SFR, and sSFR (SFR per unit stellar mass) of galaxies as marks. Additionally, we explore the effect of apparent flux limits on the correlation between small-scale clustering and galaxy properties. For this purpose, we impose various flux limits to the parent sample from the GAMA survey. Further we compare the MCFs using different marks to see how these functions are affected by the change in flux limit. We also compare the measurements in our GAMA sample with those from SDSS. We select samples from SDSS because it provides a larger number of galaxies to brighter flux limits than GAMA.

This paper is structured as follows. In Sect. 2 we describe properties of the GAMA survey and our sample selection method. Different clustering techniques and their definitions are described in Sect. 3. Clustering measurements are presented in Sect. 4. The results are discussed and compared with other works in Sect. 5 and finally concluded in Sect. 6.

Throughout the paper, a flat ΛCDM cosmological model with ΩM = 0.3 and ΩΛ = 0.7 is adopted and the Hubble constant is parametrised via h = H0/100 km s−1 Mpc−1. All galaxy properties are measured using h = 0.7. The distances are expressed in comoving coordinates and are in units of h−1 Mpc.

2. Data

2.1. Galaxy and Mass Assembly

The GAMA survey is a spectroscopic and multiwavelength galaxy survey that aims to test the ΛCDM model of structure formation and to study the galaxy evolution through the latest generation of ground-based and space-borne, wide-field survey facilities (Driver et al. 2011; Liske et al. 2015). It provides a comprehensive survey of galaxy populations by bringing together data from eight ground-based surveys and four space missions. The GAMA covers three equatorial regions named G09, G12, and G15 and two southern regions G02 and G23. Detailed descriptions of the GAMA survey are provided in Driver et al. (2009, 2011), Robotham et al. (2010), and Liske et al. (2015); below we briefly describe the survey details important in the context of our work.

We exploit the main r-band limited data from GAMA II equatorial regions with targets drawn primarily from SDSS DR7 (Abazajian et al. 2009). For extinction-corrected r-band Petrosian magnitudes (Petrosian 1976) limited at rpetro < 19.8, GAMA provides high spatial completeness and an overall redshift completeness of 98.48% in the equatorial regions. This excellent completeness of GAMA is achieved by repeated surveying of the same field (Robotham et al. 2010), thereby making GAMA an ideal survey for clustering measurements.

In this study we select galaxies from the GAMA II main survey (SURVEY_CLASS ≥ 4) in the equatorial regions with spectroscopic redshifts in the range 0.1 < z < 0.16. The redshifts of GAMA objects were measured using the software AUTOZ, as described in Baldry et al. (2014), and are corrected for the local flow via the Tonry et al. (2000) model, tapered smoothly to the cosmic microwave background rest frame for z ≥ 0.03. Liske et al. (2015) provide a detailed assessment of the quality and reliability of these redshifts. We only use secure redshifts with quality flag nQ ≥ 3, which assures that the redshift has > 90% chance of being correct. In addition to this redshift quality cut, we only select objects with VIS_CLASS = 0, VIS_CLASS = 1, or VIS_CLASS = 255, so that we avoid sources that are visually classified to be deblends of star or parts of other galaxies (Baldry et al. 2010).

As GAMA combines photometric data from several surveys, it has a very wide wavelength range, from X-ray to radio. In this work, we make use of STELLARMASSESLAMBDARV20 DMU (Taylor et al. 2011; Wright et al. 2016). The stellar masses are based on the methods of Taylor et al. (2011) applied to the LAMBDAR photometry of Wright et al. (2016). These methods are based on the stellar population synthesis (SPS) modelling of broadband photometry using stellar evolution models by Bruzual & Charlot (2003), assuming a Chabrier (2003) initial mass function and Calzetti et al. (2000) dust law. The absolute magnitudes are inferred from the SPS fits and are corrected for internal dust extinction and k-corrected to z = 0. As the fits are constrained to the rest-frame wavelength range 300−1100 nm, the J and K-band absolute magnitudes we use are extrapolations of the fit to data. The stellar masses and absolute magnitudes are fluxscale corrected as described in Taylor et al. (2011) to account for the difference in aperture matched and Sérsic photometry. Galaxies with physically unrealistic fluxscale values are not considered for our analysis. The SFRs and sSFRs are taken from the DMU MAGPHYSV06 and are estimated using the energy balance spectral energy distribution (SED)-fitting code MAGPHYS (da Cunha et al. 2008). All the quantities are derived for the concordance (ΩM, ΩΛ, h) = (0.3, 0.7, 0.7) cosmology. For clustering measurements, we use the GAMA random galaxy catalogue (RANDOMSV02 DMU) by Farrow et al. (2015). In the random catalogue, we assign stellar mass to each random galaxy by matching on the CATAID of the real galaxy.

2.2. Sample selection

Our full sample counts 34 455 galaxies in the redshift range 0.1 < z < 0.16 with apparent flux limit rpetro < 19.8 distributed over the three equatorial regions in the sky, for which the details such as the number of galaxies, median redshift, and total area are given in Table 1. The aforementioned redshift range is chosen to optimally select volume-limited samples that include low-mass galaxies. All these galaxies have reliable spectroscopic redshift and well-measured absolute magnitudes, stellar mass, SFR, and sSFR.

Table 1.

Properties of the GAMA equatorial regions used in this work.

To study the environmental dependence of luminosity in different bands, stellar mass, SFR, and sSFR, we define a nearly stellar-mass-complete sample by applying an additional stellar mass cut of log (M/M)min = 9.3. This is referred to as sample 𝒜1. As a representative example of the selection technique, in Fig. 1 we show the selection cut for the sample 𝒜1 with flux limit r < 19.8.

thumbnail Fig. 1.

Selection of galaxy sample 𝒜1 used in this work. The grey dots represent the GAMA galaxies with flux limit r < 19.8. Top and right histograms: distribution of redshift and stellar mass, respectively. The red lines represent the stellar mass cut and the redshift limit of the sample 𝒜1.

To further investigate how the environmental dependence of the galaxy properties varies with different flux limits, we select two stellar-mass-selected samples in the same redshift range with a stellar mass cut of log (M/M)min = 10.4 but different flux limits: sample ℬ1 with r < 19.8 and ℬ2 with r < 17.8. The stellar mass and redshift distribution of these samples are shown in Fig. 2. The definition and properties of all the selected samples are given in Tables 2 and 3 respectively. All the selected samples contain a sufficient number of galaxies for reliable clustering measurements. For all the selected samples, random samples are selected from RANDOMSV02 DMU (Farrow et al. 2015) after applying corresponding r-band apparent magnitude cut and stellar mass cut.

thumbnail Fig. 2.

Redshift and stellar mass distribution of galaxies in stellar-mass-selected samples ℬ1 (r < 19.8; brown dots) and ℬ2 (r < 17.8; green circles). Top and right histograms: distribution of redshift and stellar mass, respectively.

Table 2.

Definitions of the galaxy samples in the redshift range 0.1 < z < 0.16, as used in this study.

2.3. SDSS samples for comparison

Apart from the comparisons between ℬ1 and ℬ2, we also compare the results with the SDSS (York et al. 2000) to understand better how brighter flux limits can affect our measurements. We use the LSS catalogue and the corresponding random catalogue1 generated from SDSS III Baryon Oscillation Spectroscopic Survey (BOSS; Dawson et al. 2013) Data Release 12 (DR12; Alam et al. 2015). The SDSS BOSS DR12 catalogue encompasses massive galaxies partitioned into two non-overlapping redshift bins named as ‘LOWZ’ and ‘CMASS’ that cover galaxies in the redshift ranges z < 0.43 and z > 0.43, respectively. Reid et al. (2016) describe the methods used in the target selection of the SDSS galaxy data sets, and give details of the MKSAMPLE code used to create the LSS catalogue and random catalogue. The total sky coverage of the LOWZ DR12 sample is 8337.47 deg2. In this work, we make use of LOWZ galaxies in the North Galactic Cap in the redshift range 0.1 < z < 0.16.

All the galaxies in the LSS catalogue are assigned the r-band apparent magnitude from SPECPHOTOALL table by matching the angular position within 2″. Stellar masses are then assigned by cross-matching with the table STELLARMASSSTARFORMINGPORT using specObjID. The stellar masses are estimated from the best-fit SED obtained from the stellar population model of Maraston et al. (2009). The fits are performed on the observed ugriz-magnitudes of BOSS galaxies with the spectroscopic redshift determined using an adaptation of Hyper-Z code of Bolzonella et al. (2000). The magnitudes used are extinction-corrected model magnitudes that are scaled to the i-band cmodel magnitude. The GAMA and SDSS stellar masses are derived using different photometry. Despite the systematic differences between SDSS and GAMA photometry, we find relatively good agreement between GAMA and SDSS stellar masses of overlapping galaxies, taking into account that the common sample is rather small (388 objects). The median offset between the masses is 0.18 dex with a scatter of the order of 0.2 dex. Hence we use the same stellar mass cuts in GAMA and SDSS to define samples for comparison.

For better statistics, we fix the brightest magnitude limit in our work to be r < 16.8. For the comparison between GAMA and SDSS, we define five stellar-mass-selected samples with the same stellar mass cut of log (M/M) > 10.8, but different flux limits. This gives samples 𝒞1, 𝒞2, 𝒞3 from our parent sample in GAMA with flux limits of r < 19.8, 18.8, 17.8, respectively. From SDSS, we have samples 𝒞4, 𝒞5 with flux limit of r < 17.8, 16.8, respectively. More details of these samples are given in Tables 2 and 3, respectively, and their mass distribution is shown in Fig. 3. The SFRs, sSFRs, and K-band absolute magnitudes of SDSS samples mentioned in Table 3 are derived from the same SED fits from which SDSS stellar masses are derived. The absolute magnitudes in the u, g, and r bands of SDSS samples are taken from the PHOTOZ table (Beck et al. 2016). The angular distribution of random galaxies for each SDSS sample are taken from the LSS random catalogue and the redshift is randomly assigned from a smoothened N(z) distribution of the corresponding real galaxy sample.

thumbnail Fig. 3.

Redshift and stellar mass distribution of stellar-mass-selected samples mentioned in Table 2. Left panel: 𝒞1 (r < 19.8; black dots), 𝒞2 (r < 18.8; red circles), and 𝒞3 (r < 17.8; green squares) galaxy samples from GAMA and the right panel shows the 𝒞4 (r < 17.8; blue dots) and 𝒞5 (r < 16.8; orange squares) galaxy samples from SDSS. Top and right histograms of both the panels show the distribution of redshift and stellar mass, respectively.

Table 3.

Properties of the galaxy samples defined in Table 2.

3. Measurement methods

3.1. The galaxy two-point correlation function

The galaxy 2pCF, ξ(r), is a statistical tool used to measure the clustering of galaxies. It is defined as the excess probability above random of observing a pair of galaxies at a given spatial separation r in a volume element dV (Peebles 1980), that is

d P = n [ 1 + ξ ( r ) ] d V , $$ \begin{aligned} \mathrm{d}P = n\,[1 + \xi (r)]\,\mathrm{d}V, \end{aligned} $$(1)

where n is the number density of galaxies.

It has been observed that CF mostly follows a power law (Groth & Peebles 1977) given by

ξ ( r ) = ( r r 0 ) γ , $$ \begin{aligned} \xi (r)=\left(\frac{r}{r_0}\right)^{-\gamma }, \end{aligned} $$(2)

where r0 and γ are the correlation length and slope, respectively.

In practice, owing to the limitations of galaxy surveys, various estimators of ξ(r) have been proposed to minimise the effects related to the limited number of objects and limited survey areas (e.g. Davis & Peebles 1983; Hamilton 1993). The Landy & Szalay (1993) estimator is the most widely used because of its capability to minimise the above-mentioned problems and is defined by,

ξ ( r ) = D D ( r ) 2 D R ( r ) + R R ( r ) R R ( r ) , $$ \begin{aligned} \xi (r)= \frac{\langle DD(r) \rangle - 2 \langle DR(r) \rangle + \langle RR(r) \rangle }{\langle RR(r) \rangle } , \end{aligned} $$(3)

where DD(r) is the observed number of galaxy-galaxy pairs with the separation in the bin centred at r in the real galaxy sample, RR(r) is the expected number of such pairs from a random galaxy distribution, DR(r) is the number of cross pairs of galaxies between the real and random sample, and ⟨⟩ refers to the quantity normalised by the total number of such pairs. The random galaxy sample reflects the same sky distribution and redshift distribution of the real galaxy sample. The number of random galaxies used for the computation is set to be significantly greater (5−10 times) than the number of real galaxies to avoid shot noise on smaller scales.

To account for the distortions in CF measurements caused by galaxy peculiar velocities, the co-moving redshift space separation between the galaxies is split into two components: parallel (π) and perpendicular (rp) to the line of sight. The CF thus takes form of a two-dimensional grid ξ(rp, π). Integrating ξ(rp, π) over the line-of-sight (π) direction gives us the projected 2pCF, ωp(rp), which can be used to recover the real space CF devoid of redshift space distortions (Davis & Peebles 1983). It is defined as

ω p ( r p ) = 2 0 π max ξ ( r p , π ) d π . $$ \begin{aligned} \omega _{\rm p}(r_{\rm p}) = 2 \, \int _0^{\pi _{\rm max}} \xi (r_{\rm p},\pi ) \, \mathrm{d}\pi . \end{aligned} $$(4)

The limit of integration πmax has to be reasonable enough to include all the correlated pairs and reduce the noise in the estimator. Following Appendix B of Loveday et al. (2018), we choose the value of πmax to be 40 h−1 Mpc.

There have been many studies in GAMA using 2pCF in the past. The dependence of projected galaxy clustering on various properties was studied by Farrow et al. (2015). Loveday et al. (2018) used 2pCF to measure the pairwise velocity distribution in a set of luminosity-selected samples from GAMA. The small-scale clustering properties of star-forming galaxies were used by Gunawardhana et al. (2018) to study the interactions between galaxies. Christodoulou et al. (2012) used CF as a tool to check the robustness of their photometric redshift estimates. Jarrett et al. (2017) analysed the spatial distribution of mid-infrared Wide-field Infrared Survey Explorer (WISE) sources observed in the G12 equatorial region of GAMA. Large-scale clustering of radio galaxies in the Very Large Array Faint Images of the Radio Sky at Twenty-cm (FIRST) survey over the GAMA survey area was studied by Lindsay et al. (2014). The clustering measurements in GAMA were also used by Alam et al. (2021) to model the redshift space distortions. The clustering properties of low-redshift (z < 0.3) submillimeter galaxies detected at 250 μm in the Herschel-ATLAS (Eales et al. 2010) using the redshift information from GAMA was carried out by van Kampen et al. (2012). van Uitert et al. (2018) used angular CF of GAMA galaxies as one of the probes to constrain cosmological parameters. In our work, we examine how galaxy clustering depends on various galaxy properties, such as luminosities in different passbands, stellar mass, and SFR, using MCF.

3.2. Marked correlation function

The 2pCF characterises the galaxy clustering. It can, and successfully has been (as described in Sect. 1), used to study clustering dependences on various properties of galaxies. This is done by defining the galaxy samples based on the property of interest (e.g. luminosity, colour, or stellar mass). However, after the selection and further during the analysis, these properties are left unconsidered and each galaxy is weighted equally during CF measurements. On the other hand, marked statistics allows us to study the properties of galaxy clustering by taking the physical properties (called marks) of each galaxy in the sample into account. These marks can be discrete or continuous values such as luminosity, colour, stellar mass, SFR, and morphology (Sheth et al. 2005).

The MCF allows for the efficient study of the spatial distribution of galaxy properties and their correlation with the environment (Skibba et al. 2013). The two-point MCF is defined as

M ( r ) = 1 + W ( r ) 1 + ξ ( r ) , $$ \begin{aligned} M(r) = \frac{1 + W(r)}{1 + \xi (r)}, \end{aligned} $$(5)

where ξ(r) is the galaxy 2pCF defined by Eq. (3) and W(r) is the weighted CF obtained with the same estimator, but with pair counts computed by weighting each real galaxy in the pair. That is

W ( r ) = W W ( r ) 2 W R ( r ) + R R ( r ) R R ( r ) · $$ \begin{aligned} W(r)= \frac{\langle WW(r) \rangle - 2 \langle WR(r) \rangle + \langle RR(r) \rangle }{\langle RR(r) \rangle }\cdot \end{aligned} $$(6)

We adopt multiplicative scheme for pair weighting, that is

W W ( r ) = ij w i × w j , $$ \begin{aligned} WW(r) = \sum _{ij} w_i \times w_j , \end{aligned} $$(7)

where wi is the weight of the ith galaxy given by the ratio of its mark to the mean mark across the sample.

The projected two-point MCF is defined as

M p ( r p ) = 1 + W p ( r p ) / r p 1 + ω p ( r p ) / r p · $$ \begin{aligned} M_{\rm p}(r_{\rm p}) = \frac{1 + W_{\rm p}(r_{\rm p})/r_{\rm p}}{1 + \omega _{\rm p}(r_{\rm p})/r_{\rm p}}\cdot \end{aligned} $$(8)

Essentially, MCF at a scale r tells us if galaxies in pairs separated by r tend to have larger or smaller values of their mark than the mean mark in the entire sample (Sheth & Tormen 2004).

3.3. Rank-ordered marked correlation function

For a given property, a stronger MCF signal at a certain scale indicates greater probability of finding galaxy pairs for which the given property has a larger value for both the galaxies. Hence the property that is more dependent on environment is that corresponding to a larger MCF. However, comparing different MCFs obtained using different properties is not straightforward (Skibba et al. 2013). When computing the MCF in a traditional approach, the value of the physical property of a galaxy is considered as its mark and the CF is directly weighted by the ratio of the given mark to the mean mark of the sample. Hence, the amplitude of MCF depends on the distribution of the marks and the variations in their formulation (e.g. log or linear). This makes it impossible to directly compare different MCFs measured using different properties as marks if these properties have different distribution or formulation. Skibba et al. (2013) developed a solution to this problem. Each galaxy is given a rank based on the relative strength of the value of its property in the sample; that is a galaxy with the lowest value is given the lowest rank and another one with a greater value is given a higher rank. This is called rank-ordering the marks. The rank of each galaxy is then used as its mark to weight the CF. Since all the ranks have a uniform distribution on [1, N], the amplitudes of the MCFs thus obtained using different properties can be compared. However, since the weight is given by the rank rather than the property value, any information contained in the shape of the distribution of the property is lost. As we are interested in relative importance of different properties for correlation measurements, all the MCFs shown in this work are rank ordered.

3.4. Error estimates

Since the galaxies are clustered and the pair counts in different bins of rp can include the same galaxies, the values of ωp for different bins are correlated. Hence the statistical errors associated with clustering measurements are estimated using the covariance matrix obtained from various methods of internal error estimation. In our work, we use the jackknife resampling method (Norberg et al. 2009), in which we divide the entire sky region into Njk subsamples of equal area. Then Njk different jackknife copies of the parent sample are created by omitting one of these subsamples in turn. Then the CF is measured in each jackknife copy. We adopt Njk = 12, which we found to be an optimal number so that the size of each subsample is larger than the maximum scale at which we measure ωp(rp) (∼10 h−1 Mpc).

The associated covariance matrix is given by

C ij = N jk 1 N jk k = 1 N jk ( ω p k ( r i ) ω ¯ p ( r i ) ) ( ω p k ( r j ) ω ¯ p ( r j ) ) , $$ \begin{aligned} C_{ij} = \frac{N_{\rm jk} - 1}{N_{\rm jk}}\sum \limits _{k=1}^{N_{\rm jk}} \left(\omega _{\rm p}^k(r_i) - \bar{\omega }_{\rm p}(r_i) \right) \,\, \left(\omega _{\rm p}^k(r_j) - \bar{\omega }_{\rm p}(r_j) \right), \end{aligned} $$(9)

where ω p k ( r j ) $ \omega_{\rm p}^k(r_j) $ represents the measurement of ωp at rp = rj in the kth jackknife copy and ω p ¯ $ \bar{\omega_{\mathrm{p}}} $ is the average from Njk copies. The square root of the diagonal elements of Cij gives the error bar for the ωp at the corresponding bin.

We estimate r0 and γ (the power-law fit parameters of ξ(r)) from the projected function ωp(rp). Using the parametrisation given in Eq. (2), the integral in Eq. (4) can be analytically performed to give the power-law fit parameters as

ω p ( r p ) = r p ( r p r 0 ) γ Γ ( 1 2 ) Γ ( γ 1 2 ) Γ ( γ 2 ) , $$ \begin{aligned} \omega _{\rm p}(r_{\rm p}) = r_{\rm p}\,\left(\frac{r_{\rm p}}{r_0} \right)^{-\gamma }\,\frac{\Gamma \left(\frac{1}{2}\right)\,\Gamma \left(\frac{\gamma -1}{2} \right)}{\Gamma \left(\frac{\gamma }{2}\right)}, \end{aligned} $$(10)

where Γ(n) is Euler’s Gamma function (Davis & Peebles 1983).

The power-law fit parameters are usually estimated by minimising the generalised χ2 using the inverse of full covariance matrix (Fisher et al. 1994; Pollo et al. 2005). But, limited number of jackknife resamplings can introduce noise to the non-diagonal elements of covariance matrix. This can lead to unreliable power-law fit parameters. Hence, we adopt the fitting procedure previously done on GAMA data by Farrow et al. (2015). For that, we first normalise the covariance matrix to give the correlation matrix ( C ij $ \tilde{C}_{ij} $) given by

C ij = C ij σ i σ j , $$ \begin{aligned} \tilde{C}_{ij} = \frac{C_{ij}}{\sigma _i \sigma _j} , \end{aligned} $$(11)

where σ i = C ii $ \sigma_i = \sqrt{C_{ii}} $, the error bar associated with the ωp measurement at the ith bin.

The correlation matrix C $ \mathbf{\tilde{C}} $ is then transformed into a matrix D using singular value decomposition (SVD) given by C = U T D U $ \mathbf{\tilde{C}} = \mathbf{U}^{\mathrm{T}}\,\mathbf{D}\,\mathbf{U} $. The diagonal matrix D has λ ij 2 δ ij $ \lambda^2_{ij}\,\delta_{ij} $ as elements, where δij is the Kronecker delta. Columns of the matrix U are the eigen modes of the correlation matrix, where λ ij 2 δ ij $ \lambda^2_{ij} \delta_{ij} $ are their corresponding eigen values. The inverse of the diagonal matrix D is given by D ij 1 = ( 1 / λ ij 2 ) δ ij $ D^{-1}_{ij} = (1/\lambda^2_{ij})\,\delta_{ij} $. Then we transform the correlation matrix back using C 1 = U D 1 U T $ \mathbf{\tilde{C}^{-1}} = \mathbf{U}\,\mathbf{D}^{-1}\,\mathbf{U}^{\mathrm{T}} $ and while doing this calculation, we set D ii 1 = 0 $ D^{-1}_{ii} = 0 $ for those eigen values λ ii 2 < 2 / N jk $ \lambda_{ii}^2 < \sqrt{2/N_{\mathrm{jk}}} $ (Gaztañaga & Scoccimarro 2005). This effectively removes the influence of least significant eigen modes that are more likely to suffer from noise.

The transformed inverse of the correlation matrix is then used to minimise the χ2, defined as

χ 2 = i , j ( ω p mod ( r i ) ω p ( r i ) ) σ i C ij 1 ( ω p mod ( r j ) ω p ( r j ) ) σ j , $$ \begin{aligned} \chi ^2 = \sum \limits _{i,j} \frac{\left(\omega _{\rm p}^\mathrm{mod}(r_i) - \omega _{\rm p}(r_i)\right)}{\sigma _i} \,\, \tilde{C}^{-1}_{ij} \,\, \frac{\left(\omega _{\rm p}^\mathrm{mod}(r_j) - \omega _{\rm p}(r_j)\right)}{\sigma _j}, \end{aligned} $$(12)

where ωp(ri) is the measured value of CF at rp = ri and ω p mod $ \omega_{\mathrm{p}}^{\mathrm{mod}} $ is the power-law model value given by Eq. (10). The uncertainties in the power-law parameters r0 and γ presented in this work are defined by the 68.3% joint confidence levels (Chapt. 15.6 of Press et al. 1992). Gaztañaga & Scoccimarro (2005) and Marín et al. (2013) provide detailed descriptions of the fitting procedure using SVD to estimate the power-law parameters.

As the errors in Wp(rp) and ωp(rp) are strongly correlated, simply summing these in quadrature gives an overestimate for the error in Mp(rp). A much better approximation of the uncertainty is obtained by randomly scrambling the marks among the galaxies and remeasuring Mp. This is repeated ∼100 times and the standard deviation around the mean gives the uncertainty in Mp (Skibba et al. 2006).

4. Results

In this section, we present our results of the environmental dependence of galaxy luminosity, stellar mass, SFR, and sSFR. We measure 2pCFs and rank-ordered MCFs for various stellar-mass-selected samples described in Table 2 selected from GAMA in the redshift range 0.1 < z < 0.16. Each 2pCF is fitted with a power-law model and the MCFs are measured using eight different properties as marks: absolute magnitudes in u, g, r, J, and K bands, stellar mass, SFR, and sSFR. Details are described in Sect. 3, especially Eqs. (4) and (8). In all the samples, we could reliably measure the CFs in the range 0.1 < rp < 10 h−1 Mpc. The errors in ωp(rp) are obtained from 12 jackknife realisations and the errors in Mp(rp) are obtained by randomising the marks (100 times) as described in Sect. 3.4. The best-fitting power-law parameters for 2pCFs in all the samples are given in Table 4.

Table 4.

Best-fitting power-law parameters for all the galaxy samples used in this work.

4.1. Two-point and marked correlation functions

In Fig. 4 we show the 2pCF and MCFs obtained for galaxies with a flux limit of r < 19.8 (sample 𝒜1). The left panel shows the projected 2pCF ωp(rp), which at first approximation obeys a power-law model. The best-fit parameters are correlation length r0 = 5.06 ± 0.70 h−1 Mpc and slope γ = 1.78 ± 0.05. This can be compared with the results of Farrow et al. (2015), although their samples vary from ours. To be compared with our sample 𝒜1, the most appropriate sample of theirs is the sample with mass limits 10 < log M/Mh−2 < 10.5 in the redshift range 0.14 < z < 0.18. The parameter values for that sample were measured to be r0 = 5.94 ± 0.46 h−1 Mpc and γ = 1.86 ± 0.05. The correlation length varies by a factor of 1.1σ from our result. This small deviation in correlation length is expected as their stellar mass cuts vary from ours.

thumbnail Fig. 4.

Projected 2pCF ωp(rp) with a power-law fit (filled markers; left panel) and rank-ordered projected MCFs Mp(rp) (unfilled markers; right panel) for sample 𝒜1 described in Sect. 2.2. In the left panel, error bars of ωp(rp) are the square root of the diagonals of the covariance matrix obtained from the jackknife resampling method. The inset in the left panel shows the power-law fit parameters (filled star) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). In the right panel the different symbols represent measurements with different marks (as labelled), and the error bars are obtained by random scrambling of the marks. The error bars of Mp(rp) are too small to be visible.

The right panel of Fig. 4 presents rank-ordered MCFs Mp(rp) obtained for different galaxy properties used as marks (as described in the legend). The presented MCFs strongly deviate from unity on small scales (rp < 1 h−1 Mpc) for all luminosity, stellar mass, SFR, and sSFR marks. This deviation then decreases, but still remains at larger scales. In general, a stronger deviation of MCF from unity, means stronger correlation of the corresponding galaxy property with the environment (as described in Sect. 3.3). As shown in Fig. 4, stellar mass MCF deviation from unity is greater than any of the luminosity MCFs. This means that stellar mass catches the galaxy over-density in small scales better than other properties used in this work. Thus, stellar mass can be considered as a more direct tracer of environment than luminosity and SFR.

Among the MCFs measured using luminosities in different passbands, the K-band MCF has the highest amplitude and the u-band MCF has the lowest. This means that the K-band luminosity traces the environment better than any other band of those used in this work. Moreover, the K-band MCF shows similar, although not exactly the same, behaviour as the stellar mass MCF; the K-band luminosity and stellar mass are therefore correlated with environment in a similar fashion. This, in turn, confirms that K-band luminosity can be used as a tracer (or a proxy) of stellar mass when we lack the measurements of stellar mass (Kochanek et al. 2001; Baldry et al. 2012).

On the other hand, SFR and sSFR MCFs behave opposite to the stellar mass MCF; this presents low values (high deviation from unity) on small scales and high values (close to unity) on large scales. This means that there is only a small number of close pairs of galaxies with strong star formation activity. That is the densest regions of the local universe are mainly populated by old or quiescent galaxies. Similar, although weaker, behaviour is presented by the u-band MCF. This suggests that the u-band luminosity traces galaxy SFRs (or at least serves as a proxy; Madau & Dickinson 2014).

4.2. Dependence on flux limit

To study the impact of survey flux limit on the correlations observed in Sect. 4.1, we select two distinct samples with the same stellar mass cut log (M/M) > 10.4, but different flux limits: sample ℬ1 with r < 19.8 and sample ℬ2 with r < 17.8 (see Table 2 for details). In Fig. 5, we show the measurements of MCFs for both of these galaxy samples. A direct comparison of results is shown in Fig. 6, where we present a 2pCFs (left panel) and the ratio between the MCFs (right panel) measured for sample ℬ2 and sample ℬ1 using different galaxy properties as marks.

thumbnail Fig. 5.

Rank-ordered projected MCFs Mp(rp) for samples ℬ1 (left) and ℬ2 (right) with different flux limits, as labelled. The different symbols represent different marks considered for the Mp(rp) measurement (as labelled) and the error bars are obtained by random scrambling of the marks.

thumbnail Fig. 6.

Projected 2pCFs ωp(rp) for samples ℬ1 and ℬ2 (filled markers; left panel) and the ratio Mℬ2/Mℬ1 between MCFs obtained for these two samples (unfilled markers; right panel). The different symbols in the right panel represent the ratio between different MCFs measured with corresponding galaxy property chosen as a mark. Small offsets along x-axis have been added for clarity. The error bars for ωp(rp) are obtained from jackknife resampling method. The inset in the left panel shows the power-law fit parameters (filled stars) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). In the right panel, the errors are calculated in quadrature.

For the same flux limit (r < 19.8), sample ℬ1 exhibits a larger correlation length than 𝒜1, although within 1σ (see Table 4). By definition, 𝒜1 contains many more less massive galaxies than ℬ1. Our observation that ℬ1 shows slightly stronger clustering than 𝒜1 therefore agrees with the common observation that more massive galaxies exhibit stronger clustering than less massive ones (e.g. Skibba et al. 2015; Durkalec et al. 2018). We also observe in the left panel of Fig. 5 that the stellar mass and g, r, J, and K-band MCFs of ℬ1 sample exhibit lower amplitude than those in 𝒜1, with g, r, J, and K-band MCFs falling below unity.

When it comes to sample ℬ2, most of the MCFs show a higher amplitude relative to ℬ1 (as shown in the Fig. 5). This can also be observed in the right panel of Fig. 6 in which the ratio of MCFs of ℬ2 to that of ℬ1 are above unity on most of the scales. This can be associated with the stellar mass incompleteness effect and Sect. 5.5 deals with a discussion on this effect. The stellar mass can still be used as a good indicator of galaxy environment; MCF, with the stellar mass used as mark, exhibits stronger deviation from unity than MCF with any of the luminosity marks.

However, the change in the flux limit does not affect our 2pCF measurements. The correlation lengths of samples ℬ1 and ℬ2 agree within 1σ. The best-fit power-law parameters for the sample ℬ1 take the values of r0 = 6.10 ± 0.86 h−1 Mpc and γ = 1.86 ± 0.06, whereas the same parameters for sample ℬ2 are r0 = 5.64 ± 0.70 h−1 Mpc and γ = 1.80 ± 0.07.

4.3. Comparison with the SDSS

To further check the effect of survey flux limit on clustering measurements, we extend our studies to even lower magnitude cuts. For that we use data from the SDSS survey. We select a total number of five samples: three from GAMA (𝒞1, 𝒞2, 𝒞3) and two from SDSS (𝒞4, 𝒞5). Each of these samples have the same stellar mass cut (log (M/M) > 10.8), but different flux limits. The details are given in Table 2. Figure 7 shows the results of 2pCF and MCF (with stellar mass used as mark) measurements for each of these samples. Additionally, the best-fit power-law parameters for each 2pCF are given in Table 4.

thumbnail Fig. 7.

Projected 2pCFs (filled markers; left panel) and stellar mass MCFs (unfilled markers; right panel) in GAMA and SDSS surveys. The different symbols represent the measurements in different samples as labelled. The error bars for ωp(rp) are obtained from jackknife resampling method and those for Mp(rp) are obtained by random scrambling of the marks. The inset in the left panel shows the power-law fit parameters (filled stars) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). Small offsets along x-axis have been added for clarity.

The correlation lengths obtained for the GAMA samples 𝒞1, 𝒞2, and 𝒞3 are comparable within 1σ. In case of MCFs, the GAMA samples (𝒞1, 𝒞2 and 𝒞3) agree between each other within the error bars. However, the stellar mass MCFs of samples 𝒞3 and 𝒞4 (with the same flux limit and stellar mass limit) are noticeably different from each other on most of the scales although their correlations lengths agree with each other. We also observe significant differences between correlation lengths and MCFs of the SDSS samples 𝒞4 and 𝒞5.

5. Discussion

5.1. Environmental dependence of galaxy properties

The MCF is a useful tool to study the environmental dependence of galaxy properties. Empirically, it is a ratio of terms involving the weighted and the unweighted CF. For a given galaxy property, all galaxy pairs in a sample carry a weight, that is the product of property values for both galaxies in units of its mean value. At each spatial scale, the MCF signal depends on the value of the weights at that scale; the larger amplitude of rank-ordered MCF, measured using a particular galaxy property, implies it has a stronger correlation with environment.

As described in Sect. 4.1 and shown in Fig. 4, we observe different amplitudes for MCFs measured using different galaxy properties. All galaxy properties (luminosities, stellar mass, SFR, and sSFR) correlate with the environment, each in different way. For example, the stellar mass and luminosity (in the g, r, J, and K bands) MCFs take values higher than unity on small scales, indicating an abundance of close galaxy pairs with these property values greater than the average of the sample. This agrees with the well-known observation that the most massive and luminous galaxies (in g, r, J, and K bands) are mostly found in dense regions (Norberg et al. 2002; Coil et al. 2006; Pollo et al. 2006; Meneux et al. 2009; Bolzonella et al. 2010; Abbas & Sheth 2006; Abbas et al. 2010; Zehavi et al. 2011; Marulli et al. 2013; Skibba et al. 2014; Farrow et al. 2015; Durkalec et al. 2018; Cochrane et al. 2018). This phenomenon can be explained in the framework of the hierarchical structure formation (Press & Schechter 1974; White & Rees 1978; Mo & White 1996; Springel et al. 2005).

5.2. u-band luminosity and star formation rate dependence on the environment

We observe that the u-band marked correlation shows different behaviour in comparison to other passbands (g, r, J, and K, see Fig. 4). The u-band MCF takes values smaller than unity on scales rp < 3 h−1 Mpc, indicating low probability of finding pairs of two galaxies similarly bright in this band. This is in complete opposition to the results from other bands, where the most g, r, and K luminous galaxies were the most strongly clustered. This special behaviour of the u-band MCF has also been observed by Deng (2012) and it agrees with the semi-analytical galaxy formation models (see Fig. 2 of Sheth et al. 2005).

The u-band and UV light are thought to be primarily emitted by starburst galaxies with young stellar populations (Cram et al. 1998). The UV-selected galaxies exhibit low clustering in the local universe (Heinis et al. 2004, 2007; Milliard et al. 2007). Heinis et al. (2004) measured a correlation length of r 0 = 3 . 2 2.3 + 0.8 h 1 $ r_0= 3.2^{+0.8}_{-2.3}\,h^{-1} $ Mpc for low-redshift UV galaxies from the FOCA survey and Milliard et al. (2007) found it to be 3.7 ± 0.6 Mpc from the GALEX survey. Both these observations confirm the weak clustering of UV-selected galaxies in the local universe. This result agrees with the work done by Barsanti et al. (2018) using GAMA groups, in which they found a rise in star formation in galaxies that are located away from the group centre compared to those in the central regions.

Similarly, SFR and sSFR MCFs also have values smaller than unity on all scales. The similar behaviour of SFR and sSFR MCFs show that both these properties correlate with environment in a similar way. This means that active star-forming galaxies are rarely found in close proximity to each other. This observation can be connected to the known fact that, at low redshifts, active star formation takes place mainly in low-density environments (Lewis et al. 2002; Gómez et al. 2003): less evolved (young) galaxies that formed in less dense regions exhibit strong star formation activity. Our observations agree with the SFR MCF measurements by Sheth et al. (2005) and, to some extent, with recent GAMA studies by Gunawardhana et al. (2018); their results show however weaker deviation of SFR and sSFR MCF from unity than ours. This can be due to the apparent absence of the SFR-density relationship as a consequence of selecting the star-forming sample of galaxies (McGee et al. 2011; Wijesinghe et al. 2012). Gunawardhana et al. (2018) measurements were made for samples of actively star-forming galaxies.

The similar behaviour of u-band and SFR MCFs has a practical interpretation. The SFR of a galaxy can be estimated by applying a scaling factor to the luminosity measurements sensitive to star formation (Condon 1992; Kennicutt 1998; Madau & Dickinson 2014). Since the u-band light is dominated by starburst galaxies with young stellar populations, it is more closely correlated to SFR than to stellar mass in galaxies (Hopkins et al. 2003) and is hence considered to be an indicator of SFR. This correlation is reflected in our results; u-band MCF follows SFR MCF in tracing the galaxy environment at small scales. This suggests that u-band luminosity can be a good proxy of SFR in the context of galaxy clustering.

5.3. The most reliable tracer of galaxy environment

The amplitude of rank-ordered MCFs, computed using various marks, can be used to find the galaxy property with the strongest environmental dependence. Our observations suggest that this parameter is the stellar mass. We therefore agree with the past results that the distribution of massive galaxies is strongly correlated with dark matter over-densities (Kauffmann et al. 2004; Scodeggio et al. 2009; Davidzon et al. 2016). This dependence is expected from the hierarchical structure formation theory according to which the local density contrast is connected to the properties of the hosting halo, which is further related to the galaxy stellar mass (Moster et al. 2010; Wechsler & Tinker 2018). The strong correlation between stellar mass and environment has been taken into account in studies that explore the environmental dependence of galaxy properties (e.g. Kauffmann et al. 2004; Peng et al. 2010).

Our observations can be also interpreted in terms of galaxy evolution. Stellar mass plays an important role in shaping galaxy star formation history (Gavazzi et al. 1996; Kauffmann et al. 2003; Heavens et al. 2004), which is yet another parameter that is strongly correlated with the local environment of a galaxy (Kauffmann et al. 2004; Blanton et al. 2005).

It is important to note a feature of our studies that might influence our results, namely the galaxy sample selection. In this work we measure MCFs in stellar mass selected samples only. That is our samples are nearly complete only in stellar mass and not in other properties. Our observation – stellar mass MCF being more enhanced than other MCFs – could be an outcome of this selection. However, in our preliminary analysis (Sureshkumar et al. 2020), we measured the same MCFs, for u, g, r, and K-band luminosities, and stellar mass marks, using samples selected based on the corresponding property. We observed similar trends as here, ruling out the effect of sample selection on our observation.

5.4. K-band luminosity as a proxy for stellar mass

In our work we also studied the behaviour of MCFs measured using different photometric wavebands and we observe clear differences between MCFs marked with u, g, r, J, and K luminosities (see Fig. 4). This means that luminosity in different passbands correlates differently with environment. Most significantly, the K band (the reddest considered) shows a stronger environmental dependence than bluer bands. This means that there is higher probability to find close pairs of galaxies similarly luminous in the K band than in other bands. In other words, galaxies luminous in the K band are strongly clustered. This observation is in agreement with various clustering studies (e.g. Oliver et al. 2004; de la Torre et al. 2007). In particular, Sobral et al. (2010), using sample of Hα emitters from HiZELS survey, show a strong increase in galaxy clustering with increasing K-band luminosity, but a weak trend in case of the B band. This difference in clustering strength between various bands is reflected in the varying amplitude of MCFs in Fig. 4.

Comparing the amplitudes of stellar mass and K-band MCFs, we observe that these galaxy properties trace the environment in a similar way. This means that the K-band luminosity can be used as the second-most reliable galaxy property (among those considered here), after galaxy stellar mass, to trace the environment. In other words, a sample that is complete in K-band luminosity can be a good substitute of a stellar-mass-complete sample.

This observation agrees with the existing results. It has been shown that longer-wavelength luminosities (e.g. K band) are dominated by evolved stellar populations and are least affected by dust extinction, making them directly related to the galaxy stellar mass (Cowie et al. 1994; Gavazzi et al. 1996; Kauffmann & Charlot 1998; Kochanek et al. 2001; Baldry et al. 2012; Jarrett et al. 2013; Cluver et al. 2014). Kauffmann & Charlot (1998) pointed out that infrared light is a much more robust tracer of stellar mass than optical light out to z ∼ 1−2. They observed that galaxies of the same stellar masses have the same K-band luminosities, independent of their star formation histories. Additionally, near-IR (mainly K-band) luminosity functions were well utilised to estimate stellar mass distribution in the local universe (Cole et al. 2001; Kochanek et al. 2001; Bell et al. 2003; Drory et al. 2004; Zhu et al. 2010; Meidt et al. 2014). Using a matched GAMA-WISE catalogue, Cluver et al. (2014) explored the usability of mid-IR wavelengths (W1 and W2) for stellar mass estimation.

However, the relation between K-band luminosity and stellar mass is not entirely direct (van der Wel et al. 2006; Kannappan & Gawiser 2007). In our results, we also observe the differences between stellar mass and K-band MCFs on small scales (see Fig. 4). The close galaxy pairs (rp < 1 h−1 Mpc) show stronger signals when weighted using stellar mass than when weighted with K-band luminosity. Given the fact that the latter is proportional to stellar mass, one possible reason for the difference between both MCFs could be the environmental dependence of the correlation between stellar mass and K-band luminosity. This is in qualitative agreement with the predictions of semi-analytic galaxy formation models proposed by Sheth et al. (2005). Their measurements were made on a volume-limited sample with z = 0.2 with a limiting stellar mass of 2 × 1010h−1M that corresponds to log (M/M) > 10.45 in our cosmological model. Using the observed difference between K-band and stellar mass MCFs (see Fig. 5 of Sheth et al. 2005), they concluded that the correlation between both these properties depends on environment. The difference in K-band and stellar mass MCFs could also be the result of the stronger correlation between stellar mass and halo mass with respect to K-band luminosity (Moster et al. 2010). Hence the stellar mass MCF picks up the environmental dependence of halo mass better than K-band MCF.

There is yet another caveat in using K-band luminosity selected samples as a proxy for stellar mass selection. By using this approximation the most evolved, red galaxies in their sample are missed, especially on small scales. Cochrane et al. (2018) observed that K-band derived stellar masses are underestimated with respect to full SED stellar masses. This means that some red galaxies might wrongly end up below the applied stellar mass cut and not be selected when K-band stellar mass approximation is used. This is extremely important for high-detail clustering studies. It has been shown that red galaxies tend to occupy dense environments (see e.g. Zehavi et al. 2005, 2011; Coil et al. 2008; Palamara et al. 2013). Samples unrepresented owing to the K-band selection might, therefore, show weaker than actual clustering properties.

5.5. Mass incompleteness effect of flux-limited galaxies

Samples ℬ1 and ℬ2 have the same stellar mass limit, but different flux limits in the r band. Namely, sample ℬ1 goes fainter with r < 19.8 and sample ℬ2 is brighter with r < 17.8 (see Table 2 for details). As a result of the imposed flux limit, sample ℬ2 is partially incomplete and significantly less numerous than the complete sample ℬ1; sample ℬ2 lacks galaxies that have mass log (M/M) > 10.4 but are not luminous enough to cross the flux limit of r < 17.8 (Fig. 2).

From the measurements of 2pCFs, we observe that the correlation length does not differ significantly between ℬ1 and ℬ2 (see left panel of Fig. 6). On the other hand, the MCFs are affected by this difference. In the right panel of Fig. 6, we show the ratio of MCFs between samples ℬ2 and ℬ1. It is to be noted that correlation studies based on the u-band luminosity are influenced at small scales by the imposed r-band flux thresholds. A difference between these two samples is visible in the u-band MCF on the smallest scale (rp ∼ 0.1 h−1 Mpc) in the right panel of Fig. 6, where a brighter sample shows a weaker MCF signal. This means that close pairs of galaxies that are similarly brighter in the u band drop out from the sample with the lower magnitude limit. The same behaviour is reflected in the SFR and sSFR MCFs. This is an important observation because the u-band luminosity is a very good tracer of star formation processes. So, even though samples ℬ1 and ℬ2 have the same stellar mass limit, clustering studies on a brighter flux-limited sample (r < 17.8) lose information about starburst galaxies. This effect has to be properly corrected for, for example by using methods discussed in Meneux et al. (2008, 2009).

As mentioned in Sect. 4.2, we observe a lowering of most of the MCF amplitudes while shifting from sample 𝒜1 to ℬ1. Interestingly, we also observe a rise in the amplitudes while shifting from ℬ1 to ℬ2, particularly at the larger scales. It is to be noted that, for the flux limit r < 19.8, sample ℬ1 is more stellar mass complete than 𝒜1. Additionally, it is evident from the stellar mass histograms in Fig. 2 that sample ℬ1 is more stellar mass complete than ℬ2 for the same stellar mass limit log (M/M)min = 10.4.

This suggests that the different behaviour of MCFs we observe in samples 𝒜1 and ℬ2 is due to their stellar mass incompleteness, which in turn comes from the apparent flux limit of the survey and the applied stellar mass limit.

We observe an enhancement in the MCFs of ℬ2 sample relative to ℬ1 at larger scales (Fig. 5 and right panel of Fig. 6). One possible reason could be the influence of the flux limit on the redshift distribution. In Fig. 2, it is clear that the ℬ1 sample is dominated by galaxies at the higher redshift. But the applied flux limit of ℬ2 reduces the number of galaxies at higher redshift by a factor of almost three. This effect could have been propagated to the MCFs causing them to deviate from unity at larger scale.

Mass incompleteness effects can also be observed in MCF measurements based on two differently flux-limited SDSS samples (see Fig. 7). The stellar mass MCF of the more incomplete sample 𝒞5 shows a lower value than 𝒞4 in most of the scales; the sample 𝒞5 value varies from 𝒞4 on average by ΔMp(rp) = 0.07 ± 0.06. These differences occur even though both samples have the same stellar mass limit.

We do not observe this mass incompleteness effect in GAMA samples, where amplitudes of MCFs do not change significantly with the apparent flux limit. However, although within error bars, the 𝒞3 stellar mass MCF shows a consistently lower value than 𝒞1 and 𝒞2. This behaviour is similar to the 𝒞4 and 𝒞5 stellar mass MCFs. This can be associated with the variation in the number of galaxies that enter the flux-limited sample even with the same stellar mass limit. The SDSS sample with the brighter flux limit r < 16.8 (𝒞5) has 11 426 fewer galaxies than that with the fainter limit r < 17.8 (𝒞4). At the same time, the differences between the number of galaxies in the GAMA samples (𝒞1, 𝒞2, and 𝒞3) are very small (see Table 2). Therefore the lack of flux limit dependence of CFs in the GAMA survey can be the result of little variation in the number of galaxies between the GAMA samples 𝒞1, 𝒞2, and 𝒞3.

To study this effect further, we would have to understand the clustering properties of the galaxies missing in brighter flux-limited samples. Such studies were previously done by Meneux et al. (2008, 2009), Marulli et al. (2013); however, these studies were only based on measurements of projected CFs. To understand detailed connections between these galaxies, we would like to measure the behaviour of different MCFs. As for now, the real data GAMA samples for this kind of studies would only consist of 300−400 galaxies (which are missing in the brighter samples). So this could only be done using catalogues built from simulations coupled with semi-analytical galaxy formation models. This kind of study is beyond the scope of the present paper.

6. Summary and conclusions

In this paper, we studied the environmental dependence of galaxy properties such as luminosity (in u, g, r, J, and K bands), stellar mass, SFR, and sSFR in the redshift range 0.1 < z < 0.16 using a spectroscopic sample of galaxies from the GAMA survey. We checked which of these properties is a better tracer of the environment, and showed how the results of clustering measurements can be influenced by selecting samples using different properties. To achieve these aims, we measured the projected correlation function and MCFs in a nearly stellar-mass-complete sample with a flux limit r < 19.8 and stellar mass cut log (M/M) > 9.3. The MCFs were measured using different marks: luminosities in the u, g, r, J, and K bands, stellar mass, SFR, and sSFR. Additionally we studied the dependence of MCF on the survey flux limit, by repeating the same measurements in two samples with different flux limits (r < 19.8 and r < 17.8) but the same stellar mass cut (log (M/M) > 10.4). We also did measurements in samples with GAMA galaxies having mass log (M/M) > 10.8 and compared with SDSS galaxies with the same stellar mass cut, but different flux limits.

The summary of our main results and conclusions of our study are as follows:

  • We observed that different galaxy properties trace the environment differently in the separation scales rp < 10 h−1 Mpc. Based on the behaviour of MCFs in Fig. 4, we concluded that the close pairs of galaxies are more luminous in the g, r, J, and K bands than distant pairs. It is also more probable to find close pairs of massive galaxies (with masses above sample average) than pairs including one less massive galaxy. However, this trend is reversed if the luminosity is measured in the u band. The u-band luminous galaxies tend to occupy less dense regions and the faint u-band galaxies tend to exist in denser regions. The same is true for actively star-forming galaxies, which tend to occupy less dense environments.

  • From the comparisons of the amplitudes of rank-ordered luminosity, stellar mass, SFR, and sSFR MCFs in Fig. 4, we concluded that stellar mass is more reliable to trace galaxy environment than the other properties.

  • We showed that a sample complete in K-band luminosity can be a good substitute for a stellar-mass-complete sample. But in such a case we tend to miss closer pairs of evolved, red galaxies.

  • From the similarity in behaviour of different MCFs, we suggested the usefulness of u-band luminosity as a proxy of SFR in the context of galaxy clustering.

  • From the comparative study of MCFs in different samples with the same mass selection, but different apparent magnitude limits, we concluded that closer pairs of star-forming galaxies drop out of the sample when the survey gets shallower in terms of limiting magnitude.

Our measurements are the first of this kind in the redshift range 0.1 < z < 0.16 and with galaxies as faint as r < 19.8. These measurements can be a useful reference for clustering studies with flux-limited surveys, specially the next generation galaxy surveys such as Vera C. Rubin Observatory (LSST Science Collaboration 2009) and Euclid (Laureijs et al. 2011). We intend to extend our measurements using galaxy catalogues from simulations and high-redshift surveys, with which we hope to provide better constraints on models of galaxy formation and evolution.


Acknowledgments

We thank the anonymous referee for the useful comments and suggestions. U.S. and A.D. are supported by the Polish National Science Centre grant UMO-2015/17/D/ST9/02121. U.S. is supported by Jagiellonian University DSC grant 2019-N17/MNS/000045. U.S. and A.P. are supported by the Polish National Science Centre grant UMO-2018/30/M/ST9/00757. MB is supported by the Polish National Science Center through grants no. 2020/38/E/ST9/00395, 2018/30/E/ST9/00698, and 2018/31/G/ST9/03388. This work is supported by Polish Ministry of Science and Higher Education grant DIR/WK/2018/12. A.H.W. is supported by an European Research Council Consolidator Grant (No. 770935). GAMA is a joint European-Australasian project based around a spectroscopic campaign using the Anglo-Australian Telescope. The GAMA input catalogue is based on data taken from the Sloan Digital Sky Survey and the UKIRT Infrared Deep Sky Survey. Complementary imaging of the GAMA regions is being obtained by a number of independent survey programmes including GALEX MIS, VST KiDS, VISTA VIKING, WISE, Herschel-ATLAS, GMRT, and ASKAP providing UV to radio coverage. GAMA is funded by the STFC (UK), the ARC (Australia), the AAO, and the participating institutions. The GAMA website is http://www.gama-survey.org/. Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 179.A-2004. During this research, we made use of Tool for OPerations on Catalogues And Tables (TOPCAT; Taylor 2005) and NASA’s Astrophysics Data System Bibliographic Services. This research was supported in part by PLGrid Infrastructure and OAUJ cluster computing facility.

References

  1. Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543 [Google Scholar]
  2. Abbas, U., & Sheth, R. K. 2006, MNRAS, 372, 1749 [NASA ADS] [CrossRef] [Google Scholar]
  3. Abbas, U., de la Torre, S., Le Fèvre, O., et al. 2010, MNRAS, 406, 1306 [NASA ADS] [Google Scholar]
  4. Alam, S., Albareti, F. D., Allende Prieto, C., et al. 2015, ApJS, 219, 12 [NASA ADS] [CrossRef] [Google Scholar]
  5. Alam, S., Peacock, J. A., Farrow, D. J., Loveday, J., & Hopkins, A. M. 2021, MNRAS, 503, 59 [NASA ADS] [CrossRef] [Google Scholar]
  6. Alpaslan, M., Robotham, A. S. G., Driver, S., et al. 2014, MNRAS, 438, 177 [NASA ADS] [CrossRef] [Google Scholar]
  7. Alpaslan, M., Driver, S., Robotham, A. S. G., et al. 2015, MNRAS, 451, 3249 [NASA ADS] [CrossRef] [Google Scholar]
  8. Baldry, I. K., Robotham, A. S. G., Hill, D. T., et al. 2010, MNRAS, 404, 86 [NASA ADS] [Google Scholar]
  9. Baldry, I. K., Driver, S. P., Loveday, J., et al. 2012, MNRAS, 421, 621 [NASA ADS] [Google Scholar]
  10. Baldry, I. K., Alpaslan, M., Bauer, A. E., et al. 2014, MNRAS, 441, 2440 [NASA ADS] [CrossRef] [Google Scholar]
  11. Barsanti, S., Owers, M. S., Brough, S., et al. 2018, ApJ, 857, 71 [NASA ADS] [CrossRef] [Google Scholar]
  12. Baugh, C. M. 2006, Rep. Progr. Phys., 69, 3101 [Google Scholar]
  13. Beck, R., Dobos, L., Budavári, T., Szalay, A. S., & Csabai, I. 2016, MNRAS, 460, 1371 [NASA ADS] [CrossRef] [Google Scholar]
  14. Behroozi, P., Wechsler, R. H., Hearin, A. P., & Conroy, C. 2019, MNRAS, 488, 3143 [NASA ADS] [CrossRef] [Google Scholar]
  15. Beisbart, C., & Kerscher, M. 2000, ApJ, 545, 6 [NASA ADS] [CrossRef] [Google Scholar]
  16. Bell, E. F., McIntosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289 [Google Scholar]
  17. Benson, A. J. 2012, New Astron., 17, 175 [NASA ADS] [CrossRef] [Google Scholar]
  18. Beutler, F., Blake, C., Colless, M., et al. 2013, MNRAS, 429, 3604 [NASA ADS] [CrossRef] [Google Scholar]
  19. Blanton, M. R., & Moustakas, J. 2009, ARA&A, 47, 159 [NASA ADS] [CrossRef] [Google Scholar]
  20. Blanton, M. R., Eisenstein, D., Hogg, D. W., Schlegel, D. J., & Brinkmann, J. 2005, ApJ, 629, 143 [NASA ADS] [CrossRef] [Google Scholar]
  21. Bolzonella, M., Miralles, J. M., & Pelló, R. 2000, A&A, 363, 476 [NASA ADS] [Google Scholar]
  22. Bolzonella, M., Kovač, K., Pozzetti, L., et al. 2010, A&A, 524, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Bond, J. R., Kofman, L., & Pogosyan, D. 1996, Nature, 380, 603 [NASA ADS] [CrossRef] [Google Scholar]
  24. Bravo, M., Lagos, C. D. P., Robotham, A. S. G., Bellstedt, S., & Obreschkow, D. 2020, MNRAS, 497, 3026 [NASA ADS] [CrossRef] [Google Scholar]
  25. Brough, S., Croom, S., Sharp, R., et al. 2013, MNRAS, 435, 2903 [NASA ADS] [CrossRef] [Google Scholar]
  26. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  27. Burton, C. S., Jarvis, M. J., Smith, D. J. B., et al. 2013, MNRAS, 433, 771 [NASA ADS] [CrossRef] [Google Scholar]
  28. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  29. Chabrier, G. 2003, ApJ, 586, L133 [NASA ADS] [CrossRef] [Google Scholar]
  30. Christodoulou, L., Eminian, C., Loveday, J., et al. 2012, MNRAS, 425, 1527 [NASA ADS] [CrossRef] [Google Scholar]
  31. Cluver, M. E., Jarrett, T. H., Hopkins, A. M., et al. 2014, ApJ, 782, 90 [NASA ADS] [CrossRef] [Google Scholar]
  32. Cochrane, R. K., Best, P. N., Sobral, D., et al. 2018, MNRAS, 475, 3730 [NASA ADS] [CrossRef] [Google Scholar]
  33. Cole, S., Norberg, P., Baugh, C. M., et al. 2001, MNRAS, 326, 255 [NASA ADS] [CrossRef] [Google Scholar]
  34. Coil, A. L., Newman, J. A., Cooper, M. C., et al. 2006, ApJ, 644, 671 [NASA ADS] [CrossRef] [Google Scholar]
  35. Coil, A. L., Newman, J. A., Croton, D., et al. 2008, ApJ, 672, 153 [NASA ADS] [CrossRef] [Google Scholar]
  36. Condon, J. J. 1992, ARA&A, 30, 575 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  37. Cowie, L. L., Gardner, J. P., Hu, E. M., et al. 1994, ApJ, 434, 114 [NASA ADS] [CrossRef] [Google Scholar]
  38. Cram, L., Hopkins, A., Mobasher, B., & Rowan-Robinson, M. 1998, ApJ, 507, 155 [NASA ADS] [CrossRef] [Google Scholar]
  39. Croton, D. J., Gao, L., & White, S. D. M. 2007, MNRAS, 374, 1303 [Google Scholar]
  40. da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 [Google Scholar]
  41. Davidzon, I., Cucciati, O., Bolzonella, M., et al. 2016, A&A, 586, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Davis, M., & Peebles, P. J. E. 1983, ApJ, 267, 465 [NASA ADS] [CrossRef] [Google Scholar]
  43. Davies, L. J. M., Robotham, A. S. G., Driver, S. P., et al. 2016, MNRAS, 455, 4013 [NASA ADS] [CrossRef] [Google Scholar]
  44. Davies, L. J. M., Lagos, C. D. P., Katsianis, A., et al. 2019a, MNRAS, 483, 1881 [Google Scholar]
  45. Davies, L. J. M., Robotham, A. S. G., Lagos, C. D. P., et al. 2019b, MNRAS, 483, 5444 [Google Scholar]
  46. Dawson, K. S., Schlegel, D. J., Ahn, C. P., et al. 2013, AJ, 145, 10 [Google Scholar]
  47. de Lapparent, V., Geller, M. J., & Huchra, J. P. 1986, ApJ, 302, L1 [Google Scholar]
  48. de la Torre, S., Le Fèvre, O., Arnouts, S., et al. 2007, A&A, 475, 443 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Deng, X.-F. 2012, AJ, 143, 15 [CrossRef] [Google Scholar]
  50. Dolley, T., Brown, M. J. I., Weiner, B. J., et al. 2014, ApJ, 797, 125 [NASA ADS] [CrossRef] [Google Scholar]
  51. Driver, S. P., Norberg, P., Baldry, I. K., et al. 2009, Astron. Geophys., 50, 5.12 [NASA ADS] [CrossRef] [Google Scholar]
  52. Driver, S. P., Hill, D. T., Kelvin, L. S., et al. 2011, MNRAS, 413, 971 [Google Scholar]
  53. Drory, N., Bender, R., Feulner, G., et al. 2004, ApJ, 608, 742 [NASA ADS] [CrossRef] [Google Scholar]
  54. Durkalec, A., Le Fèvre, O., Pollo, A., et al. 2018, A&A, 612, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Eales, S., Dunne, L., Clements, D., et al. 2010, PASP, 122, 499 [NASA ADS] [CrossRef] [Google Scholar]
  56. Farrow, D. J., Cole, S., Norberg, P., et al. 2015, MNRAS, 454, 2120 [NASA ADS] [CrossRef] [Google Scholar]
  57. Fisher, K. B., Davis, M., Strauss, M. A., Yahil, A., & Huchra, J. 1994, MNRAS, 266, 50 [NASA ADS] [CrossRef] [Google Scholar]
  58. Gavazzi, G., Pierini, D., & Boselli, A. 1996, A&A, 312, 397 [NASA ADS] [Google Scholar]
  59. Gaztañaga, E., & Scoccimarro, R. 2005, MNRAS, 361, 824 [NASA ADS] [CrossRef] [Google Scholar]
  60. Gómez, P. L., Nichol, R. C., Miller, C. J., et al. 2003, ApJ, 584, 210 [NASA ADS] [CrossRef] [Google Scholar]
  61. Grootes, M. W., Tuffs, R. J., Popescu, C. C., et al. 2017, AJ, 153, 111 [NASA ADS] [CrossRef] [Google Scholar]
  62. Groth, E. J., & Peebles, P. J. E. 1977, ApJ, 217, 385 [NASA ADS] [CrossRef] [Google Scholar]
  63. Grylls, P. J., Shankar, F., Leja, J., et al. 2020, MNRAS, 491, 634 [NASA ADS] [Google Scholar]
  64. Gu, M., Conroy, C., & Behroozi, P. 2016, ApJ, 833, 2 [NASA ADS] [CrossRef] [Google Scholar]
  65. Gunawardhana, M. L. P., Norberg, P., Zehavi, I., et al. 2018, MNRAS, 479, 1433 [NASA ADS] [CrossRef] [Google Scholar]
  66. Guo, H., Zheng, Z., Zehavi, I., et al. 2015, MNRAS, 453, 4368 [NASA ADS] [CrossRef] [Google Scholar]
  67. Hamilton, A. J. S. 1993, ApJ, 417, 19 [NASA ADS] [CrossRef] [Google Scholar]
  68. Hartley, W. G., Almaini, O., Cirasuolo, M., et al. 2010, MNRAS, 407, 1212 [NASA ADS] [CrossRef] [Google Scholar]
  69. Heavens, A., Panter, B., Jimenez, R., & Dunlop, J. 2004, Nature, 428, 625 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  70. Heinis, S., Treyer, M., Arnouts, S., et al. 2004, A&A, 424, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  71. Heinis, S., Milliard, B., Arnouts, S., et al. 2007, ApJS, 173, 503 [NASA ADS] [CrossRef] [Google Scholar]
  72. Hopkins, A. M., Miller, C. J., Nichol, R. C., et al. 2003, ApJ, 599, 971 [NASA ADS] [CrossRef] [Google Scholar]
  73. Jarrett, T. H., Masci, F., Tsai, C. W., et al. 2013, AJ, 145, 6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Jarrett, T. H., Cluver, M. E., Magoulas, C., et al. 2017, ApJ, 836, 182 [Google Scholar]
  75. Kannan, R., Stinson, G. S., Macciò, A. V., et al. 2014, MNRAS, 437, 3529 [NASA ADS] [CrossRef] [Google Scholar]
  76. Kannappan, S. J., & Gawiser, E. 2007, ApJ, 657, L5 [NASA ADS] [CrossRef] [Google Scholar]
  77. Kauffmann, G., & Charlot, S. 1998, MNRAS, 297, L23 [NASA ADS] [CrossRef] [Google Scholar]
  78. Kauffmann, G., Heckman, T. M., White, S. D. M., et al. 2003, MNRAS, 341, 54 [Google Scholar]
  79. Kauffmann, G., White, S. D. M., Heckman, T. M., et al. 2004, MNRAS, 353, 713 [NASA ADS] [CrossRef] [Google Scholar]
  80. Kennicutt, R. C., Jr. 1998, ARA&A, 36, 189 [Google Scholar]
  81. Kochanek, C. S., Pahre, M. A., Falco, E. E., et al. 2001, ApJ, 560, 566 [NASA ADS] [CrossRef] [Google Scholar]
  82. Lagos, C. D. P., Tobar, R. J., Robotham, A. S. G., et al. 2018, MNRAS, 481, 3573 [CrossRef] [Google Scholar]
  83. Landy, S. D., & Szalay, A. S. 1993, ApJ, 412, 64 [Google Scholar]
  84. Laureijs, R., Amiaux, J., Arduini, S., et al. 2011, ArXiv e-prints [arXiv:1110.3193] [Google Scholar]
  85. Lewis, I., Balogh, M., De Propris, R., et al. 2002, MNRAS, 334, 673 [NASA ADS] [CrossRef] [Google Scholar]
  86. Lindsay, S. N., Jarvis, M. J., Santos, M. G., et al. 2014, MNRAS, 440, 1527 [NASA ADS] [CrossRef] [Google Scholar]
  87. Linke, L., Simon, P., Schneider, P., et al. 2020, A&A, 640, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Liske, J., Baldry, I. K., Driver, S. P., et al. 2015, MNRAS, 452, 2087 [NASA ADS] [CrossRef] [Google Scholar]
  89. Loveday, J., Christodoulou, L., Norberg, P., et al. 2018, MNRAS, 474, 3435 [NASA ADS] [CrossRef] [Google Scholar]
  90. LSST Science Collaboration (Abell, P. A., et al.) 2009, ArXiv e-prints [arXiv:0912.0201] [Google Scholar]
  91. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [NASA ADS] [CrossRef] [Google Scholar]
  92. Mao, Y.-Y., Zentner, A. R., & Wechsler, R. H. 2018, MNRAS, 474, 5143 [NASA ADS] [CrossRef] [Google Scholar]
  93. Maraston, C., Strömbäck, G., Thomas, D., Wake, D. A., & Nichol, R. C. 2009, MNRAS, 394, L107 [NASA ADS] [CrossRef] [Google Scholar]
  94. Marín, F. A., Blake, C., Poole, G. B., et al. 2013, MNRAS, 432, 2654 [NASA ADS] [CrossRef] [Google Scholar]
  95. Marulli, F., Bolzonella, M., Branchini, E., et al. 2013, A&A, 557, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. McCarthy, I. G., Font, A. S., Crain, R. A., et al. 2012, MNRAS, 420, 2245 [NASA ADS] [CrossRef] [Google Scholar]
  97. McGee, S. L., Balogh, M. L., Wilman, D. J., et al. 2011, MNRAS, 413, 996 [NASA ADS] [CrossRef] [Google Scholar]
  98. McNaught-Roberts, T., Norberg, P., Baugh, C., et al. 2014, MNRAS, 445, 2125 [NASA ADS] [CrossRef] [Google Scholar]
  99. Meidt, S. E., Schinnerer, E., van de Ven, G., et al. 2014, ApJ, 788, 144 [NASA ADS] [CrossRef] [Google Scholar]
  100. Meneux, B., Le Fèvre, O., Guzzo, L., et al. 2006, A&A, 452, 387 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Meneux, B., Guzzo, L., Garilli, B., et al. 2008, A&A, 478, 299 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  102. Meneux, B., Guzzo, L., de la Torre, S., et al. 2009, A&A, 505, 463 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Milliard, B., Heinis, S., Blaizot, J., et al. 2007, ApJS, 173, 494 [NASA ADS] [CrossRef] [Google Scholar]
  104. Mo, H. J., & White, S. D. M. 1996, MNRAS, 282, 347 [Google Scholar]
  105. More, S., van den Bosch, F. C., Cacciato, M., et al. 2009, MNRAS, 392, 801 [NASA ADS] [CrossRef] [Google Scholar]
  106. Mostek, N., Coil, A. L., Cooper, M., et al. 2013, ApJ, 767, 89 [NASA ADS] [CrossRef] [Google Scholar]
  107. Moster, B. P., Somerville, R. S., Maulbetsch, C., et al. 2010, ApJ, 710, 903 [Google Scholar]
  108. Moster, B. P., Naab, T., & White, S. D. M. 2020, MNRAS, 499, 4748 [Google Scholar]
  109. Muldrew, S. I., Croton, D. J., Skibba, R. A., et al. 2012, MNRAS, 419, 2670 [NASA ADS] [CrossRef] [Google Scholar]
  110. Naab, T., & Ostriker, J. P. 2017, ARA&A, 55, 59 [Google Scholar]
  111. Norberg, P., Baugh, C. M., Hawkins, E., et al. 2001, MNRAS, 328, 64 [NASA ADS] [CrossRef] [Google Scholar]
  112. Norberg, P., Baugh, C. M., Hawkins, E., et al. 2002, MNRAS, 332, 827 [NASA ADS] [CrossRef] [Google Scholar]
  113. Norberg, P., Baugh, C. M., Gaztañaga, E., & Croton, D. J. 2009, MNRAS, 396, 19 [NASA ADS] [CrossRef] [Google Scholar]
  114. Oliver, S., Waddington, I., Gonzalez-Solares, E., et al. 2004, ApJS, 154, 30 [NASA ADS] [CrossRef] [Google Scholar]
  115. Palamara, D. P., Brown, M. J. I., Jannuzi, B. T., et al. 2013, ApJ, 764, 31 [NASA ADS] [CrossRef] [Google Scholar]
  116. Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe (Princeton: Princeton University Press) [Google Scholar]
  117. Peng, Y.-J., Lilly, S. J., Kovač, K., et al. 2010, ApJ, 721, 193 [NASA ADS] [CrossRef] [Google Scholar]
  118. Petrosian, V. 1976, ApJ, 210, L53 [NASA ADS] [CrossRef] [Google Scholar]
  119. Pollo, A., Meneux, B., Guzzo, L., et al. 2005, A&A, 439, 887 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  120. Pollo, A., Guzzo, L., Le Fèvre, O., et al. 2006, A&A, 451, 409 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  121. Pollo, A., Takeuchi, T. T., Solarz, A., et al. 2013a, Earth Planets Space, 65, 1109 [NASA ADS] [CrossRef] [Google Scholar]
  122. Pollo, A., Takeuchi, T. T., Suzuki, T. L., & Oyabu, S. 2013b, Earth Planets Space, 65, 273 [NASA ADS] [CrossRef] [Google Scholar]
  123. Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [NASA ADS] [CrossRef] [Google Scholar]
  124. Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in FORTRAN. The Art of Scientific Computing (Cambridge: Cambridge University Press) [Google Scholar]
  125. Reid, B., Ho, S., Padmanabhan, N., et al. 2016, MNRAS, 455, 1553 [NASA ADS] [CrossRef] [Google Scholar]
  126. Riggs, S. D., Barbhuiyan, R. W. Y. M., Loveday, J., et al. 2021, MNRAS, 506, 21 [NASA ADS] [CrossRef] [Google Scholar]
  127. Robotham, A., Driver, S. P., Norberg, P., et al. 2010, PASA, 27, 76 [NASA ADS] [CrossRef] [Google Scholar]
  128. Schaefer, A. L., Croom, S. M., Allen, J. T., et al. 2017, MNRAS, 464, 121 [NASA ADS] [CrossRef] [Google Scholar]
  129. Schaefer, A. L., Croom, S. M., Scott, N., et al. 2019, MNRAS, 483, 2851 [NASA ADS] [CrossRef] [Google Scholar]
  130. Scodeggio, M., Vergani, D., Cucciati, O., et al. 2009, A&A, 501, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  131. Sheth, R. K. 2005, MNRAS, 364, 796 [NASA ADS] [CrossRef] [Google Scholar]
  132. Sheth, R. K., & Tormen, G. 2004, MNRAS, 350, 1385 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  133. Sheth, R. K., Connolly, A. J., & Skibba, R. 2005, ArXiv e-prints [arXiv:astro-ph/0511773] [Google Scholar]
  134. Sheth, R. K., Jimenez, R., Panter, B., & Heavens, A. F. 2006, ApJ, 650, L25 [NASA ADS] [CrossRef] [Google Scholar]
  135. Skibba, R., Sheth, R. K., Connolly, A. J., & Scranton, R. 2006, MNRAS, 369, 68 [NASA ADS] [CrossRef] [Google Scholar]
  136. Skibba, R. A., Sheth, R. K., Croton, D. J., et al. 2013, MNRAS, 429, 458 [NASA ADS] [CrossRef] [Google Scholar]
  137. Skibba, R. A., Smith, M. S. M., Coil, A. L., et al. 2014, ApJ, 784, 128 [NASA ADS] [CrossRef] [Google Scholar]
  138. Skibba, R. A., Coil, A. L., Mendez, A. J., et al. 2015, ApJ, 807, 152 [NASA ADS] [CrossRef] [Google Scholar]
  139. Sobral, D., Best, P. N., Geach, J. E., et al. 2010, MNRAS, 404, 1551 [NASA ADS] [Google Scholar]
  140. Somerville, R. S., & Davé, R. 2015, ARA&A, 53, 51 [Google Scholar]
  141. Springel, V., White, S. D. M., Jenkins, A., et al. 2005, Nature, 435, 629 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
  142. Stoyan, D., & Stoyan, D. H. 1994, Fractals, Random Shapes and Point Fields: Methods of Geometrical Statistics (New York: Wiley) [Google Scholar]
  143. Sureshkumar, U., Durkalec, A., Pollo, A., Bilicki, M., & GAMA Collaboration 2020, in XXXIX Polish Astronomical Society Meeting, eds. K. Małek, M. Polińska, A. Majczyna, et al., 10, 346 [NASA ADS] [Google Scholar]
  144. Taylor, M. B. 2005, in Astronomical Data Analysis Software and Systems XIV, eds. P. Shopbell, M. Britton, & R. Ebert, ASP Conf. Ser., 347, 29 [Google Scholar]
  145. Taylor, E. N., Franx, M., van Dokkum, P. G., et al. 2009, ApJ, 694, 1171 [NASA ADS] [CrossRef] [Google Scholar]
  146. Taylor, E. N., Hopkins, A. M., Baldry, I. K., et al. 2011, MNRAS, 418, 1587 [NASA ADS] [CrossRef] [Google Scholar]
  147. Tonry, J. L., Blakeslee, J. P., Ajhar, E. A., & Dressler, A. 2000, ApJ, 530, 625 [NASA ADS] [CrossRef] [Google Scholar]
  148. van der Wel, A., Franx, M., Wuyts, S., et al. 2006, ApJ, 652, 97 [NASA ADS] [CrossRef] [Google Scholar]
  149. van Dokkum, P. G., Quadri, R., Marchesini, D., et al. 2006, ApJ, 638, L59 [NASA ADS] [CrossRef] [Google Scholar]
  150. van Kampen, E., Smith, D. J. B., Maddox, S., et al. 2012, MNRAS, 426, 3455 [NASA ADS] [CrossRef] [Google Scholar]
  151. van Uitert, E., Joachimi, B., Joudaki, S., et al. 2018, MNRAS, 476, 4662 [NASA ADS] [CrossRef] [Google Scholar]
  152. Vázquez-Mata, J. A., Loveday, J., Riggs, S. D., et al. 2020, MNRAS, 499, 631 [Google Scholar]
  153. Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, MNRAS, 436, 3031 [NASA ADS] [CrossRef] [Google Scholar]
  154. Wang, L., Norberg, P., Brough, S., et al. 2018, A&A, 618, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  155. Wechsler, R. H., & Tinker, J. L. 2018, ARA&A, 56, 435 [NASA ADS] [CrossRef] [Google Scholar]
  156. White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 [NASA ADS] [CrossRef] [Google Scholar]
  157. Wijesinghe, D. B., Hopkins, A. M., Brough, S., et al. 2012, MNRAS, 423, 3679 [NASA ADS] [CrossRef] [Google Scholar]
  158. Wright, A. H., Robotham, A. S. G., Bourne, N., et al. 2016, MNRAS, 460, 765 [NASA ADS] [CrossRef] [Google Scholar]
  159. Yang, X., Mo, H. J., van den Bosch, F. C., Zhang, Y., & Han, J. 2012, ApJ, 752, 41 [NASA ADS] [CrossRef] [Google Scholar]
  160. York, D. G., Adelman, J., Anderson, J. E., Jr., et al. 2000, AJ, 120, 1579 [CrossRef] [Google Scholar]
  161. Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2005, ApJ, 630, 1 [Google Scholar]
  162. Zehavi, I., Zheng, Z., Weinberg, D. H., et al. 2011, ApJ, 736, 59 [NASA ADS] [CrossRef] [Google Scholar]
  163. Zentner, A. R., Hearin, A. P., & van den Bosch, F. C. 2014, MNRAS, 443, 3044 [NASA ADS] [CrossRef] [Google Scholar]
  164. Zheng, Z., Berlind, A. A., Weinberg, D. H., et al. 2005, ApJ, 633, 791 [NASA ADS] [CrossRef] [Google Scholar]
  165. Zhu, Y.-N., Wu, H., Li, H.-N., & Cao, C. 2010, Res. Astron. Astrophys., 10, 329 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Properties of the GAMA equatorial regions used in this work.

Table 2.

Definitions of the galaxy samples in the redshift range 0.1 < z < 0.16, as used in this study.

Table 3.

Properties of the galaxy samples defined in Table 2.

Table 4.

Best-fitting power-law parameters for all the galaxy samples used in this work.

All Figures

thumbnail Fig. 1.

Selection of galaxy sample 𝒜1 used in this work. The grey dots represent the GAMA galaxies with flux limit r < 19.8. Top and right histograms: distribution of redshift and stellar mass, respectively. The red lines represent the stellar mass cut and the redshift limit of the sample 𝒜1.

In the text
thumbnail Fig. 2.

Redshift and stellar mass distribution of galaxies in stellar-mass-selected samples ℬ1 (r < 19.8; brown dots) and ℬ2 (r < 17.8; green circles). Top and right histograms: distribution of redshift and stellar mass, respectively.

In the text
thumbnail Fig. 3.

Redshift and stellar mass distribution of stellar-mass-selected samples mentioned in Table 2. Left panel: 𝒞1 (r < 19.8; black dots), 𝒞2 (r < 18.8; red circles), and 𝒞3 (r < 17.8; green squares) galaxy samples from GAMA and the right panel shows the 𝒞4 (r < 17.8; blue dots) and 𝒞5 (r < 16.8; orange squares) galaxy samples from SDSS. Top and right histograms of both the panels show the distribution of redshift and stellar mass, respectively.

In the text
thumbnail Fig. 4.

Projected 2pCF ωp(rp) with a power-law fit (filled markers; left panel) and rank-ordered projected MCFs Mp(rp) (unfilled markers; right panel) for sample 𝒜1 described in Sect. 2.2. In the left panel, error bars of ωp(rp) are the square root of the diagonals of the covariance matrix obtained from the jackknife resampling method. The inset in the left panel shows the power-law fit parameters (filled star) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). In the right panel the different symbols represent measurements with different marks (as labelled), and the error bars are obtained by random scrambling of the marks. The error bars of Mp(rp) are too small to be visible.

In the text
thumbnail Fig. 5.

Rank-ordered projected MCFs Mp(rp) for samples ℬ1 (left) and ℬ2 (right) with different flux limits, as labelled. The different symbols represent different marks considered for the Mp(rp) measurement (as labelled) and the error bars are obtained by random scrambling of the marks.

In the text
thumbnail Fig. 6.

Projected 2pCFs ωp(rp) for samples ℬ1 and ℬ2 (filled markers; left panel) and the ratio Mℬ2/Mℬ1 between MCFs obtained for these two samples (unfilled markers; right panel). The different symbols in the right panel represent the ratio between different MCFs measured with corresponding galaxy property chosen as a mark. Small offsets along x-axis have been added for clarity. The error bars for ωp(rp) are obtained from jackknife resampling method. The inset in the left panel shows the power-law fit parameters (filled stars) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). In the right panel, the errors are calculated in quadrature.

In the text
thumbnail Fig. 7.

Projected 2pCFs (filled markers; left panel) and stellar mass MCFs (unfilled markers; right panel) in GAMA and SDSS surveys. The different symbols represent the measurements in different samples as labelled. The error bars for ωp(rp) are obtained from jackknife resampling method and those for Mp(rp) are obtained by random scrambling of the marks. The inset in the left panel shows the power-law fit parameters (filled stars) and their 1σ, 2σ, and 3σ error contours (solid, dashed, and dotted, respectively). Small offsets along x-axis have been added for clarity.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.