Soil Freeze-thaw State Research Articles

Comparing Three Freeze-Thaw Schemes Using C-Band Radar Data in Southeastern New Hampshire, USA

Soil freeze-thaw (FT) cycles over agricultural lands are of great importance due to their vital role in controlling soil moisture distribution, nutrient availability, health of microbial communities, and water partitioning during flood events. Active microwave sensors such as C-band Sentinel-1 synthetic aperture radar (SAR) can serve as powerful tools to detect field-scale soil FT state. Using Sentinel-1 SAR observations, this study compares the performance of two FT detection approaches, a commonly used seasonal threshold approach (STA) and a computationally inexpensive general threshold approach (GTA) at an agricultural field in New Hampshire, US. It also explores the applicability of an interferometric coherence approach (ICA) for FT detection. STA and GTA achieved 85% and 78% accuracy, respectively, using VH polarization. We find a marginal degradation in the performance of STA (82%) and GTA (76%) when employing VV-polarized data. While there was approximately a 6 percentage point difference between STA’s and GTA‘s overall accuracy, we recommend GTA for FT detection using SAR images at sub-field-scale over extended regions because of its higher computational efficiency. Our analysis shows that interferometric coherence is not suitable for detecting FT transitions under mild and highly dynamic winter conditions. We hypothesize that the relatively mild winter conditions and therefore the subtle FT transitions are not able to significantly reduce the correlation between the phase values. Also, the ephemeral nature of snowpack in our study area, further compounded by frequent rainfall, could cause decorrelation of SAR images even in the absence of a FT transition. We conclude that despite Sentinel-1’s ~80% mapping accuracy at a mid-latitude site, understanding the cause of misclassification remains challenging, even when detailed ground data are readily available and employed in error attribution efforts.

Spatial and Temporal Changes in Soil Freeze-Thaw State and Freezing Depth of Northeast China and Their Driving Factors

It is necessary to further investigate the spatial considerations, temporal characteristics, and drivers of change affecting the beginning and end of soil freezing and thawing, including the maximum depth of the seasonal freezing (MDSF) and the active layer thickness (ALT) in Northeast China. Hourly soil temperature, among other data, from 1983–2022 were investigated, showing a delay of about 6 days in freezing. In contrast, thawing and complete thawing advanced by about 26 and 20 d, respectively. The freezing period and total freeze-thaw days decreased by about 29 and 23 days, respectively. The number of complete thawing period days increased by about 22 days, while the MDSF decreased by about 25 cm. The ALT increased by about 22 cm. Land Surface Temperature (LST) is the main factor influencing the beginning and end of soil freezing and thawing, MDSF and ALT changes in Northeast China; air temperature, surface net solar radiation, and volumetric soil water content followed. The influence of the interacting factors was greater than the single factors, and the interactive explanatory power of the LST and surface net solar radiation was highest when the soil started to freeze (0.858). The effect of the LST and the air temperature was highest when the soil was completely thawed (0.795). LST and the volumetric soil water content interacted to have the first explanatory power for MDSF (0.866) and ALT (0.85). The results of this study can provide scientific reference for fields such as permafrost degradation, cold zone ecological environments, and agricultural production in Northeast China.

A Spatiotemporal Enhanced SMAP Freeze/Thaw Product (1980–2020) over China and Its Preliminary Analyses

