Non-random Co-occurrence of Juvenile White Sharks (Carcharodon carcharias) at Seasonal Aggregation Sites in Southern California

Many terrestrial and aquatic taxa are known to form periodic aggregations, whether across life history or solely during specific life stages, that are generally governed by the availability and distribution of resources. Associations between individuals during such aggregation events are considered random and not driven by social attraction or underlying community structure. White sharks (Carcharodon carcharias) have been described as a species that exhibits resource-driven aggregative behaviors across ontogenetic stages and juvenile white sharks are known to form aggregations at specific nursery sites where individuals may remain for extended periods of time in the presence of other individuals. We hypothesized juvenile white sharks form distinct communities during these critical early phases of ontogeny and discuss how a tendency to co-occur across life stages may be seeded by the formation of these communities in early ontogeny. We present results from a series of social network analyses of 86 juvenile white sharks derived from 6 years of passive acoustic telemetry data in southern California, demonstrating the likelihood of association of tagged juvenile white sharks is greater when sharks are of similar size-classes. Individuals in observed networks exhibited behaviors that best approximated fission-fusion dynamics with spatiotemporally unstable group membership. These results provide evidence of possible non-resource driven co-occurrence and community structure in juvenile white sharks during early life stages.


INTRODUCTION
Sociality, or the propensity of individuals to form social groups (Merriam-Webster., 2021), is a trait exhibited across taxa and varies between loose temporary aggregations to life-long associations (Alexander, 1974;Sabol et al., 2020). Social interactions have been shown to reduce predation risk and increase foraging efficiency, navigational capability, and reproductive opportunities Berdahl et al., 2018;Campbell et al., 2018;Diaz-Aguirre et al., 2019). Thus, animal social groups usually form when the associated benefits of group behavior outweigh the costs, while the influence of such associated costs upon the decision to engage in social behaviors is driven by extrinsic (e.g., prey availability) and intrinsic (e.g., competition) variability (Krause and Ruxton, 2002;Silk et al., 2014). Animal social groups are thus characterized by motivated cohesion and can be distinguished from aggregative behaviors, which are temporary assemblages of individuals in response to non-social forcing factors, such as seasonal resource availability (Jacoby et al., 2011;Meese and Lowe, 2019;Grueter et al., 2020).
White sharks (Carcharodon carcharias) are circumglobally distributed in both temperate and tropical waters (Compagno, 2002). In the sub-adult to adult ontogenetic stages, they are known to occur as solitary individuals but are also found in seasonal aggregations at feeding grounds such as pinniped haul outs (Bruce et al., 2006;Robbins, 2007;Domeier and Nasby-Lucas, 2008;Jorgensen et al., 2010;Duffy et al., 2012;Kock et al., 2013Kock et al., , 2018Kanive et al., 2021). To date, the extent to which socially motivated factors play a role in co-occurrences of individual white sharks remains unresolved, with only two published studies on the subject arriving at different conclusions (Findlay et al., 2016;Schilds et al., 2019). Specifically, Findlay et al. (2016) reported that associations between white sharks in their study were random, although they exhibited weak structuring by sex and body size, whereas Schilds et al. (2019) reported nonrandom, sex-dependent associations with temporal variability. These studies drew upon observations of sub-adult and adult sharks co-occurring under baited and chummed (burleyed) conditions, at locations proximal to pinniped haul-outs, with very different time thresholds constituting co-occurrence. Juvenile white sharks use nearshore and beach habitat, where they have been observed to exhibit high degrees of residency and aggregation site fidelity (Weng et al., 2007;Bruce and Bradford, 2008;Werry et al., 2012;Dicken and Booth, 2013;Harasti et al., 2017;Oñate-González et al., 2017;Curtis et al., 2018;Tamburin et al., 2019;White et al., 2019;Spaet J. L. Y. et al., 2020;Anderson et al., 2021a,b). Such site fidelity is common among elasmobranchs in early ontogeny, as these habitats typically provide opportunities for optimal growth and foraging, and reduced predation risk (Simpfendorfer and Milward, 1993;Heupel et al., 2007;Chapman et al., 2009). Preliminary high-spatial-resolution movement data derived from hyperbolic positioning via a high density acoustic array off Carpinteria, California, indicate individual juvenile white sharks form aggregations (5-30 individuals of multiple size-classes) and use overlapping, spatially restricted areas (∼8 km 2 ) continuously, for periods of weeks to months (Spurgeon et al., unpublished data; Anderson et al., 2021b).
It is unknown whether co-occurrence of individual sharks is random, driven by resource availability and/or environmental factors, or is a function of socially mediated behaviors. Cooccurrence may be a function of overlap in home ranges and exploited niches, or may be due to preferred associations between individuals of phenotypic characteristics. In network theory, this propensity for individuals with similar traits to co-occur is termed assortativity (Newman, 2002;Noldus and Van Mieghem, 2015). Co-occurrence of individual sharks may afford the opportunity for the transfer of information and the development of association preferences. Thus, seasonal resourcedriven aggregations may seed the development of social groups and assortative interactions (Jacoby et al., 2011).
In this study, we applied a network analysis to passive acoustic telemetry data gathered from tagged juvenile white sharks [ < 150-350 cm Total Length (TL)] in southern California over a 6-year period (2014-2019). This approach was designed to gather insight as to whether tagged shark co-occurrence at acoustically monitored sites was random, an artifact of resource-driven aggregation behavior, or may be driven by social preferences. We aimed to examine (1) whether apparent structure could be identified within generated networks, (2) the extent to which tagged sharks formed associations with other tagged individuals, (3) the extent of association stability across time, and (4) to characterize assortativity in association preferences.

