scholarly journals Impact of Land Surface Initialization Approach on Subseasonal Forecast Skill: A Regional Analysis in the Southern Hemisphere

2014 ◽  
Vol 15 (1) ◽  
pp. 300-319 ◽  
Author(s):  
Annette L. Hirsch ◽  
Jatin Kala ◽  
Andy J. Pitman ◽  
Claire Carouge ◽  
Jason P. Evans ◽  
...  

Abstract The authors use a sophisticated coupled land–atmosphere modeling system for a Southern Hemisphere subdomain centered over southeastern Australia to evaluate differences in simulation skill from two different land surface initialization approaches. The first approach uses equilibrated land surface states obtained from offline simulations of the land surface model, and the second uses land surface states obtained from reanalyses. The authors find that land surface initialization using prior offline simulations contribute to relative gains in subseasonal forecast skill. In particular, relative gains in forecast skill for temperature of 10%–20% within the first 30 days of the forecast can be attributed to the land surface initialization method using offline states. For precipitation there is no distinct preference for the land surface initialization method, with limited gains in forecast skill irrespective of the lead time. The authors evaluated the asymmetry between maximum and minimum temperatures and found that maximum temperatures had the largest gains in relative forecast skill, exceeding 20% in some regions. These results were statistically significant at the 98% confidence level at up to 60 days into the forecast period. For minimum temperature, using reanalyses to initialize the land surface contributed to relative gains in forecast skill, reaching 40% in parts of the domain that were statistically significant at the 98% confidence level. The contrasting impact of the land surface initialization method between maximum and minimum temperature was associated with different soil moisture coupling mechanisms. Therefore, land surface initialization from prior offline simulations does improve predictability for temperature, particularly maximum temperature, but with less obvious improvements for precipitation and minimum temperature over southeastern Australia.

2009 ◽  
Vol 48 (10) ◽  
pp. 2181-2196 ◽  
Author(s):  
R. Hamdi ◽  
A. Deckmyn ◽  
P. Termonia ◽  
G. R. Demarée ◽  
P. Baguis ◽  
...  

Abstract The authors examine the local impact of change in impervious surfaces in the Brussels capital region (BCR), Belgium, on trends in maximum, minimum, and mean temperatures between 1960 and 1999. Specifically, data are combined from remote sensing imagery and a land surface model including state-of-the-art urban parameterization—the Town Energy Balance scheme. To (i) isolate effects of urban growth on near-surface temperature independent of atmospheric circulations and (ii) be able to run the model over a very long period without any computational cost restrictions, the land surface model is run in a stand-alone mode coupled to downscaled 40-yr ECMWF reanalysis data. BCR was considered a lumped urban volume and the rate of urbanization was assessed by estimating the percentage of impervious surfaces from Landsat images acquired for various years. Model simulations show that (i) the annual mean urban bias (AMUB) on minimum temperature is rising at a higher rate (almost 3 times more) than on maximum temperature, with a linear trend of 0.14° and 0.05°C (10 yr)−1, respectively, (ii) the 40-yr AMUB on mean temperature is estimated to be 0.62°C, (iii) 45% of the overall warming trend is attributed to intensifying urban heat island effects rather than to changes in local–regional climate, and (iv) during summertime, a stronger dependence between the increase of urban bias on minimum temperature and the change in percentage of impervious surfaces is found.


2016 ◽  
Vol 18 (1) ◽  
pp. 85-108 ◽  
Author(s):  
Paul A. Dirmeyer ◽  
Subhadeep Halder

Abstract Retrospective forecasts from CFSv2 are evaluated in terms of three elements of land–atmosphere coupling at subseasonal to seasonal time scales: sensitivity of the atmosphere to variations in land surface states, the magnitude of variability of land states and fluxes, and the memory or persistence of land surface anomalies. The Northern Hemisphere spring and summer seasons are considered for the period 1982–2009. Ensembles are constructed from all available pairings of initial land and atmosphere/ocean states taken from the Climate Forecast System Reanalysis at the start of April, May, and June among the 28 years, so that the effect of initial land states on the evolving forecasts can be assessed. Finally, improvement and continuance of forecast skill derived from accurate land surface initialization is related to the three coupling elements. It is found that soil moisture memory is the most broadly important element for significant improvement of realistic land initialization on forecast skill. However, coupling strength manifested through the elements of sensitivity and variability are necessary to realize the potential predictability provided by memory of initial land surface anomalies. Even though there is clear responsiveness of surface heat fluxes, near-surface temperature, humidity, and daytime boundary layer development to variations in soil moisture over much of the globe, precipitation in CFSv2 is unresponsive. Failure to realize potential predictability from land surface states could be due to unfavorable atmospheric stability or circulation states; poor quality of what is considered realistic soil moisture analyses; and errors in the land surface model, atmospheric model, or their coupled interaction.