The soil freeze/thaw (FT) state has emerged as a critical role in the ecosystem, hydrological, and biogeochemical processes, but obtaining representative soil FT state datasets with a long time sequence, fine spatial resolution, and high accuracy remains challenging. Therefore, we propose a decision-level spatiotemporal data fusion algorithm based on Convolutional Long Short-Term Memory networks (ConvLSTM) to expand the SMAP-enhanced L3 landscape freeze/thaw product (SMAP_E_FT) temporally. In the algorithm, the Freeze/Thaw Earth System Data Record product (ESDR_FT) is sucked in the ConvLSTM and fused with SMAP_E_FT at the decision level. Eight predictor datasets, i.e., soil temperature, snow depth, soil moisture, precipitation, terrain complexity index, area of open water data, latitude and longitude, are used to train the ConvLSTM. Direct validation using six dense observation networks located in the Genhe, Maqu, Naqu, Pali, Saihanba, and Shandian river shows that the fusion product (ConvLSTM_FT) effectively absorbs the high accuracy characteristics of ESDR_FT and expands SMAP_E_FT with an overall average improvement of 2.44% relative to SMAP_E_FT, especially in frozen seasons (averagely improved by 7.03%). The result from indirect validation based on categorical triple collocation also shows that ConvLSTM_FT performs stable regardless of land cover types, climate types, and terrain complexity. The findings, drawn from preliminary analyses on ConvLSTM_FT from 1980 to 2020 over China, suggest that with global warming, most parts of China suffer from different degrees of shortening of the frozen period. Moreover, in the Qinghai–Tibet region, the higher the permafrost thermal stability, the faster the degradation rate.

Effects of Arctic Wetland Dynamics on Tower-Based GNSS Reflectometry Observations

A tower-based global navigation satellite system reflectometry (GNSS-R) experiment is set up in an Arctic wetland environment for investigating the possibility of monitoring wetland inundation and freeze/thaw (FT) dynamics which are additionally impacted by snow on the ground. Effects of inundation, snow cover, and soil FT state on observed GNSS-R signal-to-noise ratios (SNRs) are analyzed for horizontal (H) and vertical (V) polarizations. A simple classification approach is suggested to detect the inundated, frozen, or thawed soil state. A simple forward reflectivity model is formulated to evaluate the influence of snow cover, overlying frozen, or thawed soil, on the reflected GNSS signals. Reflectivity time series are simulated in H- and V-polarizations using <italic xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">in situ</i> observations of the Arctic wetland site. The simulations are used to verify the tower-based observations, which show a significant impact of wet snow on reflectivity during melting conditions in spring. The observed SNR is strongly correlated with the Sentinel-1 backscatter coefficient. Generally, soil states detected by GNSS-R are in high agreement with ground truth soil states, especially for inundated and frozen soils. Wet snow conditions, however, complicate the correct timing estimation of soil thawing by inducing reflectivities of a similar order as thawing soil. It is recommended that GNSS-R land application models and retrieval algorithms consider snow cover effects to reduce false classification, especially in FT detection. Overall, the outcome of this study is relevant to the upcoming ESA HydroGNSS mission.

Spatiotemporal characteristics of the soil freeze-thaw state and its variation under different land use types - A case study in Northeast China

With global warming, the soil freeze-thaw state has greatly changed, which has profoundly impacted the balance between water resources, ecological security and regional agriculture and forestry production in frozen regions. Based on daily soil freezing/thawing depth data in Heilongjiang Province over the last 60 years, Mann-Kendall trend analysis, Sen's slope and Hurst exponent were performed to analyze the spatiotemporal characteristics of the soil freeze-thaw state, and the effects of nonclimatic factors "land use types" and snow cover on the change in the soil freeze-thaw state were discussed. The results indicated that (1) the soil freeze-thaw state significantly changed over the past 60 years. The maximum seasonally frozen depth decreased approximately by 48 cm, the first freezing date was delayed by approximately 10 d, the first thawing date and last date advanced approximately by 6 and 27 d, respectively, the freezing and thawing periods decreased by approximately 17 and 21 d, respectively, and the duration and actual number of freezing days decreased by 38 and 34 d, respectively. (2) The soil freeze-thaw state changed abruptly in 1988 and exhibited a significant trend after the 1990s. The Hurst exponent of the soil freeze-thaw state variables at approximately 2/3 of the stations reached the strong persistence level. (3) In forestland and grassland areas, vegetation coverage reduced the climate warming impact on the soil freeze-thaw state. In addition, the maximum snow depth was significantly correlated with the soil freezing depth and first thawing date, and the change in snow depth also offset the effects of temperature rise on soil freezing depth and first thawing date to a certain extent. Changes in the soil freeze-thaw state have both advantages and disadvantages on agricultural production. Therefore, it is necessary to strengthen agricultural monitoring and use territorial knowledge to reasonably regulate agricultural production.

