Satellite Remote Sensing in Shark and Ray Ecology, Conservation and Management

have compiled a list of popular SRS data sources and sensors for common environmental variables in marine science. Our review of 71 papers (55 published in the last 10 years) identiﬁed ten SRS-derived environmental variables that have been used in elasmobranch studies, from a range of satellite sensors and data sources. Sea surface temperature and ocean productivity were the most frequently used variables. Articles primarily analyzed variables individually or in pairs, with few studies looking at a suite of interacting variables. Here, we present a summary of the current state of knowledge on the application of SRS, current gaps and limitations, and discuss some of the potential future directions in which we envisage this ﬁeld developing. Threatened elasmobranch populations inhabit some of the world’s most remote marine ecosystems. With often global coverage, SRS presents an opportunity to analyze the important environmental drivers of elasmobranch ecology to aid management decisions for the conservation of declining and threatened populations.


INTRODUCTION
Global populations of elasmobranchs (cartilaginous fish, including sharks, rays, and chimeras) have declined dramatically over the past 50 years (Baum et al., 2003;Burgess et al., 2005;Myers et al., 2007;Dulvy et al., 2014). Particularly susceptible to extinction (Musick et al., 2000), their K-selected life histories, slow growth rate, large body size, long gestation periods, late sexual maturity, and low fecundity make them acutely vulnerable to overexploitation (Reynolds et al., 2001;Dulvy et al., 2004), as well as less adaptable to environmental change through demographic responses (Perry et al., 2005;Rosa et al., 2014). Recent estimates have calculated that only 37% of elasmobranchs are considered safe (of least concern) by the IUCN Red List, the lowest proportion of non-threatened species among vertebrate groups studied to date (Molina and Cooke, 2012;Dulvy et al., 2014). The true extent, to which elasmobranch species are threatened, however, is still largely unknown. Many species are data deficient (Molina and Cooke, 2012;Dulvy et al., 2014), fished above recommended levels (Baum et al., 2003;Burgess et al., 2005;Simpfendorfer et al., 2011) and fisheries data on elasmobranchs is often poorly documented (Dulvy et al., 2008;Heupel and Simpfendorfer, 2010;Bornatowski et al., 2014). These factors combined, may mask their true population trajectories.
Reversing these declines through effective conservation and management strategies is important as elasmobranchs play a valuable role in marine ecosystems. As upper trophic level predators, many elasmobranchs exert top down effects by influencing prey populations through mortality and anti-predator behavior (Creel and Christianson, 2008;Heithaus et al., 2008;Ferretti et al., 2010). Some elasmobranch species also provide bottom up effects through the provision of nutrients into reef ecosystems from pelagic environments (Williams et al., 2018). The most appropriate conservation strategy can be difficult to define in many species due to their broad geographic distribution and migratory nature (Knip et al., 2012). Successful conservation of elasmobranch populations requires detailed knowledge of their breeding and feeding habitat, and their routes, timings and rates of movement between these areas Speed et al., 2010;Knip et al., 2012). As such, understanding how and why elasmobranchs move, migrate and aggregate is essential for their conservation and management at multiple spatial scales (Bonfil, 1997;Jacoby et al., 2012). Increased knowledge of movement patterns not only aid conservation strategies by optimally designing protected areas, but it also increases understanding of the risk posed by threats such as industrial fishing, habitat degradation, pollution and climate change (Chin et al., 2010;Espinoza et al., 2014;Schlaff et al., 2014).
Until recently, research has mainly focused on the "where" rather than the "why" aspects of elasmobranch movement ecology; where individuals or populations move and migrate, but not always why they do so (Hammerschlag et al., 2011). How the physical environment influences movement, is one of the key questions for understanding the movement ecology of marine megafauna (Speed et al., 2010;Hays et al., 2016;Jacoby and Freeman, 2016). Biotic factors, predator avoidance and prey availability, play a significant role in the movement ecology of elasmobranch species (Heithaus, 2001;Heithaus et al., 2002;Torres et al., 2006;Collins et al., 2007). However, permanent and transitory abiotic factors, for example bathymetry, sea surface temperature (SST), salinity and ocean currents, are also thought to strongly influence marine megafauna movements (Block et al., 2011). Schlaff et al. (2014) reviewed the influence of different environmental variables on the movement of elasmobranchs and identified strong environmental correlates of movement in some species. SST, for example, influenced the site fidelity of zebra sharks, Stegostoma fasciatum (Dudgeon et al., 2013), and the aggregation behavior of several shark species in north Western Australia (Speed et al., 2011). Furthermore, light intensity can affect daily movement patterns of elasmobranchs, acting as a trigger for foraging (Nelson et al., 1997;Cartamil et al., 2003;Andrews et al., 2009) and resting behavior (Papastamatiou and Lowe, 2012;Dudgeon et al., 2013).
Correlates of elasmobranch movement with other environmental variables are less clear. For example, salinity can alter the distribution and local abundance of elasmobranch species (Hopkins and Cech, 2003;Ubeda et al., 2009), but as most elasmobranch species are stenohaline, movement in response to salinity changes is most likely a response to avoid physiological stress (Schlaff et al., 2014). Productivity, often measured as the distribution and density of chlorophyll-a, has been linked to occurrence, movement and aggregation behavior in large filter feeding elasmobranchs, such as whale sharks (Rhincodon typus) and manta ray species (Heyman et al., 2001;Sleeman et al., 2007;Dewar et al., 2008;Couturier et al., 2011). In contrast, the relationship between primary productivity and movement is less clear in high trophic level species (Gende and Sigler, 2006). Most investigations have focused on how a small number of abiotic factors influence movement (Schlaff et al., 2014). However, abiotic factors vary in their level of impact between species, sex and life-history stage, and seldom act in isolation. It is often difficult to establish the relative contribution of different variables to movement and spatial patterns (Sleeman et al., 2007). As such, how the physical environment interacts with individual behavioral motivation to influence the movement of elasmobranchs is still not fully understood.
Data for the marine environment can be derived from both in situ or remotely sensed sources. In situ data sources include localized sensors and profiling floats [(e.g., the Argo program 1 , (Roemmich et al., 2009)]. These provide highly accurate point measurements but can have limited temporal and geographical coverage, and data is not always relayed in real time (Guinehut et al., 2004;Roemmich et al., 2009). Alternatively, remote sensing acquires information about the ocean surface without being in physical contact with the medium being measured (Martin, 2004). This can be done using sensors mounted on ships or aircraft, however, the most effective and common method is via satellite (Martin, 2004). The use of satellite sensors allows wide ranging observations of oceanic processes over various spatiotemporal scales (Yang et al., 2013). Satellite sensors can be passive, measuring natural sources of reflected or emitted electromagnetic radiation, or active, measuring backscattered electromagnetic radiation that they have produced [for more detail see Pettorelli et al. (2018)]. Although it provides a more indirect measurement than in situ methods, numerous environmental variables are being detected by satellite remote sensing (SRS) (Chassot et al., 2011), with the accuracy and resolution of several variables now comparable to in situ instrumentation (Zieger et al., 2009;Young et al., 2011;Lee, 2016). In addition, SRS does not have the spatial limitations of in situ environmental data (Zieger et al., 2009) and many satellites provide global coverage with open data policies, such as the European Space Agency's Sentinel missions and the NASA/USGS Landsat Program (Malenovský et al., 2012;Roy et al., 2014). In recent years, SRS has become widely used for oceanic environmental data collection (Guinehut et al., 2004;Zieger et al., 2009;Yang et al., 2013).
Integrating SRS with information on elasmobranch distributions and movements allows behavior to be directly linked to spatially extensive measurements of ocean environmental variables. This can aid the prediction and mapping of likely future space use, distribution and behavior (Handcock et al., 2009), further improving our understanding of elasmobranch ecology. Biotelemetry techniques are frequently adopted to explore the movements and behavior of wide-ranging elasmobranchs Donaldson et al., 2014). The integration of SRS data with animal-borne tag data is a relatively new field, largely due to the spatial and temporal disparities between the two data types Rose et al., 2015). Recently, however, these disparities have begun to be resolved through modern technological developments in tagging, greater processing ability and increased access to data (Santos, 2000).
Here we review the application of SRS for elasmobranch ecology, conservation and management, with the following aims: (1) examine current uses of SRS in elasmobranch research and its integration with tagging technologies, (2) investigate which environmental variables and associated metrics are used most frequently, (3) provide an overview of the myriad of available SRS data sources and sensors, and (4) highlight gaps in knowledge and potential future directions in this emerging field of research.

