Improvements and Evaluation of the FLake Model in Dagze Co, Central Tibetan Plateau

: FLake has been one of the most extensively used lake models in many studies for lake thermal structure simulations. However, due to the scarcity of lake temperature observations, its applicability and performance on lakes over the Tibetan Plateau are still poorly investigated, especially in small-to medium-sized lakes. In this study, based on water proﬁle observations in Dagze Co, a medium-sized lake on the central Tibetan Plateau, the sensitivity of lake thermal features to three key parameters in FLake was investigated. The performance of FLake in reproducing the lake thermal features was evaluated and improved by optimizing these key parameters. The results showed that the FLake model with default parameter settings can generally reproduce the thermal features of Dagze Co, but there are still signiﬁcant deviations compared to observation. The sensitive experiments demonstrated that the thermal structure of the lake obviously responds to the change in the water extinction coefﬁcient ( K d ), friction velocity (u*), and ice albedo ( α ice ). Based on previous studies and sensitive experiments, the three key parameters were set to the optimized value, which substantially improved the performance of FLake. The values of bias and RMSE of simulated lake surface water temperature decreased from 3.08 ◦ C and 3.62 ◦ C to 2.0 ◦ C and 2.48 ◦ C after parameter optimization. The integration of a simple salinity scheme further improved the ability of FLake to reproduce the observed thermal features of Dagze Co. These results will improve our understanding of thermal processes in lakes on the Tibetan Plateau, as well as the applicability of lake models.


Introduction
There are more than 1400 lakes with a surface area larger than 1 km 2 in the Tibetan Plateau (TP), contributing to over half of the lake area in China with a total area of 5 × 10 4 km 2 in 2018 [1,2]. Lakes are regarded as sentinels of climate change, as water temperature, surface areas, water levels, clarity, physical and biogeochemical properties, mixing regimes, and the ice phenology of lakes are very sensitive and rapidly respond to climate change [3][4][5][6][7][8][9]. One of the most important and direct responses of a lake to climate warming is the rising lake surface water temperature (LSWT), which may further affect the internal thermal processes of the lake by changing the physical characteristics of the lake [10][11][12].
Characterized by lower albedo, smaller roughness, larger heat capacity, and primary moisture sources, lakes can significantly affect the weather and climate at local to regional scales by altering the surface mass and energy exchange [13][14][15][16]. In the TP region, with an average altitude of 4000 m above mean sea level, the climatic environment is extraordinarily Dagze Co (31 • 49 -31 • 59 N, 87 • 25 -87 • 39 E, 4450 m a.m.s.l.) is a brackish lake located in the central TP ( Figure 1). The lake has a surface area of 245 km 2 , which is mainly fed by precipitation and the Bogcarg Zangbo River. The maximum depth of the lake is 38 m, roughly divided into several layers in terms of salinity: a layer from the lake surface to approximately 24 m with a slowly increasing salinity from 14.7 g L −1 to about 16.4 g L −1 , a halocline at 25-29 m with a rapidly increasing salinity from 16.4 g L −1 to about 21 g L −1 , and a layer with almost unchanged salinity of 21.4 g L −1 from 30 m to the bottom of the lake [54]. The thermocline appears at a depth of 16-23 m in summer (observed on 18 August 2012), while the ice-covered period lasts from November to April. The mean annual precipitation and air temperature observed by XainZa Meteorological Station, which is 150 km southeast of the lake, were 316 mm and 0.55 • C, respectively [55]. A multiyear averaged air pressure of 553 hPa and wind speed of 4.9 m s −1 were obtained by Shuanghu station 180 km away in the northeast direction of the lake.
Water 2023, 15,3135 3 of 20 to three key model parameters was investigated by conducting a series of sensitivity experiments. Based on the sensitivity experiments and observations, we further optimized the performance of the FLake model in experiments with tuned key parameters. This study aimed to improve our understanding of the suitability of FLake for lakes on the TP.

Study Area
Dagze Co (31°49′-31°59′ N, 87°25′-87°39′ E, 4450 m a.m.s.l.) is a brackish lake located in the central TP ( Figure 1). The lake has a surface area of 245 km 2 , which is mainly fed by precipitation and the Bogcarg Zangbo River. The maximum depth of the lake is 38 m, roughly divided into several layers in terms of salinity: a layer from the lake surface to approximately 24 m with a slowly increasing salinity from 14.7 g L −1 to about 16.4 g L −1 , a halocline at 25-29 m with a rapidly increasing salinity from 16.4 g L −1 to about 21 g L −1 , and a layer with almost unchanged salinity of 21.4 g L −1 from 30 m to the bottom of the lake [54]. The thermocline appears at a depth of 16-23 m in summer (observed on 18 August 2012), while the ice-covered period lasts from November to April. The mean annual precipitation and air temperature observed by XainZa Meteorological Station, which is 150 km southeast of the lake, were 316 mm and 0.55 °C, respectively [55]. A multiyear averaged air pressure of 553 hPa and wind speed of 4.9 m s −1 were obtained by Shuanghu station 180 km away in the northeast direction of the lake.

GloboLakes Surface Water Temperature
The Global Observatory of Lake Responses to Environmental Changes (GloboLakes) data product v4.0 was employed in this study to verify the simulated lake surface water temperature (LSWT). The GloboLakes provides daily values of LSWT for 2000 lakes globally distributed from 1995 to 2019 with a horizontal spatial resolution of 0.05° × 0.05°. The temperatures were retrieved from different orbit instruments, including the Along Track Scanning Radiometer (ATSR), Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Sea and Land Surface Temperature Radiometer (SLSTR), with the same algorithm and harmonized between sensors using overlap periods to ensure consistency [56]. The data represent the current stateof-the-art for LSWT record production and perform well in comparison with in situ measurements [57].

