Conceptual Problems in Land Surface Data Assimilation

A land data assimilation system (LDAS) merges observations (or satellite retrievals) of land surface hydrological conditions, including soil moisture, snow, and terrestrial water storage (TWS), into a numerical model of land surface processes. In theory, the output from such a system is superior to estimates based on the observations or the model alone, thereby enhancing our ability to understand, monitor, and predict key elements of the terrestrial water cycle. In practice, however, several conceptual problems can interfere with realizing the potential improvements from data assimilation. Of particular concern is the frequent mismatch between the assimilated observations and the land surface model variables of interest. The seminar will discuss recent research with the ensemble-based NASA GEOS-S LDAS to address various aspects of this mismatch. These aspects include (i) the assimilation of coarse-scale observations into higher-resolution land surface models, (ii) the partitioning of satellite observations (such as TWS retrievals) into their constituent water cycle components, (iii) the forward modeling of microwave brightness temperatures over land for radiance-based land surface data aSSimilation, and (iv) the selection of the most relevant types of observations for the analysis of a specific water cycle variable (such as root zone soil moisture). At its core, the solution to the above challenges involves the careful construction of an observation operator that maps from the land surface model variables of interest to the space of the assimilated observations

Assessment of MERRA-2 Land Surface Energy Flux Estimates

In MERRA-2, observed precipitation is inserted in place of model-generated precipitation at the land surface. The use of observed precipitation was originally developed for MERRA-Land(a land-only replay of MERRA with model-generated precipitation replaced with observations).Previously shown that the land hydrology in MERRA-2 and MERRA-Land is better than MERRA. We test whether the improved land surface hydrology in MERRA-2 leads to the expected improvements in the land surface energy fluxes and 2 m air temperatures (T2m)

Spring Hydrology Determines Summer Net Carbon Uptake in Northern Ecosystems

Increased photosynthetic activity and enhanced seasonal CO2 exchange of northern ecosystems have been observed from a variety of sources including satellite vegetation indices (such as the Normalized Difference Vegetation Index; NDVI) and atmospheric CO2 measurements. Most of these changes have been attributed to strong warming trends in the northern high latitudes (greater than or equal to 50N). Here we analyze the interannual variation of summer net carbon uptake derived from atmospheric CO2 measurements and satellite NDVI in relation to surface meteorology from regional observational records. We find that increases in spring precipitation and snow pack promote summer net carbon uptake of northern ecosystems independent of air temperature effects. However, satellite NDVI measurements still show an overall benefit of summer photosynthetic activity from regional warming and limited impact of spring precipitation. This discrepancy is attributed to a similar response of photosynthesis and respiration to warming and thus reduced sensitivity of net ecosystem carbon uptake to temperature. Further analysis of boreal tower eddy covariance CO2 flux measurements indicates that summer net carbon uptake is positively correlated with early growing-season surface soil moisture, which is also strongly affected by spring precipitation and snow pack based on analysis of satellite soil moisture retrievals. This is attributed to strong regulation of spring hydrology on soil respiration in relatively wet boreal and arctic ecosystems. These results document the important role of spring hydrology in determining summer net carbon uptake and contrast with prevailing assumptions of dominant cold temperature limitations to high-latitude ecosystems. Our results indicate potentially stronger coupling of boreal/arctic water and carbon cycles with continued regional warming trends

Observation Error Specification for Near-Surface Soil Moisture Assimilation: Does It Matter?

The observation error variances assumed in the data assimilation (DA) of remotely sensed near-surface soil moisture observations are known to have unrealistic spatial patterns. The temporal variability in the errors is also typically neglected. Soil moisture errors cannot be directly calculated at large scales, and the globally distributed errors are unknown. We apply several approaches to estimating the observation errors (R matrix) for near-surface soil moisture DA experiments, to establish whether the assumed errors are important to the DA performance, and hence whether uncertainty in the observation error specification is a likely source of uncertainty in soil moisture DA output

Observation-Corrected Precipitation Estimates in GEOS-5

