Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameterization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

Sea State Based Estimation of White Cap Fraction: Implications for Primary Marine Aerosol Fluxes

Oceanic whitecaps (hereafter, W) or the characteristic whiteness of the sea foam is an important feature for predicting exchange of gases, sea spray aerosols (SSAs), heat and momentum transfer between the ocean and the atmosphere at the air-sea interface. Due to its increased surface emission and brightness temperature, whitecaps are critical for satellite retrievals of ocean albedo, ocean color, ocean surface wind vectors from satellite borne radiometer and microwave instruments. Most of the existing models predict W using wind speed and sea surface temperature (SST). However, numerous publications have pointed out that there are large uncertainties in the predicted W and using parameterizations based on wind-wave state can improve the precision of the predicted W. Here, we integrate the University of Miami Wave Model - 2.0 (UMWM) in Goddard Earth Observing System (GEOS) and use wave diagnostics to predict W. We choose the year 2006 for our global UMWM/GEOS runs because of the availability of W dataset from satellite observations. We run UMWM/GEOS at 0.5o x 0.5o by replaying to MERRA2 meteorology and evaluate the wave diagnostics using measurements from fixed buoys and satellite altimeters. We use three different parameterizations for W based on: 1) Reynolds number, 2) wave dissipation energy, and 3) volume of air entrained by breaking waves. We compare our results of W with previous studies and also with the satellite based observational dataset. Predicting W is important for understanding the processes at the air-sea interface. Therefore, this work is a step further in improving the uncertainties in the aerosol and atmospheric chemistry modules of the global models.

Oceanic Whitecaps↗

The Use of MERRA-2 Near Surface Meteorology to Understand the Behavior of Planetary Boundary Layer Heights Derived from Wind Profiler Data over the US Great Plains

The atmospheric general circulation model (GCM) that underlies the MERRA-2 reanalysis includes a suite of physical parameterizations that describe the processes that occur in the planetary boundary layer (PBL). The data assimilation system assures that the atmospheric state variables used as input to these parameterizations are constrained to the best fit to all of the available observations. Many studies, however, have shown that the GCM-based estimates of MERRA-2 PBL heights are biased high, and so are not reliable for boundary layer studies.A 20-year record of PBL heights was derived from Wind Profiler (WP) backscatter data measured at a wide network of stations throughout the US Great Plains and has been validated against independent estimates. The behavior of these PBL heights shows geographical and temporal variations that are difficult to attribute to particular physical processes without additional information that are not part of the observational record.In the present study, we use information on physical processes from MERRA-2 to understand the behavior of the WP derived PBL heights. The behavior of the annual cycle of both MERRA-2 and WP PBL heights shows four classes of behavior: (i) canonical, characterized by a monthly progression in PBL height that follows the solar insolation, (ii) double peak, characterized by canonical behavior that is interrupted by a minimum in July, (iii) late peak, characterized by a suppressed heights in May and June, and return to canonical in July and August, and (iv) early peak where the PBL height rises with solar insolation but is suppressed later in the summer. The explanation for these behaviors and the relationship to local precipitation, temperature, sensible and latent heat fluxes, net radiation and aerosol load is articulated using information from MERRA-2.

Salmun, Haydee↗

Global Cloud-Resolving Models

Global cloud-resolving models (GCRMs) are a new category of atmospheric global models designed to solve different flavors of the nonhydrostatic equations through the use of kilometer-scale global meshes. GCRMs make it possible to explicitly simulate deep convection, thereby avoiding the need for cumulus parameterization and allowing for clouds to be “resolved” by microphysical models responding to grid-scale forcing. GCRMs require high-resolution discretization over the globe, for which a variety of mesh structures have been proposed and employed. The first GCRM was constructed 15 years ago, and in recent years, other groups have also begun adopting this approach, enabling the first intercomparison studies of such models. Because conventional general circulation models (GCMs) suffer from large biases associated with cumulus parameterization, GCRMs are attractive tools for researchers studying global weather and climate. In this review, GCRMs are described, with some emphasis on their historical development and the associated literature documenting their use. The advantages of GCRMs are presented, and currently existing GCRMs are listed and described. Future prospects for GCRMs are also presented in the final section.

Satoh, Masaki↗

Extratropical Cyclone Clouds in the GFDL Climate Model: Diagnosing Biases and the Associated Causes

