Next Article in Journal
Security Enhancement Mechanism Based on Contextual Authentication and Role Analysis for 2G-RFID Systems
Next Article in Special Issue
Mapping Wide Row Crops with Video Sequences Acquired from a Tractor Moving at Treatment Speed
Previous Article in Journal
DNA Strand Patterns on Aluminium Thin Films
Previous Article in Special Issue
Geosensors to Support Crop Production: Current Applications and User Requirements
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of Sensor Imagery Data for Surface Boundary Conditions in Regional Climate Modeling

Department of Civil Engineering, Yeungnam University, 214-1 Dae-dong, Gyeongsan-si, Gyeongbuk-do, 712-749, Korea
Submission received: 25 March 2011 / Revised: 22 June 2011 / Accepted: 27 June 2011 / Published: 28 June 2011
(This article belongs to the Special Issue Sensors in Agriculture and Forestry)

Abstract

: Mesoscale climate and hydrology modeling studies have increased in sophistication and are being run at increasingly higher resolutions. Data resolution sufficiently finer than that of the computational model is required not only to support sophisticated linkages and process interactions at small scales but to assess their cumulative impact at larger scales. The global distributions at fine spatial and temporal scales can be described by means of various senor imagery data collected through remote sensing techniques, sensor image and photo programs, scanning and digitizing skills for existing maps, etc. The availability of global sensor imagery maps facilitates assimilation in land surface models to account for terrestrial dynamics. This study focuses on the use of global imagery data for development and construction of surface boundary conditions (SBCs) specifically designed for mesoscale regional climate model (RCM) applications. The several SBCs are currently presented in a RCM domain for the continent of Asia at 30-km spacing by using sensor imagery data. Geographic Information System (GIS) software application tools are mainly used to convert data information from various raw data onto RCM-specific grids. The raw data sources and processing procedures are elaborated in detail, by which the SBCs can be readily constructed for any specific RCM domain anywhere in the world.
Keywords:
sensor imagery; RCM; SBC; GIS

1. Introduction