Several GEOS-5 applications, including the GEOS-5 seasonal forecasting system and the MERRA-Land data product, rely on global precipitation data that have been corrected with satellite and or gauge-based precipitation observations. This document describes the methodology used to generate the corrected precipitation estimates and their use in GEOS-5 applications. The corrected precipitation estimates are derived by disaggregating publicly available, observationally based, global precipitation products from daily or pentad totals to hourly accumulations using background precipitation estimates from the GEOS-5 atmospheric data assimilation system. Depending on the specific combination of the observational precipitation product and the GEOS-5 background estimates, the observational product may also be downscaled in space. The resulting corrected precipitation data product is at the finer temporal and spatial resolution of the GEOS-5 background and matches the observed precipitation at the coarser scale of the observational product, separately for each day (or pentad) and each grid cell

Contribution of Soil Moisture Information to Streamflow Prediction in the Snowmelt Season: A Continental-Scale Analysis

In areas dominated by winter snowcover, the prediction of streamflow during the snowmelt season may benefit from three pieces of information: (i) the accurate prediction of weather variability (precipitation, etc.) leading up to and during the snowmelt season, (ii) estimates of the amount of snow present during the winter season, and (iii) estimates of the amount of soil moisture underlying the snowpack during the winter season. The importance of accurate meteorological predictions and wintertime snow estimates is obvious. The contribution of soil moisture to streamflow prediction is more subtle yet potentially very important. If the soil is dry below the snowpack, a significant fraction of the snowmelt may be lost to streamflow and potential reservoir storage, since it may infiltrate the soil instead for later evaporation. Such evaporative losses are presumably smaller if the soil below the snowpack is wet. In this paper, we use a state-of-the-art land surface model to quantify the contribution of wintertime snow and soil moisture information -- both together and separately -- to skill in forecasting springtime streamflow. We find that soil moisture information indeed contributes significantly to streamflow prediction skill

Multiscale assimilation of Advanced Microwave Scanning Radiometer-EOS snow water equivalent and Moderate Resolution Imaging Spectroradiometer snow cover fraction observations in northern Colorado

Eight years (2002–2010) of Advanced Microwave Scanning Radiometer–EOS (AMSR-E) snow water equivalent (SWE) retrievals and Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) observations are assimilated separately or jointly into the Noah land surface model over a domain in Northern Colorado. A multiscale ensemble Kalman filter (EnKF) is used, supplemented with a rule-based update. The satellite data are either left unscaled or are scaled for anomaly assimilation. The results are validated against in situ observations at 14 high-elevation Snowpack Telemetry (SNOTEL) sites with typically deep snow and at 4 lower-elevation Cooperative Observer Program (COOP) sites. Assimilation of coarse-scale AMSR-E SWE and fine-scale MODIS SCF observations both result in realistic spatial SWE patterns. At COOP sites with shallow snowpacks, AMSR-E SWE and MODIS SCF data assimilation are beneficial separately, and joint SWE and SCF assimilation yields significantly improved root-mean-square error and correlation values for scaled and unscaled data assimilation. In areas of deep snow where the SNOTEL sites are located, however, AMSR-E retrievals are typically biased low and assimilation without prior scaling leads to degraded SWE estimates. Anomaly SWE assimilation could not improve the interannual SWE variations in the assimilation results because the AMSR-E retrievals lack realistic interannual variability in deep snowpacks. SCF assimilation has only a marginal impact at the SNOTEL locations because these sites experience extended periods of near-complete snow cover. Across all sites, SCF assimilation improves the timing of the onset of the snow season but without a net improvement of SWE amounts

Inferring Land Surface Model Parameters for the Assimilation of Satellite-Based L-Band Brightness Temperature Observations into a Soil Moisture Analysis System

