Autonomous Aerial Vehicles (AAVs) as a Tool for Improving the Spatial Resolution of Snow Albedo Measurements in Mountainous Regions

We present technical advances and methods to measure effective broadband physical albedo in snowy mountain headwaters using a prototype dual-sensor pyranometer mounted on an Autonomous Aerial Vehicle (an AAV). Our test flights over snowy meadows and forested areas performed well during both clear sky and snowy/windy conditions at an elevation of ~2650 m above mean sea level (MSL). Our AAV-pyranometer platform provided high spatial (m) and temporal resolution (sec) measurements of effective broadband (310–2700 nm) surface albedo. The AAV-based measurements reveal spatially explicit changes in landscape albedo that are not present in concurrent satellite measurements from Landsat and MODIS due to a higher spatial resolution. This AAV capability is needed for validation of satellite snow albedo products, especially over variable montane landscapes at spatial scales of critical importance to hydrological applications. Effectively measuring albedo is important, as annually the seasonal accumulation and melt of mountain snowpack represent a dramatic transformation of Earth’s albedo, which directly affects headwaters’ water and energy cycles.


Introduction
Albedo is the proportion of incoming incident radiation that is reflected back to the atmosphere, and is a critical component of the Earth's energy and water balance [1]. This simple ratio governs many of the calculations of land-atmosphere interactions and climate models [2,3], quantifies the energy reflected from the surface back to the atmosphere [4], and varies considerably across time and space [5].
In mountainous regions that are characterized by seasonal snowpack, the spatio-temporal variability of albedo is particularly relevant as snow accumulation and melt cycles drive some of the most profound changes on the surface of the Earth [6,7]. This seasonal phenomenon dramatically alters the energy balance of the planet's mountain headwaters. In snow-free conditions lower albedo vegetation, soils, and rock can absorb~90% of short-wave radiation. Comparatively, higher albedo snowpack reflects up to 90% of visible radiation [8]. Additionally, the albedo of snow is not fixed, and decays as a function of time as snow crystals' physical properties change, particularly during

Opportunities to Improve Measurements of Albedo
Autonomous Aerial Vehicles (AAVs) bridge the scaling and measurement challenges associated with fixed position, in situ data providing near-Earth measurements horizontally and at multiple altitudes. This flexibility allows AAVs to support measurements for field campaigns that range from ha to km 2 . By design, AAVs are flexible and not fixed to a particular field location, allowing one craft to measure multiple sites with minimal marginal costs. The nascent AAV industry and associated novel sensors are projected to have a compound annual growth rate of 15.8% between 2019 and 2026 [27], resulting in improved aircraft and sensors at lower prices [28]. Combined, these decreased costs, improved hardware, and sensors create new opportunities to measure the Earth under a range of seasonal conditions. These data can help bridge the existing gap between ground-truthing field data and satellites [21], measure dramatic alterations of the Earth's surface [18], or be integrated into models as forcing or assimilation data.
Previous AAV-albedo focused work demonstrated the ability of a fixed-wing AAV to measure albedo over the Greenland ice sheet [29]. This study used small silicon pyranometers (a sensor that measures solar irradiance over an area (Wm −2 )) with the ability to capture a spectral signal of 300-1100 nm. This limited spectral range excluded some solar energy contributions to the surface energy balance and was tested on ice sheets with albedos ranging from 0.11-0.64 [29]. This range (0.11-0.64) captures on the lower scope of snow albedo, which ranges from 0.95 for fresh snow to 0.4 for shallow dirty snow [9]. Levy et al. [21] measured surface albedo across relatively level, forested stands in the northeastern United States using more advanced pyranometers (Kipp & Zonen CMP3 with a spectral range: 310-2700 nm and a mass of 350 g [30]). Reflected irradiance was measured from the AAV, and incoming irradiance was measured from a fixed position (a single pole). This provided an advance in the application of AAVs but did not measure albedo spatially across snowy landscapes. Also this study did not resolve the issue of a single fixed-point measurement of incoming radiation. Collecting fixed point data requires additional equipment (e.g., mounting pole/tower) and are not collocated with the AAV measurements. Webster & Jonas [12] furthered AAV-based measurements by employing upward-and downward-facing CMP3 sensors to measure albedo over a forested stand in Switzerland. Their study focused on forest snow albedo over a single~0.75 ha study site with flight paths of~170 m. Their work demonstrated that measurements of incoming and outgoing irradiance from a single AAV extends the ability to measure albedo variability over forested areas, at a range of altitudes, and at detailed spatial resolutions. All these data extend well beyond fixed-location measurements and at a much higher spatial resolution than satellite data.
The spatial flexibility of AAVs in the x, y, z plane begins to address the scaling issues (from tens of meters to several km 2 ) that have commonly hampered the calibration and validation of satellite albedo products such as MODIS and Landsat, and have not been addressed in previous research. Additionally, AAVs provide the ability to assess the albedo of the terrain under a range of conditions and in locations that are inaccessible or potentially hazardous for human safety (e.g., avalanche-prone areas).
To address these issues, we present viable materials and methodology to improve the measurement of effective broadband physical albedo from AAVs based upon the results of four AAV flights in an alpine catchment. Results document changes in albedo at the scale of tens of meters that directly address Hydrology 2020, 7, 41 4 of 16 the known albedo scaling issues [14,18,22]. This research focuses on effective broadband physical albedo (310-2700 nm) across snowy landscapes comprised of pure snow and coniferous forest.