GloboLakes Surface Water Temperature
The Global Observatory of Lake Responses to Environmental Changes (GloboLakes) data product v4.0 was employed in this study to verify the simulated lake surface water temperature (LSWT). The GloboLakes provides daily values of LSWT for 2000 lakes globally distributed from 1995 to 2019 with a horizontal spatial resolution of 0.05 • × 0.05 • . The temperatures were retrieved from different orbit instruments, including the Along Track Scanning Radiometer (ATSR), Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS), and Sea and Land Surface Temperature Radiometer (SLSTR), with the same algorithm and harmonized between sensors using overlap periods to ensure consistency [56]. The data represent the current state-of-the-art for LSWT record production and perform well in comparison with in situ measurements [57] There were 18 temperature loggers deployed along the water column, with  1-m intervals above 11 m depth, 2-m intervals from 11 to 23 m, and 3-m intervals below  23 m. The loggers commenced operation at 00:00 on 19 August 2012 and were retrieved on 20 August 2013, during which the water temperature was recorded at one-hourly intervals. The uppermost logger was placed at a~4-m depth to prevent the influence of ice freezing and melting [54].

In situ Water Clarity Dataset
The water clarity data of Dagze Co was extracted from the in situ water quality parameters of the lakes on the Tibetan Plateau (2009-2020), which provides the in situ lake water parameters of 124 lakes in the TP [58]. These in situ parameters include water temperature, salinity, pH, chlorophyll-a concentration, blue-green algae concentration, turbidity, dissolved oxygen, fluorescent dissolved organic matter, and water clarity measured by Secchi Disk Depth (SDD).

China Meteorological Forcing Dataset (CMFD)
The forcing data of FLake were extracted from the China Meteorological Forcing Dataset (CMFD), which is a high spatial-temporal resolution gridded meteorological dataset for land surface processes simulation in China. The dataset is derived from a fusion of remote sensing products (GEWEX-SRB and TRMM), a reanalysis dataset (Princeton and GLDAS), and CMA in situ observation data, and provides 7 near-surface meteorological elements, including surface pressure, 2 m air temperature and specific humidity, 10 m wind speed, downward shortwave radiation, downward longwave radiation and precipitation rate with a temporal resolution of 3 h, and a spatial resolution of 0.1 • from 1979 (currently up to 2018) [59]. It has become one of the most widely used datasets in China due to its better accuracy than the internationally available reanalysis data [18,50,52].

FLake Description and Configuration
FLake is a one-dimensional (1D) lake model capable of predicting the vertical temperature and mixing conditions in lakes of various depths on time scales from a few hours to many years [45,60]. FLake separates the water column into two layers, an upper mixed layer and an underlying thermocline layer, to parametrically represent the evolving temperature profile. The water temperature in the upper mixed layer is assumed to be vertically uniform. The thermocline layer between the upper mixed layer and the lake bottom is described based on the concept of self-similarity of the temperature-depth curve. The temperature profile of the thermocline can be fairly accurately parameterized through the following function of dimensionless depth: where t is time (s), z is the depth (m), θ s (t) is the upper mixed layer water temperature (K), h(t) is the mixed layer depth, ∆θ(t) = θ s (t) − θ b (t) is the temperature difference across the thermocline of depth ∆h(t), and θ b (t) is the temperature at the bottom of the thermocline. Φ θ (ζ) is a dimensionless "universal" function of the relative depth ζ = [z − h(t)]/∆h(t), which satisfies two boundary conditions, Φ θ (0) = 0 and Φ θ (1) = 1. It follows that the two-layer parameterization of the vertical temperature profile at time t is given by: where the dimensionless shape function Φ θ (ζ) is calculated by: where C θ is the shape factor and can be computed by: where t rc is the relaxation time scale (s), which indicates the time of the temperature profile in the thermocline evolving from one curve to the other, following the change of the sign in dt . The minimum and maximum values of the shape factor are empirically set to C min θ = 0.5 and C max θ = 0.8. The concept described above is also applied to describe the thermal structure of one ice layer, one snow layer on the lake ice, and two thermally active layers of lake sediments, which are included in the FLake.
The ice albedo in FLake is assumed to be represented by an empirical formulation based on the surface temperature: where α max is the albedo of white ice (set to 0.6 by default), α min is the albedo of blue ice (set to 0.1 by default), T f is the freezing temperature (K), and T s is the ice surface temperature (K). The solar radiation transfers in water and ice/snow are described by a one-band exponential approximation of the Beer-Lambert decay law, for which an extinction coefficient of 3 m −1 in water and 1.0 × 10 −7 in ice/snow is used in the default setting. Lake momentum fluxes, sensible heat fluxes, and latent heat fluxes were estimated using a parametric scheme that integrates the relationships between roughness length, potential temperature, specific humidity, and wind speed. The function used to calculate friction velocity is: where u(z) is the wind velocity at height z, u s is the wind velocity at the lake surface, u * is the friction velocity, k is the von Karman constant, z 0u is roughness length of wind velocity, and Ψ u is Monin-Obukhov dimensionless functions in the similarity theory. FLake intends to be coupled to different numerical weather prediction models (NWPs) and regional climate models (RCMs) as a lake parameterization scheme [23,28,61,62], but it can also be used as a stand-alone lake model. FLake has been well-tested by many studies worldwide, including a few lakes in the TP, and showed comparable or even better performance to many other 1D lake models in reproducing lake water temperature [39,43,46,63]. The simple two-layer stratification gives FLake the prominent merits of being computationally efficient and a small number of parameters have to be specified [18], but also limits its performance for lakes deeper than 50 m due to the absence of the hypolimnion layer between the thermocline and the lake bottom in deep lakes. However, this fact will not have much impact in our case as the water depth of Dagze Co is less than 50 m.