The Soil Moisture and Ocean Salinity (SMOS) satellite mission provides global measurements of L-band brightness temperatures at horizontal and vertical polarization and a variety of incidence angles that are sensitive to moisture and temperature conditions in the top few centimeters of the soil. These L-band observations can therefore be assimilated into a land surface model to obtain surface and root zone soil moisture estimates. As part of the observation operator, such an assimilation system requires a radiative transfer model (RTM) that converts geophysical fields (including soil moisture and soil temperature) into modeled L-band brightness temperatures. At the global scale, the RTM parameters and the climatological soil moisture conditions are still poorly known. Using look-up tables from the literature to estimate the RTM parameters usually results in modeled L-band brightness temperatures that are strongly biased against the SMOS observations, with biases varying regionally and seasonally. Such biases must be addressed within the land data assimilation system. In this presentation, the estimation of the RTM parameters is discussed for the NASA GEOS-5 land data assimilation system, which is based on the ensemble Kalman filter (EnKF) and the Catchment land surface model. In the GEOS-5 land data assimilation system, soil moisture and brightness temperature biases are addressed in three stages. First, the global soil properties and soil hydraulic parameters that are used in the Catchment model were revised to minimize the bias in the modeled soil moisture, as verified against available in situ soil moisture measurements. Second, key parameters of the "tau-omega" RTM were calibrated prior to data assimilation using an objective function that minimizes the climatological differences between the modeled L-band brightness temperatures and the corresponding SMOS observations. Calibrated parameters include soil roughness parameters, vegetation structure parameters, and the single scattering albedo. After this climatological calibration, the modeling system can provide L-band brightness temperatures with a global mean absolute bias of less than 10K against SMOS observations, across multiple incidence angles and for horizontal and vertical polarization. Third, seasonal and regional variations in the residual biases are addressed by estimating the vegetation optical depth through state augmentation during the assimilation of the L-band brightness temperatures. This strategy, tested here with SMOS data, is part of the baseline approach for the Level 4 Surface and Root Zone Soil Moisture data product from the planned Soil Moisture Active Passive (SMAP) satellite mission

Prediction of Hydrological Drought: What Can We Learn From Continental-Scale Offline Simulations?

Land surface model experiments are used to quantify, across the coterminous United States, the contributions (isolated and combined) of soil moisture and snowpack initialization to the skill of seasonal streamflow forecasts at multiple leads and for different start dates. Forecasted streamflows are compared to naturalized streamflow observations where available and to synthetic (model-generated) streamflow data elsewhere. We find that snow initialization has a major impact on skill in the mountainous western U.S. and in a portion of the northern Great Plains; a mid-winter (January 1) initialization of snow in these areas leads to significant skill in the spring melting season. Soil moisture initialization also contributes to skill, and although the maximum contributions are not as large as those seen for snow initialization, the soil moisture contributions extend across a much broader geographical area. Soil moisture initialization can contribute to skill at long leads (up to 5 or 6 months), particularly for forecasts issued during winter

A Dynamic Approach to Addressing Observation-Minus-Forecast Mean Differences in a Land Surface Skin Temperature Data Assimilation System

In land data assimilation, bias in the observation-minus-forecast (O-F) residuals is typically removed from the observations prior to assimilation by rescaling the observations to have the same long-term mean (and higher-order moments) as the corresponding model forecasts. Such observation rescaling approaches require a long record of observed and forecast estimates, and an assumption that the O-F mean differences are stationary. A two-stage observation bias and state estimation filter is presented, as an alternative to observation rescaling that does not require a long data record or assume stationary O-F mean differences. The two-stage filter removes dynamic (nonstationary) estimates of the seasonal scale O-F mean difference from the assimilated observations, allowing the assimilation to correct the model for synoptic-scale errors without adverse effects from observation biases. The two-stage filter is demonstrated by assimilating geostationary skin temperature (Tsk) observations into the Catchment land surface model. Global maps of the O-F mean differences are presented, and the two-stage filter is evaluated for one year over the Americas. The two-stage filter effectively removed the Tsk O-F mean differences, for example the GOES-West O-F mean difference at 21:00 UTC was reduced from 5.1 K for a bias-blind assimilation to 0.3 K. Compared to independent in situ and remotely sensed Tsk observations, the two-stage assimilation reduced the unbiased Root Mean Square Difference (ubRMSD) of the modeled Tsk by 10 of the open-loop values