Pyranometers and AAV
We tested the efficacy of Kipp & Zonen PR1 pyranometers (calibrated in Nov 2017), a retrofitted sensor built for this study with Kipp & Zonen's SMP3 hardware and encased in a smaller, lighter shell. These digital sensors have a spectral range of 310-2700 nm, a response time (95%) of 0.2 s, an effective field of view of 180 • , a mass of 145 g, and connect to a Meteon 2.0 digital data logger [31]. The PR1 provides a more efficient sensing platform for AAVs than the silicon pyranometers and CMP3 pyranometers used in previous studies [12,21,29]. CMP3 sensors employ passive correction, whereas the SMP3/PR1 hardware provides active digital temperature correction ranging from −20 • C to 50 • C. As a result, the PR1 provides higher accuracy measurements. The PR1 is also 59% lighter than the CMP3 (350 g lighter per sensor) [30], allowing for increased flight times and better handling for this AAV and payload configuration.
Upward and downward-facing PR1 pyranometers were mounted on fixed-position gimbals aboard a DJI Matrice 210 V2 (M210, Figure 1). This is a quadcopter design with a total wingspan of about 88 cm (35 inches) and without payload weighs approximately 4.53 kg (10 lbs). The upward-facing sensor measured the incoming irradiance ( K ↓) , and the downward-facing sensor measured the reflected irradiance from the surface ( K ↑) . Atmospheric corrections of albedo were assumed to be negligible at the low altitude of the AAV flights. For these test flights we assume that the ratio between incoming and reflected irradiance represents the effective broadband surface albedo (∝) of the snow surface [9] as measured by the pyranometers (Equation (1)).