The Fourth Assessment Report of the United Nations Intergovernmental Panel on Climate Change (IPCC AR4) [1] has brought to the fore the severity of impacts of global climate changes. In Asia, at its current rate climate change is expected to put close to 50 million extra people at risk of hunger by 2020, rising to an additional 130 million in 2050. Glaciers in the Himalayas could, at a similar rate, disappear altogether by 2035. Further, the IPCC AR4 [1] notes that Europe’s vast reaches of low lying coast are vulnerable to rising sea levels likely to threaten up to 2.5 million people. As global warming and extreme weather pose a severe threat to the safety of life and property over the whole world, climate information from regional climate models (RCMs) has become regarded as a primary tool to address climate and weather variability, changes, and impacts at local and regional scales. Since numerous RCMs have been developed and applied, the next-generation Weather Research and Forecasting (WRF) model ( http://www.wrf-model.org/) was developed by a broad community of government and university researchers [24]. The WRF was based upon the most advanced supercomputing technologies and promises greater efficiency in computation and flexibility in new module incorporation. The climate extension of the WRF (CWRF) has been developed by the Illinois State Water Survey in collaboration with the WRF Working Groups for incorporating inclusively all WRF functionalities for numerical weather prediction (NWP) while enhancing the capability for climate applications [5].

For all RCMs, one essential component is the representation of surface-atmosphere interactions, which generally requires specification of surface boundary conditions (SBCs). A comprehensive set of SBCs based on best observational data is desired for RCM general applications to all effective, dynamically coupled or uncoupled, combinations of the surface modules as well as portable to any specific region of the world. A critical requirement in constructing the SBCs for the RCM use is that each field must be globally defined with no missing value and physical consistency must be maintained across all relevant parameters. Missing data, if any, must be appropriately filled. The mesoscale weather and climate models, both global and regional, have increased in sophistication and are being run at increasingly higher resolutions. Hence, the raw data should be available at the finest possible resolution and global observation for assimilation in land surface models to improve predictability. The global distributions at fine spatial and temporal scales can be available by various sensor imagery data such as remote sensing observations, sensor photography images, scanned and digitized maps, and so on. This will facilitate a more realistic representation of surface heterogeneity effects. When the data resolution is sufficiently finer than the RCM grid, the subgrid effects can be further incorporated using composite, mosaic or statistical-dynamical approaches [610]. With a large volume of global imagery data of the Earth’s terrestrial surface becoming available, precisely monitoring the dynamics of the land surface state variables becomes possible [11].

There is no universal, complete set of SBCs because different modules may require more or less surface parameters to be specified. The input parameter requirement generally depends on the formulation complexity of the surface modules. The CWRF have incorporated into the Common Land Model (CLM) [12], a state-of-the-art model for Soil-Vegetation-Atmosphere Transfer (SVAT). For construction of the primary SBCs in CLM that play an important role in surface-atmosphere interaction, this study has collected many raw datasets at high resolution with global sensor imagery observations freely available. Among the static SBCs, fundamental input fields, independent or defining other derivatives, are orography and vegetation parameters. The surface orography parameters that examine subgrid topography effects on momentum and radiation dynamics include the mean and standard deviation of the terrain elevation at both subgrid-scale and resolved-scale (SOAVE, SOSTD, ROAVE, and ROSTD), and the anisotropy parameter and the angle of the maximum gradient of resolved-scale orography (ROANI and ROANG). The vegetation fields consist of the land cover category (LCC), the fractional vegetation cover (FVC), and the leaf area index (LAI) to determine contribution partitioning between bare soil and vegetation for fluxes crucial to land-atmosphere interactions.

The existing sensor observational databases have various resolutions, a wide range of map projections, and different data formats, and often contain missing values or inconsistencies between variables. This presents significant challenges and requires labor-intensive efforts to process the data onto the RCM-specific grid mesh and input data format. This study employs the Geographic Information System (GIS) software application tools, ArcInfo and ArcMap, from Environmental Systems Research Institute, Inc., particularly to determine the geographic conversion information from a specific map of each raw data to the identical RCM grid system. The information includes location indices, geometric distances, and fractional areas of all input cells contributing to each RCM grid. The grid representative mean values are usually determined by a bilinear interpolation method in terms of the geometric distances and a mass conservative approach as weighted by the fractional areas. For LCC, a categorical field, this study first calculates the total fractional area of each distinct surface category contributing to a given RCM grid and then chooses the dominant one that occupies the largest fraction of the grid.

2. General Information

The SBCs data quality, value representation, and visual display largely depend on the RCM computational domain and grid resolution. For climate applications of the Asia region, the domain is centered at 30.0°N and 110.0°W using the Lambert Conformal Conic map projection and 30-km horizontal grid spacing, with total grid numbers of 301 (west-east) × 251 (south–north). Figure 1(a) shows the country map overlaid with latitude and longitude lines projected on the Lambert Conformal Conic map, and Figure 1(b) denotes the 301 × 251 dimensional 30-km spacing grids over the Asia RCM domain. The domain covers the most parts of Asia and represents the regional climate that results from interactions between the planetary circulation and surface processes, including orography and vegetation fields. In this study, the fundamental SBCs using sensor imagery observations are constructed and displayed on this RCM domain, suitable for Asia applications.

For convenience, the geographic location of a point is hereafter referred as a “pixel” for raw data and a “grid” for the RCM result. A given pixel or grid value represents the area surrounding the point as defined by its respective horizontal spacing. The following section elaborates in detail on raw data sources and processing procedures used to construct any specific RCM domain over the globe. Most procedures use ArcInfo and ArcMap commands. In particular, IMAGEGRID and GRIDPOLY convert input data from the sensor image to the ArcGIS raster grid and to the polygon coverage formats, respectively; PROJECT remaps the raw input data onto the RCM grid projection; UNION and CLIP geometrically intersect polygon features of input data with the RCM grid mesh and extract the fractional area of each pixel contributing to the grid; GRID DOCELL and IF statements conditionally merge, replace, or adjust different input datasets for an improved product. Even relatively finer resolution (8-km) input data, LAI, has missing value pixels due to cloud contamination and atmospheric effects, which are filled by the spatial average over the nearby data pixels having the same land cover category within a certain radius around a missing point. The number of pixels and the range of radius used for filling depend on the resolution of the raw input data. Since the 30-km LAI for the RCM are generated from the two different resolution data of the 1-km FVC and the 8-km LAI, the 1-km data are integrated onto the 8-km map and then a smoothing filter is introduced to remove abnormal values due to inconsistency of the two data at individual pixels (see Section 3.4 for details).

3. Surface Boundary Conditions

The various surface orography parameters (SOAVE, SOSTD, ROAVE, ROSTD, ROANI, and ROANG) have been introduced for use of analysis on orographic turbulence effects under stable atmospheric conditions by using the terrain elevation model through the multi sensor imagery dataset. The vegetation characteristic fields (LCC, FVC, and LAI) for use of prediction on fluxes crucial to land-atmosphere interactions have been constructed by using satellite remote sensing data. The details about the raw data sources and processing procedures for each SBCs for RCM uses in the Asia domain are discussed below.

3.1. Surface Orography Parameters

These fields are constructed from the U.S. Geological Survey (USGS) HYDRO1k Digital Elevation Model (DEM) with a 1-km nominal cell size ( http://eros.usgs.gov/#/Find_Data/Products_and_Data_Available/HYDRO1K), which is developed at the U.S. Geological Survey’s Earth Resources Observation and Science (EROS) Data Center to provide to users, on a continent by continent basis, hydrologically correct DEMs for use in continental and regional scale modeling and analyses. It is based on the Global 30-arc-second elevation dataset (GTOPO30) derived from several different raster and vector sources of topographic information. The GTOPO30 data were produced jointly by several national and international organizations. National Imagery and Mapping Agency (NIMA) produced Digital Terrain Elevation Data (DTED) which was used for the source of most parts of Eurasia in GTOPO30. DTED is a raster topographic data collected through remote sensing techniques, aerial image sensor photography, digitization of interpolated contour lines on existing maps, and actual on-site surveying measurements. The DEM data is available in the band interleaved by line (BIL) image format on the Lambert Azimuthal equal area projection. For construction of surface orography parameters, the raw imagery data are converted into ArcGIS raster grid format and then remapped onto the RCM projection. Subsequently, the ArcInfo/GRID commands, ZONEALMEAN and ZONALSTD, are used to calculate the mean and standard deviation of the subgrid-scale elevations, SOAVE and SOSTD, respectively, within each RCM grid. The mean and deviation from the centroid of each grid are picked up by ArcInfo’s Arc Macro Language (AML) program Gridspot70 for SOAVE and SOSTD.

Given a grid spacing Δx, the subgrid variability effects of topography are currently not accounted in most RCMs. Subgrid orography effects need to be incorporated in RCMs through certain parameterizations especially in mesoscale global models. The values of the surface orography parameters depend crucially on the resolution of the raw topography data derived for a target model grid that has a lower resolution than the source data pixel. It is necessary to filter the subgrid-scale orography variables in a proper scale to avoid numerical noise for the model integrations. For the numerical stability the terrain height is filtered to remove 4-grid waves, and the subgrid effects are denoted based on the mean terrain averaged over an area of (4Δx)2 surrounding each grid center. Thus the resolved mesoscale orographic parameterizations are here calculated from the terrain elevation (hi in meters) by the HYDRO1k DEM with a 1-km nominal cell size. Following Rontu [13], the resolved resolved-scale mean elevation, ROAVE (H4Δx in meters), is calculated for 4 grid-lengths (4Δx) using double smoothing method:

H 4 Δ x = 1 N 4 Δ x i = 1 N 4 Δ x h 9 i
where h 9 i = 1 9 i = 1 9 h i (in meters) is 9-point averaged elevation including the point and its 8 neighbor points based on 1-km DEM (hi), Δx is the RCM grid-length which is 30 km for the current model, and N4Δx is the number of 1-km DEM pixels within a 4Δx-resolution grid. As defined in Rontu [13], the standard deviation of resolved-scale elevations, ROSTD (σm in meters), for the difference between h9i and H4Δx is calculated as:
σ m = 1 N Δ x 1 i = 1 N Δ x ζ 2
where ζ = h9iH4Δx is h9i departure from H4Δx, and NΔx is the number of 1-km data points for each Δx-resolution grid. Figure 2 compares geographic distributions of mean terrain elevations, and standard deviations of terrain elevations at subgrid-scale and resolved-scale.

Other parameters for resolved-scale orography are anisotropy, ROANI (α), characterizing the anisotropy of orography, and direction angle, ROANG (θ), representing the mean slope within each grid. The parameters ROANI (α) and ROANG (θ) are based on the orographic gradient correlation tensor by Lott and Miller [14]:

H i j = ζ x i ζ x j ¯
where xi and xj are the two principle coordinates of the model grids. Note that the parameters with overbar (hereafter) are model grid average values. Following Rontu et al. [15], the two intermediate parameters are defined respectively as:
T = ζ x i ζ x i ¯ + ζ x j ζ x j ¯
D = ( ζ x i ζ x i ¯ ζ x j ζ x j ¯ ) 2 + 4 ( ζ x i ζ x j ¯ ) 2
Thus, the anisotropy ROANI (α) is finally calculated as:
α = | T D | | T + D |

The anisotropy is zero for a two-dimensional ridge and increases for circular-shaped mountains. The direction angle ROANG (θ in radian) between the maximum gradient of resolved-scale orography and the xi -axis of the model grid is defined as:

θ = arctan ( ζ x j ζ x j ¯ ζ x i ζ x i ¯ + D 2 ζ x i ζ x j ¯ )

Figure 3 illustrates the geographic distributions of anisotropy and direction angle over the Asia RCM domain.

3.2. Land Cover Category

The RCM uses the 24-category USGS land cover classification (see Table 1) developed from the April 1992–March 1993 Advanced Very High Resolution Radiometer (AVHRR) satellite-derived Normalized Difference Vegetation Index (NDVI) composites. This data is based on a flexible data base structure and seasonal land cover regions concepts. The regions are composed of relatively homogeneous land cover associations which exhibit distinctive phenology, and have common levels of primary production. The raw data are available at 1-km spacing on the geographic coordinate system in BIL image format ( http://edc2.usgs.gov/glcc/glcc.php), converted into the ArcGIS raster grid and polygon coverage, and remapped onto the RCM projection. The fractional area of each pixel contributing to the grid is extracted after the result is intersected with the RCM grid mesh. The contributing area for each of the 24 LCCs is summed over all pixels of the same category within each RCM grid. The category contributing the largest area is chosen as the LCC for the grid. When the fractional area of water bodies is less than 0.5 but dominates the grid, the category chosen is the one contributing the second largest area.

Figure 4 illustrates the LCC geographic distribution over the RCM domain. Note that the USGS raw data do not contain LCC types 4 (mixed dryland/irrigated cropland & pasture) and 20 (herbaceous tundra) over the globe, and additionally LCC type 23 (bare ground tundra) is not chosen for LCC majority category in the Asia RCM domain. Therefore, the final LCC includes only 21 categories of LCCs over the present RCM domain.

3.3. Fractional Vegetation Cover

The FVC is one ecological parameter that determines contribution partitioning between bare soil and vegetation for surface evapotranspiration, photosynthesis, albedo, and other fluxes crucial to land-atmosphere interactions. Following Zeng et al. [16,17], the time-invariant FVC is derived from the 10-day April 1992–March 1993 composites of the global 1-km AVHRR NDVI product. The annual maximum NDVI (Np,max) for each LCC are chosen to minimize the effect of cloud contamination on data quality. For each pixel, the vegetation cover is computed by:

C v = N p , max N s N c , v N s
where Nc,v is the NDVI value for a complete coverage of a specific USGS LCC over the pixel and Ns for bare soil.

Using a commercial imagery database, Zeng et al. [16] determined Nc,v by examining percentiles of the Np,max histogram for each LCC of the International Geosphere Biosphere Programme (IGBP) classification [18,19]. Liang et al. [20] calculated the Nc,v values for the 24 USGS LCCs from those of the 17 IGBP categories by intersecting the USGS and IGBP land cover maps and computing the fractional areas of individual IGBP categories contributing to each USGS category. The final Nc,v is the average of all contributing IGBP values weighted by their corresponding fractional areas. Corresponding Nc,v values and occupation rates for the USGS and IGBP categories are listed in Table 1.

The resultant Cv point data at 1-km spacing are converted to polygon coverage data, remapped onto the RCM projection, and intersected with the RCM grid mesh. The fractional area of each pixel contributing to the grid is extracted. The final FVC is obtained by the area-weighted averaging of Cv values for all pixels within each RCM grid. Figure 5 illustrates the FVC geographic distributions derived from the AVHRR NDVI over the Asia RCM domain.

3.4. Leaf Area Index

The LAI is defined as the total one-sided area of all green canopy elements over vegetated ground area, which are constructed from the global monthly mean distributions of green vegetation leaf area index, based on the July 1981–December 1999 AVHRR NDVI data at 8-km spacing on the Interrupted Goode Homolosine projection provided by Boston University [21,22]. LAI has missing value pixels due to cloud contamination and atmospheric effects, which are filled by the spatial average over nearby data pixels having the same LCC within a certain radius starting from 16 km (24 pixels) around a missing point and increasing until a 3-pixel minimum is obtained. Filled data are converted into the raster grid, then the polygon coverage, and remapped onto the RCM projection. After this result is further adjusted to be confined by the USGS LCC (see Section 3.2) for a consistent representation of water bodies, the adjusted raw data (Lraw) with respect to unit ground area is divided by local vegetation cover Cv to define the green leaf area index (Lgv) with respect to vegetated area only following Zeng et al. [17]. There is inconsistency between Cv and Lraw data at individual pixels mainly because Cv was derived based on the 24 USGS LCCs at 1-km spacing, but Lraw in terms of six alternative biomes with distinct vegetation structures at an 8-km interval. For the RCM, the 1-km Cv data are integrated onto the 8-km Lraw map to compute Lgv guess values and then a smoothing filter is introduced to remove abnormal values due to inconsistency between Cv and Lraw data at individual pixels. Liang et al. [20] designed the filter through trial and error by examining the frequency distribution of abnormal Lgv values and considering the canopy displacement height in the CLM for each USGS LCC. The point value that exceeds the filter threshold listed in Table 2 is filled by the average over nearby data pixels having the same LCC within a certain radius starting from 16 km (24 pixels) around the point and increasing until a 3-pixel minimum is obtained.

In addition, Lgv data contain large uncertainties in winter due to cloud contamination, especially for the USGS LCC types 13 and 14 (evergreen broadleaf and needleleaf forests). Following Zeng et al. [17], Lgv values in winter months for these two categories are adjusted by:

L g v = max ( L g v , c L g v , max )
where correction coefficient c is 0.8 and 0.7 for LCC types 13 and 14 respectively, and Lgv,max is the maximum Lgv. For the climatology, the maximum can be determined from all monthly values during the entire period, while for interannual variations it is taken in three consecutive years.

After extreme value removal and winter adjustment at each 8-km pixel, the new Lgv is multiplied by its respective Cv and then intersected with the RCM grid mesh. The fractional area of each pixel contributing to the grid is extracted. The area-weighted averaging of all pixels within each RCM grid results in the new LAI per unit ground, which will be divided by local FVC (see Section 3.3) to produce the final LAI. Figure 6 depicts January, April, July, and October mean LAI distributions of the AVHRR climatologies over the RCM domain.

4. Conclusions

The increase in resolution of regional climate and mesoscale atmospheric models can be facilitated by availability of sensor imagery observations in global distributions at the finest spatial and temporal scales. A large volume of global sensor imagery data of the Earth’s terrestrial surface can improve model predictability by supporting better model parameterizations for dynamics of the land surface-atmosphere interactions and increasingly sophisticated assimilation schemes in land surface models. This study focuses on the construction of fundamental SBCs desired for general RCM applications to all effective, dynamically coupled or uncoupled, combinations of the surface modules, as well as portability to any specific region of the world. The new SBCs development motivated by the limitations and inconsistencies of existing SBCs can be readily incorporated into any RCM suitable for climate and hydrology modeling studies. The primary SBCs constructed by using sensor imagery data include surface orography parameters such as means and deviations of terrain elevation at both subgrid-scale and resolved-scale, along with anisotropy and direction angle, and vegetation parameters such as land cover category, fractional vegetation cover, and leaf area index. A critical requirement in constructing the SBCs for RCM use is that each variable must be defined globally with physical consistency across all relevant parameters. This study tried to appropriately manage and rectify existing databases that have various resolutions, a wide range of map projections, different data formats, and often contain missing values or inconsistencies between variables. The GIS application tools such as ArcInfo and ArcMap are mainly used to process vast amount of raw data and to determine the geographic conversion information from a specific map projection of raw data to the identical RCM grid system. Given that data quality and value representation depend on the RCM computational domain and grid resolution, all the SBCs are constructed onto the 30-km RCM domain suitable for the continent of Asia. The raw data sources and processing procedures for the SBCs can be readily used for any specific RCM domain in the world. Even for choosing the best-available data quality, comprehensive processing procedures, and consistency between alternatives, the SBCs so constructed carry over uncertainties inherent in the raw data. Future studies will be required to assess impacts of the SBCs treatments, and an upcoming paper will address the RCM climate sensitivity to these SBCs.

Acknowledgments

This research was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2010-0003954).

References

  1. Intergovernmental Panel on Climate Change. In Fourth Assessment Report, Climate Change 2007—Summary for Policymakers; IPCC: Geneva, Switzerland, 2007.
  2. Chen, F; Dudhia, J. Coupling an advanced land surface-hydrology model with the Penn State—NCAR MM5 modeling system. Part I: Model implementation and sensitivity. Mon. Wea. Rev 2001, 129, 569–585. [Google Scholar]
  3. Klemp, J; Skamarock, W; Dudhia, J. Conservative Split-Explicit Time Integration Methods for the Compressible Nonhydrostatic Equations. WRF Eulerian Prototype Model Equations on Height and Mass Coordinates. 2000. Available online: http://www.mmm.ucar.edu/wrf/users/wrf-dyn-num.html (accessed on 22 June 2011). [Google Scholar]
  4. Michalakes, J. Weather Research and Forecast Model 1.0: Software Design and Implementation. 2000. Available online: http://www.mmm.ucar.edu/wrf/users/wrf-architect.html (accessed on 22 June 2011). [Google Scholar]
  5. Liang, X-Z; Li, L; Dai, A; Kunkel, KE. Regional climate model simulation of summer precipitation diurnal cycle over the United States. Geophys Res Lett 2004, 31, L24208. [Google Scholar] [CrossRef]
  6. Avissar, R; Pielke, RA. A parameterization of heterogeneous land surfaces for atmospheric numerical models and its impact on regional meteorology. Mon. Wea. Rev 1989, 117, 2113–2136. [Google Scholar]
  7. Dickinson, RE; Henderson-Sellers, A; Kennedy, PJ; Wilson, MF. Biosphere Atmosphere Transfer Scheme (BATS) Version 1e as Coupled for Community Climate Model; NCAR Technical Note; Note NCAR/TN–378+STR; Climate and Global Dynamics Division, NCAR: Boulder, CO, USA, 1993; p. 72. [Google Scholar]
  8. Giorgi, F. An approach for the representation of surface heterogeneity in land surface models. Part I: Theoretical framework. Mon. Wea. Rev 1997, 125, 1885–1899. [Google Scholar]
  9. Koster, RD; Suarez, MJ. A comparative analysis of two land surface heterogeneity representations. J. Climate 1992, 5, 1379–1390. [Google Scholar]
  10. Leung, LR; Ghan, SJ. Parameterizing subgrid orographic precipitation and surface cover in climate models. Mon. Wea. Rev 1998, 126, 3271–3291. [Google Scholar]
  11. Zhan, XW; Gao, W; Qi, JG; Houser, PR; Slusser, JR; Pan, XL; Gao, ZQ; Ma, YJ. Remote sensing and modeling the dynamics of soil moisture and vegetative cover of arid and semiarid areas. Proc. SPIE 2004, 5153, 51–60. [Google Scholar]
  12. Dai, Y; Zeng, X; Dickinson, RE; Baker, I; Bonan, GB; Bosilovich, MG; Denning, AS; Dirmeyer, PA; Houser, PR; Niu, G; Oleson, KW; Schlosser, CA; Yang, Z-L. The common land model. Bull. Amer. Meteor. Soc 2003, 84, 1013–1023. [Google Scholar]
  13. Rontu, L. A study on parameterization of orography-related momentum fluxes in a synoptic-scale NWP model. Tellus A 2006, 58, 69–81. [Google Scholar]
  14. Lott, F; Miller, M. A new subgrid-scale orographic drag parameterization: Its formulation and testing. Quart. J. Roy. Met. Soc 1997, 123, 101–127. [Google Scholar]
  15. Rontu, L; Sattler, K; Sigg, R. Parameterization of Subgrid-Scale Orography Effects in HIRLAM, Technical Report 56; HIRLAM, 2002. Available online: http://hirlam.knmi.nl (accessed on 22 June 2011).
  16. Zeng, X; Dickinson, RE; Walker, A; Shaikh, M; DeFries, RS; Qi, J. Derivation and evaluation of global 1–km fractional vegetation cover data for land modeling. J. Appl. Meteor 2000, 39, 826–839. [Google Scholar]
  17. Zeng, X; Shaikh, M; Dai, Y; Dickinson, RE; Myneni, R. Coupling of the common land model to the NCAR community climate model. J. Climate 2002, 15, 1832–1854. [Google Scholar]
  18. Belward, AS. The IGBP–DIS Global 1 km Land Cover Data Set (DISCover): Proposal and Implementation Plans. IGBP–DIS Working Paper No. 13; IGBP–DIS Office: Toulouse, France, 1996; p. 61. [Google Scholar]
  19. Loveland, TR; Reed, BC; Brown, JF; Ohlen, DO; Zhu, Z; Yang, L; Merchant, JW. Development of a global land cover characteristics database and IGBP DISCover from 1 km AVHRR data. Int. J. Remote Sens 2000, 21, 1303–1330. [Google Scholar]
  20. Liang, X-Z; Choi, HL; Kunkel, KE; Dai, Y; Joseph, E; Wang, JXL; Kumar, P. Surface boundary conditions for mesoscale regional climate models. Earth Interactions 2005, 9, 1–28. [Google Scholar]
  21. Buermann, W; Wang, Y; Dong, J; Zhou, L; Zeng, X; Dickinson, RE; Potter, CS; Myneni, RB. Analysis of a multiyear global vegetation leaf area index data set. J Geophys Res 2002, 107, 4646. [Google Scholar] [CrossRef]
  22. Zhou, L; Tucker, CJ; Kaufmann, RK; Slayback, D; Shabanov, NV; Myneni, RB. Variations in northern vegetation activity inferred from satellite data of vegetation index during 1981 to 1999. J. Geophys. Res 2001, 106, 20069–20083. [Google Scholar]
Figure 1. Asia RCM domain overlaid with (a) latitude and longitude and (b) the 301 × 251 dimensional 30-km spacing grids.
Figure 1. Asia RCM domain overlaid with (a) latitude and longitude and (b) the 301 × 251 dimensional 30-km spacing grids.
Sensors 11 06728f1 1024
Figure 2. The comparison of geographic distributions of surface orography parameters over the Asia RCM domain. (a) and (b) denote mean terrain elevations at subgrid-scale and resolved-scale, respectively; (c) and (d) denote standard deviations of terrain elevations at subgrid-scale and resolved-scale, respectively.
Figure 2. The comparison of geographic distributions of surface orography parameters over the Asia RCM domain. (a) and (b) denote mean terrain elevations at subgrid-scale and resolved-scale, respectively; (c) and (d) denote standard deviations of terrain elevations at subgrid-scale and resolved-scale, respectively.
Sensors 11 06728f2 1024
Figure 3. The geographic distributions of (a) anisotropy and (b) direction angle at resolved-scale over the Asia RCM domain.
Figure 3. The geographic distributions of (a) anisotropy and (b) direction angle at resolved-scale over the Asia RCM domain.
Sensors 11 06728f3 1024
Figure 4. The geographic distribution of 21 LCC categories over the Asia RCM domain.
Figure 4. The geographic distribution of 21 LCC categories over the Asia RCM domain.
Sensors 11 06728f4 1024
Figure 5. The geographic distributions of FVC derived from the April 1992–March 1993 AVHRR over the Asia RCM domain.
Figure 5. The geographic distributions of FVC derived from the April 1992–March 1993 AVHRR over the Asia RCM domain.
Sensors 11 06728f5 1024
Figure 6. The geographic distributions of mean LAI based on the original 1981–1999 AVHRR climatology data in January, April, July, and October over the Asia RCM domain.
Figure 6. The geographic distributions of mean LAI based on the original 1981–1999 AVHRR climatology data in January, April, July, and October over the Asia RCM domain.
Sensors 11 06728f6 1024
Table 1. Comparison of Nc,v between USGS and IGBP Land Cover Legends.
Table 1. Comparison of Nc,v between USGS and IGBP Land Cover Legends.
USGS land use/land cover legendIGBP land cover legend
TypeDescriptionNc,vTypeDescriptionNc,vOccupation ratio (%)
1Urban and Built-Up Land0.6213Urban and Built-Up0.62100
2Dryland Cropland and Pasture0.6112Croplands0.61100
3Irrigated Cropland and Pasture0.6112Croplands0.6194.41
14Cropland/Natural Vegetation Mosaic0.655.59
5Cropland/Grassland Mosaic0.6512Croplands0.611.33
14Cropland/Natural Vegetation Mosaic0.6598.67
6Cropland/Woodland Mosaic0.6514Cropland/Natural Vegetation Mosaic0.65100
7Grassland0.4910Grasslands0.49100
8Shrubland0.606Closed Shrublands0.6014.81
7Open Shrublands0.6081.00
8Woody Savannas0.624.19
9Mixed Shrubland/Grassland0.596Closed Shrublands0.6015.28
7Open Shrublands0.6076.53
10Grasslands0.498.19
10Savanna0.608Woody Savannas0.6250.48
9Savanna0.5849.52
11Deciduous Broadleaf Forest0.702Evergreen Broadleaf Forest0.6920.41
4Deciduous Broadleaf Forest0.7068.26
5Mixed Forest0.6811.33
12Deciduous Needleleaf Forest0.634Deciduous Needleleaf Forest0.63100
13Evergreen Broadleaf Forest0.692Evergreen Broadleaf Forest0.69100
14Evergreen Needleleaf Forest0.631Evergreen Needleleaf Forest0.63100
15Mixed Forest0.685Mixed Forest0.68100
16Water Bodies17Water Bodies100
17Herbaceous Wetland0.5611Permanent Wetlands0.56100
18Wooded Wetland0.5611Permanent Wetlands0.56100
19Barren or Sparsely Vegetated0.6016Barren or Sparsely Vegetated0.60100
21Wooded Tundra0.617Open Shrublands0.6068.11
8Woody Savannas0.6231.89
22Mixed Tundra0.6016Barren or Sparsely Vegetated0.60100
23Bare Ground Tundra0.6016Barren or Sparsely Vegetated0.60100
24Snow or Ice15Snow and Ice100

Note: Land cover types 4 (Mixed Dryland/Irrigated Cropland & Pasture) and 23 (Herbaceous Tundra) do not exist in the global dataset.

Table 2. Parameters in deriving LAI for each USGS Land Cover.
Table 2. Parameters in deriving LAI for each USGS Land Cover.
TypeUSGS land cover legendFVC (1 km)Displacement height (m)Lgv filter threshold
1Urban and Built-Up Land0.7350.6677
2Dryland Cropland and Pasture0.8750.6677
3Irrigated Cropland and Pasture0.8040.6677
5Cropland/Grassland Mosaic0.7290.6677
6Cropland/Woodland Mosaic0.8690.6677
7Grassland0.7110.6676
8Shrubland0.3810.3335
9Mixed Shrubland/Grassland0.3910.3335
10Savanna0.8480.6677
11Deciduous Broadleaf Forest0.87113.3338
12Deciduous Needleleaf Forest0.92013.3338
13Evergreen Broadleaf Forest0.95323.3338
14Evergreen Needleleaf Forest0.89513.3338
15Mixed Forest0.87513.3338
16Water Bodies0.667
17Herbaceous Wetland0.94713.3336
18Wooded Wetland0.8350.6678
19Barren or Sparsely Vegetated0.0730.3334
21Wooded Tundra0.7140.6676
22Mixed Tundra0.3230.3336
23Bare Ground Tundra0.0180.3336
24Snow or Ice0.667

Note: Land cover types 4 (Mixed Dryland/Irrigated Cropland & Pasture) and 23 (Herbaceous Tundra) do not exist in the global dataset.

Share and Cite

MDPI and ACS Style

Choi, H.I. Use of Sensor Imagery Data for Surface Boundary Conditions in Regional Climate Modeling. Sensors 2011, 11, 6728-6742. https://0-doi-org.brum.beds.ac.uk/10.3390/s110706728

AMA Style

Choi HI. Use of Sensor Imagery Data for Surface Boundary Conditions in Regional Climate Modeling. Sensors. 2011; 11(7):6728-6742. https://0-doi-org.brum.beds.ac.uk/10.3390/s110706728

Chicago/Turabian Style

Choi, Hyun Il. 2011. "Use of Sensor Imagery Data for Surface Boundary Conditions in Regional Climate Modeling" Sensors 11, no. 7: 6728-6742. https://0-doi-org.brum.beds.ac.uk/10.3390/s110706728

Article Metrics

Back to TopTop