Forest cover is a crucial factor that influences the performance of optical satellite-based snow cover monitoring algorithms. However, evaluation of such algorithms in forested landscapes is rare due to lack of reliable in situ data in such regions. In this investigation, we assessed the performance of the operational snow detection (SCA) and fractional snow cover estimation (FSC) algorithms employed by the Copernicus Land Monitoring Service for High-Resolution Snow & Ice Monitoring (HRSI) with a combination of Sentinel-2 and Landsat-7/8 satellite scenes, lidar-based, and in situ datasets. These algorithms were evaluated over test sites located in the forested mountainous landscape of the Pyrenees in Spain and the Sierra Nevada in the USA. Over the Pyrenees site, the effectiveness of snow cover detection was evaluated with respect to a time-series of in situ snow depth measurements logged over test plots with different aspects, canopy cover, and solar irradiance. Over the Sierra Nevada site, the impact of ground vegetation was assessed over the under canopy fractional snow cover retrievals using airborne lidar-derived fractional vegetation cover information. The analyses over the Pyrenees indicated a good accuracy of snow detection with the exception of plots with either dense canopy cover or insufficient solar exposure (shaded forested slope), or both. The operational HRSI algorithm yielded similar performances (25-30% RMSE) as the computationally intensive spectral unmixing approach while retrieving the subcanopy ground FSC over the Sierra Nevada site. It was observed that a more accurate lidar-derived tree cover density map did not improve the subcanopy FSC retrievals.

Owing to their important roles in biogeochemical cycles, phytoplankton functional types (PFTs) have been the aim of an increasing number of ocean color algorithms. Yet, none of the existing methods are based on phytoplankton carbon (C) biomass, which is a fundamental biogeochemical and ecological variable and the “unit of accounting” in Earth system models. We present a novel bio-optical algorithm to retrieve size-partitioned phytoplankton carbon from ocean color satellite data. The algorithm is based on existing methods to estimate particle volume from a power-law particle size distribution (PSD). Volume is converted to carbon concentrations using a compilation of allometric relationships. We quantify absolute and fractional biomass in three PFTs based on size – picophytoplankton (0.5–2 µm in diameter), nanophytoplankton (2–20 µm) and microphytoplankton (20–50 µm). The mean spatial distributions of total phytoplankton C biomass and individual PFTs, derived from global SeaWiFS monthly ocean color data, are consistent with current understanding of oceanic ecosystems, i.e., oligotrophic regions are characterized by low biomass and dominance of picoplankton, whereas eutrophic regions have high biomass to which nanoplankton and microplankton contribute relatively larger fractions. Global climatological, spatially integrated phytoplankton carbon biomass standing stock estimates using our PSD-based approach yield ∼ 0.25 Gt of C, consistent with analogous estimates from two other ocean color algorithms and several state-of-the-art Earth system models. Satisfactory in situ closure observed between PSD and POC measurements lends support to the theoretical basis of the PSD-based algorithm. Uncertainty budget analyses indicate that absolute carbon concentration uncertainties are driven by the PSD parameter N_o which determines particle number concentration to first order, while uncertainties in PFTs' fractional contributions to total C biomass are mostly due to the allometric coefficients. The C algorithm presented here, which is not empirically constrained a priori, partitions biomass in size classes and introduces improvement over the assumptions of the other approaches. However, the range of phytoplankton C biomass spatial variability globally is larger than estimated by any other models considered here, which suggests an empirical correction to the N_o parameter is needed, based on PSD validation statistics. These corrected absolute carbon biomass concentrations validate well against in situ POC observations.

Snow disappearance date (SDD) affects the ecohydrological dynamics of montane forests, by altering water availability, forest fire regime, and the land surface energy budget. The forest canopy modulates SDD through competing processes; dense canopy intercepts snowfall and enhances longwave radiation while shading snowpack from shortwave radiation and sheltering it from the wind. Limited ground‐based observations of snow presence and absence have restricted our ability to unravel the dominant processes affecting SDD in montane forests. We apply a lidar‐derived method to estimate fractional snow cover area (fSCA) at two relatively warm sites in the Sierra Nevada and two colder sites in the Rocky Mountains, which we link to SDD. With the exception of late season snowpack and low fSCA, snow retention is longer under low vegetation density than under high vegetation density in both warm and cold sites. Warm forests consistently have longer snow retention in open areas compared to dense under canopy areas, particularly on south‐facing slopes. Cold forests tend to have longer snow retention under lower density canopy compared to open areas, particularly on north‐facing slopes. We use this empirical analysis to make process inferences and develop an initial framework to predict SDD that incorporates the role of topography and vegetation structure. Building on our framework will be necessary to provide better forest management recommendations for snowpack retention across complex terrain and heterogenous canopy structure.