Pyranometers and AAV
We tested the efficacy of Kipp & Zonen PR1 pyranometers (calibrated in Nov 2017), a retrofitted sensor built for this study with Kipp & Zonen's SMP3 hardware and encased in a smaller, lighter shell. These digital sensors have a spectral range of 310-2700 nm, a response time (95%) of 0.2 s, an effective field of view of 180°, a mass of 145 g, and connect to a Meteon 2.0 digital data logger [31]. The PR1 provides a more efficient sensing platform for AAVs than the silicon pyranometers and CMP3 pyranometers used in previous studies [12,21,29]. CMP3 sensors employ passive correction, whereas the SMP3/PR1 hardware provides active digital temperature correction ranging from −20 °C to 50 °C. As a result, the PR1 provides higher accuracy measurements. The PR1 is also 59% lighter than the CMP3 (350 g lighter per sensor) [30], allowing for increased flight times and better handling for this AAV and payload configuration.
Upward and downward-facing PR1 pyranometers were mounted on fixed-position gimbals aboard a DJI Matrice 210 V2 (M210, Figure 1). This is a quadcopter design with a total wingspan of about 88 cm (35 inches) and without payload weighs approximately 4.53 kg (10 lbs). The upwardfacing sensor measured the incoming irradiance ( ↓), and the downward-facing sensor measured the reflected irradiance from the surface ( ↑). Atmospheric corrections of albedo were assumed to be negligible at the low altitude of the AAV flights. For these test flights we assume that the ratio between incoming and reflected irradiance represents the effective broadband surface albedo (∝) of the snow surface [9] as measured by the pyranometers (Equation (1)). To ensure that the sensors were equally standardized, both PR1s were removed from the AAV and over the course of ~50 h recorded incoming irradiance, capturing a range of solar zenith angles on a rooftop with minimal topography or buildings in the field of view. Of concern was the influence that the AAV's non-retractable landing gear would have on the downward facing irradiance measurements ( Figure 1). A separate calibration was performed to quantify this measurement interference. For reference conditions, the upward-facing sensor was set up with an unobstructed field of view. For AAV-obstructed conditions the downward-facing sensor was mounted on the underside of the M210, which was then flipped upside down providing sky view conditions with the same obstructions that would occur during flight. This calibration was conducted over ~6 h on multiple days, across zenith conditions (10 h-18 h).
Because the PR1 sensors are prototypes, one of the primary challenges we encountered was securely fastening the pyranometers to the gimbals on the AAV. Three-dimensional (3D) printing facilitated the design and fabrication of mounts and allowed our group to develop, prototype, and test the mounting brackets and hardware at a relatively minimal cost. Mounts were based on measurements from standard engineering calipers, designed in AutoDesk Fusion360 (CAD software), and printed with a Flashforge Inventor II 3D printer using PLA filaments. The first test To ensure that the sensors were equally standardized, both PR1s were removed from the AAV and over the course of~50 h recorded incoming irradiance, capturing a range of solar zenith angles on a rooftop with minimal topography or buildings in the field of view. Of concern was the influence that the AAV's non-retractable landing gear would have on the downward facing irradiance measurements ( Figure 1). A separate calibration was performed to quantify this measurement interference. For reference conditions, the upward-facing sensor was set up with an unobstructed field of view. For AAV-obstructed conditions the downward-facing sensor was mounted on the underside of the M210, which was then flipped upside down providing sky view conditions with the same obstructions that would occur during flight. This calibration was conducted over~6 h on multiple days, across zenith conditions (10 h-18 h).
Because the PR1 sensors are prototypes, one of the primary challenges we encountered was securely fastening the pyranometers to the gimbals on the AAV. Three-dimensional (3D) printing facilitated the design and fabrication of mounts and allowed our group to develop, prototype, and test the mounting brackets and hardware at a relatively minimal cost. Mounts were based on measurements from standard engineering calipers, designed in AutoDesk Fusion360 (CAD software), and printed with a Flashforge Inventor II 3D printer using PLA filaments. The first test flights were conducted in a nearby agricultural field with fabricated PR1 replicas of the same dimensions and mass, which were mounted on the AAV to test for weight distribution issues and evaluate the integrity of the mounting brackets. These flights also provided us with an estimate of maximum flight duration. The Meteon Data Logger was attached to the M210 leg on a mount that was designed and 3D printed for this project. To enable implementation of these methods by other groups, all 3D CAD files used in the project are freely available for download at https://github.com/GEOSWRL/UAV_Albedo.

Flight Planning
The M210 provided a stable, robust platform to mount the sensors and data loggers. One of the challenges associated with this new-to-market (in 2019) AAV was instability with regards to several different flight planning software programs. We tested multiple software products (e.g., DJI, other proprietary, and open source software) during practice flights in a nearby agricultural field. All software underperformed. The firmware on the new-to-market AAV was not backwards compatible with existing flight planning software. A down-looking camera would greatly facilitate data analysis and interpretation. However, the camera for the AAV was backordered and unavailable, thus no high-resolution images were captured during our flights. The inherent hardware and software limitations of the AAV at the time of testing required that all paths and plans were flown manually, using line of sight and the field-of-view camera installed within the M210.
The M210 s internal GPS and Inertial Measurement Unit (IMU) recorded geographic location and altitude in addition to the pitch, roll, and yaw of the AAV. The PR1 sampled irradiance values every second, which were recorded on the data logger. For tables detailing the characteristics of the M210 and the PR1 sensors please refer to Appendices A and B.

