Fine-scale variation in microhabitat conditions influence physiology and metabolism in an Antarctic insect

Microhabitats with distinct biotic and abiotic properties exist within landscapes, and this microhabitat variation can have dramatic impacts on the phenology and physiology of the organisms occupying them. The Antarctic midge Belgica antarctica inhabits diverse microhabitats along the Western Antarctic Peninsula that vary in macrophyte composition, hygric qualities, nutrient input, and thermal patterns. Here, we compare seasonal physiological changes in ve populations of B. antarctica living in close proximity but in different microhabitats in the vicinity of Palmer Station, Antarctica. Thermal regimes among our sample locations differed in both mean temperature and thermal stability. Between the warmest and coldest sites, seasonal mean temperatures differed by 2.6˚C and degree day accumulations above freezing differed by a factor of 1.7. Larval metabolic and growth rates varied among the sites, and adult emergence occurred at different times. Distinct microhabitats also corresponded with differences in body composition, as lipid and carbohydrate content of larvae differed across sites. Further, seasonal changes in carbohydrate and protein content were dependent on site, indicating ne-scale variation in the biochemical composition of larvae as they prepare for winter. Together, these results demonstrate that variation in microhabitat properties inuences the ontogeny, phenology, physiology, and biochemical makeup of midge populations living in close proximity. These results have implications for predicting responses of Antarctic ecosystems to environmental change.