Spring photosynthetic onset and net CO2 uptake in Alaska triggered by landscape thawing.

The springtime transition to regional-scale onset of photosynthesis and net ecosystem carbon uptake in boreal and tundra ecosystems are linked to the soil freeze-thaw state. We present evidence from diagnostic and inversion models constrained by satellite fluorescence and airborne CO2 from 2012 to 2014 indicating the timing and magnitude of spring carbon uptake in Alaska correlates with landscape thaw and ecoregion. Landscape thaw in boreal forests typically occurs in late April (DOY 111±7) with a 29±6day lag until photosynthetic onset. North Slope tundra thaws 3weeks later (DOY 133±5) but experiences only a 20±5day lag until photosynthetic onset. These time lag differences reflect efficient cold season adaptation in tundra shrub and the longer dehardening period for boreal evergreens. Despite the short transition from thaw to photosynthetic onset in tundra, synchrony of tundra respiration with snow melt and landscape thaw delays the transition from net carbon loss (at photosynthetic onset) to net uptake by 13±7days, thus reducing the tundra net carbon uptake period. Two global CO2 inversions using a CASA-GFED model prior estimate earlier northern high latitude net carbon uptake compared to our regional inversion, which we attribute to (i) early photosynthetic-onset model prior bias, (ii) inverse method (scaling factor+optimization window), and (iii) sparsity of available Alaskan CO2 observations. Another global inversion with zero prior estimates the same timing for net carbon uptake as the regional model but smaller seasonal amplitude. The analysis of Alaskan eddy covariance observations confirms regional scale findings for tundra, but indicates that photosynthesis and net carbon uptake occur up to 1month earlier in evergreens than captured by models or CO2 inversions, with better correlation to above-freezing air temperature than date of primary thaw. Further collection and analysis of boreal evergreen species over multiple years and at additional subarctic flux towers are critically needed.

Parameterization of the freeze/thaw discriminant function algorithm using dense in-situ observation network data

ABSTRACTThe near-surface soil freeze–thaw (FT) transition is an important factor affecting land-atmosphere exchanges, hydrology and carbon cycles. Thus, effectively monitoring the temporal–spatial changes of soil FT processes is crucial to climate change and environment research. Several approaches have been developed to detect the soil FT state from satellite observations. The discriminant function algorithm (DFA) uses temperature and emissivity information from Advanced Microwave Scanning Radiometer Enhanced (AMSR-E) passive microwave satellite observations. Although it is well validated, it was shown to be insufficiently robust for all land conditions. In this study, we use in-situ observed soil temperature and AMSR-E brightness temperature to parameterize the DFA for soil FT state detection. We use the in-situ soil temperature records at 5 cm selected from available dense networks in the Northern Hemisphere as a reference. Considering the distinction between ascending and descending orbits, two different sets of parameters were acquired for each frequency pair. The validation results indicate that the overall discriminant accuracy of the new function can reach 90%. We further compared the Advanced Microwave Scanning Radiometer 2 discriminant results using the new function to the Soil Moisture Active Passive freeze/thaw product, and a reasonable consistency between them was found.

Characterizing Snowpack and the Freeze–Thaw State of Underlying Soil via Assimilation of Multifrequency Passive/Active Microwave Data: A Case Study (NASA CLPX 2003)