2004 ◽  
Vol 5 (6) ◽  
pp. 1049-1063 ◽  
Author(s):  
Randal D. Koster ◽  
Max J. Suarez ◽  
Ping Liu ◽  
Urszula Jambor ◽  
Aaron Berg ◽  
...  

Abstract Forcing a land surface model (LSM) offline with realistic global fields of precipitation, radiation, and near-surface meteorology produces realistic fields (within the context of the LSM) of soil moisture, temperature, and other land surface states. These fields can be used as initial conditions for precipitation and temperature forecasts with an atmospheric general circulation model (AGCM). Their usefulness is tested in this regard by performing retrospective 1-month forecasts (for May through September, 1979–93) with the NASA Global Modeling and Assimilation Office (GMAO) seasonal prediction system. The 75 separate forecasts provide an adequate statistical basis for quantifying improvements in forecast skill associated with land initialization. Evaluation of skill is focused on the Great Plains of North America, a region with both a reliable land initialization and an ability of soil moisture conditions to overwhelm atmospheric chaos in the evolution of the meteorological fields. The land initialization does cause a small but statistically significant improvement in precipitation and air temperature forecasts in this region. For precipitation, the increases in forecast skill appear strongest in May through July, whereas for air temperature, they are largest in August and September. The joint initialization of land and atmospheric variables is considered in a supplemental series of ensemble monthly forecasts. Potential predictability from atmospheric initialization dominates over that from land initialization during the first 2 weeks of the forecast, whereas during the final 2 weeks, the relative contributions from the two sources are of the same order. Both land and atmospheric initialization contribute independently to the actual skill of the monthly temperature forecast, with the greatest skill derived from the initialization of both. Land initialization appears to contribute the most to monthly precipitation forecast skill.


2008 ◽  
Vol 136 (6) ◽  
pp. 1971-1989 ◽  
Author(s):  
Keith M. Hines ◽  
David H. Bromwich

Abstract A polar-optimized version of the fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model (MM5) was developed to fill climate and synoptic needs of the polar science community and to achieve an improved regional performance. To continue the goal of enhanced polar mesoscale modeling, polar optimization should now be applied toward the state-of-the-art Weather Research and Forecasting (WRF) Model. Evaluations and optimizations are especially needed for the boundary layer parameterization, cloud physics, snow surface physics, and sea ice treatment. Testing and development work for Polar WRF begins with simulations for ice sheet surface conditions using a Greenland-area domain with 24-km resolution. The winter month December 2002 and the summer month June 2001 are simulated with WRF, version 2.1.1, in a series of 48-h integrations initialized daily at 0000 UTC. The results motivated several improvements to Polar WRF, especially to the Noah land surface model (LSM) and the snowpack treatment. Different physics packages for WRF are evaluated with December 2002 simulations that show variable forecast skill when verified with the automatic weather station observations. The WRF simulation with the combination of the modified Noah LSM, the Mellor–Yamada–Janjić boundary layer parameterization, and the WRF single-moment microphysics produced results that reach or exceed the success standards of a Polar MM5 simulation for December 2002. For summer simulations of June 2001, WRF simulates an improved surface energy balance, and shows forecast skill nearly equal to that of Polar MM5.


2014 ◽  
Vol 15 (1) ◽  
pp. 69-88 ◽  
Author(s):  
Randal D. Koster ◽  
Gregory K. Walker ◽  
Sarith P. P. Mahanama ◽  
Rolf H. Reichle