Introduction
Landscapes are heterogenous, and even a seemingly uniform landscape can have a vast array of microhabitats that differ in abiotic and biotic characteristics (Risser, 1987). Small ectotherms like arthropods can experience considerable variation in environmental conditions across a landscape, either as a consequence of moving among microhabitats or being isolated within a microhabitat. For example, leaves directly exposed to sunlight might be several degrees warmer than leaves in the shade, exposing insects like leaf miners to signi cantly different thermal conditions even if they are living on the same tree (Pincebourde et al., 2007). Thus, it is important to consider ne-scale variation in microhabitat conditions when assessing phenology, physiology, behavior and potential responses to climate change (see reviews by Terblanche, 2006, andWoods, 2020). During winter and in polar regions, even slight variation in microhabitat conditions can have signi cant impacts. For example, microclimates can vary considerably depending on whether hibernacula are above or below the snowpack, and these slight differences in microenvironment signi cantly in uence overwintering success and fecundity (Irwin and Lee, 2003;Marshall and Sinclair, 2012). For the Antarctic springtail Cryptopygus antarcticus, microhabitat has a notable effect on cold acclimatization, as individuals directly exposed to air are more cold hardy than individuals embedded a few centimeters into moss (Hawes et al., 2008). Despite the importance of microhabitat for overwintering arthropods, this ne-scale variation is often overlooked when assessing local adaptation and responses to climate change.
In Antarctica, a small percentage of land is ice-free during the austral summer, and these ice-free areas contain most of the terrestrial biodiversity on the continent (Convey and Stevens, 2007;Convey, 2013;Convey et al., 2014). These habitats are exceptionally heterogeneous in composition, yielding an array of microhabitats that present unique challenges to their inhabitants (Holdgate, 1977). Abiotic factors that contribute to this microhabitat diversity include freshwater availability, topography, proximity to the ocean, and exposure to prevailing winds, among others (Convey, 1997;Convey et al., 2014). Local patterns of biodiversity can also in uence microhabitat conditions. For example, localized nitrogen deposition by marine vertebrates (i.e., penguins and seals) creates nutrient-rich areas that support more diverse plant and invertebrate communities (Bokhorst et al., 2019). These differences in plant communities in turn in uence soil temperature, as moss-dominated habitats tend to have cooler active layers than habitats dominated by vascular plants (Guglielmin et al., 2008). Together, variation in habitat leads to a diverse thermal pro le across terrestrial Antarctic habitats, such that invertebrates experience differences in mean temperature, thermal variability, and timing of freeze-thaw events (Convey et al., 2018).
However, while the consequences of microhabitat variation have been addressed for insects at temperate latitudes, particularly in forested ecosystems (Pincebourde and Woods, 2020), similar investigations for Antarctic insects are lacking.
The Antarctic midge Belgica antarctica, the only insect endemic to Antarctica, has a 2 + year life cycle spent mainly in the larval stage; adults live only 10-14 days and emerge only during the brief summer period (Lee and Denlinger, 2015). Larvae have limited mobility, and although the wingless adults have been found airborne in moderate wind (> 20km/h) and rafting on the surface of meltwater streams (Edwards and Baust, 1981), the midge has an otherwise limited ability to disperse (Peckham, 1971; Atchley and Davis, 1979 Prior to the onset of winter, larvae undergo seasonal metabolic depression that appears to be genetically programmed (Spacht et al., 2020), but whether this metabolic depression is affected by microhabitat conditions has not been assessed.
Populations of the midge are scattered along the Western Antarctic Peninsula and adjacent islands, and distribution within islands is often patchy, with densities ranging from few to tens of thousands of individuals per square meter (Peckham, 1971;Rico and Quesada, 2013;Potts et al., 2020). Larvae are found in a wide array of microhabitats frequently associated with several types of vegetation, including the terrestrial algae Prasiola crispa, several native moss species, and in grassy areas containing the Antarctic hair grass Deschampia antarctica (Tilbrook, 1967;Peckham, 1971;Gantz et al., 2018). Larvae may also be abundant where vegetation is sparse or absent, such as at the peripheries of seal wallows, penguin rookeries, and within former or active nesting sites of several other bird species (Gressitt, 1967;Strong, 1967;Peckham, 1971). These habitats re ect signi cant variation in nutrients (carbon and nitrogen), hygric conditions, and micronutrients, indicating that this species can thrive in a variety of environments (Potts et al., 2020).
While populations of B. antarctica clearly occupy diverse microhabitats, the extent to which this habitat variation in uences life history and physiology has not been assessed. For this study, we selected ve populations of B. antarctica near Palmer Station, Antarctica that occupy notably different microhabitats. We hypothesized that this variation in microhabitat conditions requires larvae to have exible seasonal patterns for growth, development, and metabolism. At these sites, we collected continuous temperature data and conducted weekly measurements of larval mass, metabolic rates, and adult emergence patterns to assess development and phenology. Additionally, we conducted body composition analysis across the season to test the hypothesis that variation in diet and temperature would result in variation in energy storage in larvae, which may affect their ability to prepare for the long Antarctic winter.  Table 1). The four populations on Cormorant Island were within an area of ~ 1 hectare. Sampling was undertaken during the austral summer, from mid-December 2016 through the end of March 2017. At the start of the season, larvae from CorC were collected from under snow and ice each week until the site thawed in early January. The CorD collection site was inaccessible until the rst week of February due to the depth of snow and ice at the site. Larvae were collected from each location within a 2-3 m radius from where the temperature loggers were placed. On each collection day, we also recorded the presence or absence of adults to assess the phenology of adult emergence.

Insects
Adults were recorded as either present and alive, present and dead, or absent. Relative abundance of adults was recorded, though this information is not reported here because sites were not visited frequently enough to accurately quantify patterns of adult abundance, given the short lifespan of adults (~ 7-10 days) and the relatively short, synchronized period of emergence. Larvae were transported to the Palmer Station laboratory and extracted from substrate using a modi ed Berlese setup. Brie y, substrate was placed on metal screens over trays of ice water. Heat lamps placed above the screens were used to drive larvae from the substrate, through the metal screens, and into underlying trays of ice water overnight. The following day, larvae were removed from the ice water for analysis of metabolic rate and body composition.

Field Recording of Temperature
Temperature at each sampling location was recorded every 30 minutes using HOBO U23® temperature/relative humidity data loggers (ONSET® Computer Corporation, Bourne, MA, USA). Data loggers were placed just below the substrate surface at each site to record microhabitat temperatures experienced by the larvae. At the beginning of the season, the data logger at CorC was placed under the snow and ice until the site thawed naturally. For a period of three days, the temperature logger at CorC was accidentally dislodged and exposed to ambient air prior to the second week of collection; these dates were excluded from the analyses. A temperature logger was not placed at CorD until the site was accessible during the rst week of February. As a proxy for heat accumulation across these sites, we calculated cumulative degree days above 0°C.