Incidental Catch
Commercial gillnet fishers, working in collaboration with CSULB and CICESE researchers, brought incidentally captured juvenile white sharks to the nearest port in a large fish tote (1.2 × 1.2 × 1.2 m) with flowing seawater. Incidentallycaught sharks were caught offshore (>3 nm) outside of State waters. Researchers physically assessed, measured, and surgically implanted a plasma sterilized V16 transmitter into the abdominal cavity of the shark through a small incision (5 cm). The wound was closed with 2-3 interrupted sutures, before the animal was released approximately 2 km offshore. All individuals were visually monitored during release and behavior was recorded. A previous study demonstrated post-release survival rates of juvenile white sharks retrieved live from gillnets to be in the order of 93% (Lyons et al., 2013), while a more recent analysis indicated that minimizing handling and on-deck holding times may in turn minimize non-lethal post-release negative effects (Raoult et al., 2019). Approximate elapsed times from researchers taking possession of the shark to offshore release were < 30 min.

Targeted Catch
Juvenile white sharks were caught in shallow inshore waters either via a purse-seine, or via targeted quick-extraction with a gillnet. Sharks were transferred to a custom-made staging tank and ventilated, before being sexed, sized, and outfitted with a surgically implanted V16 transmitter as described above. Approximate handling times in all instances were < 20 min.

Dart Tagging
Juvenile white sharks swimming close to, or at the surface (at least 2 m visibility) were approached by either a small boat or personal watercraft from behind, with the aim of minimizing stress and flight response from the targeted animal and tagged with either a V13 or V16 acoustic tag using a 3 m modified pole-spear fitted with a tag applicator. Tags were inserted into the dorsal musculature at the base of the first dorsal fin using a titanium dart tethered to the tag. When possible, targeted sharks were visually checked for the presence of existing tags, and sex of individual was determined using a pole-mounted dip camera. Tagging effort via this method was largely directed at locations where sharks were known to seasonally aggregate. Shark size (TL cm) was estimated by comparison to an object of known size (the watercraft) from orthogonal aerial drone footage. Comparative methods such these have been shown to produce valid estimates of size (Sequeira et al., 2016;May et al., 2019). All tagged sharks were assigned a size-class based upon their measured or estimated size at time of tagging (Table 1). For sharks detected in more than one calendar year, a growth rate of 25 cm per year was assumed (Cailliet et al., 1985).