The clouds in southern hemisphere extratropical cyclones generated by the GFDL climate model are analyzed against MODIS, CloudSat and CALIPSO cloud and precipitation observations. Two model versions are used: one is a developmental version of AM4, a model GFDL will utilize for CMIP6, the other is the same model with a different parameterization of moist convection. Both model versions predict a realistic top-of-atmosphere cloud cover in the southern oceans, within 5% of the observations. However, an examination of cloud cover transects in extratropical cyclones reveals a tendency in the models to overestimate high-level clouds (by differing amounts) and underestimate cloud cover at low-levels (again by differing amounts), especially in the post-cold frontal (PCF) region, when compared to observations. Focusing on only the models, their differences in high and mid-level clouds are consistent with their differences in convective activity and relative humidity (RH), but the same is not true for the PCF region. In this region, RH is higher in the model with less cloud fraction. These seemingly contradictory cloud and RH differences can be explained by differences in the cloud parameterization tuning parameters that ensure radiative balance. In the PCF region, the model cloud differences are smaller than either of the model biases with respect to observations, suggesting other physics changes are needed to address the bias. The process-oriented analysis used to assess these model differences will soon be automated and shared.

Catherine M Naud↗

Extending XBAER Algorithm to Aerosol and Cloud Condition

The retrieval of cloud optical properties for aerosol contaminated water cloud is challenging because of the complexity of physical processes in such situations. Conventionally, cloud optical data products are typically derived, ignoring the aerosol impacts on radiative transfer in the retrieval process. This is potentially a significant source of error. In this paper, the eXtensible Bremen Aerosol Retrieval (XBAER) algorithm has been optimized for the retrieval of aerosol/cloud properties for the aerosol contaminated cloud (ACC) scenarios. This version of XBAER delivers cloud optical thickness (COT) and cloud effective radius (CER) for ACCs and simultaneously retrieves aerosol optical thickness (AOT). The surface parameterization and aerosol types used in the standard XBAER algorithm have been adapted in this retrieval to account for the ACC conditions. Aerosol types in XBAER for the retrieval of ACC scenarios have been parameterized to comprise weak and strong absorptions. The comparisons of COT, CER, and AOT retrieved using this adapted XBAER algorithm and two new NASA algorithms show good agreements, especially for biomass burning aerosol. The correlation coefficients are >0.9 for AOT, ~0.8 for CER, and ~0.7 for COT. There is also a good agreement for dust plume contaminated cloud scene. The XBAER derived AOT values for dust aerosols are systematically smaller than the NASA retrieval, but both products have the same spatial distribution patterns. The comparison of COT retrieved using the adapted XBAER algorithm and that retrieved from ground-based microwave radiometer (MWR) measurements shows much better agreement for ACC conditions with high AOT.

satellite↗

The Diurnal Cycle of Precipitation and Organized Convection in a Set of Global Mesoscale Simulations with the NASA GEOS AGCM

The NASA Global Modeling and Assimilation Office (GMAO) has conducted a series of 40-day nonhydrostatic global simulations with horizontal grid spacing ranging from 200 km to 3 km, as part of the DYAMOND model intercomparison project. The Goddard Earth Observing System (GEOS) model was run with the Grell-Freitas scale-aware convection scheme, which smoothly reduces parameterized deep convection with increasing resolution. Here we evaluate the diurnal cycle and other statistics of precipitation and organized convection as a function of resolution. For validation we use the 0.1 degree IMERG precipitation and 4 km Merged IR brightness temperature datasets, focusing on four regions: the continental United States, Amazonia, the equatorial Indian ocean, and the Maritime Continent. Early results indicate good phase agreement but excessive magnitude of the continental diurnal cycle of precipitation at coarser resolutions, with improved magnitude as resolution increases and the role of parameterization is reduced. Convective cloud clusters are identified with a brightness temperature threshold, and we find realistic numbers of the largest clusters (>10^4 km sq) at all resolutions, while the number of smaller clusters increases with resolution, approaching observations when dx=3 km. The observed diurnal cycle in the cluster size distribution is also reproduced, with realistic magnitude in the highest resolution runs. Precipitation characteristics across cluster sizes are also examined. The results show the potential for global mesoscale simulations as a community science resource, and we invite collaboration to explore these runs in greater detail.

Arnold, Nathan↗

Utilization of Airborne Observations to Assess Model Parametrizations of Critical RH Profiles in the Arctic Ocean