Metabolic Rate
Larvae collected weekly from each site were held in the laboratory for fewer than 24 h prior to assessing oxygen consumption as a proxy for metabolic rate. Oxygen consumption was measured using closed system respirometry according to Lee and Baust (1982). Brie y, 10 µl microcapillary tubes were glued into tips of 1 ml plastic syringes a xed with metal washers as weights. For each site, 8 replicates of 10 larvae were placed in a small space near the tips of the syringes with moistened cotton and sealed with plungers. Syringes were then partially submerged in a temperature-controlled circulating bath set to 10˚C. Once in the bath, a 10% solution of KOH was placed into the end of the microcapillary tubes, and the respirometers were allowed to equilibrate for approximately 30 minutes prior to taking measurements every 30 minutes over the course of the following hour (for a total of 3 measurements). The KOH solution absorbs CO 2 released by larvae, creating negative pressure in the syringes, thus causing the KOH solution to move along the microcapillary tube as the larvae respire. Distances the KOH plugs travelled were measured, and oxygen consumption rate was calculated from these values. Larvae (still in groups of 10) were removed from the respirometer and dried to a constant mass at 65˚C, and oxygen consumption was calculated per mg dry mass.

Nutrient Analysis
To assess the impact of microhabitat and collection date on biochemical composition of larvae, we measured lipid, carbohydrate, and protein levels. Larvae that were previously dried following respirometry experiments (see above) were stored at -80°C until measurement.
Lipids were quanti ed using vanillin-phosphoric acid reagent, according to Teets et al. (2011) with some modi cations. For each sampling date and location, three groups of two larvae were homogenized in 500 µl 1:1 chloroform:methanol in a bead homogenizer, and tubes were shaken for 15 min at 300 rpm. After centrifugation, the supernatant was transferred to a clean glass test tube and heated at 100°C for 20 min to evaporate the solvent. 100 µl of sulfuric acid was added to the dried extract and heated for 10 min at 90°C, after which 1 ml of vanillin-phosphoric acid reagent (1.2 mg vanillin /ml in 80% phosphoric acid) was added to each tube for color development. After 20 min at room temperature, three aliquots of each sample were added to a 96-well microtiter plate, and absorbance was measured at 525 nm on a microplate spectrophotometer (CLARIOStar, BMG LABTECH, Cary, NC). Absorbance values were compared to an 8-point standard curve containing 0-300 µg vegetable oil.
For carbohydrate measurements, we used the anthrone assay according to Mercer

Statistical Analysis
Analyses were carried out using R Version 3.6.3 (R Core Team, 2020), implemented in RStudio Version 1.2.5042 (RStudio Team, 2020) and methods following Thomas (2017). For each variable (dry mass, metabolic rate, carbohydrate, lipid, and protein), a generalized additive model (GAM) in the mgcv package (Wood, 2019) was used. For these models, the response variable was t as a function of 'week', which was smoothed by site to account for non-linear changes over time. 'Site' was included in the model to test for differences in each variable across site, and 'mass' was included as a covariate (excluding the analysis where 'mass' was the independent variable) (Table S1). In each model, the smoothing parameter estimation method was GCV.Cp (General Cross Validation). For nutrient analyses, 'block' (i.e., the plate used for analysis of each sample) was included as a random effect to account for plate-to-plate variation in the biochemical assays. Using the dataset collected with HOBO dataloggers, mean temperature during the previous week at each site was also included as a covariate in each model. We initially tested four models to determine the best performance for models including temperature: Model 1) temperature included as a smoothed variable, Model 2) temperature included as a smoothed variable, separated by site, Model 3) temperature as a non-smoothed variable and Model 4) temperature as a non-smoothed variable interacting with Site. Our model selection criteria (see details below) indicated that model 3 gave the best model performance, and hence Temperature was included as a sole response variable (Table S1). In addition to models with site-speci c smoothing terms, we also t a GAM with a single