Experimental Design
The FLake model was driven by CMFD from 20 April 1990 to 31 December 2018 with a time step of 10 min. The running from 1990 to 1994 was for spin-up, while the result from 1995 to 2018 was used for analysis to avoid uncertainty of results in a single year. The simulation started on the ice-off date to simplify the initial condition. The lake depth of Dagze Co, which is one of the few parameters that need to be specified in FLake, was set to 37 m, equal to the water depth at the temperature profile monitoring site.
To evaluate the performance of the FLake model with different parameter settings in reproducing the LSWT and thermal structure of Dagze Co, a control experiment (CTRL, Table 1) with the default model configuration was first carried out. In the CTRL experiment, the extinction coefficient is set to 3 m −1 , the friction velocity was calculated by Equation (6), and the ice-albedo was obtained from Equation (5). In addition to the CTRL experiment, a series of sensitive experiments (SenExps in Table 1) were conducted by tuning some key parameters, i.e., extinction coefficient (K d ), friction velocity (u*), and ice albedo (α ice ), to explore the sensitivity of simulated water temperature to these key parameters. The extinction coefficient (K d ), which indicates the attenuation of solar radiation within the water column, is one of the most important parameters in controlling the water temperature profile. To investigate how the K d affects water temperature, we performed a series of sensitivity experiments (SenExp_K d ), gradually changing K d from 0.15 to 0.5 with a 0.05 interval, and from 0.6 to 2.8 with a 0.2 interval. The friction velocity (u*) in FLake, according to a previous study, also has a significant impact on water temperature [64]. Hence, experiments (SenExp_u*) with u* multiplied by a series of scale factors (varies from 1.2 to 2.0 with a 0.2 interval) were conducted to analyze the sensitivity of water temperature to u*. Previous studies have shown that the white ice albedo (α max ) is a key parameter that determines lake ice phenology [65,66]. Thus, a series of sensitive experiments with changed white ice albedo (from 0.2 to 0.6 with a 0.05 interval) were conducted to explore its influence on ice phenology and water temperature.
To further improve the performance of FLake, three experiments were conducted by gradually optimizing each parameter mentioned above (OptExps in Table 1). According to previous studies, the K d is closely correlated with water clarity (SDD) and can be approximated by the empirical equation [67]: The K d in the CTRL experiment (3 m −1 , corresponding to~0.5 m SDD) seems too high for Dagze Co, which has a mean in situ SDD of~7 m. Hence, the K d was tuned to~0.37 m −1 according to Equation (7) in the parameter optimization experiments (OptExp_1). Then, the friction velocity was also tuned by multiplying an optimized scale factor to minimize the deviation of the simulated LSWT (OptExp_2). Based on the sensitive experiments, the white ice albedo (α max ) in FLake has generally been overestimated by many lake models before in TP lakes, which are barely covered by any snow in winter, and was also tuned to an optimized value to further minimize the deviation of the reproduced LSWT. The tuning process of these parameters was both efficient and physically reasonable for the representation of lake properties. In addition to the above parameter tuning process, a salinity parameterization scheme was added to FLake to consider the effect of salinity (S) on the maximum density temperature (T m ) and water freezing point (T f ). The equations are as follows [68]: where S is the salinity (g L −1 ) of the lake water.

Methodology
The statistics including the mean bias error (BIAS), correlation (R), and root-meansquare-error (RMSE) are adopted to evaluate and improve the model performance. The formulas for these statistics are shown as follows: where N is the number of points (in time or space).
is the mean value of the simulation (observation) over N sample points. Both lower BIAS and RMSE suggest a better performance of the simulation, while higher R indicates a larger similarity in the temporal variation between the simulation and observation.

Lake Surface Water Temperature
To reveal the ability of the FLake model to simulate the structure of the Dagze Co with default parameter settings, we first compare the LSWT from the CTRL experiment with the observation. In the same way as observations [54], the CTRL experiment shows that Dagze Co is frozen from late November or early December to April each year. During the ice-covered period, the simulated lake surface temperature sharply fluctuates because the lake surface is frozen, which reduces the specific heat capacity of the lake surface and blocks the heat exchange with the deep lake water, so the lake surface temperature is greatly influenced by the air temperature. Figure 2 exhibits that although the seasonal variation of LSWT and ice on/off characteristics of Dagcze Co can be roughly simulated by FLake, there is still a significant deviation between simulation and observation. In particular, the LSWT is underestimated in winter (December-February) with a large systematic negative bias, which may be caused by an overestimation of FLake in lake ice albedo for alpine lakes on the TP [69]. While the LSWT is generally overestimated in summer (June-August) with a large systematic positive bias.

