Impacts of Soil Heating Condition on Precipitation Simulations in the Weather Research and Forecasting Model

2009 ◽  
Vol 137 (7) ◽  
pp. 2263-2285 ◽  
Author(s):  
Xingang Fan

Soil temperature is a major variable in land surface models, representing soil energy status, storage, and transfer. It serves as an important factor indicating the underlying surface heating condition for weather and climate forecasts. This study utilizes the Weather Research and Forecasting (WRF) model to study the impacts of changes to the surface heating condition, derived from soil temperature observations, on regional weather simulations. Large cold biases are found in the 40-yr European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis project (ERA-40) soil temperatures as compared to observations. At the same time, a warm bias is found in the lower boundary assumption adopted by the Noah land surface model. In six heavy rain cases studied herein, observed soil temperatures are used to initialize the land surface model and to provide a lower boundary condition at the bottom of the model soil layer. By analyzing the impacts from the incorporation of observed soil temperatures, the following major conclusions are drawn: 1) A consistent increase in the ground heat flux is found during the day, when the observed soil temperatures are used to correct the cold bias present in ERA-40. Soil temperature changes introduced at the initial time maintain positive values but gradually decrease in magnitude with time. Sensible and latent heat fluxes and the moisture flux experience an increase during the first 6 h. 2) An increase in soil temperature impacts the air temperature through surface exchange, and near-surface moisture through evaporation. During the first two days, an increase in air temperature is seen across the region from the surface up to about 800 hPa (∼1450 m). The maximum near-surface air temperature increase is found to be, averaged over all cases, 0.5 K on the first day and 0.3 K on the second day. 3) The strength of the low-level jet is affected by the changes described above and also by the consequent changes in horizontal gradients of pressure and thermal fields. Thus, the three-dimensional circulation is affected, in addition to changes seen in the humidity and thermal fields and the locations and intensities of precipitating systems. 4) Overall results indicate that the incorporation of observed soil temperatures introduces a persistent soil heating condition that is favorable to convective development and, consequently, improves the simulation of precipitation.

2020 ◽  
Vol 14 (8) ◽  
pp. 2581-2595 ◽  
Author(s):  
Bin Cao ◽  
Stephan Gruber ◽  
Donghai Zheng ◽  
Xin Li

Abstract. ERA5-Land (ERA5L) is a reanalysis product derived by running the land component of ERA5 at increased resolution. This study evaluates ERA5L soil temperature in permafrost regions based on observations and published permafrost products. We find that ERA5L overestimates soil temperature in northern Canada and Alaska but underestimates it in mid–low latitudes, leading to an average bias of −0.08 ∘C. The warm bias of ERA5L soil is stronger in winter than in other seasons. As calculated from its soil temperature, ERA5L overestimates active-layer thickness and underestimates near-surface (<1.89 m) permafrost area. This is thought to be due in part to the shallow soil column and coarse vertical discretization of the land surface model and to warmer simulated soil. The soil temperature bias in permafrost regions correlates well with the bias in air temperature and with maximum snow height. A review of the ERA5L snow parameterization and a simulation example both point to a low bias in ERA5L snow density as a possible cause for the warm bias in soil temperature. The apparent disagreement of station-based and areal evaluation techniques highlights challenges in our ability to test permafrost simulation models. While global reanalyses are important drivers for permafrost simulation, we conclude that ERA5L soil data are not well suited for informing permafrost research and decision making directly. To address this, future soil temperature products in reanalyses will require permafrost-specific alterations to their land surface models.


2009 ◽  
Vol 48 (7) ◽  
pp. 1362-1376 ◽  
Author(s):  
Jonathan E. Pleim ◽  
Robert Gilliam