Abstract Offline simulations over the conterminous United States (CONUS) with a land surface model are used to address two issues relevant to the forecasting of large-scale seasonal streamflow: (i) the extent to which errors in soil moisture initialization degrade streamflow forecasts, and (ii) the extent to which a realistic increase in the spatial resolution of forecasted precipitation would improve streamflow forecasts. The addition of error to a soil moisture initialization field is found to lead to a nearly proportional reduction in large-scale seasonal streamflow forecast skill. The linearity of the response allows the determination of a lower bound for the increase in streamflow forecast skill achievable through improved soil moisture estimation, for example, through the assimilation of satellite-based soil moisture measurements. An increase in the resolution of precipitation is found to have an impact on large-scale seasonal streamflow forecasts only when evaporation variance is significant relative to precipitation variance. This condition is met only in the western half of the CONUS domain. Taken together, the two studies demonstrate the utility of a continental-scale land surface–modeling system as a tool for addressing the science of hydrological prediction.


2005 ◽  
Vol 6 (2) ◽  
pp. 146-155 ◽  
Author(s):  
M. Rodell ◽  
P. R. Houser ◽  
A. A. Berg ◽  
J. S. Famiglietti

Abstract Improper initialization of numerical models can cause spurious trends in the output, inviting erroneous interpretations of the earth system processes that one wishes to study. In particular, soil moisture memory is considerable, so that accurate initialization of this variable in land surface models (LSMs) is critical. The most commonly employed method for initializing an LSM is to spin up by looping through a single year repeatedly until a predefined equilibrium is achieved. The downside to this technique, when applied to continental- to global-scale simulations, is that regional annual anomalies in the meteorological forcing accumulate as artificial anomalies in the land surface states, including soil moisture. Nine alternative approaches were tested and compared using the Mosaic LSM and 15 yr of global meteorological forcing. Results indicate that the most efficient way to initialize an LSM, if possible and given that multiple years of preceding forcing are not available, is to use climatological average states from the same model for the precise time of year of initialization. Three other approaches were also determined to be preferable to the single-year spinup method. In addition, low-resolution spinup scenarios were devised and tested, and based on the results, an effective yet computationally economical technique is proposed.


2010 ◽  
Vol 11 (6) ◽  
pp. 1358-1372 ◽  
Author(s):  
Theodore J. Bohn ◽  
Mergia Y. Sonessa ◽  
Dennis P. Lettenmaier

Abstract Multimodel techniques have proven useful in improving forecast skill in many applications, including hydrology. Seasonal hydrologic forecasting in large basins represents a special case of hydrologic modeling, in which postprocessing techniques such as temporal aggregation and time-varying bias correction are often employed to improve forecast skill. To investigate the effects that these techniques have on the performance of multimodel averaging, the performance of three hydrological models [Variable Infiltration Capacity, Sacramento/Snow-17, and the Noah land surface model] and two multimodel averages [simple model average (SMA) and multiple linear regression (MLR) with monthly varying model weights] are examined in three snowmelt-dominated basins in the western United States. These evaluations were performed for both simulating and forecasting [using the Ensemble Streamflow Prediction (ESP) method] monthly discharge, with and without monthly bias corrections. The single best bias-corrected model outperformed the multimodel averages of raw models in both retrospective simulations and ensemble mean forecasts in terms of RMSE. Forming an MLR multimodel average from bias-corrected models added only slight improvements over the best bias-corrected model. Differences in performance among all bias-corrected models and multimodel averages were small. For ESP forecasts, both bias correction and multimodel averaging generally reduced the RMSE of the ESP ensemble means at lead times of up to 6 months in months when flow is dominated by snowmelt, with the reduction increasing as lead time decreased. The primary reason for this is that aggregating simulated streamflows from daily to monthly time scales increases model cross correlation, which in turn reduces the effectiveness of multimodel averaging in reducing those components of model error that bias correction cannot address. This effect may be stronger in snowmelt-dominated basins because the interannual variability of winter precipitation is a common input to all models. It was also found that both bias correcting and multimodel averaging using monthly varying parameters yielded much greater error reductions than methods using time-invariant parameters.


Author(s):  
Shaini Naha ◽  
Praveen K. Thakur ◽  
S. P. Aggarwal