Ground-based passive microwave observations at 18.7- and 36.5-GHz frequencies and active microwave observations in L- [1.4 GHz] and Ku- [15.5 GHz] bands are used within an ensemble-based data assimilation (DA) framework to characterize the snow water equivalent (SWE) and the underlying soil freeze-thaw state (including soil surface temperature and both soil ice/liquid water content). The proposed framework is tested at the local-scale observation site of the National Aeronautics and Space Administration (NASA) Cold Land Processes Experiment field campaign during the third intensive observation period (February 18-26, 2003) for which the best set of collocated ground-based passive/active microwave observations, SWE, soil surface temperature, and moisture measurements are available. The DA approach effectively merges an a priori estimate of the soil freeze-thaw state and SWE generated by a land surface model (LSM) with information contained in passive/active microwave observations in order to overcome errors in the forcing data of LSM. Results indicate that the root-mean-square errors of SWE, soil surface temperature, and soil ice+liquid water content after the assimilation of passive (active) observations respectively decrease to 25.4 mm (22.8 mm), 0.61 K (0.52 K), and 0.063 (0.057) from 90.55 mm, 2.17 K, and 0.13 before assimilation, resulting in improvements of 75% (77%), 72% (76%), and 51% (56%). Also, it is found that the simultaneous assimilation of passive and active measurements further improves the estimates of SWE and soil temperature as well as soil ice/liquid water content, suggesting that there is an advantage offered by the synergistic use of passive and active measurements. Overall, the findings show that future studies can take advantage of remotely sensed microwave passive and active measurements from present and upcoming satellites such as Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E), Soil Moisture Active Passive, and COld REgion Hydrology High-resolution Observatory (CoReH2O) for monitoring SWE and the underlying soil freeze-thaw state.

Feasibility of Characterizing Snowpack and the Freeze–Thaw State of Underlying Soil Using Multifrequency Active/Passive Microwave Data

An ensemble-based data assimilation approach is developed to characterize the snow water equivalent (SWE) and underlying soil freeze-thaw state (including the soil surface temperature and both soil ice and liquid water content) using multifrequency passive and active microwave remote-sensing measurements. Its feasibility was examined using a synthetic test where passive microwave (1.4, 18.7, and 36.5 GHz) and active microwave [L-band (1.4 GHz), C-band (5.4 GHz), and Ku-band (12 GHz)] measurements at the point scale were individually and simultaneously assimilated to estimate the SWE and soil freeze-thaw state using an Ensemble Batch Smoother framework. The contribution of each channel in retrieving the true SWE, soil surface temperature, soil liquid water and ice content was investigated at the local-scale observation site of the National Aeronautics and Space Administration Cold Land Processes Experiments Field Campaign in northern Colorado during both the snow accumulation (Fall 2002-Winter 2003) and melt (Spring 2003) periods. All of the utilized passive and active measurements were found to contain valuable and complementary information for characterizing the SWE and freeze-thaw state of the underlying soil. L-band measurements were most effective for soil freeze-thaw state estimation, whereas higher frequencies were more effective at SWE characterization. In addition, results from the simultaneous assimilation of passive and active microwave data were compared to those from a modeling approach without assimilating microwave data (open loop). It was found that assimilating both passive and active microwave data decreased the errors that are associated with the open-loop approach. Finally, passive and active measurements were undersampled as expected from the overpasses of current and future satellite platforms. It was observed that the developed method can reliably estimate the soil freeze-thaw state and SWE, even with measurement sequences anticipated from the temporal frequency of existing and future satellites such as the Special Sensor Microwave/Imager, Soil Moisture Active Passive Mission, and Cold Regions Hydrology High-Resolution Observatory.

Variable infiltration capacity cold land process model updates

The Variable Infiltration Capacity (VIC) macroscale hydrologic model is distinguished from other Soil–Vegetation–Atmosphere Transfer schemes (SVATS) by its focus on runoff processes. These are represented via the variable infiltration curve, a parameterization of the effects of subgrid variability in soil moisture holding capacity, from which the model takes its name, and a representation of nonlinear baseflow. Recent upgrades to the model have improved its representation of cold land processes, and the effects of surface storage in lakes and wetlands. Specific improvements described in this paper include the following: (1) explicit representation of the canopy energy balance separate from the land surface when snow is intercepted in the canopy; (2) parameterization of the effects of spatial variability in soil freeze–thaw state and snow distribution on moisture and energy fluxes; and (3) effects of advection on snowmelt under conditions of partial snow cover. The effects of these model updates are demonstrated using data from the PILPS Phase 2(e) validation catchments within the Torne-Kalix River basin, Sweden.