Abstract The Pleim–Xiu land surface model (PX LSM) has been improved by the addition of a second indirect data assimilation scheme. The first, which was described previously, is a technique in which soil moisture is nudged according to the biases in 2-m air temperature and relative humidity between the model- and observation-based analyses. The new technique involves nudging the deep soil temperature in the soil temperature force–restore (FR) model according to model bias in 2-m air temperature only during nighttime. While the FR technique is computationally efficient and very accurate for the special conditions for which it was derived, it is very dependent on the deep soil temperature that drives the restoration term of the surface soil temperature equation. Thus, adjustment of the deep soil temperature to optimize the 2-m air temperature during the night, when surface forcing is minimal, provides significant advantages over other methods of deep soil moisture initialization. Simulations of the Weather Research and Forecasting Model (WRF) using the PX LSM with and without the new deep soil temperature nudging scheme demonstrate substantial benefits of the new scheme for reducing error and bias of the 2-m air temperature. The effects of the new nudging scheme are most pronounced in the winter (January 2006) during which the model’s cold bias is greatly reduced. Air temperature error and bias are also reduced in a summer simulation (August 2006) with the greatest benefits in less vegetated and more arid regions. Thus, the deep temperature nudging scheme complements the soil moisture nudging scheme because it is most effective for conditions in which the soil moisture scheme is least effective, that is, when evapotranspiration is not important (winter and arid climates).


2013 ◽  
Vol 26 (15) ◽  
pp. 5608-5623 ◽  
Author(s):  
Andrew G. Slater ◽  
David M. Lawrence

Abstract Permafrost is a characteristic aspect of the terrestrial Arctic and the fate of near-surface permafrost over the next century is likely to exert strong controls on Arctic hydrology and biogeochemistry. Using output from the fifth phase of the Coupled Model Intercomparison Project (CMIP5), the authors assess its ability to simulate present-day and future permafrost. Permafrost extent diagnosed directly from each climate model's soil temperature is a function of the modeled surface climate as well as the ability of the land surface model to represent permafrost physics. For each CMIP5 model these two effects are separated by using indirect estimators of permafrost driven by climatic indices and compared to permafrost extent directly diagnosed via soil temperatures. Several robust conclusions can be drawn from this analysis. Significant air temperature and snow depth biases exist in some model's climates, which degrade both directly and indirectly diagnosed permafrost conditions. The range of directly calculated present-day (1986–2005) permafrost area is extremely large (~4–25 × 106 km2). Several land models contain structural weaknesses that limit their skill in simulating cold region subsurface processes. The sensitivity of future permafrost extent to temperature change over the present-day observed permafrost region averages (1.67 ± 0.7) × 106 km2 °C−1 but is a function of the spatial and temporal distribution of climate change. Because of sizable differences in future climates for the representative concentration pathway (RCP) emission scenarios, a wide variety of future permafrost states is predicted by 2100. Conservatively, the models suggest that for RCP4.5, permafrost will retreat from the present-day discontinuous zone. Under RCP8.5, sustainable permafrost will be most probable only in the Canadian Archipelago, Russian Arctic coast, and east Siberian uplands.


2021 ◽  
Author(s):  
Bin Cao ◽  
Stephan Gruber ◽  
Donghai Zheng ◽  
Xin Li

&lt;div&gt; &lt;p&gt;ERA5 is the latest generation atmospheric reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ECMWF). ERA5-Land (ERA5L) is derived by running the land component of ERA5, Tiled ECMWF Scheme for Surface Exchanges over Land with a revised land surface hydrology (HTESSEL), at an increased resolution of 0.1&amp;#176;. This study evaluates ERA5L soil temperature in permafrost regions based on observations and published permafrost products. We find that ERA5L overestimates soil temperature in northern Canada and Alaska but underestimates it in mid&amp;#8211;low latitudes, leading to a near-zero overall bias (&amp;#8722;0.08 &amp;#730;C). The warm bias of ERA5L soil is more pronounced in winter than in other seasons. As calculated from its soil temperature, ERA5L overestimates active-layer thickness and underestimates near-surface (&lt; 1.89 m) permafrost area.This is thought to be due in part to the shallow soil column and coarse vertical discretization of the land surface model and to warmer simulated soil.&lt;/p&gt; &lt;p&gt;The soil temperature bias in permafrost regions correlates well with the bias in air temperature and with snow height. A review of the ERA5L snow parameterization in the code and a simulation example comparison with permafrost-specific processes rich model (GEOtop) both point to an error in snow metamorphism of HTESSEL leading to a low bias in ERA5L snow density as a possible cause for the warm bias in soil temperature. The apparent disagreement of station-based and areal evaluation techniques highlights challenges in our ability to test permafrost simulation models. While global reanalyses are important drivers for permafrost simulation, we conclude that ERA5L soil data are not well suited for informing permafrost research and decision making directly. To address this, future soil temperature products in reanalyses will require permafrost-specific alterations to their land surface models.&lt;/p&gt; &lt;p&gt;&lt;/p&gt;&lt;/div&gt;


