• We present four methods to calculate LAI on a daily basis from PAR. • Each method shows high linear correlation to MODIS and LAI-2000 datasets. • All methods provide a precise indication of start and end of the growing season. • PAR based LAI has broad potential to reveal phenological response to global change. Leaf area index (LAI) is a critical biophysical indicator that describes foliage abundance in ecosystems. An accurate and continuous estimation of LAI is therefore desirable to quantify ecosystem status and function (e.g. carbon and water exchange between the land surface and the atmosphere). However, deriving accurate LAI measurements at regular temporal intervals remains challenging, requiring either destructive sampling or manual collection of canopy gap fraction measurements at discrete time intervals. In this study, we present four methods to obtain continuous LAI data, simply derived from above and below canopy measurements of photosynthetically active radiation (PAR) at the Borden Forest Research Station from 1999 to 2018. We compared LAI derived using the four PAR-based methods to independent measurements of LAI from optical methods and the MODIS satellite LAI product. LAI derived from all four PAR-based methods captured the seasonal changes in observed and remotely sensed LAI and showed a close linear correspondence with one another (R 2 of 0.55 to 0.76 compared to MODIS LAI, and R 2 of 0.78 to 0.84 compared to LAI-2000 measurements). A PAR-based method using Miller's Integral theorem showed the strongest linear relationship with LAI-2000 measurements (R 2 =0.84, p<0.001, SE=0.40). In many years MODIS LAI indicated an earlier start of season and earlier end of season than the daily PAR-based LAI datasets showing systematic biases in the MODIS assessment of growing season. The four PAR-based LAI methods outlined in this study provide an LAI dataset of unprecedented temporal resolution. These methods will allow precise determination of phenological events, improve leaf to canopy scaling in process-based models, and provide valuable insight into dynamic vegetation responses to global climate change.

Abstract Twenty-seven models participated in the Earth System Model–Snow Model Intercomparison Project (ESM-SnowMIP), the most data-rich MIP dedicated to snow modeling. Our findings do not support the hypothesis advanced by previous snow MIPs: evaluating models against more variables and providing evaluation datasets extended temporally and spatially does not facilitate identification of key new processes requiring improvement to model snow mass and energy budgets, even at point scales. In fact, the same modeling issues identified by previous snow MIPs arose: albedo is a major source of uncertainty, surface exchange parameterizations are problematic, and individual model performance is inconsistent. This lack of progress is attributed partly to the large number of human errors that led to anomalous model behavior and to numerous resubmissions. It is unclear how widespread such errors are in our field and others; dedicated time and resources will be needed to tackle this issue to prevent highly sophisticated models and their research outputs from being vulnerable because of avoidable human mistakes. The design of and the data available to successive snow MIPs were also questioned. Evaluation of models against bulk snow properties was found to be sufficient for some but inappropriate for more complex snow models whose skills at simulating internal snow properties remained untested. Discussions between the authors of this paper on the purpose of MIPs revealed varied, and sometimes contradictory, motivations behind their participation. These findings started a collaborative effort to adapt future snow MIPs to respond to the diverse needs of the community.

Abstract. Thirty-year simulations of seasonal snow cover in 22 physically based models driven with bias-corrected meteorological reanalyses are examined at four sites with long records of snow observations. Annual snow cover durations differ widely between models but interannual variations are strongly correlated because of the common driving data. No significant trends are observed in starting dates for seasonal snow cover, but there are significant trends towards snow cover ending earlier at two of the sites in observations and most of the models. A simplified model with just two parameters controlling solar radiation and sensible heat contributions to snowmelt spans the ranges of snow cover durations and trends. This model predicts that sites where snow persists beyond annual peaks in solar radiation and air temperature will experience rapid decreases in snow cover duration with warming as snow begins to melt earlier and at times of year with more energy available for melting.

Abstract. The 30-year simulations of seasonal snow cover in 22 physically based models driven with bias-corrected meteorological reanalyses are examined at four sites with long records of snow observations. Annual snow cover durations differ widely between models, but interannual variations are strongly correlated because of the common driving data. No significant trends are observed in starting dates for seasonal snow cover, but there are significant trends towards snow cover ending earlier at two of the sites in observations and most of the models. A simplified model with just two parameters controlling solar radiation and sensible heat contributions to snowmelt spans the ranges of snow cover durations and trends. This model predicts that sites where snow persists beyond annual peaks in solar radiation and air temperature will experience rapid decreases in snow cover duration with warming as snow begins to melt earlier and at times of year with more energy available for melting.

Abstract. This paper describes ESM-SnowMIP, an international coordinated modelling effort to evaluate current snow schemes, including snow schemes that are included in Earth system models, in a wide variety of settings against local and global observations. The project aims to identify crucial processes and characteristics that need to be improved in snow models in the context of local- and global-scale modelling. A further objective of ESM-SnowMIP is to better quantify snow-related feedbacks in the Earth system. Although it is not part of the sixth phase of the Coupled Model Intercomparison Project (CMIP6), ESM-SnowMIP is tightly linked to the CMIP6-endorsed Land Surface, Snow and Soil Moisture Model Intercomparison (LS3MIP).