Post-Processing of Data
Post-flight processing of the M210 data used DatCon [33], an open-source product, to convert the raw AAV data to a .csv file. The Meteon 2.0 data logger required accompanying free, proprietary software to download the data. To ensure the quality of the radiometer data any records were removed that corresponded to prior/during liftoff, during the initial ascent, and until the radiometer measurements stabilized. Merging the two disparate time series (radiometer data logger and the internal GPS data of the AAV) into a single composite was completed using the "Synchronize" function in MATLAB [34]. Data with temporal mismatches between the AAV and the data logger were discarded, so that every measurement corresponded to a measured lat/lon location.
were conducted between 11:29 and 12:30 in the day, and each lasted from 10-15 min, the safe limits of the battery power in cold conditions. The fourth flight was 2 min due to battery limitations. The solar zenith angles at noon on these dates were 43° and 38°, and the solar noon was 13:25 and 13:22 respectively. Both dates had meteorological characteristics distinct from one another. 15 April was windy and snowy (12 cm in the previous 48 h), and had an air temperature of −2 °C under cloudy skies. 1 May had clear sky conditions, an air temperature of 1 °C, with 5 cm of snow falling 48 h prior. Each pyranometer was mounted on a gimbal mount that ensures stability but does not ensure nadir. Thus to reduce error associated with the AAV's pitch and roll, a tilt limit of 5 • (based upon the IMU data from the M210) was applied [35], and data records with a pitch or roll >5 • were removed. This initial pass resulted in filtering 62 of 1381 points, or 4% of the data collected indicating stable flights even in windy conditions. The landscape where the flights were conducted was classified as Meadow (an open snow-covered area), Trees (stands of trees, comprised entirely of conifers), or Edge (a 10 m buffer around the stands of trees) using recent aerial imagery at 1m resolution [32] and standard GIS functions. Using the latitude and longitude locations from the merged AAV and radiometer dataset, each point was assigned one of the three landcover classifications based upon geographic location of the measurement.

Reducing Oversampling of AAV Measurements for Comparison to Space-Borne Data
The AAV collected albedo data at higher spatial resolutions than individual 30 m Landsat grid cells. To address potential oversampling biases of the AAV-based measurements, the central tendency of the data (mean, geometric mean, and median) were calculated over the co-located Landsat cell. These calculations were completed using the pandas and geopandas libraries in Python. For example, if there were eight AAV measurements collected within the spatial domain of an individual Landsat cell, the central tendency of the eight measurements was assigned to the Landsat cell for comparison. The intent of calculating the central tendency was to develop a better quantitative comparison of the AAV and Landsat data at the same spatial resolution, and to address the technical limitations of all flight paths and plans being flown manually.

Space-Borne-Data
The space-borne observations used for comparative analysis of albedo consisted of the daily MODIS 10A1 (MOD10A1-Snow Albedo) product [16] and calculations of broadband albedo using Landsat 7 and Landsat 8 Analysis Ready Data (ARD) [36,37]. The MOD10A1 data combines measurements from both Terra and Aqua satellites to retrieve directional hemispherical reflectance at local solar noon as both spectral (seven narrow spectral bands [16]) and three broadband (300-700 nm, 0.7-5.0 µm, and 0.3-5.0 µm) quantities. The MOD10A1 Terra and Aqua flights crossed the study area at approximately 12:25 and 2:00 local time, respectively, on both field days. The solar zenith angle is parameterized in the MOD10A1 albedo calculations [38]. The algorithm's coefficients for broadband albedo calculations were applied for Landsat 7 [39] and Landsat 8 [40]. The Landsat 8 flight paths crossed the study area on 14 April (13:12 local time, zenith angle of 43 • [36]) and 30 April (13:00 h local time, zenith angle of 24 • [36]), one day before each of our AAV missions. Landsat 7 was co-located with the study area on 1 May. Because of the clear sky conditions on this day and resultant bright snow over the study area, the Landsat 7 images were oversaturated and not usable for analysis.