MATERIALS AND METHODS
This review took a systematic approach to searching the literature, enabling greater reproducibility, transparency and rigor (Tranfield et al., 2003;Armitage and Keeble-Allen, 2008;Moher et al., 2015). Several academic search tools were used in conjunction to conduct the bibliographic search, as recommended in the systematic review literature (Bramer et al., 2013(Bramer et al., , 2016Giustini and Boulos, 2013). These included ISI Web of Science 2 , Scopus 3 , and Google Scholar 4 . As the first study on remote sensing in elasmobranch ecology was undertaken in the early 1980s (Priede, 1984), the search included articles published between 1 January 1980 and a cut-off date of 1 May 2018. The term "shark" was used as it obtained greater search results than the term "elasmobranch." Due to the ambiguity of the term "ray, " however, "ray" and "elasmobranch" were used together. Additionally, the terms "movement, " "habitat, " "ecology" and "environment * " were added as Boolean terms to our searches. How these terms were applied varied depending on the search engine (see Supplementary Table S1). Two separate searches were conducted using Google Scholar. The first used traditional Boolean search terms. The "shark" search produced 20,500 results and the first 1% (205) were filtered for relevance. The "elasmobranch AND ray" search produced 1,310 results. To match the "shark" search, the first 205 results were analyzed (see Supplementary Figure S1 for numbers produced by each search, per search term). There are biases toward highly cited articles and the occurrence of a search term in an article's title in Google Scholar's ranking algorithm (Beel and Gipp, 2009). As such, standard Google Scholar search results will contain highly cited articles that provide a broad overview of the topic (Giustini and Boulos, 2013). Specific or relevant literature may be missed or ranked so low that they will not be found. To account for these biases, a second search of Google Scholar (named "Smart Scholar") was used, directly searching in the subject and title for our search terms (see Supplementary Table S1).
The initial searches yielded 412 articles for "shark" and 215 for "elasmobranch AND ray." From this initial list, articles were selected for inclusion in the review following a process adapted from Moher et al. (2015). First, articles were excluded if their titles did not contain the terms "shark" and/or "remote sensing" (or any associated terms). Second, articles were excluded if they did not contain the terms "shark" or "ray" (or correlates) in the abstract, and satellite remote sensing was not used as the source of environmental data in the methods section. Finally, the full text and the data sources of the articles were screened, and articles were excluded if; elasmobranch data were not included as a data source for the study; the article was a review or from a non-peer-reviewed/respected source; or the full text could not be obtained. Additionally, 54 duplicate articles found between the collective searches were also removed. In total, 71 articles were retained for inclusion in the review. The overall process used for identification and screening of articles for the systematic review is outlined in Supplementary Figure S1.
Key information and metadata (as summarized in Supplementary Table S2) were extracted from the final 71 articles and consolidated into a database. Study applications were identified for each article (i.e., the predominant research question for which SRS data were used). Geographic coordinates for study locations were obtained either directly from the study if reported, or by obtaining latitude and longitude of study locations 5 . If data were obtained over a survey area, rather than a single site, then the average geographic coordinates for this area were used. Within the literature, the terms used to describe environmental data are highly changeable. In this review we use the term "variable" to describe a single measure of the environment, such as "SST, " or "salinity." The term "metric" is used to describe a specific measure for quantifying a variable. For example, multiple metrics related to a single variable might be used in one article, or the same variable might be referred to differently in different articles (e.g., SST, thermal anomaly or change in SST). As such, we consolidated similar metrics into 10 different variables (Table 1).
To provide a detailed list of data sources and sensors available to researchers, those mentioned in the 71 reviewed

REVIEW OF CURRENT KNOWLEDGE
The systematic literature search identified 71 articles for review. Further details of the bibliographic search results can be found in Supplementary Table S5 and Supplementary Figure S2.
Reviewed articles were published from 1984 to May 2018. The number of publications, however, increased considerably between January 2000 to May 2018 (n = 68) compared to the previous two decades (n = 3). Journal articles made up the majority (n = 63) of the final reviewed articles. Other sources include conference proceedings (n = 2), book chapters (n = 2), Ph.D. theses" (n = 3) and reports (n = 1). This section presents an overview of the current use of satellite remote sensing (SRS) within elasmobranch research. It includes a summary of the elasmobranch data used in the studies, and the geographical and taxonomic focus of the articles. Details of the SRS-derived environmental variables used within the reviewed articles are also presented.

Elasmobranch Data
Use of tagging technologies was common, with 48% of reviewed articles using animal-borne tags as a source of elasmobranch data, compared to 25% using fisheries data, and 18% using scientific surveys. Baited remote underwater video systems and citizen science each made up 4% of elasmobranch data sources. A single article used previously published data as the source of elasmobranch data for the investigation. Out of the 71 reviewed articles, 32 combined satellite-based telemetry (e.g., SPOT tags and PAT/PSAT tags) with SRS. For example, Sleeman et al. (2010) used satellite telemetry and SRS data to show that whale sharks move independently of geostrophic currents in the Indian Ocean. The use of SRS data with acoustic telemetry was not common (2/71 articles).

Geography, Taxonomy and Threat Status
There was a broad geographic spread of study locations covering the Atlantic, Pacific and Indian Oceans (Figure 1). There was a notable absence of research in the north-western Pacific Ocean, the Artic and the Antarctic. A total of 37 publications (52%) conducted studies in the Atlantic Ocean compared to the Pacific Ocean (27%), and Indian Ocean (20%). A single study by Sequeira et al. (2014) used data taken from the Atlantic, Pacific and Indian Oceans to assess ocean current, and predict future distributions of whale sharks. Studies were more common in temperate waters (63%) than in tropical regions (37%) and more common in coastal waters, within 200 nautical miles from a major land mass (79%) than open ocean (21%). These results may reflect the location of primary research institution (see Supplementary Figure S3) and the financial and time constraints of working in oceanic, opposed to coastal, environments. It should be noted that the geographic distribution is likely to be biased by the fact that this review was only able to consider publications written in English.
Of the 12 extant elasmobranch orders (Pavan-Kumar et al., 2018), eight were represented in the reviewed articles (Figure 1), representing of 37 different species. Unsurprisingly, there was a taxonomic bias toward highly mobile species with Carcharhiniformes and Lamniformes dominating, both featuring in 34% of the reviewed articles. Orectolobiformes (28%) were the next most frequent elasmobranch Order. The other remaining Orders included two ray Orders (Myliobatiformes, and Rajiformes, 10 and 1%, respectively), and three additional shark Orders (Heterodontiformes, Squaliformes and Squatiniformes, 1% each). A further 1% of articles did not specify the species the study was undertaken on. Of the 37 different model species, only one, the Brazilian guitar fish, Pseudobatos horkelii, was listed as Critically Endangered

SRS Environmental Variables
Two remotely sensed environmental variables dominated the review articles; SST (n = 64, 90%) and ocean productivity (n = 48, 68%), with sea surface height, ocean current metrics and metrics of frontal data present in 11, 6 and 6% of the articles, respectively (Figure 2). SST metrics were utilized for a variety of research applications including investigations into habitat and distribution (Laurs et al., 2006;Curtis et al., 2014;Lea et al., 2015;Miller et al., 2015;Reynolds et al., 2017), foraging ecology (Sims and Quayle, 1998;Weng et al., 2008;Weeks et al., 2015), elasmobranch fisheries (Bigelow et al., 1999;Mitchell et al., 2014;Selles et al., 2014) and track reconstruction Queiroz et al., 2010;Ramírez-Macías et al., 2017). Metrics of ocean productivity, such as primary productivity and chlorophyll-a, were used primarily in investigations into habitat and distribution but also foraging ecology. For example, Jaine et al. (2012) showed that primary productivity was a principal driver of reef manta ray (Manta alfredi) aggregations and Humphries et al. (2010) used metrics of ocean productivity to demonstrate that predatory fish, including several species of elasmobranch, conduct behavioral switching of search strategies in response to variation in productivity and thus food availability.  Most articles focused on one or two variables (n = 27 and n = 24, respectively), with investigations analyzing three or more variables less frequently (three variables, n = 9; ≥ four variables, n = 11).