Vertical Temperature Profile
In Figure 3, the time-depth distribution of the daily water temperature of the Dagze Co from observation and simulation is shown. According to the observation (Figure 3a), the water temperature during the ice-covered period (December-April) is generally under 2 • C in the upper layer of the lake (depth of 0-20 m), and gradually increases to~4 • C at the bottom of the lake, except that after March the water temperature near the surface gradually increased to~4-7 • C. However, after the ice break-up, the temperature near the surface layer suddenly dropped to 2-3 • C during a short period in late April and then began to increase from May to August. The temperature in the epilimnion of the lake

Vertical Temperature Profile
In Figure 3, the time-depth distribution of the daily water temperature of the Dagze Co from observation and simulation is shown. According to the observation (Figure 3a), the water temperature during the ice-covered period (December-April) is generally under 2 °C in the upper layer of the lake (depth of 0-20 m), and gradually increases to ~4 °C at the bottom of the lake, except that after March the water temperature near the surface gradually increased to ~4-7 °C. However, after the ice break-up, the temperature near the surface layer suddenly dropped to 2-3 °C during a short period in late April and then began to increase from May to August. The temperature in the epilimnion of the lake eventually exceeded 10 °C in August, while the temperature in the hypolimnion below 20 m was essentially maintained at about 4 °C.   Figure 3b shows the overall vertical thermal structure of the lake reproduced by the FLake in the CTRL experiment. Compared to the observation, the simulated water temperature near the surface layer (about 0-5 m) during the ice-covered period is also lower than 2 °C, but there is a significant rise in temperature from 5-20 m, gradually increasing

Vertical Temperature Profile
In Figure 3, the time-depth distribution of the daily water temperature of the Dagze Co from observation and simulation is shown. According to the observation (Figure 3a), the water temperature during the ice-covered period (December-April) is generally under 2 °C in the upper layer of the lake (depth of 0-20 m), and gradually increases to ~4 °C at the bottom of the lake, except that after March the water temperature near the surface gradually increased to ~4-7 °C. However, after the ice break-up, the temperature near the surface layer suddenly dropped to 2-3 °C during a short period in late April and then began to increase from May to August. The temperature in the epilimnion of the lake eventually exceeded 10 °C in August, while the temperature in the hypolimnion below 20 m was essentially maintained at about 4 °C.   Figure 3b shows the overall vertical thermal structure of the lake reproduced by the FLake in the CTRL experiment. Compared to the observation, the simulated water temperature near the surface layer (about 0-5 m) during the ice-covered period is also lower than 2 °C, but there is a significant rise in temperature from 5-20 m, gradually increasing  Figure 3b shows the overall vertical thermal structure of the lake reproduced by the FLake in the CTRL experiment. Compared to the observation, the simulated water temperature near the surface layer (about 0-5 m) during the ice-covered period is also lower than 2 • C, but there is a significant rise in temperature from 5-20 m, gradually increasing to 4 • C, as well as water temperature for the bottom layer of the lake. From June to early August, the temperature in the upper layer of the lake exceeded 10 • C, especially for the epilimnion with a temperature near 20 • C. An obvious temperature gradient appears in metalimnion with a depth ranging from 5 m to 20 m, where the water temperature dropped to about 8 • C. For the hypolimnion with a depth of 20-37 m, the water temperatures were higher compared to the observation.
In conclusion, the FLake model with default parameter settings can reproduce the overall structure of the observed water temperature, with a relatively smaller vertical temperature gradient during the ice-covered period and larger vertical temperature gradients during the stratified season. However, there are still obvious deviations in the simulation compared to the observation. FLake produced a shorter ice-covered period, mainly due to a later ice-on date and an earlier ice-off date. Additionally, the water temperature in the mixed layer during the stratified season was overestimated, with a relatively shallower mixed layer depth from September to October 2012.

Extinction Coefficient
To improve the applicability of the FLake model in Dagze Co, it is necessary to systematically investigate the sensitivity of some key parameters in the FLake on lake thermal features first. By setting the water extinction coefficient (K d ) of Dagze Co to different values, the changes in simulated water temperature for both the surface and interior of the lake with different transparency conditions were investigated. According to the observation, the Secchi Disk Depth (SDD) for Dagze Co is~7 m, corresponding to a K d of~0.37 m −1 , which is much smaller than the default value of 3 m −1 in FLake. Therefore, the extinction coefficient was gradually changed from 2.8 m −1 to 0.6 m −1 with an interval of 0.2 m −1 , and from 0.5 m −1 to 0.15 m −1 with an interval of 0.05 m −1 . As shown in Figure 4 (only some of the results are shown), the curves of LSWT, mean water column temperature (MWCT), lake bottom water temperature (LBWT), and mixed layer depth (MLD) apparently changed with the extinction coefficient. For turbid lake water with larger K d , solar radiation is absorbed by the particles in the water, leading to a small depth penetrated by solar radiation. Thus, the heating of solar radiation mainly operates in a shallow layer of the lake. As shown in Figure 4, larger K d corresponds to a warmer LSWT, colder LBWT, and smaller MLD, while lower K d produces the opposite result. The higher LSWT of the turbid lake leads to an increase in the latent and sensible heat fluxes at the lake surface with accelerated heat exchange between the lake and the atmosphere, which makes the LSWT respond more quickly to the solar radiation and air temperature. This explains the fluctuations in the LSWT with higher K d and smoother LSWT with lower K d , as well as the fact that the LSWT of the higher K d is warmer in the spring and colder in the autumn, as well as lower MWCT, which is just the opposite for the result with small K d .     Figure 5 illustrates the sensitivity of LSWT to K d , representing the variation in LSWT corresponding to a 1 m −1 change in K d . Generally, LSWT is not very sensitive to changes in K d when the K d is greater than 0.5 m −1 , and there is no evident correlation between changes in LSWT and K d . However, once K d is less than 0.5 m −1 , LSWT becomes quite sensitive to changes in K d with pretty obvious patterns between the two of them. That is, LSWT generally decreases with decreasing K d from May to September, while increases with decreasing K d from October to April occur in the following year. With a decreasing K d , the most pronounced cooling in LSWT occurred in June and July, especially when the K d is equal to 0.15, in which case each 1 m −1 reduction in K d will result in a maximum decrease of 16.22 • C for LSWT. In contrast, LSWT increased the most in January and December as K d reduced, which also reaches the maximum value of 15.15 • C when K d is equal to 0.15 m −1 . Figure 5 illustrates the sensitivity of LSWT to Kd, representing the variation in LSWT corresponding to a 1 m −1 change in Kd. Generally, LSWT is not very sensitive to changes in Kd when the Kd is greater than 0.5 m −1 , and there is no evident correlation between changes in LSWT and Kd. However, once Kd is less than 0.5 m −1 , LSWT becomes quite sensitive to changes in Kd with pretty obvious patterns between the two of them. That is, LSWT generally decreases with decreasing Kd from May to September, while increases with decreasing Kd from October to April occur in the following year. With a decreasing Kd, the most pronounced cooling in LSWT occurred in June and July, especially when the Kd is equal to 0.15, in which case each 1 m −1 reduction in Kd will result in a maximum decrease of 16.22 °C for LSWT. In contrast, LSWT increased the most in January and December as Kd reduced, which also reaches the maximum value of 15.15 °C when Kd is equal to 0.15 m −1 .