Microhabitat Temperature
Microhabitat temperatures differed considerably across sites (Fig. 3a-f; Table 1). Mean temperature at the warmest and most thermally stable site, CorA, was 5.6 ± 1.6˚C, while temperature at the coldest and most variable site, CorC, was 3.4 ± 2.8˚C during our sampling period. The temperature at CorD, 3.0 ± 1.5˚C, was slightly lower and more consistent than CorC, but we have incomplete data for this site because it was covered with snow and ice during the early portion of the season. However, when we only include the dates for which CorD was accessible, it remained the coldest of the ve sites. Degree days, a measure of accumulated heat units over a period of time, also varied among sites (Table 1). Among sites for which we have a complete record, CorA accumulated the most degree days above 0° (539) and CorC accumulated the fewest (322).

Development: Dry Mass
Larval dry mass, metabolic rate, and adult presence/absence at collecting sites were used as proxies for developmental progression. At all sites, mass tended to increase over time ( Fig. 4a-b), with weight increasing 106% on average over the season. However, changes over time were site-speci c, as determined by AIC-comparison of GAM models with and without site-speci c smoothing terms (DAIC = 334.3). We observed a signi cant negative correlation between temperature in the preceding week and mass (t= -2.75, p < 0.001). Absolute weights also varied by site (p-value for effect of site < 2E-16), with CorC and CorD having the largest larvae, and CorA and CorB having the smallest larvae. Absolute weights also varied by site (p-value for effect of site < 2E-16), with CorC and CorD having the largest larvae and CorA and CorB having the smallest larvae. Complete results for the GAM analysis are presented in Table S2, mass for each site is separately plotted in Figure S1, and prediction plots for dry mass are plotted in Figure S2a.

Development: Metabolic Rate
At all sites, metabolic rate decreased over time, with an average decrease of 66% from the beginning to end of the sampling period ( Fig. 3c-d). Like body size, changes over time were site-speci c, as determined by AIC comparison of GAM models (DAIC = 134.5). We observed no correlation between mean temperature in the preceding week and metabolic rate (t-value = 1.381, p = 0.168). Metabolic rates also varied strongly by site (p-value for effect of site = 6.11E-10), with Humble having a greater metabolic rate than CorA (t = 3.752, p = < 0.001). No signi cant differences were found between CorA and CorB (t = 0.553, p = 0.58), and CorC (t = -0.083, p = 0.193). Mass-speci c metabolic rate was also highly dependent on size, with larger larvae having lower mass-speci c metabolic rates (t-value = -9.058, p < 2.2E-18). The complete results for the GAM analysis are presented in Table S3, metabolic rates for each site are separately plotted in Figure S3, and prediction plots of metabolic rates for each site are separately plotted in Figure S2b.

Development: Adult Presence/Absence
The timing of adult emergence varied among sites (Table 2). Earliest adult emergences were recorded at Hum, CorA, and CorB, and adults were present during each week we sampled in December. Given that these sites had already thawed prior to the start of our sampling, it is likely that we missed the earliest adult emergence events at each of these sites. Adults were observed at CorC for all but one week in January, and there was little or no overlap with the other sites. At CorD, adults were rst observed in the nal week of January and were present every week thereafter until the rst week of March. No adults were observed at any location during the nal three weeks of sampling.
Carbohydrate content varied over time, with most sites decreasing early in the season and increasing at the end of the season (Fig. 4a-b).
We observed no signi cant correlation between mean temperature in the preceding week and carbohydrate content  Table S4, carbohydrate contents for each site are separately plotted in Figure S4, and prediction plots of carbohydrate content for each site are separately plotted in Figure S2c.