SRS Data Sources and Sensors
In total 16 different data repositories were used to obtain remotely sensed environmental variables and 72% of variables used in the reviewed articles had a specified data repository (Supplementary Figure S4). The National Oceanic and Atmospheric Administration (NOAA) (42%) and the National Aeronautics and Space Administration (NASA) (24%) were the primary sources for data acquisition. European satellite products, such as those from Copernicus Marine Environment Monitoring Service, were also used for data acquisition but were far less common (5%). Of the 71 reviewed articles, only one (Solanki et al., 2003) used commercial SRS data, to investigate potential fishing zones for an Indian elasmobranch fishery. All other reported SRS data came from free and open sources. Outlays for commercial SRS can be high (Turner et al., 2003(Turner et al., , 2015. Despite falling costs, conservation and ecology remains underfunded (McCarthy et al., 2012) and researchers will only use commercial SRS data if they can afford them (Turner et al., 2015). In addition, the increased resolution of free and open SRS data is sufficient for many research questions related to the marine environment, precluding the need for expensive high-resolution commercial data. Furthermore, policies from institutions such as the United States Geological Survey (USGS) and the European Space Agency (ESA) to provide SRS data without cost (Woodcock et al., 2008;Donlon et al., 2012) has led to improved access to comprehensive global environmental variables that, due to the expense, were previously unavailable to a large part of the research community. Currently, there are few repositories available in the literature for researchers to obtain information on the large number of different SRS variables. Where collated data do exist, [e.g., Toth and Jóźków (2016); Pettorelli et al. (2018)] there is a heavy focus on the terrestrial realm and often reference to satellite systems rather than specific environmental variables. Chassot et al. (2011) and Kachelriess et al. (2014) both summarize remote sensing variables for marine realms, however, these are not exhaustive, and a host of new sources and metrics have been made available since their publication. In this review, we provide a collection of different web platforms available to obtain SRS marine variables at the global scale ( Table 2), as well as collating information on the various sensors which collect environmental data ( Table 3). These lists are non-exhaustive but will add to and update previous collated data from reviews. SRS data used in the studies came from multiple satellite sensors (Figure 2). A full glossary of the sensor acronyms can be found in Supplementary Table S3. The sensor used to collect the environmental data were specified for 64% of variables (Supplementary Figure S4), with the most commonly used sensors being AVHRR (n = 42), MODIS (n = 42), and SeaWiFS (n = 27), together providing 90% of the SRS variable data. The sensors used in the reviewed articles were primarily passive optical sensors, with bands at visible and infrared wavelengths. Most sensors used were onboard polar orbiting satellite platforms, delivering frequent global coverage (approximately every day) at spatial resolutions typically of the order of a few kilometers. One sensor onboard a geostationary orbiting satellite was used for SST retrieval (in combination with SST data from other polar orbiting satellites) (Domeier et al., 2012), the GOES (Geostationary Operational Environmental Satellite) Imager, a four-channel visible and infrared radiometer capable of imaging the western hemisphere every 30 min. In addition, two microwave sensors were used; the AMSR-E passive microwave radiometer was used for SST observations (n = 4) (Domeier et al., 2012) and the SeaWinds active microwave scatterometer provided wind data in one study (Wickham, 2011).
Importantly, the spatial and temporal resolution of the SRS data products used varied considerably between studies: spatial resolutions ranged from 1 km to over 100 km and temporal resolutions from daily to monthly. The spatial and temporal resolutions of SRS-derived variables used in the reviewed articles are summarized in Table 4.

Habitat Use and Distribution
Most studies identified in this review used SRS to investigate habitat use and the distribution of elasmobranch populations (n = 48, 67%). Only by understanding where animals reside, and how they move, can the impact of their threats be evaluated. Reynolds et al. (2017) combined remotely sensed SST and chlorophyll-a with movement data from satellite telemetry to show that large areas of suitable habitat for whale sharks were not covered by the marine protected area (MPA) network in the south-eastern Indian Ocean. Oh et al. (2017) used SST from the SeaWiFS sensor in combination with baited remote underwater video systems to map species specific habitat requirements for juvenile sharks, the management of which is important for the survival and persistence of coastal shark species (Cortés, 2002). Furthermore, Lea et al. (2015) found that SST may trigger summer migrations in tiger sharks (Galeocerdo cuvier) in the northwest Atlantic. This baseline information on habitat use is of critical importance in elasmobranch conservation and has been the main focus of studies using SRS, reflecting both the relative novelty of this technology as well as the lack of baseline ecological data for many elasmobranch species globally (Dulvy et al., 2014). This and the other main applications of SRS in elasmobranch research are summarized in Figure 3.

Fisheries
A total of nine studies (13%) integrated SRS variables with fishery data. When SRS was applied to elasmobranch fisheries, SST (primarily from AVHRR sensor) and chlorophyll-a (primarily from MODIS Aqua/Terra or SeaWiFS sensors) were principally utilized. In the reviewed articles, SST and chlorophyll-a were used to inform interactions for stock assessment models, explore spatial distributions, investigate factors affecting catch rates, and for developing models for predicting fishery catches of commercial elasmobranch species (Bigelow et al., 1999;Iglesias et al., 2004;Carvalho et al., 2011;Novianto et al., 2016). For example, Carvalho et al. (2011), Mitchell et al. (2014 and Novianto et al. (2016) all used SRS data to find optimal abiotic factors for peak catch per unit effort of blue sharks, Prionace glauca, in the southwest Atlantic Ocean, English Channel and Indonesia, respectively, whilst Solanki et al. (2003) found that using SST and chlorophyll-a to predict fishing zones increased catch per unit effort in elasmobranch fisheries in India.

Track Reconstruction
A total of seven studies (10%) used SRS data to reconstruct movement tracks from electronic tag datasets. SST is used in post-processing models to correct raw geolocations from tagging technologies, such as PSAT tags, to increase the horizontal accuracy of animal movement tracks (Ferreira et al., 2018). For example, Queiroz et al. (2010) matched tagderived SST values with SRS-derived SST values to acquire more precise latitudinal estimates for the movements of blue sharks in the north-eastern Atlantic Ocean. PSAT tags require a substantial amount of extrapolation and FIGURE 2 | Frequency of SRS variables split by the sensor used to obtain the variable. Multiple metrics of an SRS variable were counted (see Table 1). Number of articles using a specific variable are in parentheses. A list of sensor acronyms can be found in Supplementary Table S3. interpolation to reconstruct the horizontal movements of elasmobranchs and SRS data can be essential in this post-processing.

Foraging Ecology
A small number of articles (n = 7, 10%) used SRS to study elasmobranch foraging ecology. For example, Weeks et al. (2015) used SST, chlorophyll-a and photic depth to show that anomalous river discharge levels can trigger feeding frenzies in reef manta rays and several elasmobranch species have been found to shift foraging strategies when entering areas of high primary productivity (Sims and Quayle, 1998;Weng et al., 2008;Humphries et al., 2010). For studies exploring foraging ecology, the course spatial resolution of SRS data [which may be 50 km 2 or more (see Table 3)] compared to the area of the foraging study (typically measured in the order of seconds of latitude) has led to researchers using in situ observations, such as those from the Argos program (Roemmich et al., 2009), conductivity, temperature and depth profilers (CTDs) (Hayes et al., 1984), or sensors on tagging technologies (Wilson et al., 2005;Teo et al., 2009), rather than SRS data. However, increasing spatial resolution and reduced costs of SRS data mean that SRS can play an important role in future investigations of elasmobranch feeding and foraging.

Limitations of SRS
Despite the benefits offered by SRS in multiple areas of elasmobranch research, satellite observations also have their limitations. The relatively coarse spatial resolution of SRS data can affect their application to many research topics, as spatial and temporal variability in environmental conditions can affect all types of animal movement at any scale (Dodge et al., 2013). High variability over a research area can pose a particular issue for those studies investigating large scale movement (Dodge et al., 2013). Restrictions in temporal resolution, either from mission design or cloud cover, may reduce the capabilities of a study to capture events or changes in environmental dynamics that occur over short periods of time . Constraints on spatial resolution of SRS data can also be a limiting factor. Position data of elasmobranch occurrence for research into behavior and ecology is available from a variety of sources, such as satellite or acoustic telemetry, fishery data, and boat or aerial-based surveys (Cartamil, 2009;Nadon et al., 2012;Sequeira et al., 2014;Lea et al., 2015;Weeks et al., 2015). The spatial error of these data sources is generally much smaller than the spatial resolution for most SRS-derived data (Espinoza et al., 2011;Hammerschlag et al., 2011), which for marine environmental variables is rarely less than 300 m and can be up to 50 km (Donlon et al., 2012; Table 3). As such, there is a still degree of mismatch in resolutions between position/movement data and SRS data. This, however, is becoming less of a problem, as recent and proposed satellite launches are greatly increasing the spatial resolution of available environmental data. For example, commercial satellite sensors now have spatial resolutions < 1 m and free and open data is available with pixel sizes as low as 10-60 m, depending on the band ; all of which are within the error associated with marine tagging and movement technologies (Palmer et al., 2015). Due to technical limitations, SRS data is only sensitive to surface environmental conditions (Takano et al., 2009;Chassot et al., 2011). Ocean environmental variables are seldom uniform, and elasmobranchs rarely inhabit a narrow, shallow range in the water column. Analyzing changes of environmental variables on a vertical scale, and how they affect elasmobranch behavior, is a relatively novel, but important, area of research, particularly given the number of species that make diel vertical migrations (Takano et al., 2009). The integration of SRS data with models can help overcome this limitation and tell us what is in deeper layers of the ocean (Brodie et al., 2018;Kvile et al., 2018). For details see "Future Directions and Opportunities." In many of the reviewed articles, SRS was used as the source of environmental data to investigate the influence of abiotic factors on elasmobranch movement (Weng et al., 2008;Humphries et al., 2010;Sims et al., 2012). The use of SRS for movement ecology is dependent on the type of tag used for positional data collection. SRS-derived data is regularly used with satellite tagging studies [e.g., SST (Lea et al., 2015), sea surface height (Vandeperre et al., 2016) and primary productivity (Vaudo et al., 2017)], but for certain variables (such as temperature and salinity), tag-derived data and in situ measurements are preferred (Stokesbury et al., 2005;Wilson et al., 2006;Stevens et al., 2010;Hammerschlag et al., 2011). Studies that use acoustic telemetry for positional data primarily use in situ measurements due to the coarse spatial resolution of SRS data. However, many tagging technologies that record detailed environmental data, such as conductivity, temperature and depth satellite delay relay loggers (Boehme et al., 2009;Roquet et al., 2014), are still not appropriate for elasmobranchs, and the use of in situ measurements for species that move large distances is unrealistic (Dodge et al., 2013). As such, increased integration of SRS-derived environmental variables with both satellite and acoustic telemetry is expected to increase in the future as satellite data becomes more readily available and accessible, and as sensor capabilities continue to improve .
An important finding of the systematic review was that the inclusion of SRS metadata in research articles was extremely variable (Supplementary Figure S4). Notably, many articles did not include key metadata, making it difficult to both interpret and replicate these studies. Furthermore, for those studies that did present metadata information, the presentation was neither consistent nor standardized, further complicating critical assessment of the different methods. Details such as spatial and temporal resolution are important as the required characteristics of SRS data depend on its application. In coastal zones, processes typically occur at small spatial scales and waters are often optically complex. They require greater spatial resolution and an increased number of spectral bands compared to the open ocean, where changes occur over much larger scales and sensors with low spatial resolution and fewer spectral bands are sufficient (Lee and Carder, 2002;Lee et al., 2007). Without the spatial and temporal resolution of the data the conclusions of an investigation cannot be fully appraised. A standardized approach to reporting SRS information in research articles should be adopted to enable transparency and reproducibility, and to facilitate more efficient exchange of information between SRS experts and marine research communities (Cooper and Gavin, 2005;Balestro et al., 2013).

GAPS IDENTIFIED
Our review highlighted the following research gaps in the application of SRS within elasmobranch ecology, conservation and management:

Limited Set of Environmental Variables Used
The majority of studies used SST and chlorophyll-a, and only a small number used three or more environmental variables in concert (Weng et al., 2008;Hacohen-Domené et al., 2015;Vandeperre et al., 2016;Lopez et al., 2017). Multiple other variables available from SRS, such as salinity, wind speed and ocean current, frequently influence elasmobranch behavior (Schlaff et al., 2014), yet were rarely used in elasmobranch research (Figure 2). In particular, despite the importance of salinity in driving elasmobranch movements (Hopkins and Cech, 2003;Ubeda et al., 2009;Schlaff et al., 2014), only two studies used SRS-derived salinity data.

Limited Use of European SRS Data
Despite its utility for studying the ecology of other wide-ranging marine predators, e.g., Pardo et al. (2015), Bailey et al. (2018), ESA data were used infrequently in elasmobranch research and data from ESA's Sentinel satellites were not used in any of the studies in our review. This may reflect the fact that historically satellite data from the United States has been easier to access than European data. Until recently ESA products, such as those from the Medium Resolution Imaging Spectrometer (MERIS), were available only after an application process (Verstraete et al., 1999;Bezy et al., 2000), which may have limited their use. In addition, Sentinels 1, 2 and 3 were only relatively recently launched, and this is likely to cause a delay in the publication of studies using these satellites. The Sentinel satellites offer huge potential for investigations of marine megafauna, as detailed in the "Future Directions and Opportunities" section.

Limited Use of Google Earth Engine
Interestingly, one major repository of SRS data not found in the reviewed articles was Google Earth Engine (GEE). GEE is a cloud-based geoprocessing platform, launched in 2010, that enables access to a wide variety of SRS datasets, including the Landsat series, Sentinels-1, 2 and 3, and MODIS, for analysis either within Earth Engine or to download and analyze externally with other software (Gorelick et al., 2017). Although GEE has had some use in research on terrestrial fauna (Joshi et al., 2016;Lewis et al., 2017), use of GEE in wildlife research is still infrequent compared to other data sources, in both marine and terrestrial realms. One reason for the slow uptake may be because GEE requires knowledge of either JavaScript or Python coding languages. The time and steep learning curve associated with learning this new skill, may put some researchers off of utilizing GEE, despite its clear benefits.

SRS Data and Acoustic Telemetry Are Rarely Combined
Despite the frequent use of acoustic telemetry in elasmobranch research (Heupel et al., 2018), the use of SRS data with acoustic telemetry was not common. Although in situ data is currently the preferred source of environmental variables for studies involving acoustic telemetry, in remote oceanic regions, where many elasmobranch species reside, obtaining and managing in situ environmental data may be difficult and expensive. With the recent improvement in the spatial resolution of open access SRS data, we expect the number of studies integrating SRS and acoustic telemetry to increase markedly, particularly in remote oceanic regions.

Limited Use of SRS Data in Conservation and Management
Although SRS has become a powerful tool for the conservation of terrestrial wildlife it is less frequently used in the marine environment (Kachelriess et al., 2014). SRS can establish baselines of quality and extent of habitats, detect loss and recovery (Nagendra et al., 2013;Pettorelli et al., 2018), and monitor population sizes and anthropogenic threats to endangered species and biodiversity (Donohue and Foley, 2007;Fretwell et al., 2012;Pettorelli et al., 2018).

Minimal Representation of Threatened Taxa
Of the seven most threatened elasmobranch families (Dulvy et al., 2014) only two were represented in our review, Alopiidae and Squatinidae. No freshwater elasmobranchs, some of the most threatened elasmobranch species (Dulvy et al., 2014), featured in our selected articles. This is perhaps unsurprising as freshwater elasmobranchs have received considerably less attention than marine species (Martin, 2005). SRS of freshwaters poses its own unique challenges, as these environments are optically complex (Brezonik et al., 2015) and require higher spatial resolution data for analyses.

FUTURE DIRECTIONS AND OPPORTUNITIES
Based on the gaps identified in the systematic review, this section summarizes the potential advances in SRS sensors and methodologies for furthering elasmobranch research. It serves to highlight some of the future opportunities for integrating SRS with elasmobranch positional data, in addition to more direct applications for elasmobranch fisheries, habitat degradation and marine protected area (MPA) management.

Sentinel Satellites
ESA's Sentinel satellites 7 provide a wealth of data for marine research. The first Sentinel satellite, 1-A, was launched in 2014, carrying a C-band synthetic aperture radar (SAR). SAR is a form of radar that creates two dimensional images and has an advantage over optical sensors as it is not impeded by cloud cover or low light levels. SAR can be used to monitor oceanic waves, which provide data on wind speed, wind direction and wave height (Tomiyasu, 1978). Sentinel 1-A was followed by Sentinel 2 in 2015, which carries a high spatial resolution optical sensor primarily used for land and coastal monitoring, useful for studying change in important elasmobranch habitats, such as coastal zones and coral reefs (Malenovský et al., 2012). Sentinel 3, launched in 2016, measures ocean and land surface temperature, ocean and land surface color, and sea-surface topography. The Sentinel missions have a fundamental focus on providing ocean data at a higher resolution than earlier satellites and sensors, such as MODIS and SeaWiFS (Donlon 7 https://sentinel.esa.int/web/sentinel/missions Malenovský et al., 2012). Data from the Sentinel series has already been a useful tool for research into other marine predators, for example predicting abundance and densities of cetacean species (Pardo et al., 2015;Bailey et al., 2018). As such, data from Sentinel satellites will provide an invaluable source of information for investigating elasmobranch ecology, conservation and management in the coming years.

Derived Products
Primary measurements obtained directly from SRS data, such as SST, chlorophyll-a, and wind speed, are typically used directly in research. However, these primary variables are actually proxies for other variables, such as oceanic fronts or upwelling regions, which are better representations of oceanic habitats (Hobday and Hartog, 2014). Although not as common as primary SRS products, derived SRS products have been used to investigate foraging ecology of seabirds and cetaceans (Tew Kai et al., 2009;Scales et al., 2014), habitat hotspots for pelagic sharks (Queiroz et al., 2016), movement physiology in marine mammals (Abrahms et al., 2018), and regions of marine biodiversity (Miller and Christodoulou, 2014). Using these derived products, as opposed to primary SRS products, may lead to improvements of predictive models and a better understanding of the behavioral drivers of marine fauna (Hobday and Hartog, 2014).

Predicting Vertical Profiles
Satellite products can be assimilated into models, such as the Hybrid Coordinate Ocean Model (HYCOM) (Chassignet et al., 2007), the Regional Ocean Modeling System (ROMS) (Haidvogel et al., 2008), or the Bluelink Ocean Data Assimilation System (BODAS) (Oke et al., 2008), that develop three-dimensional ocean models taking vertical differences into consideration. HYCOM assimilates satellite altimetry (sea surface height) and satellite SST, along with various in situ datasets, and generates predictions of environmental variables, such as water velocity, SST, and salinity, at depths from 0 to 5000 m (Cummings and Smedstad, 2013). Investigating drivers of vertical migrations will improve our ability to predict species occurrence, investigate foraging strategies, and assess the effect of anthropogenic disturbance on elasmobranch species (Sims et al., 2006;Comfort and Weng, 2015). Further developments in the incorporation of SRS data into models that predict vertical profiles will impact our current understanding of many species' ecology. For example, vertical data from the Regional Ocean Modeling System (ROMS) increased the suitability of species distribution models for four pelagic species off the Californian coast (Brodie et al., 2018), and enabled modeling of the transport of zooplankton Calanus finmarchicus, an important food resource for commercially important fish species (Opdal and Vikebø, 2015).

Salinity
The integration of additional environmental variables into elasmobranch research, beyond the frequently used SST and chlorophyll-a, is likely, especially as new satellite products become available. There is potential for SRS-derived salinity data to be integrated into studies of mobile marine megafauna. Only with the launches of the European Space Agency's (ESA) Soil Moisture and Ocean Salinity (SMOS) mission (Kerr et al., 2010) and NASAs Aquarius mission in the early 2010s has SRS-derived salinity been available to researchers. Retrospective analyses using new SRS data have been used previously, for example for determining levels of suspended matter (Dekker et al., 2002). Thus, using salinity data from both the Aquarius and SMOS missions, retrospective analyses on past data could expand our knowledge of marine fauna, such as elasmobranchs.

Investigating Multiple Environmental Drivers Simultaneously
The wide range of SRS products and sensors currently available offers greater scope for investigating multiple environmental drivers simultaneously. Future investigations using three or more environment variables will be important to not only help disentangle the mechanistic drivers of elasmobranch movement but also to fully understand how elasmobranch species may respond to human mediated changes in their environment (Schlaff et al., 2014). Integrating environmental variables with animal movement is challenging due to the size and frequency of data collected and at different temporal and spatial scales (Dodge et al., 2013;Ferreira et al., 2018). This is now being facilitated by tools such as Bio-ORACLE (Tyberghein et al., 2012), MARSPEC (Sbrocco and Barber, 2013), and R packages rsMove (Remelgado et al., 2018) and rerddapXtracto (Mendelssohn, 2017), which have been developed to aid this integration process more manageably. Increased access to new variables, the rise of cloud-based processing platforms such as Google Earth Engine or Amazon Web Services and the advances in multivariate statistical approaches will continue to facilitate this process and analyses further. In addition, new artificial intelligence initiatives, such as Microsoft's "AI for Earth" and Google's "AutoML" have led to an explosion in machine learning algorithms suitable for processing SRS-derived "big data" that will further progress this field.

Acoustic Telemetry
Most articles considered in this review used SRS-derived environmental data with satellite, rather than acoustic, telemetry as the source of elasmobranch positional data. Acoustic telemetry is particularly useful for investigating short term fine-scale, and long term meso-scale patterns on individual movement and behavior (Heupel et al., 2006;Meyer et al., 2010). The longer battery life of acoustic tags and often internal implantation allows individuals to be tracked over longer time scales (years rather than months) compared to satellite tags (Hays et al., 2007;Heupel et al., 2018). Acoustic arrays may vary in size from regional through to national and global networks (Heupel et al., 2018). Researchers are now generating shared repositories for acoustic data for multiple marine species [e.g., the Australian Integrated Marine Observing System (IMOS), (Hill et al., 2010); the Ocean Tracking Network, (Cooke et al., 2011); and Movebank, (Kranstauber et al., 2011)] that combine multiple "independent" arrays and provide acoustic tracking data at spatial scales more suitable for integration with most SRS products. Data from these integrated arrays, combined with SRS data, will present opportunities to measure movement responses to large spatial and temporal scale environmental and oceanographic processes (Hazen et al., 2012). Simultaneously, the increasing spatial resolution of SRS data is beginning to align with scales more appropriate for regional acoustic array integration.

eDNA
We envisage SRS playing an important role in elasmobranch conservation, particularly in conjunction with environmental DNA (eDNA). eDNA is an efficient and non-invasive way of sampling biodiversity, particularly useful for animals with restricted accessibility, such as cryptic or endangered species [for more information see Thomsen and Willerslev (2015) and Barnes and Turner (2016)]. Combining SRS environmental data with eDNA offers an opportunity to create species distribution and habitat maps to expand the biological knowledge and aid conservation management of threatened elasmobranch species (Boothroyd et al., 2016).

Population Monitoring
High spatial resolution satellite imagery (as opposed to satellite-derived variables), with resolutions now down to less than 1 m, can be used to identify and count individual animals and estimate abundance. Examples include; southern right whales (Eubalaena australis) (Fretwell et al., 2014); emperor penguins (Mirounga leonina) (Fretwell et al., 2012); southern elephant seals (Mirounga leonina) ; Weddell seals (Leptonychotes weddellii) (LaRue et al., 2011); and southern hairy-nosed wombats (Lasiorhinus latifrons) (Swinbourne et al., 2018). LaRue et al. (2017) identified three main characteristics of species to be successfully monitored by satellite imagery; open landscape; organism color that contrasts with the landscape; and sufficient organism size. Certain elasmobranch species may therefore be suitable for monitoring using satellite imagery. Given calm conditions, the ocean surface is a suitable landscape for detection, and certain species such as whale sharks, basking sharks (Cetorhinus maximus) and manta rays have both the requisite color contrasts with the landscape and sufficient size to be successfully identified using very high-resolution satellite imagery. Although potentially limited by time at the surface, as well as the high-cost involved, satellite imagery could be used as a viable non-invasive tool for population monitoring of elasmobranch species, far exceeding the coverage obtained through drone technologies currently being employed for this purpose (Rieucau et al., 2018).

Fisheries
There is great potential in integrating environmental data from SRS with elasmobranch catch data to improve the quotas within elasmobranch targeted fisheries, both enhancing efficiency and reducing overfishing (Chassot et al., 2010). For example, blue shark catch is thought to be constrained by SST, frontal data, sea surface height, and chlorophyll-a concentration (Bigelow et al., 1999;Queiroz et al., 2012;Mitchell et al., 2014;Selles et al., 2014). By monitoring chlorophyll-a, fishing effort can be streamlined, preventing overfishing by decreasing effort in periods of low primary productivity and increasing back to appropriate levels during periods of high productivity. As well as overfishing, bycatch of elasmobranchs in non-target industrial fisheries is a primary driver of elasmobranch declines (Molina and Cooke, 2012;Dulvy et al., 2014). Current bycatch reduction methods are primarily based on regulation of fishing zones (Molina and Cooke, 2012). However, data on diel and seasonal variability in bycatch species is needed for efficient implementation, which is not always available (Molina and Cooke, 2012). SRS might be used to investigate baseline data on environmental drivers of elasmobranch bycatch, as well as subsequently informing spatial and temporal closures of fishing zones. To our knowledge this has not currently been undertaken in elasmobranch species, however, it has been used effectively to limit bycatch of loggerhead (Caretta caretta) and leatherback (Dermochelys coriacea) turtles in Hawaiian pelagic longline fisheries (Howell et al., 2008(Howell et al., , 2015. In recent years, dynamic ocean management, where management changes spatially and temporally depending on observed biological, oceanographic and socio-economic data, has increased (Maxwell et al., 2015;Dunn et al., 2016). Dynamic ocean management can greatly improve the efficiency of management of fishery catches and closures (Dunn et al., 2016). Using near-real-time SRS data, such as ocean color data from Sentinel-3, would improve the efficiency of dynamic management, helping avoid overfishing of already depleted stocks as well as preventing food security and economic uncertainty in fishery-dependent countries (Chassot et al., 2010). Increased availability and resolution, and decreased costs of SRS data might therefore improve dynamic ocean management by reducing bycatch rates and enhancing the efficiency of targeted elasmobranch fisheries (Dunn et al., 2016).
Allocation of fishing efforts, as well as control of fishing mortality, can aid in preventing overfishing (Maunder, 2002). High-resolution SRS imagery for vessel presence, used in conjunction with the Automatic Identification System (AIS), can calculate vessel density, fishing effort, rate of unreported vessels in specific fisheries, monitor vessel responses to management changes, identify illegal vessel numbers in closed regions, and reduce costs of current surveillance and control measures (Kourti et al., 2005;Corbane et al., 2010;Brusch et al., 2011;Bradley et al., 2018). Use of SRS for this purpose to date has been limited due to insufficient spatial coverage of very high-resolution imagery (Kourti et al., 2005) and because AIS systems are not a pre-requisite on all fishing vessels . However, significant progress is envisaged in the next few years due to increased satellite numbers and coverage, and mandatory vessel monitoring systems. SRS has potential for monitoring vessels where AIS may not be present, such as artisanal and illegal, unreported and unregulated IUU) fishing vessels. Tools such as the Search for Unidentified Maritime Objects (SUMO) ship detector (Greidanus et al., 2017), the Integrated System for Surveillance of Illegal, Unlicensed and Unreported Fishing (INSURE) (Kurekin et al., 2019), and VIIRS boat detection data (Elvidge et al., 2018;Oozeki et al., 2018), have enabled real-time detection of IUU vessels, and the identification of vessels down to 1 m and fishing vessels using lights to attract catch, all of which can enable detection of illegal fishing by smaller artisanal vessels further aiding enforcement decisions. These systems could make SRS a powerful tool for surveillance and control, particularly in remote areas where little or no enforcement is currently undertaken (Corbane et al., 2010;Kachelriess et al., 2014).

Habitat Degradation
Land use change damages marine as well as terrestrial biodiversity (Kachelriess et al., 2014), with fertilizer use, land clearing and sediment run off, significantly impacting marine environments (Fabricius, 2005;Chérubin et al., 2008). In addition, sea level rise may severely impact important habitats, such as nesting areas, of marine fauna (Baker et al., 2006;Fuentes et al., 2010). The spectral signatures from SRS data can be used to monitor changes in the seas and oceans, such as turbidity, depth, temperature, salinity, sea surface height, and flow, all of which can impact elasmobranch species (Doxaran et al., 2002;Knip et al., 2010;Kachelriess et al., 2014). In the last decade SRS has been developed sufficiently to detect oil spills (Brekke and Solberg, 2005;Klemas, 2010), allowing the early detection, tracking and prediction of spills to be monitored, which is important for clean-up and effective counter measures (Kachelriess et al., 2014). Furthermore, there is potential for the development of hyperspectral sensors to monitor ocean plastic pollution (Kachelriess et al., 2014), and although full scale research is yet to be undertaken in this area, a proof of concept for this methodology has already been established (Garaba et al., 2018). SRS variables, such as ocean current, wind and SST, have been used to guide surveys to monitor levels of plastic debris in the Pacific Ocean (Pichel et al., 2007) and to assess the seasonality in risk of marine debris to Hawaiian monk seal (Monachus schauinslandi) populations, aiding their management (Donohue and Foley, 2007). Thus, SRS shows huge potential for monitoring dynamic threats to marine animals, including elasmobranchs, in future.
SRS data and high-resolution imagery is also likely to play a role in monitoring important elasmobranch habitats, such as mangroves and coral reefs, that act as nursery areas or refugia by evaluating the distribution, health status, productivity, resilience and resistance to change over time (Kuenzer et al., 2011;Duncan et al., 2018). SRS-derived SST and altimetry data can also monitor events, such as bleaching and sea level rise that may harm coral reef populations (Hoegh-Guldberg, 1999;Hughes et al., 2003). By identifying key factors influencing bleaching or disease events, such as coral cover and coral type, SRS might help generate metrics for evaluating and predicting the health of coral reef ecosystems (Hochberg et al., 2003;Xu and Zhao, 2014). Although mangroves and coral reef systems have been monitored using SRS, further development, and integration with elasmobranch ecological data is required. Global mapping of both mangrove and coral systems is limited or non-existent, and unified classification systems for monitoring mangroves and coral reefs are not available (Heumann, 2011;Xu and Zhao, 2014). Development in these areas will greatly aid elasmobranch monitoring, conservation and management.

Marine Protected Areas
Outside of fishery management, a primary model for elasmobranch conservation has been the establishment of marine protected areas (MPAs) (Davidson, 2012). Although the effectiveness of MPAs for elasmobranch conservation is still under debate (Davidson, 2012;Dulvy, 2013;Hilborn, 2015), MPAs have been shown to have measurable benefits to elasmobranch populations (Garla et al., 2005;Knip et al., 2012;White et al., 2017). The success and efficacy of MPAs are determined by multiple factors, with habitat identification being a critical feature (Edgar et al., 2014). SRS will likely become increasingly useful for monitoring environmental correlates of marine biodiversity, such as primary productivity, temperature and ocean currents, which may in turn be used for the identification of appropriate habitats and biological hotspots to feed into MPA designation (Palacios et al., 2006;Platt and Sathyendranath, 2008;Kachelriess et al., 2014). Further, SRS data may also play a role in the direct assessment of current MPAs for the conservation of elasmobranch species (White et al., 2017). Considering the dynamic, changing and unpredictable nature of certain marine zones, dynamic MPAs, developed with near-real-time SRS data, that vary spatio-temporally have been proposed (Hobday and Hartmann, 2006;Game et al., 2009). In addition, enforcement plays an important role in the success of MPAs (Edgar et al., 2014). As discussed, both high-resolution optical satellites, synthetic aperture radars and AIS provide options to aid and enhance enforcement efforts in MPAs by monitoring vessel presence and movement via imagery (Kourti et al., 2005;Corbane et al., 2010). In addition, SRS-derived data could be used to identify habitat hot spots within MPAs for more directed enforcement effort (Ross et al., 2011;Irigoyen et al., 2018). Consequently, we feel SRS will play a valuable role for the selection and enforcement of MPAs in future.

CONCLUSION
To our knowledge, this is the first systematic review of the use of satellite remote sensing in shark and ray ecology, conservation and management. Substantial progress in this field has been made in the last 10 years with remotely sensed environmental data facilitating significant advances in our understanding of species distributions, habitat use and movement, particularly for highly-mobile elasmobranchs. SRS data are now regularly integrated with satellite telemetry to address questions about elasmobranch ecology and behavior. Principally from free and open data sources, researchers now have access to a suite of environmental variables for analyses at varying resolutions, from multiple satellite sensors. Therefore, to support future research in this field, we have compiled inventories of online repositories for satellite data and current sensors available for the marine realm.
This review has identified several significant research gaps where satellite remote sensing has yet to be fully utilized. Most studies focused on SST and chlorophyll-a, and the wide range of satellite variables accessible to researchers are rarely combined to explore collective influence. Future studies would benefit from investigating multiple variables simultaneously and incorporating additional SRS variables that influence elasmobranch movements, such as salinity and ocean currents. Despite being a valuable source of marine data, use of European satellites has so far been limited in elasmobranch studies. In addition, the use of SRS has been restricted to marine elasmobranchs and, although presenting its own challenges, SRS for the conservation of freshwater species should also be explored. Application of SRS is limited by the coarse spatial resolution of many current satellite products but sensor capabilities are continually improving. Furthermore, we identified a need for more consistent reporting of SRS metadata for variables used in analyses to facilitate reproducibility.
There are substantial opportunities for satellite remote sensing to further advance our understanding of elasmobranchs and other marine fauna in the future. The improved availability and resolution of free and open satellite data sources, such as ESA's Sentinel missions, and the increased use of derived and modeled SRS ocean products, will enhance the accuracy of predictive models and improve our understanding of behavioral and ecological processes. The increasing number of networks of receivers with national and international programs, and collaborations within the research community, combined with increased spatial resolution of SRS data, will improve the application of SRS with acoustic telemetry. These developments will aid in future applications for conservation and management, particularly in the areas of population monitoring, fisheries, habitat degradation and MPA design/management. This will have a huge impact particularly for species occupying remote ocean regions, where resources for science are scarce, where in situ data is often lacking and where little protection is afforded to them.

DATA AVAILABILITY
All datasets generated for this study are included in the manuscript and/or the Supplementary Files.

AUTHOR CONTRIBUTIONS
MW, ET, and DJ conceived the review. MW conducted the literature survey methodology, analyzed the data, led the writing of the paper and prepared figures and tables. ET, DJ, and TD provided guidance on the survey design and analysis. All authors contributed to editing the manuscript.

ACKNOWLEDGMENTS
This review was undertaken as part of the Ph.D. thesis of MW, at the Department of Geography, King's College London and the Institute of Zoology, London. The authors would like to thank the two reviewers for their comments to help improve the manuscript.