Observational paucity is a reality in the Arctic Ocean. This is especially true for near-surface variables such as temperature, moisture, heat fluxes and BL clouds. As a result, modeling has become one of the major avenues for understanding current and future Arctic trends. Reanalyses are frequently used to force global ocean circulation and sea-ice models. But in northern high latitudes, model integrations and reanalyses are known to have large uncertainties in temperature and humidity profiles, and in boundary layer cloudiness. These are common sources of error in the surface radiative budget terms.An important way to diagnose these biases spatially and temporally is by using satellite remote sensing data. However, remotely-sensed observations also have large uncertainties, especially in near-surface temperature and relative humidity profiles.In situ observational studies are important in bridging our knowledge gap in regions such as the Arctic Ocean. Here, we utilize airborne and ship observations from the ARISE, ACME-IV, and ASCOS campaigns to construct critical relative humidity (RH) profiles over the Beaufort Sea. Such profiles are used as parameterization inputs in the NASA GOES global model to derive the total water condensate in a model grid-box, which determines the cloud fraction. Currently, the critical RH profiles are derived by global AIRS data, relaying mostly on mid-latitude regions, which are not necessarily relevant to the Arctic.We derive campaign-wide mean, standard deviation, and critical RH values, for grid size of 50x50 km and altitude bins between 50 to 400 m, covering both open-ocean and sea-ice covered regions. We compare profiles over open ocean and sea-ice, and look at correlations between the observed critical RH values and water condensate (by cloud number concentration) from observations versus the modeled ones. We then input our calculated values of minimal critical RH values into a set of GEOS single column model (SCM) simulations over the ARISE and ASCOS regions and compare the differences between the predicted values of cloud liquid water path (LWP), ice water path (IWP) and surface fluxes with the observed ones under the range of input parameterizations. Finally, we discuss the implications on surface radiative budget predictions in this region.

Segal-Rozenhaimer, Michal↗

The Present and Future of Secondary Organic Aerosol Direct Forcing on Climate

Secondary organic aerosols (SOA), a subset of organic aerosols that are chemically produced in the atmosphere, are included in climate modeling calculations using very simple parameterizations. Estimates on their shortwave forcing on climate span almost two orders of magnitude, being potentially comparable to sulfate direct forcing. In the longwave, a neglected part of the spectrum when it comes to SOA, the direct SOA forcing could exceed that of sulfate and black carbon, although in absolute values it is much weaker than the shortwave forcing. Critical for these estimates is the vertical distribution of the climate active agents, pointing to SOA temperature-dependent volatility. Over the last few years, research also revealed the highly oxidized character of organic aerosol and its chemical aging in the atmosphere that partially leads to the formation of brown carbon, an absorbing form of organic aerosol. This review summarizes critical advances in the understanding of SOA behavior and properties relevant to direct climate forcing and puts them in perspective with regard to primary organic aerosol and brown carbon. These findings also demonstrate an emerging dynamic picture of organic aerosol that has not yet been integrated in climate modeling. The challenges for the coming years in order to reduce uncertainties in the direct organic aerosol climate impact are discussed. High priority for future model development should be given to the dynamic link between “white” and “brown” organic aerosol and between primary and secondary organic aerosol. The SOA temperature-dependent volatility parameterizations and wavelength-dependent refractive index should be also included.

Kostas Tsigaridis↗

GISS Model E2.2: A Climate Model Optimized for the Middle Atmosphere. Part 1: Model Structure, Climatology, Variability and Climate Sensitivity

We introduce a new climate model (GISS E2.2) that has been specially optimized for the Middle Atmosphere, and whose output is being contributed to the CMIP6 archive. The top of the model is at a geopotential altitude of 89 km, and parameterizations of moist convection and various forms of gravity wave drag based on tropospheric processes are chosen specifically for this optimization. We first evaluate the model in its configuration as a coupled atmosphere-chemistry model with respect to its simulation of the mean state of the Middle Atmosphere, from the mesosphere down through the Upper Troposphere/Lower Stratosphere. Then we assess its use as a coupled atmosphere-ocean climate model by exploring its mean ocean climatology. To evaluate its variability, we report on its simulation of the primary modes in the troposphere, stratosphere and ocean. Two climate change simulations are presented, the responses to instantaneous increases of 2xCO2 and 4xCO2, run with two different ocean models. Sensitivity studies are performed to illustrate the effect of parameterizations on the model results. We compare these results to the lower vertical resolution/top GISS Model E2.1, whose output has also been submitted to CMIP6. The different choices made for these models are explored. It is shown that important improvements in the circulation above and below the tropopause can be obtained when attention is paid to representation of Middle Atmosphere processes in climate model development.

D. Rind↗

TPSAS-NF1676L-22800-DND

