Please cite this article as: Narancic, B., Wolfe, B.B., Pienitz, R., Meyer, H., Lamhonwah, D., Landscape-gradient assessment of thermokarst lake hydrology using water isotope tracers, Journal of Hydrology (2016), doi: http:// dx.doi.org/10. 1016/j.jhydrol.2016.11.028 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.1 Landscape-gradient assessment of thermokarst lake hydrology using water isotope -Water isotopes are used to characterize thermokarst lake hydrology in Nunavik, Canada.-Rainfall and/or permafrost meltwater are the main lake water sources.-Maritime climate enhances the hydrological resiliency of thermokarst lakes.-Terrestrial carbon inputs from surface runoff are likely to increase in the future.-Thermokarst lakes will persist and methane emissions are likely to increase.3 ABSTRACT Thermokarst lakes are widespread in arctic and subarctic regions. In subarctic Québec (Nunavik), they have grown in number and size since the mid-20 th century.Recent studies have identified that these lakes are important sources of greenhouse gases. This is mainly due to the supply of catchment-derived dissolved organic carbon that generates anoxic conditions leading to methane production. To assess the potential role of climate-driven changes in hydrological processes to influence greenhouse-gas emissions, we utilized water isotope tracers to characterize the water balance of thermokarst lakes in Nunavik during three consecutive mid-to late summer seasons (2012)(2013)(2014). Lake distribution stretches from shrub-tundra overlying discontinuous permafrost in the north to spruce-lichen woodland with sporadic permafrost in the south.Calculation of lake-specific input water isotope compositions ( I ) and lake-specific evaporation-to-inflow (E/I) ratios based on an isotope-mass balance model reveal a narrow hydrological gradient regardless of diversity in regional landscape characteristics. Nearly all lakes sampled were predominantly fed by rainfall and/or permafrost meltwater, which suppressed the effects of evaporative loss. Only a few lakes in one of the southern sampling locations, which overly highly degraded sporadic permafrost terrain, appear to be susceptible to evaporative lake-level drawdown. We attribute this lake hydrological resiliency to the strong maritime climate in coastal regions of Nunavik. Predicted climate-driven increases in precipitation and permafrost degradation will likely contribute to persistence and expansion of thermokarst lakes throughout the region. If coupled with an increase in terrestrial carbon inputs to thermokarst lakes from surface runoff, conditions favorable for mineralizatio...