Friction Velocity
A recent study in Ngoring Lake found that the friction velocity is underestimated by the FLake compared to the observed value, resulting in a significant positive deviation in LSWT produced by FLake [64]. Hence, a series of sensitivity experiments with the friction velocity u* increased by a factor of 1.2, 1.4, 1.6, 1.8, and 2 were conducted to investigate the influence of different friction velocities on the simulated LSWT in Dagze Co. As Figure 6 shows, compared with the CTRL experiment with default parameter setting, the simulated LSWT, MWCT, LBWT, and MLD apparently changed with different friction velocities. As the friction velocity increases, LSWT generally decreases and fluctuates more dramatically. The ice-on date is advanced while the ice-off date is delayed, resulting in a longer ice-covered period. The effect of friction velocity on lake thermal structure is not as pronounced as K d , but still considerable. In addition to a decrease in the LSWT, the MWCT of the lake also slightly decreased with increasing friction velocity (Figure 6b). Although the LSWT and MWCT decreased, there was no significant change in the LBWT (Figure 6c). Variation in the MLD demonstrated that the larger the friction velocity, the deeper the MLD, as well as more obvious fluctuations. The reason for the decrease in LSWT and MCWT is that the increase in friction velocity increases the air-lake heat exchange and vertical mixing of the lake.

Ice Albedo
The lake ice albedo can affect the amount of solar radiation absorbed by the ice during the ice-covered period of the lake, further affecting the ice-off date. The climate on TP is dry and with strong solar radiation in winter, resulting in sparse snow over the lake ice and thus a generally lower lake surface albedo compared to low-lying regions. Observations indicate that the actual ice albedo for lakes on the TP is much smaller than the value used in the lake models [66]. Therefore, to verify the effect of the ice albedo on simulated lake thermal conditions, the white ice albedo in FLake is gradually reduced from the default value of 0.6 to 0.2 with an interval of 0.05. As shown in Figure 7, during the ice-covered period, LSWT slightly increases with the decrease of lake ice albedo, while the variation is not obvious in other months. As the ice albedo decreases, the simulated ice-covered period generally shortens, as well as an increase in LSWT just before ice-on and after ice-off. This is because the smaller the ice albedo is, the more solar radiation the ice absorbs, and thus the faster the ice melts. Lake ice can prevent solar radiation from penetrating the lake water, and an early ice-off will lead to more absorbed solar radiation and rapidly rising LSWT after the ice-off. The absorbed energy will be stored in the lake and released until the cooling stage of the lake. In such a case, the lake will maintain a higher temperature before the ice-on, resulting in a delayed ice-on date. Therefore, the lake ice albedo can not only affect the duration of the ice-covered period but also the LSWT before and after the ice-covered period. The ice albedo can also affect the internal thermal feature of the lake, the MWCT slightly increases from July to November as the ice albedo decreases, but not significantly in other months. The change in albedo has no obvious effect on LBWT, which slightly increases from June to December with a decrease in albedo. The MLD generally shows a slightly deepening trend, especially in the turnover period. In general, the decrease in ice albedo not only affects the ice-covered duration of the lake but also indirectly affects the solar radiation absorbed by the lake and the lake-air energy exchange, further affecting the surface and internal thermal features of the lake.