This lecture will provide an overview of an advanced turbulence closure scheme used in cloud-resolving models (CRMs) and climate modeling effort with the superparameterization approach. The CRM is embedded into each atmospheric column of a climate model as a super parameterization. It replaces the traditional subgrid-scale physical parameterizations in a conventional climate model. Two sensitivity simulations with doubling of CO2 and an increase of sea surface temperature, respectively, are performed and compared to similar simulations without this advanced turbulence closure in CRM to better understand the role of low-level clouds in climate sensitivity.

Kuan-Man Xu↗

Assessing the Impact of Soil Layer Depth Specification on the Observability of Modeled Soil Moisture and Brightness Temperature

The utility of hydrologic land surface models (LSMs) can be enhanced by using information from observational platforms, but mismatches between the two are common.This study assesses the degree to which model agreement with observations (observability) is affected by two mechanisms in particular: 1) physical incongruities between the support volumes being characterized and 2) inadequateor inconsistent parameterizations of physical processes. The Noah and Noah-MP LSMs by default characterize surface soil moisture (SSM) in the top 10 cm of the soil column. This depth may be reasonable when comparing against soil moisture from in situ probes centered at 5 cm,but it is notably different from the 5 cm (or less) sensing depth of NASA’s Soil Moisture Active Passive (SMAP) satellite mission. These depth inconsistencies are examined by using thinner model layers in the Noah and Noah-MP LSMs and comparing resultant simulations to in situ and SMAP soil moisture. In addition, a forward radiative transfer model to simulate microwave brightness temperatures (Tbs) is used to facilitate direct comparisons of LSM-based and SMAP-based L-band Tb retrievals. Observability is quantified using Kolmogorov-Smirnov distance values, calculated from empirical cumulative distribution functions of SSM and Tb time series. Experiment results depend on the particular subspace being analyzed (SSM or Tb). This study concludes that therole of increasedsoil layer discretizationon LSM observability is secondary to the influence of component parameterizations, the effects of which dominate systematic differences with observations

SMAP↗

Evaluating Subgrid-Scale Properties of Low Clouds over the Beaufort Sea in Arctic System Reanalysis using ARISE Airborne In Situ Observations

Arctic low clouds influence the evolution of the Arctic system through their effects on radiative fluxes, boundary layer mixing, stability, turbulence, humidity, and precipitation. Unfortunately, atmospheric models and retrospective analysis (reanalysis) products struggle to accurately simulate the occurrence and properties of low clouds in the Arctic. To better understand the magnitude and possible causes of this problem, we analyze the disagreements between the Arctic System Reanalysis version 2 (ASR) and data taken during the September 2014 Arctic Radiation-IceBridge Sea and Ice Experiment (ARISE) airborne campaign conducted over the Beaufort Sea. Previous analysis of grid scale-averaged meteorological variables revealed a lack of cloud water produced in ASR, related to warm/dry biases, a failure to produce sufficient cloud water in conditions that are observed to be favorable for large cloud water values, and sampling issues with the ARISE flight paths. However, ARISE collected meteorological observations at one second intervals, offering spatial data resolution comparable to the subgrid scale of ASR. Subgrid scale properties and processes must be parameterized in conventional atmospheric models and reanalyses such as ASR, and often rely on assumed distributions of meteorological variables about the grid box mean state. To enable an examination of the realism of ASR’s subgrid scale properties, we use high resolution data from ARISE determine the “subgrid scale” distributions of cloud water and thermodynamic variables in the real atmosphere. Preliminary results show that the subgrid scale distribution of thermodynamic variables about the grid scale mean resembles the triangular distribution assumed in many bulk microphysical parameterization schemes, but also with tails resembling Gaussian distributions. The main goal is to understand the partitioning of the total available atmospheric water into water vapor and cloud water. Models typically covert some percentage of total water in a grid box to cloud water based on the estimated amount of supersaturation from the assumed vapor distribution. The ARISE measurements can reveal the relationship between subgrid scale supersaturation and the partition between cloud water and water vapor.

Jason Brant Dodson↗

Clouds and Convective Self-Aggregation in a Multi-Model Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective56equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climatecsensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique amongst intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than un-aggregated simulations. Plain Language Summary This study investigates tropical clouds and climate using results from more than 30 different numerical models set up in a simplified framework. The dataset of model simulations is unique in that it includes a wide range of model types configured in a consistent manner. We address some of the biggest open questions in climate science, including how cloud properties change with warming and the role that the tendency of clouds to form clusters plays in determining the average climate and how climate changes. While there are large differences in how the different models simulate average temperature, humidity, and cloudiness, in a majority of models, the amount of high clouds decreases as climate warms. Nearly all models simulate a tendency for clouds to cluster together. There is agreement that when the clouds are clustered, the atmosphere is drier with fewer clouds overall. We don’t find a conclusive result for how cloud clustering changes as the climate warms.