Lipid Content
Unlike other variables, lipid content was relatively constant over time ( Fig. 4c-d), with levels uctuating early in the season but converging by the end of the sampling period. Also, changes over time were not site-speci c, as the GAM model with a universal smoothing term was slightly better supported than a complex model with site-speci c smoothing terms DAIC = 2.618). The universal smoothing term also  Table S5, lipid content for each site is separately plotted in Figure  S5, and prediction plots of lipid content for each site and of lipid content per mg fresh mass are separately plotted in Figures S2d and S2e, respectively.

Protein Content
Protein levels varied over time at most sites, with levels increasing early in the season and decreasing at the end of the season (Fig. 4e-f). Changes over time were site-speci c, as determined by comparison of GAM models (DAIC = 53.2). Protein levels also varied strongly by site (F = 7.638, p-value for effect of site = 1.19E-5), with larvae at CorB having the highest protein levels over time. Larvae at CorC and Hum had lower protein levels than larvae at CorA, whilst CorB presented higher protein levels (Table S6). Mass was also a strong predictor of protein content, with larger larvae having lower mass-speci c protein content (t-value = -7.672, p = 1.64E-12). We observed no signi cant correlation between mean temperature in the preceding week and protein content (t-value = -0.696, p = 0.487). The complete results for the GAM analysis are presented in Table S6, protein content for each site is separately plotted in Figure S6, and prediction plots of protein content for each site are separately plotted in Figure S2f.