Testing of AAV and Calibration of Sensors
Test flights with the replica PR1 sensors assessed the stability of the M210. These flights in April 2019 occurred in windy conditions, (up to~4 ms −1 ), and the AAV provided stable flights and responsive navigation with the replicas and data loggers mounted. Similarly, flights with the PR1 sensors, accompanying cables, and data logger in nearby agricultural fields provided stable flights and sample data from the sensors.
Ground-based instrument standardization of the PR1 sensors provided a high level of agreement without the landing gear present (R 2 = 0.998, Mean Absolute Error = 11.25 Wm −2 , Mean Absolute Error % = 2.89%) and with the landing gear present in the downward facing sensor's field of view (R 2 = 0.995, Mean Absolute Error = 24.0 Wm −2 , Mean Absolute Error % = 3.15%). The effects of the visible landing gear on the downward facing pyranometer were minimal (<3%), similar to ground-based calibrations of Webster & Jonas [12]. Field measurements from the downward facing pyranometer were calibrated by a factor of 1.0287 for the presence of landing gear. This bias correction and a high level of agreement (R 2 = 0.998) between the two pyranometers provide a high level of confidence that calculations of physical broadband albedo (a ratio) are accurate.

Metrics From Flight Data
The length and time of the four flights recorded by the AAV are provided in Table 1. While the flights were flown manually, every attempt to fly at a constant velocity and altitude during each flight was made. The flights on the first day were flown at roughly half the velocity as the second day of flights. Flights 2 and 3 were intentionally flown at a higher median altitude than flights 1 and 4. The differences between the calculated central tendencies of the data to address oversampling of albedo data over Landsat cells showed minimal variability (Figure 3). In Figure 3, there is little to no difference between the mean and geometric mean (Gmean-blue plot). The difference between the mean and median is nominal and evenly distributed (orange plot) with minimal skew or bias. Thus all data is displayed as the mean for each Landsat cell (Figure 4).
The differences between the calculated central tendencies of the data to address oversampling of albedo data over Landsat cells showed minimal variability (Figure 3). In Figure 3, there is little to no difference between the mean and geometric mean (Gmean-blue plot). The difference between the mean and median is nominal and evenly distributed (orange plot) with minimal skew or bias. Thus all data is displayed as the mean for each Landsat cell (Figure 4).  The shape of the plot widens as a function of data distribution and its probability density-where wider reflects more probable.