Allison A. Wing↗

A comparison of global estimates of marine primary production from ocean color

The third primary production algorithm round robin (PPARR3) compares output from 24 models that estimate depth-integrated primary production from satellite measurements of ocean color, as well as seven general circulation models (GCMs) coupled with ecosystem or biogeochemical models. Here we compare the global primary production fields corresponding to eight months of 1998 and 1999 as estimated from common input fields of photosynthetically-available radiation (PAR), sea-surface temperature (SST), mixed-layer depth, and chlorophyll concentration. We also quantify the sensitivity of the ocean-color-based models to perturbations in their input variables. The pair-wise correlation between ocean-color models was used to cluster them into groups or related output, which reflect the regions and environmental conditions under which they respond differently. The groups do not follow model complexity with regards to wavelength or depth dependence, though they are related to the manner in which temperature is used to parameterize photosynthesis. Global average PP varies by a factor of two between models. The models diverged the most for the Southern Ocean, SST under 10 C, and chlorophyll concentration exceeding 1mg Chlm-3. Based on the conditions under which the model results diverge most, we conclude that current ocean-color-based models are challenged by high-nutrient low-chlorophyll conditions, and extreme temperatures or chlorophyll concentrations. The GCM-based models predict comparable primary production to those based on ocean color: they estimate higher values in the Southern Ocean, at low SST, and in the equatorial band, while they estimate lower values in eutrophic regions (probably because the area of high chlorophyll concentrations is smaller in the GCMs). Further progress in primary production modeling requires improved understanding of the effect of temperature on photosynthesis and better parameterization of the maximum photosynthetic rate.

Yamanaka, Yasuhiro↗

Uncertainty Quantification of Global Net Methane Emissions From Terrestrial Ecosystems Using a Mechanistically Based Biogeochemistry Model

Quantification of methane (CH4) emissions from wetlands and its sinks from uplands is still fraught with large uncertainties. Here, a methane biogeochemistry model was revised, parameterized, and verified for various wetland ecosystems across the globe. The model was then extrapolated to the global scale to quantify the uncertainty induced from four different types of uncertainty sources including parameterization, wetland type distribution, wetland area distribution, and meteorological input. We found that global wetland emissions are 212 ± 62 and 212 ± 32 Tg CH4 year−1 (1Tg = 1012 g) due to uncertain parameters and wetland type distribution, respectively, during 2000–2012. Using two wetland distribution data sets and three sets of climate data, the model simulations indicated that the global wetland emissions range from 186 to 212 CH4 year−1 for the same period. The parameters were the most significant uncertainty source. After combining the global methane consumption in the range of −34 to −46 Tg CH4 year−1, we estimated that the global net land methane emissions are 149–176 Tg CH4 year−1 due to uncertain wetland distribution and meteorological input. Spatially, the northeast United States and Amazon were two hotspots of methane emission, while consumption hotspots were in the Eastern United States and eastern China. During 1950–2016, both wetland emissions and upland consumption increased during El Niño events and decreased during La Niña events. This study highlights the need for more in situ methane flux data, more accurate wetland type, and area distribution information to better constrain the model uncertainty.

wetland methane emission↗

Machine Learning the COSMO Model for Predicting Thermodynamics of Electrolyte Mixtures

Bottom-up design of electrolyte mixtures for battery systems requires predicting macro thermodynamic properties from molecular constituents. For instance, molten salt electrolyte batteries require conditions far above room temperature to operate. Therefore, discovering mixtures with increasingly lower eutectic melting points is desirable. A model that can approximate chemical activity is a valuable tool to search through the vast compositional design space. Machine learning can predict properties of materials such as vibrational free energies, electronic energy gaps, and thermal conductivities. Moreover, they can learn physical models such as interatomic potentials. The COSMO-SAC model uses theory and empirical parameterization to predict liquid-vapor and liquid-solid properties using first-principles calculations. However, obtaining activity coefficients required for parameterizing the COSMO-SAC model is costly and limited to a select chemical space. In this work, we explored if machine learning methods could improve the COSMO-SAC model and bridge density functional theory calculations to liquid phase thermodynamic properties. Our data-driven approach uses existing databases for sigma-profiles of organic solvents and reconciles their methodological differences via ensemble averaging. First, an optimal machine learning model is constructed for each dataset. Our machine learning algorithms use the sigma-profile as an input feature to predict binary mixtures' activity coefficients using multi-output regression. Each dataset uses different choices of functionals, methods, and basis sets. Therefore, our ensemble model attempts to predict corrected activity coefficients given the combination of all the model outputs. The activity coefficients used for training are generated using the COSMO-SAC model. This approach enables the extraction of meaningful information from the existing datasets to improve the COSMO-SAC model for obtaining thermodynamic properties of electrolyte mixtures. With the liquid phase activities, we can identify electrolyte mixtures that meet desired phase equilibria conditions.