Friction Velocity
A recent study in Ngoring Lake found that the friction velocity is underestimated by the FLake compared to the observed value, resulting in a significant positive deviation in LSWT produced by FLake [64]. Hence, a series of sensitivity experiments with the friction velocity u* increased by a factor of 1.2, 1.4, 1.6, 1.8, and 2 were conducted to investigate the influence of different friction velocities on the simulated LSWT in Dagze Co. As Figure  6 shows, compared with the CTRL experiment with default parameter setting, the simulated LSWT, MWCT, LBWT, and MLD apparently changed with different friction velocities. As the friction velocity increases, LSWT generally decreases and fluctuates more dramatically. The ice-on date is advanced while the ice-off date is delayed, resulting in a longer ice-covered period. The effect of friction velocity on lake thermal structure is not as pronounced as Kd, but still considerable. In addition to a decrease in the LSWT, the MWCT of the lake also slightly decreased with increasing friction velocity (Figure 6b). Although the LSWT and MWCT decreased, there was no significant change in the LBWT ( Figure  6c). Variation in the MLD demonstrated that the larger the friction velocity, the deeper the MLD, as well as more obvious fluctuations. The reason for the decrease in LSWT and MCWT is that the increase in friction velocity increases the air-lake heat exchange and vertical mixing of the lake.

Ice Albedo
The lake ice albedo can affect the amount of solar radiation absorbed by the ice during the ice-covered period of the lake, further affecting the ice-off date. The climate on TP is dry and with strong solar radiation in winter, resulting in sparse snow over the lake ice and thus a generally lower lake surface albedo compared to low-lying regions. Observations indicate that the actual ice albedo for lakes on the TP is much smaller than the value used in the lake models [66]. Therefore, to verify the effect of the ice albedo on simulated lake thermal conditions, the white ice albedo in FLake is gradually reduced from the default value of 0.6 to 0.2 with an interval of 0.05. As shown in Figure 7, during the icecovered period, LSWT slightly increases with the decrease of lake ice albedo, while the variation is not obvious in other months. As the ice albedo decreases, the simulated icecovered period generally shortens, as well as an increase in LSWT just before ice-on and after ice-off. This is because the smaller the ice albedo is, the more solar radiation the ice

Optimization of Key Parameters
Based on the previous studies and results from sensitive experiments conducted above, three sets of experiments with optimized parameters were conducted to gradually improve the performance of the FLake model for Dagze Co (OptExp in Table 1). The CTRL experiment with default parameters settings was used as a reference. In the OptExp_1 experiment, the K d was modified from 3.0 m −1 to 0.37 m −1 , which is calculated by Equation (7) with an SDD of~7 m in Dagze Co. In addition to optimized K d , the OptExp_2 experiment enlarged the friction velocity (u*) with a factor of 2. In the OptExp_3 experiment, the white ice albedo was further changed from 0.6 to 0.25. Through the above experiments, three key parameters were gradually optimized to enhance the performance of the FLake model. of the lake, the MWCT slightly increases from July to November as the ice albedo decreases, but not significantly in other months. The change in albedo has no obvious effect on LBWT, which slightly increases from June to December with a decrease in albedo. The MLD generally shows a slightly deepening trend, especially in the turnover period. In general, the decrease in ice albedo not only affects the ice-covered duration of the lake but also indirectly affects the solar radiation absorbed by the lake and the lake-air energy exchange, further affecting the surface and internal thermal features of the lake.

Optimization of Key Parameters
Based on the previous studies and results from sensitive experiments conducted above, three sets of experiments with optimized parameters were conducted to gradually improve the performance of the FLake model for Dagze Co (OptExp in Table 1). The CTRL experiment with default parameters settings was used as a reference. In the OptExp_1 experiment, the Kd was modified from 3.0 m −1 to 0.37 m −1 , which is calculated by Equation (7) with an SDD of ~7 m in Dagze Co. In addition to optimized Kd, the OptExp_2 experiment enlarged the friction velocity (u*) with a factor of 2. In the OptExp_3 experiment, the white ice albedo was further changed from 0.6 to 0.25. Through the above experiments, three key parameters were gradually optimized to enhance the performance of the FLake model. Figure 8 shows the LSWT of Dagze Co from GloboLakes observation and FLake simulation in different experiments with optimized parameters. In the CTRL experiment, the FLake overestimated LSWT in summer and autumn, and the lake warmed more rapidly in spring and also cooled more rapidly in autumn compared to observation. After the Kd was changed to 0.35 m −1 in the OptExp_1, although the overestimation of LSWT still exists, the seasonal variation in LSWT was more consistent with observations. In the OptExp_2,  Figure 8 shows the LSWT of Dagze Co from GloboLakes observation and FLake simulation in different experiments with optimized parameters. In the CTRL experiment, the FLake overestimated LSWT in summer and autumn, and the lake warmed more rapidly in spring and also cooled more rapidly in autumn compared to observation. After the K d was changed to 0.35 m −1 in the OptExp_1, although the overestimation of LSWT still exists, the seasonal variation in LSWT was more consistent with observations. In the OptExp_2, with friction velocity enlarged by a factor of 2, the overestimation of LSWT in summer and autumn was significantly reduced, but the ice-off date is still later than the observation. Hence, based on the above experiments, the white ice albedo in FLake was further changed to 0.25, resulting in a significant improvement in the LSWT simulation compared to the observation. with friction velocity enlarged by a factor of 2, the overestimation of LSWT in summer and autumn was significantly reduced, but the ice-off date is still later than the observation. Hence, based on the above experiments, the white ice albedo in FLake was further changed to 0.25, resulting in a significant improvement in the LSWT simulation compared to the observation. As the statistics results show in Table 2, compared with observations, the CTRL experiment has a bias of 3.08 °C, an RMSE of 3.62 °C, and a correlation coefficient of 0.91. After changing the Kd in OptExp_1, although the bias, RMSE, and correlation coefficient of simulated LSWT have almost not been reduced or are even slightly worse, the phase of annual variation of LSWT is more consistent with observation. With the friction velocity tuned, the bias and RMSE are significantly reduced to 2.20 °C and 2.72 °C, with a correlation coefficient of 0.86. In OptExp_3, following optimized Kd and friction velocity, the white ice albedo was also modified to 0.25, the bias and RMSE of LSWT further reduced to 2.00 °C and 2.48 °C, respectively, while the correlation coefficient was re-raised to 0.90. This result indicates that the capability of FLake in reproducing the LSWT of Dagze Co was substantially improved by the optimization of key parameters. As the statistics results show in Table 2, compared with observations, the CTRL experiment has a bias of 3.08 • C, an RMSE of 3.62 • C, and a correlation coefficient of 0.91. After changing the K d in OptExp_1, although the bias, RMSE, and correlation coefficient of simulated LSWT have almost not been reduced or are even slightly worse, the phase of annual variation of LSWT is more consistent with observation. With the friction velocity tuned, the bias and RMSE are significantly reduced to 2.20 • C and 2.72 • C, with a correlation coefficient of 0.86. In OptExp_3, following optimized K d and friction velocity, the white ice albedo was also modified to 0.25, the bias and RMSE of LSWT further reduced to 2.00 • C and 2.48 • C, respectively, while the correlation coefficient was re-raised to 0.90. This result indicates that the capability of FLake in reproducing the LSWT of Dagze Co was substantially improved by the optimization of key parameters. In addition to an improvement in the LSWT simulation, the optimization of the three key parameters also enhances the capability of FLake to reproduce the internal thermal structure of the lake. As Figure 9 shows, compared with the CTRL experiment (shown in Figure 3), the OptExp_3 experiment produced a warmer metalimnion and increased MLD, mainly due to a significant increase in absorbed solar radiation for water in the metalimnion by reducing the K d . Meanwhile, the overestimated MLWT in the epilimnion of the lake by the CTRL experiment was also decreased in the OptExp_3 experiment, which ascribed to a reduced K d as well as increased friction velocity that accelerated the release of heat flux from the lake into the atmosphere. In addition, the lake ice melted earlier in the OptExp_2 experiment compared to the CTRL experiment due to the reduction in ice albedo, leading to more absorption of solar radiation after ice-off, which increased the LBWT as well as delaying the ice-on date. With the optimized key parameters, the simulation result by OptExp_3 is closer to the observations, indicating the performance of FLake in reproducing the thermal structure of Dagze Co has been significantly improved.