2016 ◽  
Vol 144 (5) ◽  
pp. 1851-1865 ◽  
Author(s):  
Tatiana G. Smirnova ◽  
John M. Brown ◽  
Stanley G. Benjamin ◽  
Jaymes S. Kenyon

The land surface model (LSM) described in this manuscript was originally developed as part of the NOAA Rapid Update Cycle (RUC) model development effort; with ongoing modifications, it is now used as an option for the WRF community model. The RUC model and its WRF-based NOAA successor, the Rapid Refresh (RAP), are hourly updated and have an emphasis on short-range, near-surface forecasts including aviation-impact variables and preconvective environment. Therefore, coupling to this LSM (hereafter the RUC LSM) has been critical to provide more accurate lower boundary conditions. This paper describes changes made to the RUC LSM since earlier descriptions, including extension from six to nine levels, improved snow treatment, and new land-use data from MODIS. The RUC LSM became operational at the NOAA/National Centers for Environmental Prediction (NCEP) as part of the RUC from 1998–2012 and as part of the RAP from 2012 through the present. The simple treatments of basic land surface processes in the RUC LSM have proven to be physically robust and capable of realistically representing the evolution of soil moisture, soil temperature, and snow in cycled models. Extension of the RAP domain to encompass all of North America and adjacent high-latitude ocean areas necessitated further development of the RUC LSM for application in the tundra permafrost regions and over Arctic sea ice. Other modifications include refinements in the snow model and a more accurate specification of albedo, roughness length, and other surface properties. These recent modifications in the RUC LSM are described and evaluated in this paper.


2018 ◽  
Vol 11 (2) ◽  
pp. 541-560 ◽  
Author(s):  
Przemyslaw Zelazowski ◽  
Chris Huntingford ◽  
Lina M. Mercado ◽  
Nathalie Schaller

Abstract. Global circulation models (GCMs) are the best tool to understand climate change, as they attempt to represent all the important Earth system processes, including anthropogenic perturbation through fossil fuel burning. However, GCMs are computationally very expensive, which limits the number of simulations that can be made. Pattern scaling is an emulation technique that takes advantage of the fact that local and seasonal changes in surface climate are often approximately linear in the rate of warming over land and across the globe. This allows interpolation away from a limited number of available GCM simulations, to assess alternative future emissions scenarios. In this paper, we present a climate pattern-scaling set consisting of spatial climate change patterns along with parameters for an energy-balance model that calculates the amount of global warming. The set, available for download, is derived from 22 GCMs of the WCRP CMIP3 database, setting the basis for similar eventual pattern development for the CMIP5 and forthcoming CMIP6 ensemble. Critically, it extends the use of the IMOGEN (Integrated Model Of Global Effects of climatic aNomalies) framework to enable scanning across full uncertainty in GCMs for impact studies. Across models, the presented climate patterns represent consistent global mean trends, with a maximum of 4 (out of 22) GCMs exhibiting the opposite sign to the global trend per variable (relative humidity). The described new climate regimes are generally warmer, wetter (but with less snowfall), cloudier and windier, and have decreased relative humidity. Overall, when averaging individual performance across all variables, and without considering co-variance, the patterns explain one-third of regional change in decadal averages (mean percentage variance explained, PVE, 34.25±5.21), but the signal in some models exhibits much more linearity (e.g. MIROC3.2(hires): 41.53) than in others (GISS_ER: 22.67). The two most often considered variables, near-surface temperature and precipitation, have a PVE of 85.44±4.37 and 14.98±4.61, respectively. We also provide an example assessment of a terrestrial impact (changes in mean runoff) and compare projections by the IMOGEN system, which has one land surface model, against direct GCM outputs, which all have alternative representations of land functioning. The latter is noted as an additional source of uncertainty. Finally, current and potential future applications of the IMOGEN version 2.0 modelling system in the areas of ecosystem modelling and climate change impact assessment are presented and discussed.