Thermodynamics↗

Comparison of high-resolution airborne soil moisture retrievals to SMAP soil moisture during the SMAP validation experiment 2016 (SMAPVEX16)

The NASA's Soil Moisture Active Passive (SMAP) mission conducted a field experiment with its partners over two 40-km agricultural domains in Iowa and Manitoba in the summer of 2016 to address concerns observed in SMAP soil moisture (SM) retrievals over agricultural areas. The experiment featured airborne Passive Active L-band System (PALS) flights over each domain with intensive ground measurements and dense networks of SM monitoring stations. With two intensive observation periods separated in time (May 28–June 20 and July 14–August 16), the flights captured both early-season/low vegetation and later-season/high-vegetation conditions. The comparison of the PALS brightness temperature (TB) measurements to the SMAP TB observed over the sites resulted in root mean square difference (RMSD) of 2.8 K and 4.0 K for vertical and horizontal polarizations, respectively. The subsequent SM analysis rescaled the PALS TB with the SMAP TB to allow equitable comparisons between the SM retrievals from the two instruments. The PALS SM retrieval algorithm used the SM sampled by the ground teams during the overpass days for tuning, and was parameterized by a high-resolution vegetation water content product calibrated using vegetation samples collected during the experiment. The tuning process was not able to find a satisfactory result with a temporally constant set of parameters in the single channel algorithm for the two intensive observation periods of the experiment. This result indicated that the rapid change in the vegetation structure during the growth stages and likely variation in the surface roughness conditions were not compatible with rigid parameterization over the entire period. However, using seasonally variable parameters we found that it was possible to retrieve soil moisture with satisfactory accuracy. Comparative analysis with the SMAP SM product included aggregation of the PALS SM to the SMAP pixel-scale. The RMSD between the PALS SM and the aggregated manual field samples was <0.04 m3/m3 with Pearson correlation >0.85 for both sites. The comparison between different in situ sources indicated that the soil moisture network measurements were not the source of the large biases observed for SMAP over the sites reported in earlier studies. Therefore, the results suggested the rapidly growing vegetation and the early-season surface condition changes not captured by the SMAP algorithm caused the SMAP retrieval errors. In addition, the significant deviations of the vegetation water content used by the SMAP product from the calibrated vegetation water content obtained during the experiment compounds the problem.

Andreas Colliander↗

Plant Phenology Evaluation of CRESCENDO Land Surface Models–Part 1: Start and End of the Growing Season

Plant phenology plays a fundamental role in land–atmosphere interactions, and its variability and variations are an indicator of climate and environmental changes. For this reason, current land surface models include phenology parameterizations and related biophysical and biogeochemical processes. In this work, the climatology of the beginning and end of the growing season, simulated by the land component of seven state-of-the-art European Earth system models participating in the CMIP6, is evaluated globally against satellite observations. The assessment is performed using the vegetation metric leaf area index and a recently developed approach, named four growing season types. On average, the land surface models show a 0.6-month delay in the growing season start, while they are about 0.5 months earlier in the growing season end. The difference with observation tends to be higher in the Southern Hemisphere compared to the Northern Hemisphere. High agreement between land surface models and observations is exhibited in areas dominated by broadleaf deciduous trees, while high variability is noted in regions dominated by broadleaf deciduous shrubs. Generally, the timing of the growing season end is accurately simulated in about 25 % of global land grid points versus 16 % in the timing of growing season start. The refinement of phenology parameterization can lead to better representation of vegetation-related energy, water, and carbon cycles in land surface models, but plant phenology is also affected by plant physiology and soil hydrology processes. Consequently, phenology representation and, in general, vegetation modelling is a complex task, which still needs further improvement, evaluation, and multi-model comparison.

Phenology↗