Salinity Parameterization
As a saline lake, the salinity in Dagze Co can affect its freezing temperature as well as its maximum density temperature, which exerts an influence on the mixing regime, and thus, the thermal structure of the lake. However, the original FLake does not consider the influence of salinity. To investigate the influence of salinity on the thermal features of Dagze Co, a simplified salinity parameterization scheme (Equations (8) and (9)) was integrated into FLake in addition to the optimization of three key parameters. As shown in

Salinity Parameterization
As a saline lake, the salinity in Dagze Co can affect its freezing temperature as well as its maximum density temperature, which exerts an influence on the mixing regime, and thus, the thermal structure of the lake. However, the original FLake does not consider the influence of salinity. To investigate the influence of salinity on the thermal features of Dagze Co, a simplified salinity parameterization scheme (Equations (8) and (9)) was integrated into FLake in addition to the optimization of three key parameters. As shown in Figure 9, compared to the original FLake simulation results, the OptExp_3 experiment greatly improved the accuracy of simulated LSWT by optimizing K d , friction velocity, and ice albedo in FLake. As shown in Figure 10, with the integration of the salinity parameterization scheme, the simulated LSWT during the open-water period shows no significant change, but the ice duration period and LSWT before and after the ice duration period slightly change. The ice-on date was slightly delayed, while the ice-off date was slightly advanced, which is closer to the observation. The bias, RMSE, and correlation coefficient of LSWT produced by the OptExp_3s experiment were further enhanced from 2.0 • C, 2.48 • C, and 0.9 to 1.84 • C, 2.28 • C, and 0.91, respectively.

Salinity Parameterization
As a saline lake, the salinity in Dagze Co can affect its freezing temperature as well as its maximum density temperature, which exerts an influence on the mixing regime, and thus, the thermal structure of the lake. However, the original FLake does not consider the influence of salinity. To investigate the influence of salinity on the thermal features of Dagze Co, a simplified salinity parameterization scheme (Equations (8) and (9)) was integrated into FLake in addition to the optimization of three key parameters. As shown in Figure 9, compared to the original FLake simulation results, the OptExp_3 experiment greatly improved the accuracy of simulated LSWT by optimizing Kd, friction velocity, and ice albedo in FLake. As shown in Figure 10, with the integration of the salinity parameterization scheme, the simulated LSWT during the open-water period shows no significant change, but the ice duration period and LSWT before and after the ice duration period slightly change. The ice-on date was slightly delayed, while the ice-off date was slightly advanced, which is closer to the observation. The bias, RMSE, and correlation coefficient of LSWT produced by the OptExp_3s experiment were further enhanced from 2.0 °C, 2.48 °C, and 0.9 to 1.84 °C, 2.28 °C, and 0.91, respectively.