Albedo Measurements Across a Snowy Landscape
All four flights provided spatially explicit albedo measurements that clearly identified the different landcover types across the varying landscape, and across Landsat and MODIS grid cells (Figures 4 and 5). In Figure 4, the sub-figures in the left column (Figure 4a This flight moved at a faster velocity (~2 ms −1 ) but was still able to capture the transitions into forested areas (semi-circles Z). Flight 3 also highlighted that the albedo of snow is not uniform across open, snowy areas. This same heterogeneity is evident in Flight 4, which focused solely on snowy, open meadows but had a range of albedo values of 0.58-0.85. Of interest are the differences in the spatially averaged albedo values in between Flights 3 and 4 (Figure 4i,l), specifically in the southern central portion of the study area. In this sub-region, flight 3 shows lower albedo values than Flight 4, which was flown approximately 30 min later and at a lower altitude.
The resolution of the measurements from the AAV data is also telling. The AAV was able to provide measurements of albedo across the landscape, with a much higher degree of spatial detail as compared to the Landsat and MODIS data. different landcover types across the varying landscape, and across Landsat and MODIS grid cells (Figures 4,5). In Figure 4, the sub-figures in the left column (Figure 4a,d,g,j) present the albedo data as a function of flight distance, symbolized by the underlying landscape as defined in Section 2.4. The sub-figures in the center column (Figure 4b,e,h,k) show the individual data points and the pattern of each flight. The right column (Figure 4c,f,i,l) present the same data spatially averaged to the individual Landsat grid cell as described in Section 2.5.

AAV-Based Albedo Measurements Compared to Landsat and MODIS Data
The AAV provided a much more spatially explicit representation of albedo that is not encapsulated in measurements from the MODIS (MOD10A1-Terra and Aqua) and Landsat instruments. The ∝ values represent the entire AAV flightpath that is comprised of open meadows, transitional areas and stands of trees. These differences are not readily distinguished in the space-borne measurements, which would potentially provide albedo measurements that are markedly higher or lower than AAV-based measurements. In the cases of flights 2 and 3 (over a mix of meadow and trees) the MODIS ∝ is much closer to the AAV ∝ than the Landsat ∝, whereas in flights 1 and 4 (almost entirely over open meadow), Landsat data has a higher level of agreement with the AAV-based measurements. It is presumed that the high bias and much narrower range of Landsat albedo measurements is a result of the sensor not adequately capturing the influence of trees. This is in part due to Landsat's status as a Level-2 product that requires field-based data for validation and calibration [2,17,18]. Understandably, MODIS data and AAV data are distinct from one another. The MODIS product is estimated from narrow band measurements, while the AAV are broadband (310-2700 nm).

Discussion
The intended goal of this study was to test the ability of the AAV to fly with the prototype pyranometer and datalogger and take irradiance measurements across the landscape. In turn these AAV-based data are used to calculate broadband effective albedo (310-2700 nm) across a snowy mountain landscape. The methods and results we present focus on the testing and mounting of the pyranometers aboard the AAV and provide a case study of a viable and tangible method to collect albedo data across a mountainous catchment at these spatial scales. These data are complemented by the comparison with satellite data products.

Identifying a Path Forward for AAV-Based Measurements of Albedo
The results of this study suggest that measurements of irradiance from AAVs can improve the spatial resolution and detail of effective broadband albedo across the landscape. Tower-based measurements over forests or weather station-based measurements over meadows provide a single measurement of albedo, a deficiency that AAVs readily address. The prototype PR1 digital

Discussion
The intended goal of this study was to test the ability of the AAV to fly with the prototype pyranometer and datalogger and take irradiance measurements across the landscape. In turn these AAV-based data are used to calculate broadband effective albedo (310-2700 nm) across a snowy mountain landscape. The methods and results we present focus on the testing and mounting of the pyranometers aboard the AAV and provide a case study of a viable and tangible method to collect albedo data across a mountainous catchment at these spatial scales. These data are complemented by the comparison with satellite data products.

Identifying a Path Forward for AAV-Based Measurements of Albedo
The results of this study suggest that measurements of irradiance from AAVs can improve the spatial resolution and detail of effective broadband albedo across the landscape. Tower-based measurements over forests or weather station-based measurements over meadows provide a single measurement of albedo, a deficiency that AAVs readily address. The prototype PR1 digital pyranometers present considerable technological advances from previous studies. They have a spectral range of 310-2700 nm, are temperature corrected, have a faster response time and are much lighter than sensors used in previous research [12,21,29]. The PR1 pyranometers tested in this study are now commercially available and provide a means for an economical, lightweight sensor that is readily deployed in the field under a range of temperatures.
The AAV (M210) is readily available and performs well at high elevations in variable meteorological conditions. The design and fabrication of custom mounts using an off-the-shelf 3D printer allowed pyranometers to be securely mounted to the AAV, and we provide the design files, readily available for download (https://github.com/GEOSWRL/UAV_Albedo). The entire hardware framework requires a budget of around $15,000 USD in 2019 (~$11,000 for AAV, batteries, and accessories,~$3000 for sensors, and~$1000 for printer and filament) and can be deployed at multiple study sites that would readily support multiple studies with the same equipment. The aerial footprint of our longest flight covered over 8 ha in 14 min, and every indication implies that updated flight planning software will facilitate longer, more efficient flight paths in the future.
As mentioned previously, all flights were conducted manually and did not use flight planning software because of compatibility issues with the new-to-market AAV. During the spring of 2019, the Matrice 210 V2 did not have functioning flight planning software available. Efforts to manually complete an overlap pattern were attempted during flight 3. The goal of the flights was to test the ability and practicality of this AAV and pyranometer, and should be interpreted as an intermediate step that tests the viability of the platform. Future implementation of flight planning software that has subsequently become available will allow for more efficient and consistent flight paths flown at constant velocities and altitudes, which will provide replicable parameters for datasets in the x, y, and z plane. This study focused on the flexibility of an AAV to assess the variability of albedo in the x-y plane. An analysis of albedo variability in the z-plane was not in the scope of this research and presents opportunities for subsequent research.
As introduced in the results, there are differences in the spatially averaged albedo values of flights 3 and 4 (Figure 4f,h) in the southern central portion of the study area. Collecting more detailed data along the z-axis will allow subsequent flights to quantify the scale dependent controls on albedo. For example, as the AAV ascends over a snowy meadow how do measurements of albedo change? Is there a threshold or optimal altitude at which AAV flights should be deployed? The insights gained from these data would readily facilitate the calibration and validation of albedo calculations from the MODIS and Landsat satellites, and would offer the distinct possibility for improving validation and calibration measurements for current and subsequent satellite calculations of albedo.
Work by Bair et al. [14] highlight the influence of topography on albedo measurements. Our study area was over rolling terrain, with minimal topographic and snow surface relief (~20 m) that did not create shade during our flight times. Subsequent research in areas with higher relief will need to account for the effect of sensor angle with respect to surface topography. Potential methods include a bias correction in post-processing or physically mounting the sensors on gimbals that adjust to the slope and topography of the landscape. Similarly, the effect of solar zenith angles on AAV measurements and Landsat 8 calculations of physical albedo were beyond the scope of our study. Subsequent, more detailed comparisons and potential validation of satellite measurements will require the incorporation of solar zenith angles.
The effects of the landing gear on measurements of reflected irradiance were characterized (quantified and calibrated). There are AAV platforms with retractable landing gear, which are larger and heavier than the M210. One of the strengths of the M210 platform is its portability, requiring only moderate effort for a two-person team to ski into remote locations for measurements. We have tested a larger DJI M600 at another field location, but its size and weight limits its deployment to sites that are accessible by vehicle. Better quantifying the influence of the landing gear will provide enhanced bias correction of the reflected shortwave radiation measurements.
Our broadband pyranometers do not address albedo's variability across spectral wavelengths. It is of note that gimbal-mounted hyperspectral sensors are commercially available (i.e., Headwall or Resonon). These products are designed primarily for commercial agricultural applications, have a low temperature threshold of 0 • C, are heavier (requiring a bigger AAV), and are more expensive than the PR1 sensors presented in this research.
Our AAV framework operates on a smaller AAV (packable in size), at a considerably reduced cost, and is temperature corrected down to −20 • C. We do not presume that AAVs will replace space-borne observations that provide global coverage over an extended period of record. Instead we suggest that AAV-based measurements from a representative portion of a mountain headwater throughout the accumulation and melt season provide a means to augment existing EO and the ability to improve the validation and calibration of satellite data. The methods and resulting data presented in this paper provide a tangible way forward to produce data at relevant spatial scales to improve snow melt hydrological models and provide insights into the scaling issues associated with space-borne measurements of albedo.
Understanding these discrepancies in albedo values would have a considerable effect on energy and water balance calculations in this study area, especially computations focused on the spring melt and runoff [41], the effects of which propagate downstream through ecosystems and society. Improved albedo measurements of snowpack would also improve drought prediction and monitoring in snow dependent regions [42].

Conclusions
The advances we present in measuring albedo from an AAV improve scientific capacity, and subsequent efforts in the lab and field will improve the robustness of the platform and its measurements. Our measurements readily encapsulate the spatial variability of albedo across landscapes that is not discernable in coarser resolution satellite data. Our efforts should be interpreted as a first stage to directly address the timely, complex challenge of improved measurements and models of albedo [13]. These improvements will facilitate workflows to measure albedo at the hillslope and small watershed scale, greatly expanding from the current status quo of point-based measurements and broad scale remote sensing. These insights will provide better data for more geographically focused energy and water balance simulations, critical for improved hydrological modelling. Additionally, subsequent data from AAVs present great opportunities to augment space-borne calibration and validation measurements. These new data will not be tethered to a fixed location and have the inherent capacity to measure across rugged and hazardous terrain. While this study focuses on 11 ha total, this approach could easily be scaled and flown in other snowy mountains across the Earth.
All of the data presented in this manuscript have been published online (https://doi.org/10.5281/ zenodo.3386418, [43]) and are licensed under Creative Commons v.4. to the study area and for logistical support. Finally, we would like to thank the two peer reviewers for their suggestions and feedback that improved this manuscript.

Conflicts of Interest:
While the authors mention several companies and their products in this manuscript, it is provided for informational purposes only. The authors do not have a conflict of interest in providing these details.  Operating Temperature −20 • to +50 • C