Sarith P. P. Mahanama

Assessment and Enhancement of MERRA Land Surface Hydrology Estimates

AbstractThe Modern-Era Retrospective Analysis for Research and Applications (MERRA) is a state-of-the-art reanalysis that provides, in addition to atmospheric fields, global estimates of soil moisture, latent heat flux, snow, and runoff for 1979–present. This study introduces a supplemental and improved set of land surface hydrological fields (“MERRA-Land”) generated by rerunning a revised version of the land component of the MERRA system. Specifically, the MERRA-Land estimates benefit from corrections to the precipitation forcing with the Global Precipitation Climatology Project pentad product (version 2.1) and from revised parameter values in the rainfall interception model, changes that effectively correct for known limitations in the MERRA surface meteorological forcings. The skill (defined as the correlation coefficient of the anomaly time series) in land surface hydrological fields from MERRA and MERRA-Land is assessed here against observations and compared to the skill of the state-of-the-art ECMWF...

The Rhône-Aggregation Land Surface Scheme Intercomparison Project: An Overview

The Rhone-Aggregation (Rhone-AGG) Land Surface Scheme (LSS) intercomparison project is an initiative within the Global Energy and Water Cycle Experiment (GEWEX)/Global Land-Atmosphere System Study (GLASS) panel of the World Climate Research Programme (WCRP). It is a intermediate step leading up to the next phase of the Global Soil Wetness Project (GSWP) (Phase 2), for which there will be a broader investigation of the aggregation between global scales (GSWP-1) and the river scale. This project makes use of the Rhone modeling system, which was developed in recent years by the French research community in order to study the continental water cycle on a regional scale. The main goals of this study are to investigate how 15 LSSs simulate the water balance for several annual cycles compared to data from a dense observation network consisting of daily discharge from over 145 gauges and daily snow depth from 24 sites, and to examine the impact of changing the spatial scale on the simulations. The overall evapotranspiration, runoff, and monthly change in water storage are similarly simulated by the LSSs, however, the differing partitioning among the fluxes results in very different river discharges and soil moisture equilibrium states. Subgrid runoff is especially important for discharge at the daily timescale and for smaller-scale basins. Also, models using an explicit treatment of the snowpack compared better with the observations than simpler composite schemes. Results from a series of scaling experiments are examined for which the spatial resolution of the computational grid is decreased to be consistent with large-scale atmospheric models. The impact of upscaling on the domain-averaged hydrological components is similar among most LSSs, with increased evaporation of water intercepted by the canopy and a decrease in surface runoff representing the most consistent inter-LSS responses. A significant finding is that the snow water equivalent is greatly reduced by upscaling in all LSSs but one that explicitly accounts for subgrid-scale orography effects on the atmospheric forcing.

Simulation of high-latitude hydrological processes in the Torne-Kalix basin: PILPS Phase 2(e) 1: Experiment description and summary intercomparisons

Abstract Twenty-one land-surface schemes (LSSs) participated in the Project for Intercomparison of Land-surface Parameterizations (PILPS) Phase 2(e) experiment, which used data from the Torne–Kalix Rivers in northern Scandinavia. Atmospheric forcing data (precipitation, air temperature, specific humidity, wind speed, downward shortwave and longwave radiation) for a 20-year period (1979–1998) were provided to the 21 participating modeling groups for 218 1/4° grid cells that represented the study domain. The first decade (1979–1988) of the period was used for model spin-up. The quality of meteorologic forcing variables is of particular concern in high-latitude experiments and the quality of the gridded dataset was assessed to the extent possible. The lack of sub-daily precipitation, underestimation of true precipitation and the necessity to estimate incoming solar radiation were the primary data concerns for this study. The results from two of the three types of runs are analyzed in this, the first of a three-part paper: (1) calibration–validation runs—calibration of model parameters using observed streamflow was allowed for two small catchments (570 and 1300 km2), and parameters were then transferred to two other catchments of roughly similar size (2600 and 1500 km2) to assess the ability of models to represent ungauged areas elsewhere; and 2) reruns—using revised forcing data (to resolve problems with apparent underestimation of solar radiation of approximately 36%, and certain other problems with surface wind in the original forcing data). Model results for the period 1989–1998 are used to evaluate the performance of the participating land-surface schemes in a context that allows exploration of their ability to capture key processes spatially. In general, the experiment demonstrated that many of the LSSs are able to capture the limitations imposed on annual latent heat by the small net radiation available in this high-latitude environment. Simulated annual average net radiation varied between 16 and 40 W/m2 for the 21 models, and latent heat varied between 18 and 36 W/m2. Among-model differences in winter latent heat due to the treatment of aerodynamic resistance appear to be at least as important as those attributable to the treatment of canopy interception. In many models, the small annual net radiation forced negative sensible heat on average, which varied among the models between −11 and 9 W/m2. Even though the largest evaporation rates occur in the summer (June, July and August), model-predicted snow sublimation in winter has proportionately more influence on differences in annual runoff volume among the models. A calibration experiment for four small sub-catchments of the Torne–Kalix basin showed that model parameters that are typically adjusted during calibration, those that control storage of moisture in the soil column or on the land surface via ponding, influence the seasonal distribution of runoff, but have relatively little impact on annual runoff ratios. Similarly, there was no relationship between annual runoff ratios and the proportion of surface and subsurface discharge for the basin as a whole.