2021 ◽  
Author(s):  
Sujeong Lim ◽  
Claudio Cassardo ◽  
Seon Ki Park

&lt;p&gt;The ensemble data assimilation system is beneficial to represent the initial uncertainties and flow-dependent background error covariance (BEC). In particular, the inevitable model uncertainties can be expressed by ensemble spread, that is the standard deviation of ensemble BEC. However, the ensemble spread generally suffers from under-estimated problems. To alleviate this problem, recent studies employed stochastic perturbation schemes to increases the ensemble spreads by adding the random forcing in the model tendencies (i.e., physical or dynamical tendencies) or parameterization schemes (i.e., PBL, convective scheme, etc.). In this study, we focus on the near-surface uncertainties which are affected by the interactions between the land and atmosphere process. The land surface model (LSM) provides various fluxes as the lower boundary condition to the atmosphere, influencing the accuracy of hourly-to-seasonal scale weather forecasting, but the surface uncertainties were not much addressed yet. In this study, we developed the stochastically perturbed parameterization (SPP) scheme for the Noah LSM. The Weather Research and Forecasting (WRF) ensemble system is used for regional weather forecasting over East Asia, especially over the Korean Peninsula. As a testbed experiment with the newly-developed Noah LSM-SPP system, we first perturbed the soil temperature &amp;#8212; a crucial variable for the near-surface forecasts by affecting sensible heat fluxes, land surface skin temperature and surface air temperature, and hence lower-tropospheric temperature. Here, the random forcing used in perturbation is made by the tuning parameters for amplitude, length scale, and time scales: they are commonly determined empirically by trial and error. In order to find optimal tuning parameter values, we applied a global optimization algorithm &amp;#8212; the micro-genetic algorithm (micro-GA) &amp;#8212; to achieve the smallest root-mean-squared errors. Our results indicate that optimization of the random forcing parameters contributes to an increase in the ensemble spread and a decrease in the ensemble mean errors in the near-surface and lower-troposphere uncertainties. Further experiments will be conducted by including soil moisture in the testbed.&lt;/p&gt;


2013 ◽  
Vol 10 (7) ◽  
pp. 4465-4479 ◽  
Author(s):  
K. L. Hanis ◽  
M. Tenuta ◽  
B. D. Amiro ◽  
T. N. Papakyriakou

Abstract. Ecosystem-scale methane (CH4) flux (FCH4) over a subarctic fen at Churchill, Manitoba, Canada was measured to understand the magnitude of emissions during spring and fall shoulder seasons, and the growing season in relation to physical and biological conditions. FCH4 was measured using eddy covariance with a closed-path analyser in four years (2008–2011). Cumulative measured annual FCH4 (shoulder plus growing seasons) ranged from 3.0 to 9.6 g CH4 m−2 yr−1 among the four study years, with a mean of 6.5 to 7.1 g CH4 m−2 yr−1 depending upon gap-filling method. Soil temperatures to depths of 50 cm and air temperature were highly correlated with FCH4, with near-surface soil temperature at 5 cm most correlated across spring, fall, and the shoulder and growing seasons. The response of FCH4 to soil temperature at the 5 cm depth and air temperature was more than double in spring to that of fall. Emission episodes were generally not observed during spring thaw. Growing season emissions also depended upon soil and air temperatures but the water table also exerted influence, with FCH4 highest when water was 2–13 cm below and lowest when it was at or above the mean peat surface.