The snow cover plays an important role in Himalayan region as it contributes a useful amount to the river discharge. So, besides estimating rainfall runoff, proper assessment of snowmelt runoff for efficient management and water resources planning is also required. A Land Surface Model, VIC (Variable Infiltration Capacity) is used at a high resolution grid size of 1 km. Beas river basin up to Thalot in North West Himalayas (NWH) have been selected as the study area. At first model setup is done and VIC has been run in its energy balance mode. The fluxes obtained from VIC has been routed to simulate the discharge for the time period of (2003-2006). Data Assimilation is done for the year 2006 and the techniques of Data Assimilation considered in this study are Direct Insertion (D.I) and Ensemble Kalman Filter (EnKF) that uses observations of snow covered area (SCA) to update hydrologic model states. The meteorological forcings were taken from 0.5 deg. resolution VIC global forcing data from 1979-2006 with daily maximum temperature, minimum temperature from Climate Research unit (CRU), rainfall from daily variability of NCEP and wind speed from NCEP-NCAR analysis as main inputs and Indian Meteorological Department (IMD) data of 0.25 °. NBSSLUP soil map and land use land cover map of ISRO-GBP project for year 2014 were used for generating the soil parameters and vegetation parameters respectively. The threshold temperature i.e. the minimum rain temperature is -0.5°C and maximum snow temperature is about +0.5°C at which VIC can generate snow fluxes. Hydrological simulations were done using both NCEP and IMD based meteorological Forcing datasets, but very few snow fluxes were obtained using IMD data met forcing, whereas NCEP based met forcing has given significantly better snow fluxes throughout the simulation years as the temperature resolution as given by IMD data is 0.5°C and rainfall resolution of 0.25°C. The simulated discharge has been validated using observed data from BBMB (Bhakra Beas Management Board) and coefficient of Correlation(R<sup>2</sup>) measured for (2003-2006) was 0.67 and 0.61 for the year 2006.But as VIC does not consider snowmelt runoff as a part of the total discharge, snowmelt runoff has been estimated for the simulation both with and without D.A. The snow fluxes as generated from VIC gives basin average estimates of Snow Cover, SWE, Snow Depth and Snow melt. It has been observed to be overestimated when model predicted snow cover is compared with MODIS SCA of 500 m resolution from MOD10A2 for each year. So MODIS 8-day snow cover area has been assimilated directly into the model state as well as by using EnKF after every 8 days for the year 2006.D.I Technique performed well as compared to EnKF. R<sup>2</sup> between Model SCA and MODIS SCA is estimated as 0.73 after D.I with Root Mean Square Error (RMSE) of +0.19. After direct Insertion of D.A, SCA has been reduced comparatively which resulted in 7% reduction of annual snowmelt contribution to total discharge.The assimilation of MODIS SCA data hence improved the snow cover area (SCA) fraction and finally updated other snow components.


2019 ◽  
Author(s):  
Clément Albergel ◽  
Yongjun Zheng ◽  
Bertrand Bonan ◽  
Emanuel Dutra ◽  
Nemesio Rodríguez-Fernández ◽  
...  

Abstract. This study demonstrates that LDAS-Monde, a global and offline Land Data Assimilation System (LDAS), that integrates satellite Earth observations into the ISBA (Interaction between Soil Biosphere and Atmosphere) Land Surface Model (LSM), is able to detect, monitor and forecast the impact of extreme weather on land surface states. LDAS-Monde jointly assimilates satellite derived Earth observations of surface soil moisture (SSM) and Leaf Area Index (LAI). It is run at global scale forced by ERA5 (LDAS_ERA5), the latest atmospheric reanalysis from the European Centre for Medium Range Weather Forecast (ECMWF) over 2010–2018 leading to a 9-yr, ~ 0.25° × 0.25° spatial resolution reanalysis of Land Surface Variables (LSVs). This reanalysis is then used to compute anomalies of land surface states, in order to (i) detect regions exposed to extreme weather such as droughts and heatwave events and (ii) address specific monitoring and forecasting requirements of LSVs for those regions. In this study, LDAS_ERA5 analysis is first successfully evaluated worldwide using several satellite-based datasets (SSM, LAI, Evapotranspiration, Gross Primary Production and Sun Induced Fluorescence), as well as in situ measurements (SSM, evapotranspiration and river discharge). The added value of assimilating the soil moisture and LAI is demonstrated with respect to a model simulation (openloop, with no assimilation). Since the global LDAS_ERA5 has relatively coarse resolution, two higher spatial resolution experiments over two areas particularly affected by heatwaves and/or droughts in 2018 were run: North Western Europe and the Murray-Darling basin in South Eastern Australia. These experiments were forced with ECMWF Integrated Forecasting System (IFS) high resolution operational analysis (LDAS_HRES, ~ 0.10° × 0.10° spatial resolution) over 2017–2018, and both openloop and analysis experiments compared once again. Since the IFS is a forecast system, it also allows LDAS-Monde to be used in forecast mode, and we demonstrate the added value of initializing 4- and 8-day LDAS-HRES forecasts of the LSVs, from the LDAS-HRES assimilation run, compared to the openloop experiments. This is particularly true for LAI that evolves on longer time space than SSM and is more sensitive to initial conditions than to atmospheric forcing, even at an 8-day lead time. This confirms that slowly evolving land initial conditions are paramount for forecasting LSVs and that LDAS-systems should jointly analyse both soil moisture and vegetation states. Finally evaluation of the modelled snowpack is presented and the perspectives for snow data assimilation in LDAS-Monde are discussed.