Spatial and Temporal Detection Analysis
Acoustic detection data were acquired from a wide-spread acoustic receiver array of up to 75 Vemco VR2 and VR2W receivers, deployed between Estero Bay (35.448, −120.952) and San Diego [32.876, −117.260 (Figure 1)], but also included offshore island monitoring sites at Santa Catalina Island (33.389,. From north to south, these included locations proximal to the following locations/landmarks: Estero Bay, Morro Bay, Pismo Beach, Santa Barbara, Ventura, Santa Catalina Island, Santa Monica Bay, Long Beach-Huntington Beach, Newport Beach-Laguna Beach, Dana Point-San Clemente, Oceanside, and San Diego (Figure 1). Receivers were largely installed and maintained by CSULB, but also included receivers of collaborators: University of California San Diego (UCSD), the Southern California Acoustic Telemetry Tracking Network (SCATTN), and the Ensenada Center for Scientific Research and Higher Education (CICESE). Acoustic receivers were deployed in shallow (∼ 3-40 m depth) coastal waters within 1,000 m of the shoreline (50% within 400 m), across a range of environments, including calm protected habitats (e.g., embayments, harbors), exposed high energy sand substrata (e.g., beach habitat), rocky reef, and kelp forest habitats. Receiver detection range (nominal range 150-700 m) varied by site, habitat type, and transmitter power output (Heupel et al., 2006;Kessel et al., 2014;Huveneers et al., 2016 To account for the possibility of false detections, raw receiver data were filtered to include only individuals that were detected two or more times per day (Simpfendorfer et al., 2015). Detection data were visually inspected to check for the possibility of double-tagged animals. Where double-tagged animals were identified, all data pertaining to the 2nd tag were removed from the analyses.

Social Network Construction and Association Definition
A gambit-of-the-group approach (Cairns and Schwager, 1987;Franks et al., 2010) was used to build proximity-based social networks using individual tagged sharks as nodes and strength of association as edges, calculated via the simple ratio index (SRI) using the R packages spatsoc (version 0.1.14; Robitaille et al., 2019) and igraph (version 1.2.5; Csardi and Nepusz, 2006). The SRI score's an individual's strength of association between 0 and 1, with 0 indicating no co-occurrence while a score close to 1 indicates a high level of co-occurrence (Aplin et al., 2013). SRI was chosen over other indices (e.g., half-weight index) as the properties and nature of acoustic detection data render them unlikely to violate the assumptions required for SRI: recorded associations are accurate, the probability of identification is independent of whether an individual is associated or not (Stehfest et al., 2013;Lilly et al., 2020). To comply with formatting requirements, detections were rounded to the nearest hour (i.e., hourly presence) and any subsequent detections of the same animal at the same receiver within that hour removed. Animal detections and associated relocations were grouped according to a temporal threshold of 1-h bins (time group), before being spatially grouped according to a threshold of 1 km radius (group). Thus, animals were considered to co-occur in space and time if they were detected at the same receiver within an hour of each other, or at separate receivers within an hour of each other, provided the two receivers were ≤ 1 km from each other. A 1 km spatial threshold was chosen as this allowed for individuals to be considered as being within the same group in instances where they were detected on separate receivers that had overlapping detection ranges, although creating such a spatial buffer introduces the possibility that sharks may be grouped together, or regarded as associating, when in-fact they could be up to 2 km from each other (assuming a 500 m detection range). Data were then subset by Julian year (2014-2019) before generating networks and associated graphs. Network communities were identified through weighted eigenvector community detection (Newman, 2006), and community modularity (Q, the extent to which communities within a network are distinct from each other) was calculated using the R package igraph. Homophilic propensity, the tendency of individuals to associate with others with similar traits, was examined by calculating assortment coefficients based on estimated shark size. Strength of associations were calculated as edge lists using a temporal based nearest neighbor approach, whereby the nearest neighbor to each individual within each time group was calculated, incorporating a distance threshold of 1 km radius.
Pre-network permutation tests based on mean values from 10,000 randomized networks generated in the R package spatsoc were used to examine whether the observed overall network structure (all years combined) differed from structure of randomized networks. Coefficients of variation (CV's) of SRI indices between vertices from observed networks were compared to those from the randomized networks using two-sample t-tests.
To further confirm the non-random nature of the observed structure, we used pre-network permutation tests (10,000) for each constructed network in each individual year of the study. CV's of SRI indices from observed and random networks were compared to examine observed and expected relationships (Bejder et al., 1998;Farine and Whitehead, 2015;, whereby tagged sharks were assumed to exhibit preferred co-occurrences where coefficients of variation from observed networks were greater than 97.5% of coefficients of variations from randomized (permuted) networks (Findlay et al., 2016). It was not possible to replicate this same process for the overall (all years combined) network as it was not possible to account for phenotypic changes in those individuals that were as identified as nodes in networks across more than 1 year.

Association Preference Examination
We used multiple linear regression to examine structural predictors of association strength (strength) in observed networks, with separate models run for each annual network. Predictor variables included estimated shark size-class (sizeclass), number of detections (detections), detection period (days), and number of acoustic receiver stations visited (statvis). The inclusion of size-class allows for identification of homophilic association preferences, while the inclusion of the number of detections, detection period, and the number of stations visited allows for identification of possible structure inherent to the data. The sex of the shark was not included as a predictor variable due to the number of animals in the study with unconfirmed sex (Supplementary Table 1). The most parsimonious model for each year was identified from a global model via AICc values using the R package glmulti (Calcagno and de Mazancourt, 2010), and confirmed by individually comparing the final model to the next two most parsimonious model iterations. Final models were also compared against null models (strength ∼ 1) to examine significance of predictors used (Supplementary Table 1). The global model was written as strength ∼ size class + statvis + detections + days + detections * days.
A multiple regression quadratic assignment procedure with double-semi-partitioning (MRQAP-DSP; Farine, 2013) was used to further quantify the influence of specific structural factors (i.e., size-class) included in multiple regression models in each year. To examine whether association strength differed for sharks grouped in the same size-class, a binary matrix was constructed that encoded size-class and tested for correlation with a binary version of the SRI matrix using Mantel tests (999 permutations) following the method described by Farine and Whitehead (2015), using the R package vegan (Oksanen et al., 2011). This process was repeated for each yearly dataset (2014-2019).

RESULTS
A total of 86 tagged juvenile white sharks were detected on at least 1 day across the study period (2014-2019) and included in the overall analyses. Of these, 29 (34%) were outfitted with internal tags, while 57 (66%) were tagged externally. Total number of days detected on acoustic receivers and sum detection period varied by shark, ranging from 1 to 424 days detected [66 ± 75 (mean ± SD)], across a sum period of 1-1,216 (265 ± 277) days (Figure 2). A total of 26 sharks were detected in 2 or more calendar years (range 2-4 years, 2.42 ± 0.58). The number of returning tagged sharks showed an increasing trend across the study (Figure 2) and was positively correlated with the cumulative number of sharks tagged (Pearson's product-moment correlation; r = 0.82, p = 0.048). Only externally tagged sharks were observed to represent nodes in networks in more than 1 year (Supplementary Table 2). Tagged juvenile white sharks exhibited both temporal and spatial overlap in detection patterns, suggesting possible association preferences (Figure 2).

Social Organization and Association Preference
Overall observed network structure (Figure 3) was compared against a null network generated from 10,000 pre-network permutations based upon the CV's of individuals, accounting for years in which individual sharks were present in a network.
Observed network strength values (mean = 0.36 ± 0.41) were significantly greater than expected values (0.27 ± 0.36; Welch Two Sample t-test; t = −2.513, p = 0.012), indicating there was apparent structure (i.e., association preferences) in juvenile white shark seasonal aggregations. Additionally, CV's of observed association strength indices vs. randomized networks for each individual year indicated that, with the exception of 2016, observed network structure was significantly stronger than expected from a randomly structured network (Supplementary Figure 1). Thus, co-occurrence of tagged juvenile white sharks at acoustically monitored locations was considered non-random. Community modularity for all detections across all years was high (O = 0.63), which reflects that although communities appear to be connected, they were spatiotemporally distinct with limited individual connections to other communities. Fifty-five of the 86 sharks (64%) included in analyses formed associations with another shark. Thirty-nine of the 86 sharks (45%) were a part of a significant network (Supplementary Table 2), with nine sharks (10%) identified as being part of a significant network in more than 1 year. The number of associations formed by individuals within a single year ranged from 1 to 10 (2.21 ± 1.65). The total number of associations per shark varied both within and between years throughout the study period but showed no discernable trend across the study period (Supplementary Table 1). Mean number of associations per shark was highest in 2015, which is reflected by the associated mean SRI value of 0.7. Lowest mean SRI values were observed in 2016 and 2019 (0.03 and 0.05, respectively), suggesting tagged sharks had low likelihoods of spatiotemporal overlap within acoustically monitored locations in those years ( Table 2).
The number of associations between sharks was not correlated with individual year (Pearson's product-moment correlation; FIGURE 3 | Sociogram showing associations of all tagged juvenile white sharks across all years. Nodes represent individual sharks, node color indicates the year in which they were detected or formed strongest associations (in the case of sharks detected over more than year). Edge thickness indicates the relative number of connections (associations) between nodes. Colors encircling multiple nodes represent identified distinct communities. Detected sharks that did not form dyadic connections to at least one other tagged shark (n = 31) are not shown. Community modularity for all detections across all years was high (Q = 0.63), which reflects that although communities appear to be connected, they were spatiotemporally distinct with limited individual connections to other communities. Network analyses based on weighted eigenvector community detection demonstrated that for all years combined (Figure 3), as well as each individual year, some aggregating sharks formed distinct, significant communities ( Table 2). In years where spatiotemporally separate aggregations formed, relatively high community modularity was observed, as indicated by high calculated O -values (Table 2, Figure 4, and Supplementary Video 1). Similarly, assortment coefficients were observed to be highest in years with higher O -values ( Table 2). Of the 39 tagged animals that were members of significant communities, 18 were detected in two or more spatially distinct locations (considered separate aggregations) across 2014 (n = 4), 2015 [n = 8 (three sharks from the significant network in 2014)], 2017 (n = 3), 2018 (n = 1), and 2019 (n = 2) (Supplementary Table 2). Sixteen sharks were identified as being members of significant communities across two or more years. Fourteen were sharks identified in communities in two separate years. Two sharks (Shark IDs 15_19 and 17_08) were part of significant communities in three separate years (2016, 2017, 2018 and 2017, 2018, 2019, respectively).

Observed Social Structure and Size Class
General linear models were run for each year in the dataset to examine potential predictors of association strength. For all years, the best fitting model was strength ∼ size class + detections * days (Figure 5 and Supplementary Table 3). The interaction term of detections * days was significant in all years for which it was included, except for 2018, where the interaction was not significant and these terms were removed, and the model rerun (Carey, 2013). Model correlation coefficients ranged from 0.19 to 0.83 (0.52 ± 0.23). The interaction term and response variable (association strength) were found to be strongly collinear across all models where included and is reflected by high associated variance inflation factors (Figure 5). Shark size-class was included in all final models, and was found to be a significant predictor of association for 3 of 6 final models [2017, 2018, and 2019 (Figure 5)]. The results from MRQAP-DSP regression indicated that shark size-class was a significant predictor of co-occurrence in 3 of 6 yearly datasets examined (2015, 2017, 2019; Table 3). Tagged shark co-occurrence was found to be significantly correlated with shark size-class in those same years [Mantel tests: (2015) r = 0.332, p = 0.001; (2017) r = 0.08, p = 0.04; (2019) r = 0.144, p = 0.014; Table 4].

DISCUSSION
Juvenile white sharks in southern California are known to seasonally aggregate in annually variable, spatially discrete nearshore locations (Lyons et al., 2013;White et al., 2019;Anderson et al., 2021a). The drivers behind these aggregation "hot spot" patterns, and the selection of specific habitat locations over an abundance of comparable available habitat remain unclear, and are not consistently explained by environmental correlates Anderson et al., 2021a). While the study included detection data from tagged juvenile white sharks at locations across southern California, a majority of annual detections and co-occurrences of tagged sharks occurred at seasonal aggregation hot spots, as indicated in Figures 2, 4. It must be acknowledge that a majority of tagged sharks included in the study (65%) were tagged externally, and the majority of these sharks were tagged at known aggregation locations (see Supplementary Table 1 for tagging locations). Thus, to some extent, detection patterns and observed network structures may be a function of tagging methodologies used. However, where this is the case, the patterns in the detection data also reflect that these sharks had already formed aggregations and may have also already formed preferential associations with individual conspecifics present within the aggregation. For instance, individuals caught offshore by commercial fishers were less likely to be found within nearshore aggregations. Individuals that spatiotemporally overlap each other are more likely to be associated with each other, due to inherently greater probability of co-occurrence. Pre-network (data stream) permutations can be used to examine whether a metric of interest (in this case observed association strength) differs to that which might occur at random. Our comparison of observed and permuted association strengths demonstrated that except for 2016 (a notably data deficient year) observed association preferences differed significantly to expected values, indicating apparent co-occurrence of individuals in space and time is non-random (Supplementary Figure 1). Observed relative association strengths in animal social networks may be a function of resource availability or site fidelity (Lusseau et al., 2006;Armansin et al., 2016). Thus, we incorporated methods that aimed to discern whether observed network structure could be explained by social factors, or were likely   driven by other forcing (e.g., environmental conditions and/or resource availability).

Assortativity and Community Structure
Animals that exhibit homophilic association preferences may do so with respect to phenotypic characteristics such as species, sex, size, and kinship (Mourier and Planes, 2021). As animal size was the only phenotypic characteristic we were able to consistently gather, we were restricted to using size-class as a metric of assortment. Although the resulting assortativity coefficients were variable, our analyses suggest the likelihood of association of tagged juvenile white sharks is greater where sharks are of similar size-classes. A comparable relationship was described in a study of juvenile white sharks in South Africa (Findlay et al., 2016). The low variance explained by our multiple regression and MRQAP analyses indicates community structure was also driven by other factors we were unable to quantify, which could include spatiotemporal overlap due to environmental factors and/or resource availability, as well as sampling effects. However, black tip reef sharks (Carcharhinus melanopterus) have been documented to exhibit preferential associations to conspecifics of similar size and sex regardless of prey availability (Mourier et al., 2012), which lends support to apparent preferential association toward conspecifics of similar size in the current study.

Shark Co-occurrence and Group Dynamics
Individuals in observed networks exhibit behaviors that best approximate fission-fusion dynamics (formation and dissolution of groups over time) with spatiotemporally unstable group membership. Fission-fusion dynamics are common in animal groups across taxa, including elasmobranchs (Haulsee et al., 2016;Perryman et al., 2019;Papastamatiou et al., 2020), and are regarded as being advantageous in the exploitation of heterogeneous environments (Ramos-Fernández and Morales, 2014;Silk et al., 2014;. Dyadic relationships seen between individual juvenile white sharks in our study are not represented by predictable or stable co-occurrence (association) as might be expected by animals that exhibit central-place foraging behavioral characteristics. Rather, the dyadic patterns observed in our study suggest individuals cooccur (fusion) punctuated by variable spatiotemporal intervals (fission). Such behaviors may be a function of exploitation of shared resources (e.g., food). A recent study of Australasian gannets (Morus serrator) determined that social associations in the study population were context dependent, and were most prevalent during foraging (constituting local enhancement). These associations were also most prevalent in regions where resources were clustered (Jones et al., 2020). Southern California juvenile white shark aggregation hot spots have been largely associated with shallow, sandy beach habitat, although across southern California, such habitat constitutes < 30% of total available shore type (Anderson et al., 2021a). Although we were not able to quantify resource homogeneity at aggregation hot spots, hot spot locations themselves were annually spatially variable. If spatial locations of aggregation hot spots are resource driven, we can expect there to be a comparatively higher degree of resource clustering at hot spot locations with respect to ostensibly similar neighboring habitat. Juvenile white sharks in the study were observed to co-occur at different receivers both within the same spatiotemporal aggregation, as well as separate, spatially discrete locations and aggregations within the same year (Supplementary Video 1). Although this behavior, which linked nodes between network communities, could simply indicate individuals at similar ontogenetic stages were responding in similar ways to the same environmental cues, it may also be attributable to following behaviors and association borne from local enhancement, where individuals are attracted to actively foraging conspecifics (Poysa, 1992). Thus, the interplay of environmental drivers and context-specific social behaviors may govern co-occurrence and apparent association of individual juvenile white sharks in the study. Although 10% of tagged sharks were identified at the same locations across more than 1 year, there was no evidence of stability in structure across multiple years, as returning sharks did not form repeated associations at monitored locations in subsequent years. Juvenile white sharks aggregating at nearshore locations, including those in this study, have been demonstrated to exhibit overlapping, restricted area use (Lyons et al., 2013;Anderson et al., 2021a). An individual's tendency toward association may be indicated by the propensity for spatiotemporal overlap and successive return to the same locations (Klimley and Holloway, 1999;Lilly et al., 2020), a trait exhibited in juvenile white sharks (Bruce et al., 2019;Spaet J. L. Y. et al., 2020;Anderson et al., 2021a). Observations of juvenile white shark interactions derived from drone surveys undertaken during the study period (Rex et al., unpublished data) indicate that individuals may often be within 10-20 m of each other while in these loose aggregations, and do not appear to be schooling in the manner seen in other aggregative shark species [e.g., scalloped hammerhead sharks (Sphyrna lewini), black tip reef sharks (Carcharhinus melanopterus)]. In addition, observations of these interactions suggest the presence of dominance hierarchies and conspecific aggression. Although tagged sharks generally exhibited low frequency of associations with other tagged conspecifics (Supplementary Table 1), the strength of those associations, as indicated by SRI indices, were relatively strong, with the exception of sharks in 2016 and 2019, which were both years with low mean and total associations ( Table 1 and Supplementary Table 1). Individuals within heterogeneous environments likely do not experience equal benefits from local enhancement and information derived from social behaviors and co-occurrence (Jones et al., 2020). Thus, fission may not be a singular synchronous event, but rather an ongoing social construct governed by intrinsic plasticity. Aggregation dynamics in juvenile white sharks may therefore be governed by the interplay of resource quality, environmental heterogeneity, and the continuous evaluation of the cost-benefit relationship of sociality. This can be likened to the hypothesis of the ecological loop that affects dispersal, proposed by Bowler and Benton (2005), whereby the size of a population (in this case at an aggregation site), which is mediated by the environment (which varies in space and time), governs interactions between individuals for resources.
Although sociality has been examined and described in a number of elasmobranch species (Sims et al., 2000;Jacoby et al., 2010Jacoby et al., , 2011Guttridge et al., 2013;Wilson et al., 2014Wilson et al., , 2015Armansin et al., 2016;Mourier et al., 2019;Perryman et al., 2019;Schilds et al., 2019;Mourier and Planes, 2021;Papastamatiou et al., 2020), the extent to which cooccurrence in white sharks as a function of social forcing factors is unknown. The difficulty of interpreting the true level of association between animals is a function of the data collection method, and the spatiotemporal thresholds used. For example, Schilds et al. (2019) used photographic identification of individuals present at a baited location within the same day as a means of identifying co-occurrence, thus arrival of connected nodes at the monitored location could potentially be hours apart. In the present study, shark presence and associations at monitored aggregation locations may be underestimates as they only reflect tagged sharks within spatially discrete locations and time thresholds (1 h), which may better represent natural co-occurrence. However, the potential for temporally co-occurring individuals to be up to 2 km apart from each other in some instances, due to the nature of acoustic telemetered data, may also represent over-estimates of spatio-temporal co-occurrence.
While the primary drivers behind co-occurrence in shark species may be extrinsic factors such as resource availability and use, additional adaptive benefits of social interactions beyond direct biological benefits (e.g., increased opportunity for reproduction) have been identified and described (e.g., Jacoby et al., 2010;Mourier et al., 2012Mourier et al., , 2019Armansin et al., 2016;Mourier and Planes, 2021;Papastamatiou et al., 2020). Papastamatiou et al. (2020) concluded that for sharks exhibiting central place foraging life history characteristics, information transfer during social foraging increased foraging efficiency and may lead to temporally stable social groups across years. Other shark species that exhibit more solitary life histories but aggregate seasonally at feeding sites may also exhibit preferential associations based on phenotypic characteristics (Findlay et al., 2016;Haulsee et al., 2016).
For naïve individuals, information transfer facilitates exploitation of patchy resources already identified by others (Aplin et al., 2013). The specific drivers governing where and when juvenile white sharks aggregate in southern California have yet to be defined but are thought to include water temperature and resource availability (White et al., 2019;Anderson et al., 2021a), and decisions by individual sharks to remain within a specific aggregation location may therefore reflect resource quality. Thus, perceived habitat quality derived from cues generated by both con and hetero-specifics, as well as transfer of information (e.g., following of "experienced" individuals by naïve individuals) within aggregation locations may lead to increased likelihood of association at a given location.

CONCLUSION
Although individual presence at aggregation sites may be driven by environmental cues and resource availability, the tendency of individuals to remain within aggregations may be the result of conscious behaviors, potentially driven by sociality, which may stem from naïve individuals acquiring information and experiencing local enhancement. To better understand the dynamics of associations between individuals, far higher resolution data is required. This could be achieved through the use of high-density acoustic arrays across much smaller spatial scales, specifically centered at aggregation hot spots [e.g., VPS (Vemco Positioning System) arrays], where the movements of individuals in three-dimensional space throughout an array can be achieved via trilateration of acoustic tag transmissions or by using proximity-based tags. Such studies would provide a more objective means of quantifying co-occurrence, fusion, and fission between individuals and groups, and elucidate the potential social dynamics that may govern space use and niche exploitation at juvenile white shark aggregation locations.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

ETHICS STATEMENT
The animal study was reviewed and approved by the CSULB Institutional Animal Care and Use Committee.

AUTHOR CONTRIBUTIONS
JA, AC, EB, and CL conceived the study. CL, EB, EM, CFW, RL, CW, PR, JM, KL, JA, AC, EG-R, and OS-N tagged the sharks and collected the data. CFW, EB, EM, BS, JA, PR, and JM organized the data. JA analyzed the data, created the main figures, and drafted the manuscript. All authors contributed equally to manuscript revisions.

FUNDING
Funding for tagging and acoustic monitoring was provided by Monterey Bay Aquarium and AltaSea. Funding for data maintenance and analysis was provided by State of California. Inkind support for tagging and monitoring was provided by Los Angeles and Orange County Lifeguards, Cities of Huntington Beach, Seal Beach, Ventura, Santa Barbara, Newport Beach, Laguna, Long Beach, and San Clemente.