Multimodel Ensemble Reconstruction of Drought over the Continental United States

Abstract Retrospectively simulated soil moisture from an ensemble of six land surface/hydrological models was used to reconstruct drought events over the continental United States for the period 1920–2003. The simulations were performed at one-half-degree spatial resolution, using a common set of atmospheric forcing data and model-specific soil and vegetation parameters. Monthly simulated soil moisture was converted to percentiles using Weibull plotting position statistics, and the percentiles were then used to represent drought severities and durations. An ensemble method, based on an inverse mapping of the average of the individual model’s soil moisture percentiles, was also used to combine all models’ simulations. Major results are 1) all models and the ensemble reconstruct the known severe drought events during the last century. The spatial extents and severities of drought are plausible for the individual models although substantial among-model disparities exist. 2) The simulations are in more agreem...

Simulation of high latitude hydrological processes in the Torne–Kalix basin: PILPS Phase 2(e): 2: Comparison of model results with observations

Model results from 21 land-surface schemes (LSSs) designed for use in numerical weather prediction and climate models are compared with each other and with observations in the context of the Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) Phase 2(e) model intercomparison experiment. This experiment focuses on simulations of land-surface water and energy fluxes in the 58,000-km2 Torne and Kalix river systems in northern Scandinavia, during the period 1989–1998. All models participating in PILPS Phase 2(e) capture the broad dynamics of snowmelt and runoff, but large differences in snow accumulation and ablation, turbulent heat fluxes, and streamflow exist. The greatest among-model differences in energy and moisture fluxes in these high-latitude environments occur during the spring snowmelt period, reflecting different model parameterizations of snow processes. Differences in net radiation are governed by differences in the simulated radiative surface temperature during the winter months and by differences in surface albedo during the spring/early summer. Differences in net radiation are smallest during the late summer when snow is absent. Although simulated snow sublimation is small for most models, a few models show annual snow sublimation of about 100 mm. These differences in snow sublimation appear to be largely dependent on differences in snow surface roughness parameterizations. The models with high sublimation generally lose their snowpacks too early compared to observations and underpredict the annual runoff. Differences in runoff parameterizations are reflected in differences in daily runoff statistics. Although most models show a greater variability in daily streamflow than the observations, the models with the greatest variability (as much as double the observed variability), produce most of their runoff through fast response, surface runoff mechanisms. As a group, those models that took advantage of an opportunity to calibrate to selected small catchments and to transfer calibration results to the basin at large had a smaller bias and root mean squared error (RMSE) in daily streamflow simulations compared with the models that did not calibrate.

Assimilation of Satellite-Derived Skin Temperature Observations into Land Surface Models