2014 ◽  
Vol 18 (5) ◽  
pp. 1761-1783 ◽  
Author(s):  
O. Branch ◽  
K. Warrach-Sagi ◽  
V. Wulfmeyer ◽  
S. Cohen

Abstract. A 10 × 10 km irrigated biomass plantation was simulated in an arid region of Israel to simulate diurnal energy balances during the summer of 2012 (JJA). The goal is to examine daytime horizontal flux gradients between plantation and desert. Simulations were carried out within the coupled WRF-NOAH atmosphere/land surface model. MODIS land surface data was adjusted by prescribing tailored land surface and soil/plant parameters, and by adding a controllable sub-surface irrigation scheme to NOAH. Two model cases studies were compared – Impact and Control. Impact simulates the irrigated plantation. Control simulates the existing land surface, where the predominant land surface is bare desert soil. Central to the study is parameter validation against land surface observations from a desert site and from a 400 ha Simmondsia chinensis (jojoba) plantation. Control was validated with desert observations, and Impact with Jojoba observations. Model evapotranspiration was validated with two Penman–Monteith estimates based on the observations. Control simulates daytime desert conditions with a maximum deviation for surface 2 m air temperatures (T2) of 0.2 °C, vapour pressure deficit (VPD) of 0.25 hPa, wind speed (U) of 0.5 m s−1, surface radiation (Rn) of 25 W m−2, soil heat flux (G) of 30 W m−2 and 5 cm soil temperatures (ST5) of 1.5 °C. Impact simulates irrigated vegetation conditions with a maximum deviation for T2 of 1–1.5 °C, VPD of 0.5 hPa, U of 0.5 m s−1, Rn of 50 W m−5, G of 40 W m−2 and ST5 of 2 °C. Latent heat curves in Impact correspond closely with Penman–Monteith estimates, and magnitudes of 160 W m−2 over the plantation are usual. Sensible heat fluxes, are around 450 W m−2 and are at least 100–110 W m−2 higher than the surrounding desert. This surplus is driven by reduced albedo and high surface resistance, and demonstrates that high evaporation rates may not occur over Jojoba if irrigation is optimized. Furthermore, increased daytime T2 over plantations highlight the need for hourly as well as daily mean statistics. Daily mean statistics alone may imply an overall cooling effect due to surplus nocturnal cooling, when in fact a daytime warming effect is observed.


2010 ◽  
Vol 2 (2) ◽  
Author(s):  
Diandong Ren

AbstractBased on a 2-layer land surface model, a rather general variational data assimilation framework for estimating model state variables is developed. The method minimizes the error of surface soil temperature predictions subject to constraints imposed by the prediction model. Retrieval experiments for soil prognostic variables are performed and the results verified against model simulated data as well as real observations for the Oklahoma Atmospheric Surface layer Instrumentation System (OASIS). The optimization scheme is robust with respect to a wide range of initial guess errors in surface soil temperature (as large as 30 K) and deep soil moisture (within the range between wilting point and saturation). When assimilating OASIS data, the scheme can reduce the initial guess error by more than 90%, while for Observing Simulation System Experiments (OSSEs), the initial guess error is usually reduced by over four orders of magnitude.Using synthetic data, the robustness of the retrieval scheme as related to information content of the data and the physical meaning of the adjoint variables and their use in sensitivity studies are investigated. Through sensitivity analysis, it is confirmed that the vegetation coverage and growth condition determine whether or not the optimally estimated initial soil moisture condition leads to an optimal estimation of the surface fluxes. This reconciles two recent studies.With the real data experiments, it is shown that observations during the daytime period are the most effective for the retrieval. Longer assimilation windows result in more accurate initial condition retrieval, underlining the importance of information quantity, especially for schemes assimilating noisy observations.


Sign in / Sign up

Export Citation Format

Share Document