Discussion
When modeling distributions and climate change responses in small animals like insects, it can be challenging to work at an appropriately small scale that is relevant for the species of interest (Pincebourde and Woods, 2021). In fact, it has long been known that microclimatic and microhabitat variation in uence arthropod physiology, behavior, and development (Uvarov, 1931;Cloudsley-Thompson, 1962;Huey, 1991). Insects often select microhabitats for characteristics such as temperature and water availability (Edney et al., 1974;Harris et al., 2015), thereby maintaining body temperature, limiting desiccation, and maintaining energy reserves when overwintering (Sinclair, 2015). For polar environments like Antarctica, growing seasons are incredibly short, so slight microclimatic variation in temperature can have major impacts on soil-dwelling organisms, especially when temperatures often uctuate around the lower limits for locomotion -B. antarctica larvae become inactive at temperatures between 0 and − 2˚C, and adults don't resume activity until ambient temperature reaches ~ 5˚C (Peckham, 1971;Convey and Block, 1996). Here, we observed considerable variation in thermal conditions among closely associated microhabitats (Fig. 3a-f; Table 1), with a 2.6˚C difference in average temperature between the warmest and coldest sites. Thermal buffering was also different among sites, with some sites such as CorC and Hum showing considerable thermal variation (SD = 2.8˚C), while CorA was much more stable (SD = 1.6˚C). Because the body temperature of small insects tends to closely track surface temperatures (Pincebourde et al., 2021), this variation in microclimate would lead to differences in body temperature, which likely explains some of the observed physiological differences between these populations.
While these results indicate considerable differences in microhabitat conditions across closely situated sites, our experiments do invoke the question of how ne-scale is ne enough? Larvae were sampled within a ~ 2-3 m radius of the logger, but factors such as solar radiation on rocks or evaporative cooling could cause variation within each sampling site. For some other small arthropods, temperature differences upwards of 10°C can be observed over a few cm, for example between sun-exposed rocks and surrounding soil, or between peaks and valleys on the surface of bark (Nicolai 1986;Pike et al., 2012). Thus, we likely failed to completely capture the thermal heterogeneity present in our sampling locations. Indeed, a recent review by Pincebourde and Woods (2020) indicates that the relevant microscale for an insect the size of B. antarctica (~ 5 mm in length) is ~ 20 linear cm, but a limited number of dataloggers and the destructive nature of our sampling prevented us from sampling at this resolution. Despite these limitations, we observed signi cant differences in several parameters across sites, indicating that local environmental conditions are likely having a direct impact on the physiology of larvae.
In general, colder sites had larger larvae with lower metabolic rates, while warmer sites had smaller larvae with higher metabolic rates.
When microhabitat temperature the week prior to sampling was included as a covariate in our models, it was negatively correlated with body size, and this is likely a result of temperature and body size showing opposite trends across the season (i.e., temperature tended to decrease across the season while larvae continued to grow). In contrast, temperature the week prior to measurement did not in uence metabolic rate, which is consistent with our previous work in the laboratory (Spacht et al., 2020). However, in general the warmer sites (i.e., CorA and CorB) tended to have higher metabolic rates than the cooler sites (CorC), suggesting that overall temperature conditions across the season may be in uencing metabolic rate as expected. It is worth noting that all metabolic rate measurements were conducted in the laboratory at a constant temperature (10°C), indicating that larvae from different environments retained different metabolic rates even when measured at a constant laboratory condition. While the hypothesis of metabolic cold adaptation would predict elevated metabolic rate in the populations from cold environments (Addo-Bediako et al., 2002), B. antarctica does not appear to have metabolic cold adaptation or the capacity for compensatory acclimation (Lee and Baust, 1982). Rather, differences in metabolic rate across sites may be primarily driven by differences in size. Much like other animals, smaller body sizes typically correspond with higher metabolic rates in insects (Brown et al., 2004), and our results are consistent with this trend, as body size was negatively correlated with metabolic rate. The extreme variation in body size also indicates asynchronous development rates across sites, and therefore ontogenetic changes in metabolism may be in uencing the variation in metabolic rates observed here. Nonetheless, our results indicate that ne-scale variation in microhabitat conditions can in uence seasonal growth and metabolism.
While there is limited information on thermal thresholds necessary for development of B. antarctica, developmental rate is strongly in uenced by ambient temperature (Harada et al., 2014). Our sites had considerable differences in degree day accumulation above freezing (Table 1), and these differences generally corresponded with the timing of adult emergence ( Table 2), rates of mass gain (Fig. 4a), and metabolic rate (Fig. 4b) throughout the season. Larvae from the two sites with the highest degree days gained the most mass throughout the season (CorA = 66.1%, CorB = 51.6%). Curiously, the lowest mass gain was at Hum (26.8%), even though it had the third highest degree days value, and the coldest sites had surprisingly high mass gain throughout the season (CorC = 46.6%, CorD = 30.1%). These results suggest that larvae may have the capacity for compensatory feeding, a behavior in which animals adjust food intake to meet speci c nutritional needs or to maintain growth rates (Clissold et al., 2013). Sites did not thaw synchronously, with CorA, CorB, and Hum thawing prior to our arrival, and two sites, CorC and CorD, thawing after our arrival, weeks apart. Thus, the available time to feed and develop is not uniform across locations, and in a previous eld season one of our sample sites, CorD, failed to thaw during the summer months. While extended frozen periods could result in larval mortality if energy stores are depleted, it is also possible that some populations of B. antarctica have life cycles extending beyond 2 years, an idea proposed by Harada et al. (2014). Temporal isolation due to microclimate variation could serve as a reproductive barrier among conspeci c populations, as seen in the stone y Leuctra hippopus (Boumans et al., 2017). While we did not evaluate the genetic relatedness of these populations, it is reasonable to hypothesize that physiological divergence, and perhaps speciation, could be the long-term effects of such temporal and spatial isolation (Rundle and Nosil, 2005).
Microhabitat differences also in uenced the biochemical composition of larvae, as re ected by differences in carbohydrate, lipid, and protein content among sites (Fig. 5a-c). Temperature is an important modulator of nutrient uptake and assimilation in insects (reviewed in Clissold and Simpson, 2015) and a number of other ectotherms, including crustaceans (Croll and Watts, 2007), reptiles (Du et al., 2007;Pafalis et al., 2007), and sh (Handeland et al., 2008;Fatma and Ahmed, 2020). However, temperature of the preceding week did not predict macronutrient content of larvae, indicating that week-to-week variation in biochemical composition across sites is not directly related to recent temperature conditions. Instead, differences in biochemical composition may be related to variation in habitat composition, as sites varied in the quantity and types of food available for larvae (Table 1). While some insects can select optimal microhabitats to meet nutritional demands (Clissold et al., 2013), the limited dispersal ability of B. antarctica likely prevents active foraging, so larvae are restricted to food sources in the vicinity of where females oviposit. B. antarctica is a generalist and can feed on a variety of detritus, microorganisms, and plant matter, but there is evidence they preferentially feed on the terrestrial algae Prasiola crispa when available (Peckham, 1971). P. crispa was readily available at only one of our sites, Hum, and either absent or present in limited quantities at the other locations. CorB had no visible vegetation and was also one of the warmest sites, and larvae from this site were the smallest and had the lowest lipid content (after accounting for the effect of body size, which had a strong negative effect on massspeci c lipid content). Thus, the combination of consistently high temperatures and absence of vegetation may have hindered growth and lipid accumulation at this site and could explain why larvae at CorB were smaller than a similarly warm site, CorA, which had lots of decaying moss and grass present. Together, our results suggest that a combination of abiotic (i.e., temperature) and biotic (i.e., vegetation composition) are contributing to the observed physiological differences, although manipulative studies are needed to tease apart these factors.
Investigations of physiological ecology frequently do not account for the in uence of microhabitat on basic aspects of an organism's physiology and life history. Indeed, local adaptation can happen on geographic scales much ner (e.g., tens of kilometers or tens of meters) than typically considered (Richardson et al., 2014), and so it is important to consider ne-scale variation in microenvironmental conditions when assessing ecological factors that in uence an organism's physiology. Given the heterogeneity of terrestrial Antarctic ecosystems and the relative immobility of B. antarctica, we expected to nd diversity in the physiology and phenology among populations. Our results demonstrate that there are considerable temperature differences among the microhabitats occupied by B. antarctica, and these microhabitat differences likely contribute to the observed physiological differences among the populations we sampled. Ongoing analyses investigating genomic variation in these populations will begin to tease apart the relative contributions of genotype and environment in explaining these observed differences across populations. The work presented here supports a growing body of literature on the importance that ne-scale heterogeneity in insect habitats (Pincebourde and Woods, 2021) may have for predicting how this species and other terrestrial organisms in Antarctic will respond to future environmental perturbation.  Temperature of each site over the course of the eld season (a-e), and a violin plot summarizing distribution and frequency of temperatures at each site (f). Closed circles in the violin plot represent mean temperature for the season, and error bars represent standard deviation. The gap in the CorC panel (c) is due to the HOBO logger being inadvertently dislodged and exposed to ambient air for three days before being replaced in the proper location.