Abstract Land surface (or “skin”) temperature (LST) lies at the heart of the surface energy balance and is a key variable in weather and climate models. In this research LST retrievals from the International Satellite Cloud Climatology Project (ISCCP) are assimilated into the Noah land surface model and Catchment land surface model (CLSM) using an ensemble-based, offline land data assimilation system. LST is described very differently in the two models. A priori scaling and dynamic bias estimation approaches are applied because satellite and model LSTs typically exhibit different mean values and variabilities. Performance is measured against 27 months of in situ measurements from the Coordinated Energy and Water Cycle Observations Project at 48 stations. LST estimates from Noah and CLSM without data assimilation (“open loop”) are comparable to each other and superior to ISCCP retrievals. For LST, the RMSE values are 4.9 K (CLSM), 5.5 K (Noah), and 7.6 K (ISCCP), and the anomaly correlation coefficients (R...

Radiative forcing and temperature response to changes in urban albedos and associated CO 2 offsets

Radiative forcing and temperature response to changes in urban albedos and associated CO 2 offsets Surabi Menon 1 , Hashem Akbari 1* , Sarith Mahanama 2 , Igor Sednev 1 and Ronnen Levinson 1 Lawrence Berkeley National Laboratory, Berkeley, CA, USA Global Modeling and Assimilation Office, NASA Goddard Space Flight Center, Greenbelt, MD, USA *Now at the Dept. of Building, Civil and Environmental Engineering, Concordia University, Montreal, Quebec, Canada. Abstract The two main forcings that can counteract to some extent the positive forcings from greenhouse gases from pre-industrial times to present-day are the aerosol and related aerosol-cloud forcings, and the radiative response to changes in surface albedo. Here, we quantify the change in radiative forcing and land surface temperature that may be obtained by increasing the albedos of roofs and pavements in urban areas in temperate and tropical regions of the globe by 0.1. Using the catchment land surface model (the land model coupled to the GEOS-5 Atmospheric General Circulation Model), we quantify the change in the total outgoing (outgoing shortwave+longwave) radiation and land surface temperature to a 0.1 increase in urban albedos for all global land areas. The global average increase in the total outgoing radiation was 0.5 Wm - , and temperature decreased by ~0.008 K for an average 0.003 increase in surface albedo. These

Soil Moisture, Snow, and Seasonal Streamflow Forecasts in the United States

AbstractLand surface model experiments are used to quantify, for a number of U.S. river basins, 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. Snow initialization has a major impact on skill during the spring melting season. Soil moisture initialization has a smaller but still statistically significant impact during this season, and in other seasons, its contribution to skill dominates. Realistic soil moisture initialization can contribute to skill at long leads (over 6 months) for certain basins and seasons. Skill levels in all seasons are found to be related to the ratio of initial total water storage (soil water plus snow) variance to the forecast period precipitation variance, allowing estimates of the potential for skill in areas outside the verification basins.

Influence of the Interannual Variability of Vegetation on the Surface Energy Balance—A Global Sensitivity Study

Abstract The degree to which the interannual variability of vegetation phenology affects hydrological fluxes over land is investigated through a series of simulations with the Mosaic land surface model, run both offline and coupled to the NASA Seasonal-to-Interannual Prediction Project (NSIPP) atmospheric general circulation model (GCM). Over a 9-yr period, from 1982 to 1990, interannual variations of global biophysical land surface parameters (i.e., vegetation density and greenness fraction) are derived from Normalized Difference Vegetation Index (NDVI) data collected by the Advanced Very High Resolution Radiometers (AVHRRs). First the sensitivity of evapotranspiration to interannual variations in vegetation properties is evaluated through offline simulations that ignore feedbacks between the land surface and the atmospheric models, and interannual precipitation variations. Evapotranspiration is shown to be highly sensitive to variations in vegetation over wet continental surfaces that are not densely ve...

An updated treatment of soil texture and associated hydraulic properties in a global land modeling system

The advent of new data sets describing soil texture and associated soil properties offers the promise of improved hydrological simulation. Here we describe the composition of a new soil texture data set and its implementation into a specific land surface modeling system, namely, the Catchment land surface model (LSM) of the NASA Goddard Earth Observing System version 5 (GEOS-5) modeling and assimilation framework. First, global soil texture composites are generated using data from the Harmonized World Soil Database version 1.21 (HWSD1.21) and the State Soil Geographic (STATSGO2) project, with explicit consideration of different levels of organic material. Then, the LSMs soil parameters are upgraded using the new texture data, with hydraulic parameters derived for the more extensive set of texture classes using pedotransfer functions. Other changes to the LSM parameters are included to further support simulations at increasingly fine resolutions. A suite of simulations with the original and new parameter versions shows modest yet significant improvements in the Catchment LSMs simulation of soil moisture and surface hydrological fluxes. The revised LSM parameters will be used for the forthcoming Soil Moisture Active Passive (SMAP) soil moisture assimilation product.