Author(s):  
Shaini Naha ◽  
Praveen K. Thakur ◽  
S. P. Aggarwal

The snow cover plays an important role in Himalayan region as it contributes a useful amount to the river discharge. So, besides estimating rainfall runoff, proper assessment of snowmelt runoff for efficient management and water resources planning is also required. A Land Surface Model, VIC (Variable Infiltration Capacity) is used at a high resolution grid size of 1 km. Beas river basin up to Thalot in North West Himalayas (NWH) have been selected as the study area. At first model setup is done and VIC has been run in its energy balance mode. The fluxes obtained from VIC has been routed to simulate the discharge for the time period of (2003-2006). Data Assimilation is done for the year 2006 and the techniques of Data Assimilation considered in this study are Direct Insertion (D.I) and Ensemble Kalman Filter (EnKF) that uses observations of snow covered area (SCA) to update hydrologic model states. The meteorological forcings were taken from 0.5 deg. resolution VIC global forcing data from 1979-2006 with daily maximum temperature, minimum temperature from Climate Research unit (CRU), rainfall from daily variability of NCEP and wind speed from NCEP-NCAR analysis as main inputs and Indian Meteorological Department (IMD) data of 0.25 °. NBSSLUP soil map and land use land cover map of ISRO-GBP project for year 2014 were used for generating the soil parameters and vegetation parameters respectively. The threshold temperature i.e. the minimum rain temperature is -0.5°C and maximum snow temperature is about +0.5°C at which VIC can generate snow fluxes. Hydrological simulations were done using both NCEP and IMD based meteorological Forcing datasets, but very few snow fluxes were obtained using IMD data met forcing, whereas NCEP based met forcing has given significantly better snow fluxes throughout the simulation years as the temperature resolution as given by IMD data is 0.5°C and rainfall resolution of 0.25°C. The simulated discharge has been validated using observed data from BBMB (Bhakra Beas Management Board) and coefficient of Correlation(R<sup>2</sup>) measured for (2003-2006) was 0.67 and 0.61 for the year 2006.But as VIC does not consider snowmelt runoff as a part of the total discharge, snowmelt runoff has been estimated for the simulation both with and without D.A. The snow fluxes as generated from VIC gives basin average estimates of Snow Cover, SWE, Snow Depth and Snow melt. It has been observed to be overestimated when model predicted snow cover is compared with MODIS SCA of 500 m resolution from MOD10A2 for each year. So MODIS 8-day snow cover area has been assimilated directly into the model state as well as by using EnKF after every 8 days for the year 2006.D.I Technique performed well as compared to EnKF. R<sup>2</sup> between Model SCA and MODIS SCA is estimated as 0.73 after D.I with Root Mean Square Error (RMSE) of +0.19. After direct Insertion of D.A, SCA has been reduced comparatively which resulted in 7% reduction of annual snowmelt contribution to total discharge.The assimilation of MODIS SCA data hence improved the snow cover area (SCA) fraction and finally updated other snow components.


Sign in / Sign up

Export Citation Format

Share Document