The particle size distribution (PSD) of suspended particles in near-surface seawater is a key property linking biogeochemical and ecosystem characteristics with optical properties that affect ocean color remote sensing. Phytoplankton size affects their physiological characteristics and ecosystem and biogeochemical roles, e.g., in the biological carbon pump, which has an important role in the global carbon cycle and thus climate. It is thus important to develop capabilities for measurement and predictive understanding of the structure and function of oceanic ecosystems, including the PSD, phytoplankton size classes (PSCs), and phytoplankton functional types (PFTs). Here, we present an ocean color satellite algorithm for the retrieval of the parameters of an assumed power-law PSD. The forward optical model considers two distinct particle populations: phytoplankton and non-algal particles (NAPs). Phytoplankton are modeled as coated spheres following the Equivalent Algal Populations (EAP) framework, and NAPs are modeled as homogeneous spheres. The forward model uses Mie and Aden–Kerker scattering computations, for homogeneous and coated spheres, respectively, to model the total particulate spectral backscattering coefficient as the sum of phytoplankton and NAP backscattering. The PSD retrieval is achieved via spectral angle mapping (SAM), which uses backscattering end-members created by the forward model. The PSD is used to retrieve size-partitioned absolute and fractional phytoplankton carbon concentrations (i.e., carbon-based PSCs), as well as particulate organic carbon (POC), using allometric coefficients. This model formulation also allows the estimation of chlorophyll a concentration via the retrieved PSD, as well as percent of backscattering due to NAPs vs. phytoplankton. The PSD algorithm is operationally applied to the merged Ocean Colour Climate Change Initiative (OC-CCI) v5.0 ocean color data set. Results of an initial validation effort are also presented using PSD, POC, and picophytoplankton carbon in situ measurements. Validation results indicate the need for an empirical tuning for the absolute phytoplankton carbon concentrations; however these results and comparison with other phytoplankton carbon algorithms are ambiguous as to the need for the tuning. The latter finding illustrates the continued need for high-quality, consistent, large global data sets of PSD, phytoplankton carbon, and related variables to facilitate future algorithm improvements.

Marine phytoplankton biomass and community structure are expected to change under global warming, with potentially significant impacts on ocean carbon, nutrient cycling, and marine food webs. Previous studies have indicated decreases of primary production and chlorophyll a concentrations and oligotrophic gyre expansions from satellite ocean‐color measurements, purportedly due to global warming. We review this topic via a reanalysis of a novel backscattering‐based phytoplankton functional type and phytoplankton biomass time series over the 1997–2010 period. Unlike previous work, we find that globally the biomass and the percent of large (small) phytoplankton increase (decrease). The oligotrophic gyres contract or expand depending on the chlorophyll a threshold definition employed. In the subtropical gyres, chlorophyll a trends are likely due to physiological changes, while the increasing biomass trends are due to winds and relevant mixing length scale increases.

Milankovitch theory postulates that periodic variability of Earth's orbital elements is a major climate forcing mechanism, causing, for example, the contemporary glacial–interglacial cycles. There are three Milankovitch orbital parameters: orbital eccentricity, precession and obliquity. The interaction of the amplitudes, periods and phases of these parameters controls the spatio-temporal patterns of incoming solar radiation (insolation) and the timing and duration of the seasons. This complexity makes Earth–Sun geometry and Milankovitch theory difficult to teach effectively. Here, we present "Earth Orbit v2.1": an astronomically precise and accurate model that offers 3-D visualizations of Earth's orbital geometry, Milankovitch parameters and the ensuing insolation forcing. The model is developed in MATLAB® as a user-friendly graphical user interface. Users are presented with a choice between the Berger (1978a) and Laskar et al. (2004) astronomical solutions for eccentricity, obliquity and precession. A "demo" mode is also available, which allows the Milankovitch parameters to be varied independently of each other, so that users can isolate the effects of each parameter on orbital geometry, the seasons, and insolation. A 3-D orbital configuration plot, as well as various surface and line plots of insolation and insolation anomalies on various time and space scales are produced. Insolation computations use the model's own orbital geometry with no additional a priori input other than the Milankovitch parameter solutions. Insolation output and the underlying solar declination computation are successfully validated against the results of Laskar et al. (2004) and Meeus (1998), respectively. The model outputs some ancillary parameters as well, e.g., Earth's radius-vector length, solar declination and day length for the chosen date and latitude. Time-series plots of the Milankovitch parameters and several relevant paleoclimatological data sets can be produced. Both research and pedagogical applications are envisioned for the model.