Discussion
In CTRL experiments, FLake with default parameter settings slightly underestimated the LSWT of Dagze Co in spring and winter, while significantly overestimated the LSWT in summer and autumn. The same conclusion was reached in another study after modeling 94 lakes larger than 100 km 2 in China using FLake [47]. Additionally, simulations in Qinghai Lake, the largest lake on the TP, also show an overestimation of LSWT by FLake in summer and autumn [50], demonstrating the commonality of this distinguishing feature of FLake. However, a significant difference appears between our study and the previous studies in the simulation of the ice-covered period. Compared to the previous simulations in Qinghai Lake, which generally produced delayed both ice-on and ice-off data [50], the CTRL simulation in this study produced a delayed ice-on date, but with an earlier ice-off date.
The extinction coefficients strongly affect the LSWT and vertical thermal structure during summer. Previous studies have tried tuning K d to improve the performance of FLake and achieved promising outcomes [39,47]. In the sensitive experiments, we further analyzed the sensitivity of LSWT to varying K d in different months, revealing an abrupt increase in sensitivity of LSWT to K d when it is less than 0.5 m −1 , as well as the inverse response of LSWT to K d between the warm and cold seasons. However, the sensitivity of LSWT to K d is still nonlinear when K d is less than 0.5 m −1 , indicating that there may be other factors exerting an influence on the sensitivity of LSWT to K d . This result is consistent with the findings in previous studies, which suggest that the response of LSWT to change in K d is nonlinear, and is more sensitive when the K d is less than 0.5 m −1 [70,71]. Therefore, it is quite reasonable for lake models to set the default K d to 3 m −1 since K d for most of the lakes in low-lying areas is greater than 0.5 m −1 . However, most of the lakes over TP are oligotrophic with fewer biochemical processes and higher clarity, resulting in a K d smaller than 0.5 m −1 [72]. In such a case, the LSWT will be very sensitive to changes in K d for lakes over TP, indicating the importance of optimizing K d before applying FLake to TP lakes. In previous studies, different parameters, e.g., K d , depth, water albedo, ice albedo, and temperature of maximum water density, have been optimized to improve the performance of FLake, while few studies have concerned the friction velocity [39]. Our studies demonstrated that the friction velocity can significantly influence the LSWT and thermal structure of the lake. By increasing the friction velocity in FLake, the positive deviation of LSWT produced by FLake in the summer was significantly reduced. In line with previous studies, the SenExp_α experiments in our study indicate that ice albedo is the main lake parameter that affects the ice-off date [47]. The overestimation of ice albedo in FLake means more solar radiation will be reflected by the lake, which affects the melting of lake ice and leads to a delay in the ice-off date. Compared to previous studies, the SenExp_α experiments in our study further demonstrated that the ice-albedo also can affect the LSWT before and after the ice-covered period and even the internal thermal structure of the lake. Compared to the previous study in Nam Co, which set the temperature of maximum water density to an empirical value [39], another important contribution of this study is the adoption of a simple salinity parameterization that considers the effect of salinity on the temperature of maximum water density and the freezing point, which further enhanced the ability of FLake in reproducing the ice-duration period.

Conclusions
In this study, forced by the CMFD dataset, the suitability of applying the 1D lake model FLake in Dagze Co on the central TP was investigated based on a one-year in situ vertical temperature observation in Dagze Co. Then, a series of sensitive experiments were conducted to reveal the influence of three key parameters of FLake on lake thermal features. According to the physical characteristics and sensitive experiments, these key parameters were optimized to improve the performance of FLake. The main conclusions are as follows.
FLake with default settings exhibits significant deviations in simulating the LSWT and thermal profile of Dagze Co. The LSWT is overestimated, mainly in summer and autumn, with a bias and RMSE of 3.08 • C and 3.62 • C. Meanwhile, the simulated warming of LSWT in spring and cooling of LSWT in autumn are both faster than the observation, and the MLD simulated by FLake is shallower than the observation.
The sensitive experiments on the K d , friction velocity, and ice albedo revealed that FLake is very sensitive to these parameters. The reduction in K d slows the warming and cooling rate of the LSWT during its rising and declining phases, respectively, as well as deepening MLD during the stratified period. The sensitivity of LSWT in K d is nonlinear, which abruptly changes when K d is less than 0.5 m −1 . The increase in friction velocity will lead to an obvious decrease in LSWT, MCWT, and a slight increase in the MLD. Meanwhile, the reduction in ice albedo would cause an earlier ice-off date and a higher LSWT just after the ice-off, resulting in a slight warming of MWCT and LBWT during the following open water period.
The optimization of three key parameters in FLake significantly improved the performance of FLake. With K d adjusted from 3.0 m −1 to 0.25 m −1 , the variation patterns of simulated LSWT became more consistent with observation. The enlarged friction velocity substantially diminished the overestimation of LSWT produced by FLake in summer. Setting the white ice albedo from 0.60 to 0.25 increased the accuracy of FLake in modeling the ice duration period, especially the ice-off date. The experiments with salinity parameterization indicate that salinity has a minor effect on the LSWT during the open water period, while slightly affecting the ice duration period, which delays the ice-on date and advances the ice-off date, leading to a better approximation of the ice-cover simulation to that of observation. After the parameter optimization and salinity parameterization, the bias and RMSE of simulated LSWT reduced to 1.84 • C and 2.28 • C, respectively.
However, there are still limitations in this work. The primary one is the deviation originating from the forcing data, which is not calibrated in our study because of the scarcity of observations. Wind speed extracted from CMFD, which is mainly provided by land observation stations, may be significantly underestimated over the lake surface. In such a case, the optimized friction velocity may be different from the value used in this study. In addition, K d was simply set to a constant value in this study, which is not realistic as K d varies with water transparency. Hence, the parameterization scheme of K d based on satellite-observed water transparency needs to be established in the future.