Figure 4
Weekly changes in dry mass and metabolic rate for each of the ve microhabitats. Weekly pro les are shown in (a) and (c), while early, mid, and late-season levels are plotted in (b) and (d) to illustrate seasonal trends. In (a) and (c), the shaded area around the lines represents standard error, while in (b) and (d), symbols represent mean ± sem. In (a) and (c), text in the bottom-left summarizes the results of GAM analyses test the effects of site and time on each variable. For "site", asterisks indicate a signi cant effect of site, with * for p<0.01, ** for p<0.01, and *** for p<0.0001. For "time," we t site-speci c smoothing terms, and the asterisks indicates that at least one of the sites reached the indicated signi cance level. The "site-speci c changes over time?" indicates the results of model-selection analysis where a GAM with site-speci c smoothing terms was compared to a GAM with a universal smoothing term. "Yes" indicates that the model with site-speci c smoothing terms had stronger support, while "No" indicates better support for a model with a single smoothing term. For "site", asterisks indicate a signi cant effect of site, with * for p<0.01, ** for p<0.01, and *** for p<0.0001. For "time," we t site-speci c smoothing terms, and the asterisks indicates that at least one of the sites reached the indicated signi cance level. The "site-speci c changes over time?" indicates the results of model-selection analysis where a GAM with site-speci c smoothing terms was compared to a GAM with a universal smoothing term. "Yes" indicates that the model with site-speci c smoothing terms had stronger support, while "No" indicates better support for a model with a single smoothing term.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download.