An Integrated Method Combining Remote Sensing Data and Local Knowledge for the Large-Scale Estimation of Seismic Loss Risks to Buildings in the Context of Rapid Socioeconomic Growth: A Case Study in Tangshan, China

Rapid socioeconomic development in earthquake-prone areas can cause rapid changes in seismic loss risks. These changes make it difficult to ensure that risk reduction strategies are realistic, practical and effective over time. To overcome this difficulty, ongoing changes in risk should be captured timely, definitively, and accurately and then specific and well-timed adjustments of the relevant strategies should be made. However, methods for rapidly characterizing such seismic disaster risks over a large area have not been sufficiently developed. By focusing on building loss risks, this paper presents the development of an integrated method that combines remote sensing data and local OPEN ACCESS Remote Sens. 2015, 7 2544 knowledge to resolve this problem. This method includes two key interdependent steps. (1) To extract the heights and footprint areas of a large number of buildings accurately and quickly from single high-resolution optical remote sensing images; (2) To estimate the floor areas, identify structural types, develop damage probability matrixes, and determine economic parameters for calculating monetary losses due to seismic damage to the buildings by reviewing building-relevant local knowledge based on these two parameters (i.e., the building heights and footprint areas). This method is demonstrated in the Tangshan area of China. Based on the integrated method, the total floor area of the residential and public office buildings in central Tangshan in 2009 was 3.99% lower than the corresponding area number obtained by a conventional earthquake loss estimation project. Our field-based verification indicated that the mean relative error of the method for estimating the floor areas of the assessed buildings was 2.99%. A simulation of the impacts of the 1976 Ms 7.8 Tangshan earthquake using this method indicated that the total damaged floor area of the residential and public office buildings and the associated direct monetary loses in the study area could have been 8.00 and 28.73 times greater, respectively, than in 1976 if this earthquake had recurred in 2009, which is a strong warning to the local people regarding the increasing challenges they may face.


Introduction
Sound knowledge of seismic disaster risk is a prerequisite for ensuring the effectiveness and efficiency of seismic hazard and disaster mitigation strategies, preparedness measures, and emergency response plans. Seismic disaster risks are constantly and significantly changing in earthquake-prone areas where rapid socioeconomic development is occurring (including many parts of China). One such area is Tangshan, China, which is located in the "Chinese Capital Region". This area experienced severe impacts from the Ms 7.8 earthquake that occurred in 1976 and significant population and economic expansion following 1976. From 1985 to 2010, the population of this administrative region increased from 1.38 to 3.07 million, with the GDP increasing from 0.61 to 67.48 billion USD, urbanization increasing from 24.6% to 53.4%, and urban built-up areas expanding from 100 to 209 km 2 [1]. These increases inevitably caused significant changes in the seismic vulnerabilities and risks from year to year [2][3][4]. The continual, rapid, and accelerating changes and/or growths in vulnerability and risk have made it difficult to ensure that the area's earthquake disaster risk reduction strategies have remained fact-based, practical and effective over time. Certain strategies may become irrelevant after a short period. To prevent such strategies from becoming irrelevant, associated changes in risks must be captured quickly and accurately to perform timely, specific, and pertinent adjustments.
In the fields of hazard and disaster research, disaster risk is commonly defined as "the potential disaster losses in lives, health status, livelihoods, assets and services, which could occur to a particular community or a society over some specified future time period" [5]. The classical expression of disaster risk is widely formulated as follows [6]: where is the disaster risk, is the hazard, and represent the elements and vulnerabilities of the elements, respectively, exposed to the hazard. The oldest and most developed approach for evaluating large scale seismic disaster risks is the conventional and mainstream "earthquake loss estimation" approach. This approach has served as the focus or base of various studies, such as ATC-13 [7], ATC-25 [8], EMS98 [9][10][11], RADIUS [12], and HAZUS [13][14][15]. In these studies, the value of in Equation (1) has typically been expressed as the intensity of the seismic ground motion (e.g., the modified Mercalli intensity (MMI)). The element of in Equation (1) primarily included buildings, lifelines, other man-made facilities, and the residents themselves, and was typically expressed in terms of the loss ratio curves (LRCs) or the damaged probability matrixes (DPMs) of the exposed elements ( ) induced by seismic ground motion. Assessing the values with certain reliable specific occurrence probabilities for large areas is extremely difficult; thus, researchers have generally applied deterministic methods (e.g., the scenario earthquake method) to realistically and practically specify, establish and determine the value of using in Equation (1). These specific previous studies were based on establishing various LRCs or DPMs. One primary advantage of these LRC or DPM-based methods is that they yield definite, accurate results that are usually expressed in terms of the specific numbers of buildings damaged, human casualties, or monetary losses. However, these methods require the preparation of large amounts of detailed data for in Equation (1) (e.g., building inventories). In most cases, acquiring these data requires significant labor, time and money. When assessing seismic loss risks over a considerably large area with rapid socioeconomic development, applying this type of method may be difficult or impossible. It is highly likely that the large amount of data collection required to target a large area of this type through conventional data preparation cannot keep pace with rapid socioeconomic development and evolving risks. This difficulty is likely greater on the Chinese mainland than elsewhere due to China's current breakneck speed of socioeconomic development.
There are other approaches that are capable of characterizing large-scale seismic disaster risk potentials, including the approaches given by Cutter et al. [16,17], Su et al. [2,3], Carreño et al. [18], and Cardona and Carreño [19]. Unlike the LRC and DPM-based methods, the specific methods of these approaches do not usually address all three elements ( , , and in Equation (1)) that contribute to risk. Specifically, these previous approaches only addressed the and or elements. However, these and or -focused studies can still reveal the seismic risk potentials of a place (to a certain extent) because the exposure of more elements at a certain place and/or the greater vulnerabilities of these elements generally correspond with more risk under conditions with the same seismic intensity impacts. Due to the extreme difficulties that can occur when obtaining seismic intensities with reliable occurrence probabilities for a large area, this general perception and/or understanding of seismic risks or risk potentials is justified and understandable. Thus, this risk perception manner has often been adopted in the field of earthquake disaster research. The first step in these approaches is to establish a set of and indexes or indexes with primary components that are derived from various public data sources (e.g., statistical yearbooks). Next, methods of this type are used to synthesize one or several man-made and or composite indexes with values that vary within a certain index range (e.g., from 0 to 1) to broadly and comprehensively characterize the disaster risk potentials of a given location. Subsequently, with updated data that can be more easily realized than updating the data required by the LRC or DPM-based methods, the temporal variations of the risk potentials at such locations during various periods can be conveniently observed. Because readily available public data are mainly used as input data, these and or composite index-focused methods possess several advantages over LRC-or DPM-based methods. Specifically, index-focused methods are more comprehensive, save more labor and money, easily cover large areas, and are easier to use when examining potential risk variations (even periodically, such as once every few years). However, conventional and mainstream LRC and DPM-based methods can yield explicit, specific results in terms of specific amounts of building damage, human deaths, and economic losses, whereas the and or composite index-focused methods can only provide general and relative results regarding statistical and inter-place contrast when considering if a location has risks and if the level of a risk at one given location may be greater than another. It is not always easy to determine what elements on earth at a given location are at risk and how serious the hidden problems on earth are when considering results represented as a highly composite index. Moreover, if the basic geographic or spatial unit of analysis for a large scale study is too small (such as the sub-county levels, e.g., the village level or town/township level), the and or composite index-focused methods frequently involve nearly the same levels of difficulty regarding data collection (e.g., the percentage of aged people, per capita income, and mean education level of each village or town/township) as the LRC and DPM-based methods. For example, a distinct lack of public data sources or similar databases exists regarding data at the sub-county scales at the contemporary Chinese mainland. Thus, it is difficult to apply this type of method on a large scale to assess seismic risk potentials at high spatial resolutions.
In addition to the two types of approaches mentioned above that are capable of characterizing seismic disaster risks or risk potentials over large areas, some other studies are also helpful for understanding or estimating seismic disaster risks of an area. Several researchers have developed GDP and population data-based models to broadly estimate the seismic vulnerability and risk of an area [20][21][22][23]. However, the risk implications of the results of these studies are not specific or detailed relative to those of conventional and mainstream LRC and DPM-based methods. From the standpoint of timeliness requirements that are crucial for capturing and reducing seismic disaster risks in the context of rapid socioeconomic development, the seismic disaster risk studies that are relevant to the field of remote sensing are particularly valuable. Because remote sensing images, especially those with high resolutions, may be the only currently practical and reliable tools that are capable of providing full, clear and detailed coverage of large areas when labor resources and work times are limited, this information is urgently needed. Recently, apart from applications in post-earthquake impact assessments [24][25][26][27][28][29][30][31], many studies have addressed the application of high-resolution remote sensing technology to pre-earthquake exposure, vulnerability, and risk. In these pre-earthquake studies, the and/or -relevant issues (see Equation (1)) of various buildings were widely discussed, including various optical imagery-based specific analyses (Table 1) [32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49][50], some attempts at fusing optical imagery and LiDAR or SAR data [51][52][53] to help extract the heights of buildings, and several comprehensive reviews [54][55][56][57]. Table 1. Brief summary of studies that applied high resolution optical RS to determine the seismic loss risks of buildings (excluding reviews).

Author Data
The Main Topics Addressed

Accuracy or Reliability
Dutta et al. [32] QuickBird (1) Developing an automatic method for estimating the number of buildings.
(2) Combining RS-derived roof material information, field data and existing database data to compile building inventories. The buildings were classified as low-, mid-, or high-rise reinforced or non-reinforced concrete buildings.
(3) Using the fragility curves of HAZUS and the scenario earthquake method to estimate the building damage numbers.
52 km 2 (suburban area of Bangkok, Thailand) Not mentioned.
Ehrlich et al. [33] QuickBird Developing automatic methods for extracting building footprint-related urban built-up areas and building heights.
48 km 2 (urban area of Sana'a, Yemen) The MRE of the extraction of relevant built-up area was 78.86%.
Geiβ et al. [34] IKONOS; Multi-temporal TM and ETM; nDSM Exploring the method and extent by which the seismic vulnerabilities of the buildings may be examined and assessed directly using various RS-derived building features and their environments. The in situ assessed vulnerabilities of the buildings were based on an expert scoring and EMS-98 scheme.

buildings for method development
Compared with expert scoring, the MAPE of the method was 10.6%. Compared with EMS-98, the overall accuracy of the method was 65.4%.
Geiβ et al. [35] IKONOS; Multi-temporal TM and ETM; nDSM (1) Examining the correlations of the seismic building structural types (SBSTs) with the buildings and environmental features derived from remote sensing data. The in situ determined SBSTs were based on the relevant logic of HAZUS.
(2) Local LRCs and the scenario intensity-based earthquake loss estimates.

IKONOS
(1) Developing the RS-based method for updating existing GIS building inventories. The focus of this method was used to automatically detect new buildings using high-resolution optical images.
(2) Using existing logic for classifying the seismic vulnerabilities of local buildings to characterize the seismic vulnerability features of the changed building stock.
(3) Using the fragility curves of local buildings and the scenario earthquake method to estimate building damage numbers.
Approximately 300 km 2 (Metro Manila, Philippines) Nearly 90% of the new mid and high-rise buildings can be detected successfully.
Mück et al. [39] QuickBird; DEM Exploring the method and extent by which the seismic vulnerabilities of buildings may be evaluated directly using the building features (size, shape, orientation, regularity, height, and accessibility) derived from RS data. The in-situ building vulnerability measurements were based on the expert scoring method.

buildings for method development
Compared with expert scoring-based vulnerability measurements, the overall accuracy of the method was 81%.
Panagiota et al. [40]; Panagiota et al. [41] Airborne high resolution data; DEM Exploring a method by which the seismic vulnerabilities of buildings may be directly evaluated using the building roof types and heights derived from RS data. The in-situ vulnerability assessments of buildings were based on EMS-98.

buildings for method development
Compared with EMS98-based vulnerability assessment, the overall accuracy of the best model of the method was 62.5%.
Pittore and Wieland [42]; Wieland et al. [43] Multi-temporal TM and MSS; High resolution optical images; OpenStreetMap; Omni-directional images (1) Developing the method using high-resolution optical RS images and OpenStreetMap to collect footprints. Developing the method using Omni-directional images to measure the building heights and floor numbers.
(2) Exploring the method by which the seismic vulnerabilities of buildings may be directly evaluated by integrating step-by-step building and block's features derived from multi-source images. These features include land-use/coverage and the ages of built-up areas of blocks, construction types, footprints, areas, shapes, heights, floor numbers, ages, and building uses. The seismic vulnerabilities of buildings were characterized using EMS-98.
Bishkek, Kyrgyzstan (case study area) Not mentioned  [44] High-resolution optical imagery Exploring this method by combining building features that were derived from RS (footprint, height, shape, roof material, and roof type) and ancillary information to estimate types of seismic building structures (SBSTs). The SBSTs were determined based on the relevant logic of HAZUS.

buildings for method development
The error of the best SBSTs estimation model that based on RS data and ancillary information was 13.46%.
Wu et al. [45] WorldView-2 Relating the heights, structural types, and uses of buildings and the local geologic conditions with the seismic damage grades of buildings under one specific seismic intensity (intensity VII). Building heights were derived from RS images using an existing method; the structural type and uses of the buildings and the local geologic conditions were obtained using in-situ surveys. (2) Referring to Yin's DPMs [49] to address the vulnerabilities of the buildings.

buildings
The MRE of the footprint assessment was less than 10%.
Note: (1) The topic(s) presented in bold text in the "The main topics addressed" column indicate the focuses of each relevant study that had multiple main topics.
(2) MAE represents the mean absolute error; MRE represents the mean relative error; RMSE is root mean square error; and MAPE is the mean absolute percentage error.
These remote sensing-based and/or studies showed that remote sensing technologies have great potential for overcoming timeliness and coverage problems relative to LRC or DPM-based approaches and and or composite index-focused methods for addressing seismic loss risks in the context of rapid socioeconomic growth. However, when accurately assessing the seismic loss risks of a large area determinately and in detail, remote sensing technologies still face several challenges that may be difficult or impossible to overcome when exclusively using this technology. (1) Currently, it is still difficult to use remote sensing data and technology alone to directly produce some types of information that are essential for quantitatively calculating specific seismic loss risks, especially core data or information, such as building structural types and DPMs. Several previous studies have broadly addressed the associations between the seismic vulnerabilities of buildings and the building heights, shapes, roof/surface materials, and other properties that can be observed, measured, and/or detected from remote sensing images [33,36,46]. However, these association discussions are far from the implications of the LRCs and DPMs of buildings, which are classical expressions of seismic building vulnerabilities and were summarized based on earthquake engineering knowledge and actual earthquake damages. Some recent studies have carefully explored how remote sensing may contribute to examining and assessing the seismic structure type and/or vulnerability of a building [34,35,[39][40][41][42][43][44][45], which are the technical cores of the mainstream quantitative approaches used to determine seismic loss risks of buildings. However, uncertainties remain that require further validation. For example, Wu et al. [45] used several tens of sample buildings in an interesting attempt to relate heights, structural types, and uses of buildings and the local geologic environments with the building damage grades under one specific seismic intensity (intensity VII). However, it is worth further exploring what or how much do these relationships imply relative to classical connotation and conventional expressions of seismic vulnerability of a building? In addition, how reliable are such relationships when they were derived from 35 buildings (for example, relative to the commonly accepted DPMs on the Chinese mainland that were summarized based on experiences from multiple actual earthquakes and innumerable affected buildings)? Geiß et al. [34], Mück et al. [39], Panagiota et al. [40,41], Wieland et al. [43], and Pittore and Wieland [42] attempted to directly link the seismic vulnerabilities of buildings with the features of buildings and the buildings' environments derived from remote sensing data. These attempts attempted to bypass the critical but frequently hard tasks (esp. when facing a large number of buildings) of conventional earthquake engineering methods for characterizing the seismic structural vulnerabilities of buildings. Thus, the methods of this kind were quick, capable of easily covering a large number of buildings, and universally transferable. However, the applicabilities of these new and enlightening alternatives for characterizing the seismic vulnerabilities of buildings (relative to conventional earthquake engineering-based methods), may require extensive and careful discussion and validation by experienced earthquake engineers and associated remote sensing-focused scientists. In addition, the diverse actualities of different places and the actual earthquake damages that were not sufficiently discussed in the context of these previous studies may also require careful considerations. Geiß et al. [34] stated that an open dialogue should be triggered. Given the needs and difficulties of large-scale estimations of seismic loss risks in the context of rapid socioeconomic growth, such dialogue is particularly worthwhile, urgent and expected; (2) Remote sensing imagery with high resolution contains important and specific information (e.g., footprint areas of buildings) that can be used to characterize specific seismic loss risks in detail. However, the techniques for extracting this information have not been sufficiently developed. When addressing a large number of objects (e.g., a large number of buildings) over a limited amount of time, it is difficult to extract highly accurate and specific information (e.g., the total floor areas of a building) from remote sensing imagery. Although Li and Zhai [36] and Zhao et al. [50] developed relatively accurate information, their studies only involved a very small number of buildings. This may be one important reason why existing remote sensing studies relevant to pre-earthquake disaster risks have mainly focused on developing various technological methods and have rarely presented applications of them in relatively or considerably large places with relatively or considerably large numbers of buildings. In summary, several large gaps exist between what remote sensing images and technology can currently provide and what conventional and mainstream large-scale seismic loss risk estimates specifically require. Bridging such gaps is important for promoting deep and actual practice-aimed applications of remote sensing data and technology in broad seismic loss risk reduction battlefields.
A method that rapidly and precisely captures seismic disaster risks and their changing scenarios over large areas with rapid socioeconomic growth has not been sufficiently developed to date. The purpose of this study was to propose an integrated or hybridized method that combines high-resolution remote sensing imagery with local knowledge to solve this issue. To accomplish this task, we introduced local knowledge reviews to bridge the gaps between the remote sensing and technical cores of conventional quantitative approach for determining seismic loss risks. We focused on building loss risks to demonstrate the development and use of this integrated method. Extensive and large-scale building development has been one of the most visible and significant aspects of China's rapid socioeconomic expansion within the last several decades [58][59][60][61][62]. No surveys of Chinese mainland buildings have been reported in the literature since 1985 when a nationwide building survey was performed. Such building-focused studies are urgently needed to cope with seismic disaster risks on the Chinese mainland.
In this building-focused study conducted in 2010, the term local knowledge was not used in the purely social science sense. Instead, the term building-relevant local knowledge (Br-LK) was used to refer to various sources of information regarding the buildings in a specific area. This knowledge included local building codes, local dwelling traditions and culture, and local building construction planning and policies. The Br-LK also included local building information contained in almanacs, statistical yearbooks, the chorography of local chronicles, building construction archives, urban and land-use planning reports, volunteered geographic information [63][64][65][66], academic books and papers, and scientific and social service research reports. The high-resolution remote sensing imagery used in this study was high-resolution optical remote sensing images (Hr-ORSIs) that were either satellite or airborne-based with visible primary image characteristics from ground objects that could be manually examined. By fully combining the advantages of these two sources (Br-LK and Hr-ORSIs), the definitiveness, accuracy and timeliness were all enhanced when estimating the earthquake-induced building loss risks across a large area with rapid socioeconomic growth.
The paper begins with the overall framework of the proposed integrated method. Next, it demonstrates how the heights and the footprint areas of a large number of buildings may be accurately and rapidly extracted from a single Hr-ORSI. Then, the bridging process is described by reviewing Br-LKs to help accurately and rapidly determine the various core parameters required to specifically calculate the building loss risks, including the building structural types, total and sub-total floor areas, DPMs, and economic parameters for calculating monetary losses due to seismic damage to buildings. Next, two methods are presented for verifying the accuracy of the method. In a specific application of this method, a simulation of the impacts of the Ms 7.8 Tangshan earthquake of 1976 is presented in the context of urban development as of 2009. In the penultimate section of the paper, a series of discussions are presented regarding the advantages, limitations, application prospects, and future extendibility or scalable directions of the method. Finally, the paper ends with several final remarks of this study.

Overall Framework of the Proposed Method
Within the last several decades, earthquake engineers in China have performed extensive conventional LRC-and/or DPM-based earthquake loss estimation work [67][68][69][70][71], often based on the methodologies developed by the Applied Technology Council of the USA (ATC). In this study, we followed the basic logic of the ATC and its Chinese variations to explore a novel integrated method for the timely, accurate, and large-scale characterization of seismic loss risks in areas of rapid socioeconomic growth. Specifically, (1) we directly established the seismic hazard parameter (H) in Equation (1) in terms of certain specific hypothetical or scenario intensity values based on the Chinese Seismic Intensity Scale [72]. As in the Modified Mercalli Intensity (MMI) scale used in the USA and the European Macroseismic Scale (EMS), this seismic intensity scale included twelve degrees designated by Roman Numerals I through XII. The ground motion and associated macro damage implications are generally equal to one another among the three scales at the same degree [73][74][75]. The empirical and statistical conversion/transform relationships between the intensities of the Chinese Seismic Intensity Scale and Peak Ground Accelerations (PGA) are shown in Table 2 [76]. In actual seismic loss assessment practices either before or after earthquakes on the Chinese mainland, currently available DPMs are still based on seismic intensity, such as the baseline DPMs for different regions developed by Yin [49]. No PGA-based DPMs are currently available that can be used to analyze seismic losses of buildings in groups on the Chinese mainland. Thus, we also used intensity as a seismic hazard parameter. Next, (2) we selected buildings as representative exposed elements for this study. After these selections, the integrated method was used to rapidly and accurately determine the required parameters for estimating the seismic loss risks of the buildings, including their structural types, DPMs, total and sub-total floor areas, and economic parameters for calculating the monetary losses due to seismic damage. (1) a max indicates the maximal acceleration (i.e., the peak ground acceleration (PGA)); (2) g refers to the gravitational acceleration, which is commonly defined as 1 g = 9.81 m/s 2 .
The overall framework of this integrated method is illustrated in Figure 1. Specifically, the method consists of the following main interdependent parts: (1) estimating the height and the footprint area of each individual building from a single Hr-ORSI accurately and quickly by improving and/or using existing application-oriented techniques; (2) using Br-LK reviews as a bridging tool to determine the total and sub-total floor areas, structural types, DPMs, and key economic parameters for calculating monetary losses associated with each type of building that are required by conventional LRC-and/or DPM-based earthquake loss estimates as quickly, reliably, and accurately as possible based on the building's height and footprint area; (3) assigning various specific hypothetical intensity values to the relevant equations to assess the building loss risks under various seismic hazard scenarios.  The use of this integrated method and its verifications were demonstrated in the central part of Tangshan, China (Figure 2), an area where our associated studies were based for several years [2][3][4]77,78]. As shown in Figure 2, the specific target region includes the entire urban district of Tangshan, its periurban fringes, and a few of its suburban and rural towns and townships. This target region is flat in terrain, the entire study area measures 456 km 2 , and the area had a population of 698,596 in 2009. The population density in this area in 2009 was 1532 persons per square kilometer, which was 12 times greater than the average population density of China in 2009. The region's GDP was 6.61 billion USD in 2009 with a GDP per capita of 9463 USD, which was approximately 3.5 times the national average of China in 2009 [79].

Developing the Parameters Required for Estimating Building Loss Risk Using Single Hr-ORSIs
In Sections 3 and 4 of the paper, we explain how the core parameters that were required for calculating the potential seismic building loss risks were determined by combining the Hr-ORSI and Br-LK data. These parameters included the building structural types, DPMs, total floor areas, and economic metrics that were necessary for calculating monetary losses. These parameters cannot be obtained directly from remote sensing images, regardless of the spatial resolutions of the images, because only the roofs and parts of the outside walls of buildings are visible (even in Hr-ORSIs). However, the heights and the footprint areas of the buildings can be directly estimated from an Hr-ORSI, and these two parameters are important for determining many of the required core parameters. We used GeoEye-1 images in 2009 to explore the ways in which these two beginning and basic parameters can be accurately and rapidly estimated.

Determining Building Heights
According to the literature, there are two primary approaches for estimating building heights using an Hr-ORSI. One method is to use stereo Hr-ORSI pairs [80][81][82][83][84][85], and the other method is to use the shadows of buildings from a single Hr-ORSI [33,37,[86][87][88]. Earlier studies have shown that using stereo Hr-ORSI pairs results in more accurate height estimates. However, it is often difficult to acquire the needed stereo Hr-ORSI pairs. There are two ways to obtain Hr-ORSI pairs. The first is to purchase existing pairs. However, few Hr-ORSI pairs exist that cover the Chinese mainland. For example, the Tangshan area is one of the few areas in China that has been extensively covered by various Hr-ORSIs. However, existing stereo pairs of Hr-ORSIs remain rare. The second method is to trust a relevant company to produce new stereo Hr-ORSI pairs. Existing Hr-ORSI data are expensive; the custom production of new Hr-ORSI pairs is even more expensive. Therefore, the financial cost of using stereo Hr-ORSI pairs may substantially increase when addressing a considerably large area.
For the purpose of broader applications, the proposed integrated method began by using single existing Hr-ORSIs to determine building heights. Researchers have developed two primary types of single Hr-ORSI-based methods for this task. Some researchers have estimated a building's height by measuring the entire length of the building's shadow [37,[87][88][89]. This type of method is theoretically and mathematically the easiest and was developed first. Other methods have been used to estimate building heights by measuring two other building imaging parameters, the building azimuth and shadow length normal to building strike/trend [47,50,86,90,91]. However, it is difficult and often impractical to apply these two types of methods to large areas containing a large number of buildings. The first type of method is not capable of resolving incomplete or overlapping shadow issues due to various interferences from neighboring buildings or other surrounding ground objects. The methods of the second type can overcome these challenges, but require many complex measurements that are labor and time intensive for areas that contain large numbers of buildings (details are presented below).
We developed an algorithm that only requires deriving the same imaging point distance (SIPD) from a single Hr-ORSI to estimate building heights. The SIPD is the distance between the imaging position of a certain point on a building's roof edge and the imaging position of that point's shadow on a single Hr-ORSI. The imaging geometry of a building and that of its SIPD are illustrated in Figure 3. In this figure, the gray bar represents a building, denotes a point on the building's roof edge, is the vertical projection of point on the ground, is the building's height, and and denote the imaging positions of point and its shadow, respectively, in a Hr-ORSI. The distance between imaging points and is the SIPD, which is denoted by in Figure 3. Using Figure 3 and referring to the basic geometry principles of triangles, we first obtained the following rudimentary expressions: where is the building's height; is the SIPD; ω and α are the elevation and azimuth angles of the satellite, respectively; and θ and β are the elevation and azimuth angles of the sun, respectively. Based on these three basic formulas, we obtained the following expression: For simplification, we introduced the coefficient φ to simplify Equation (5) as follows: The coefficient φ, as shown in Equation (7), is determined using the imaging mechanism and parameters of the Hr-ORSI. To obtain (the building's height), the only task left is to measure the SIPD, namely , in Equation (6).
The existing single Hr-ORSI-based methods explained above must be based on either measuring the entire length of a building's shadow or on measuring a building's azimuth and its shadow length normal to its strike/trend. When applied to a large number of buildings over a large area with rapid socioeconomic development, the algorithm we developed is easier, quicker, and more accurate, feasible and practical than the two types of existing methods. (1) The SIPD method does not consider the entire length of a building's shadow (i.e., it does not consider point in Figure 3). Therefore, the SIPD method allows one to avoid extreme difficulties in identifying this point, which can result from various shading effects or from possible high gray-tone similarities between this point and its surroundings; (2) The two imaging points and shown in Figure 3 originate from the same reference point (i.e., the same point of a building's roof edge). To serve as a reference, we always selected a distinct point, such as a building's roof edge corner. This selection made the identification of points and much easier. Moreover, using a point on a building's roof edge as a reference was not the only option. All of these advantages of this algorithm allow one to avoid various shading or gray-tone similarity interferences when identifying points and ; (3) Compared with using a building's azimuth and its shadow length normal to the building's strike/trend to determine a building's height, the SIPD-based method only removes one required parameter in theory. However, this new measurement saves significant time and labor by avoiding the complex process of tracking and determining each building's azimuths, normal, and strike/trends separately. Saving time and labor is critical when targeting a large area with a high number of buildings fueled by rapid socioeconomic growth.
We used ArcGIS 9.3 to calculate the building heights. During this process, two key steps were performed. First, the SIPDs were manually identified and vectorized using ArcGIS 9.3. This step was neither difficult nor particularly time or labor consuming. For example, all of the densely packed buildings (approximately 1200 individuals) in a 50-km 2 zone in the study area were vectorized within 4 to 5 h by a single skilled person. The second step was to calculate all of the SIPDs, in Equation (6) using the software's Calculate Geometry and Field Calculator modules. This calculation was simple and rapid.

Estimating Building Footprint Areas
We based the estimate of a building's footprint area on the outline information of the building's roof provided by the Hr-ORSI. To date, experts have developed several automatic and semiautomatic procedures to rapidly delineate a building's roof outline [33,39,45,[92][93][94][95][96][97][98]. However, it is still difficult to apply these technically advanced procedures to irregularly shaped buildings and to identify and eliminate shading, interferences, and high gray-tone similarity effects resulting from either the buildings themselves or from the buildings' environments. To ensure accuracy, a great deal of post-processing work must be performed. When addressing a large area with a large number of buildings, applying these automatic or semiautomatic procedures remains difficult.
To achieve the good accuracy that is required (and essential in fact) for seismic loss risk reduction practices, this study included manual vectorization to estimate the buildings' footprint areas. Using GIS, the roof outlines of buildings were manually digitized and the footprint areas were estimated accordingly. Indeed, the manual vectorization technique for determining the outlines of building roofs using ArcGIS 9.3 is easy and not very time or labor intensive. Specifically, the required time and labor essentially equaled that required to manually digitize the buildings' SIPDs.

Developing the Required Parameters for Building Loss Risk Estimations by Reviewing Various Sources of Br-LK
In China, the most established and commonly used system for classifying the seismic vulnerabilities of buildings is the system first developed by Yin [68]. In this system, buildings are assigned to the following five categories: (1) earth/wood buildings (E/W buildings); (2) single-story brick buildings (S/B buildings); (3) brick and concrete buildings (B/C buildings); (4) reinforced concrete buildings (R/C buildings); and (5) specially structured buildings. However, not all of these building categories were addressed in this study. We concentrated on the first four categories of Yin's classification system by focusing on residential buildings and public service offices (For clarity, we use the following phrases/expressions: addressed buildings, assessed buildings, buildings for residential and public office usages, or buildings of interest as needed in the associated statements/sentences of the article to indicate this concentration/focus.). This focus was selected due to China's high population density, with residential buildings representing the vast majority of its total building stock. Furthermore, during the recent decades of rapid socioeconomic development in Mainland China, the growth in residential and public service buildings has constituted the most significant change in its total building stock [60-62,99-101].

Estimating the Number of Floors
In our integrated method, the number of floors in a building must be accurately estimated because the structural types, total floor areas and DPMs of certain groups of buildings are estimated directly or indirectly from the estimated number of floors. We estimated the number of floors/stories in a residential and public office building using a correlation with the building's height because the height of a single-floor residential and public office building and the height of each floor in a multi-floor residential and public office building are considerably regular in a certain area. This regularity results from the integrated impacts of long standing dwelling traditions and architectural cultures and modern building codes and urban planning practices of/in a particular local area.
Tangshan, China, is a typical development in North China. The buildings of interest there, such as those in rural and suburban areas, typically follow the North China architectural style [102][103][104][105]. In addition, this area is industrialized and modernized. Therefore, its buildings, particularly those in urban centers, meet modern building construction standards and regulations. To develop general correlations between the heights and numbers of floors of residential and public office buildings in this area, we referred to the Br-LK, which included (but was not limited to) archives, yearbooks, building codes, and urban planning reports. Next, we carefully reviewed and confirmed these correlations during a field investigation that spanned 3 days.
For this fieldwork, we first selected sample building sites based on our reviews or observations of the general features of the assessed buildings regarding their structural types, heights or high/low-rise appearances, and spatial distributions (e.g., the differences between urban and nonurban buildings or areas). These reviews were based on either relevant Br-LKs (such as urban planning materials and building codes) or on the Hr-ORSIs. Overall, our selection resulted in 82 sample sites. These sites were assigned to two categories, urban and nonurban. Next, one or more typical buildings within each site were selected as sample buildings for analysis ( Figure 4). Overall, 235 buildings were selected based on this stratified sampling method. The height of each sample building was measured using a handheld laser rangefinder, and the number of floors and structural type of each building were recorded using a handheld GPS instrument. Based on these data, we developed a correlation between the numbers of floors and the heights of the residential and public office buildings in the study area. An analysis of the arithmetical means indicated that the mean floor height of all sample buildings was 2.94 m. The sample data showed that the total heights (from ground to roof) of single-floor sample buildings and the ground to ceiling heights of the first floor in multi-floor sample buildings were between 3 to 4 m, with most closer to 4 m because the buildings in this area are generally built approximately 1 m above the ground. The box-plots and the means of the total heights of the sample buildings with different numbers of floors showed that the maximum, minimum, and the mean total heights of the single-floor sample buildings were 4.1 m, 2.9 m, and 3.3 m, respectively, and those of the buildings with two floors were 6.8 m, 6.2 m, and 6.4 m, respectively ( Figure 5), and so forth. The regression analysis showed that a good linear relationship existed between the total heights and the numbers of floors in the sample buildings ( Figure 6). This linear correlation is more defined than that observed by Wurm et al. [106] because the buildings in their study included all types of buildings in their study area. In contrast, the buildings in our study only included those for residential and public office usages. The floor heights of the residential or public office buildings within a given area are generally similar, consistent, and stabile.
Thus, we made the following general rules to rapidly estimate the numbers of floors of residential or public use buildings based on their estimated heights in the study area: buildings with an estimated height of approximately 3 to 4 m have one floor, buildings with estimated heights of approximately 6 to 7 m have two floors, buildings with estimated heights of approximately 9 to 10 m have three floors, and so forth. Using these rules and the building heights that were estimated in Section 3.1, the numbers of floors were rapidly determined for all of the assessed buildings in the study area.
Based on the numbers of floors described above and on the estimated footprint areas as described in Section 3.2, the total floor area of a given addressed building was calculated using the field calculator function in the GIS instrument.

Assessing Structural Types
We used the numbers of floors obtained in the last step as a starting point for determining the structural types of the buildings [31,33,48,50,107]. The height or number of floors of a residential and/or public office use building in a specific area usually presented a certain type of empirical correspondence with the building's structure type. This correspondence is the combined product of various local contexts, including a long tradition of dwelling and architectural cultures, current and historical economic development levels, building code requirements, building construction policies, and land-use and urban planning practices.
(1) Because of the overall synchronized and systematic post-disaster reconstruction of Tangshan after the 1976 earthquake, earth/wood buildings (E/W buildings) nearly disappeared in this area [111,113,114]. Among all of the buildings, apart from a relatively small number of specially structured buildings (belonging to the fifth type of Yin's classification system), only three primary types of buildings were identified in the area following the 1976 earthquake, single-story brick buildings (S/B buildings), brick and concrete buildings (B/C buildings), and reinforced concrete buildings (R/C buildings).
(2) Codes for the seismic design of buildings in China, including TJ11-74 [122], TJ11-78 [123], GBJ 11-89 [124], GB 50011-2001 [125], and GB 50011-2010 [126], clearly regulate the maximum number of floors for buildings of various structural types based on certain seismic design requirements. For example, based on the associated construction requirements, the maximum height of B/C buildings in the zones where the seismic fortification intensity is VIII (on the Chinese Seismic Intensity Scale) on the Chinese mainland is six floors. Buildings with more than six floors must be of R/C type. Tangshan experienced a terrible Ms 7.8 earthquake in 1976. Intensity VIII specifies the seismic design intensity/criterion in this area following this large earthquake. Because Tangshan is a relatively economically developed area on the Chinese mainland, relevant building codes have been adhered to well since the 1976 earthquake. Consequently, the maximum number of floors (i.e., six floors) in B/C type residential and public office buildings is well defined in this area. Currently, all residential and Floor number x y public office buildings with more than six floors in Tangshan are R/C type, except for a few very old buildings.
(3) The building uses, ages, spatial distributions and aggregation patterns can serve as a reference for judging the structural types of buildings. In particular, the spatial distributions and aggregate appearances are generally clearly visible in the Hr-ORSIs. This observation is one important reason for our use of these types of images as the basis of our integrated method. Figure 7 displays examples of the 2009 GeoEye-I panchromatic images of the buildings in the Tangshan area. Specifically, Figure 7a shows a group of S/B houses, each measuring 8-9 m on a side with a yard and standing in regular rows. These S/B houses were primarily constructed between 1976 and 1986 [112,114] and are widely distributed across the rural regions of the study area. Figure 7b shows a remote sensing image of a group of slab-type and tower-like buildings. These buildings are rectangular, spaced approximately 20 to 25 m apart, and regularly clustered together. Buildings of this type all consist of two to six floors. Large sections occupied by such buildings were also constructed from 1976 to 1986 [112,114]. These buildings are essentially B/C buildings and are usually used for urban residences. Figure 7c shows various high-rise buildings, most of which are used for residential purposes with some being used as public service offices. Small portions of the first floors of these buildings are used for commercial purposes. These buildings all belong to the R/C category. Most of these R/C buildings were constructed after 2000, are typically 10 floors or taller and are situated in periurban areas or along main urban or periurban roads. A few R/C buildings were built before 2000, are seven or eight stories tall and are located in or near the urban center [112,114]. To further investigate and summarize the correlations between the numbers of floors and the structural types of the residential and public office buildings in the study area, we confirmed the structural types of all sample buildings and their correlations with the numbers of floors during the 3-day field investigation mentioned in Section 4.1 (Figure 4). The following pre-assumptions were made for this analysis: all single-floor buildings are S/B type, all of the buildings with two through six floors are B/C type, and all the buildings with seven or more floors are R/C type, and vice versa. The confusion matrix of the sample data showed that the overall accuracy of these pre-assumptions was 97.02%, with a kappa coefficient of 0.95 (Table 3), which agreed well with our preceding findings presented immediately above from carefully and comprehensively reviewing various Br-LKs and observing the associated Hr-ORSIs (e.g., based on the relevant seismic building codes, the maximum height of the B/C buildings in this area is six floors tall). In fact, it was these preceding findings (esp. those from various Br-LKs) that formed the foundation and boundary conditions of our pre-assumptions.  Based on these findings, we adopted the following general rules to rapidly identify the structural types of all of the residential and public use buildings in the study area: buildings with one estimated floor are S/B, buildings with two to six estimated floors are B/C, and buildings with seven or more estimated floors are R/C. The results indicated that the buildings with relatively good seismic performance in 2009, including the R/C and B/C buildings, were primarily distributed in the urban center, western urban outskirts, and at a few additional locations (including north of the towns of Fengnan and Dazhuangtuo). The S/B buildings were distributed across the rural and suburban regions (Figure 8a). This distribution indicated large structural vulnerability variations from place to place. The most potentially resilient area in 2009 was the urban center, where S/B buildings were rare. The most seismically vulnerable location was the town of Wali, where very few R/C buildings existed at that time ( Figure 8a).
Based on the total floor area of each building estimated in Section 4.1 and on its structural type judged here, we rapidly obtained the general total floor areas of each type of building (S/B, B/C, and R/C buildings) within the entire study area in 2009 (see bars extending to the right in Figure 8d). A comparison of the building stock in the study area in 2009 and in 1976 indicated that very significant growth and changes occurred during the previous decades. For instance, the E/W buildings completely disappeared and the total floor area of all residential and public office buildings in the study area was 8.96 times greater in 2009 than the total floor area of the same buildings in 1976 (shortly before the Ms 7.8 Tangshan earthquake). In addition, the total floor area of all residential and public office buildings in the township of Guoyuan increased by an astonishing factor of 24.82 (Figure 8a-d).

Defining Building Damage Probability Matrixes (DPMs)
The DPMs used in this study were developed in three steps. (1) Yin [49] developed a set of baseline DPMs for Chinese mainland buildings that consisted of several sub-sets based on buildings in various parts of China. In his DPM set, a sub-set was included that was primarily based on the experiences of the 1976 Tangshan earthquake and was specific to the middle region of northern China. Therefore, we first based our DPMs of seismic intensities VI through X (Unless otherwise stated, these intensities and intensity XI discussed in this paragraph and those mentioned throughout the paper all correspond with the Chinese Seismic Intensity Scale.) on this sub-set; (2) Because no matrix with an intensity of XI was included in Yin's matrix set, we used observations of the building damages in the zone subjected to an intensity of XI during the 2008 Wenchuan earthquake [127]. No discussion of the matrix of intensity XI was conducted for the Chinese mainland before this large earthquake, mainly because an insufficient number of the sample buildings that were damaged by this very high intensity; (3) To update our DPMs and make them more consistent with local current conditions, we performed several additional adjustments to the matrixes obtained in the previous two steps. These adjustments were primarily based on relevant findings from several recent earthquake loss estimation studies in cities in the middle region of northern China, including Nanyang, Henan province [128], Handan, Hebei province [129], and urban Tangshan, Hebei province [130]. Our final DPMs are presented in Table 4.
Based on previous results, including the total constructed/floor areas of each type of building for residential and public office usages in the study area in 2009 and their DPMs, we calculated the damage potential of one specific type of building in this area in 2009 under a specific intensity scenario using Equation (8) [49].
where refers to the damage area of -type buildings under a specific intensity scenario ; represents the total constructed/floor area of the -type buildings; | are the buildings' DPMs; represents the specific structural types of the buildings, including the S/B, B/C and R/C types; and denotes a building's degree of damage, which is usually assigned to slight damage, moderate damage, major damage, or total collapse in China [131]. Using the result from Equation (8), the value of in the following Equation (9) can be derived accordingly. This equation expresses the total potential damage to all of the addressed types of buildings in the study area in 2009 under certain specific intensity scenarios. Example calculations are presented in the Appendix of this paper.

Developing Parameters for Calculating Monetary Losses
The direct losses in monetary terms due to building damage are usually assigned to two categories: structural damage to the buildings themselves and damage to the contents of the buildings [7,70,132].
From this standpoint, the expression widely adopted in mainland China for calculating total direct monetary losses is where represents the total direct monetary losses that resulted from seismic damage to buildings under a specific intensity of ; is the replacement cost per unit floor area of the -type buildings expressed in RMB or USD per square meter; ε denotes the loss ratio of -type buildings with the degree of damage of ; is the average/mean monetary value of the building contents per unit floor area of the -type buildings, which is usually expressed in RMB or USD per square meter; and ε indicates the loss ratio of the building contents with a degree of damage of .
Of the parameters in Equation (10), | and were obtained in earlier steps. Therefore, estimating the direct monetary losses due to building damage requires proper determinations of other parameters in the expression, including , , ε , and ε . We determined the values of these parameters by primarily reviewing the Br-LK that was relevant to the study area in 2009. Moreover, to render these parameters more appropriate and reliable, during the 3-day field investigation (Figure 4), we performed small-scale interviews involving a number of local households, public service departments, and building construction companies.
To determine , we primarily referred to 2009 data from the local Economic Statistical, local Price, and local Construction Statistical Yearbooks and used the records of interviews of local households and public service departments. In the published data alone (e.g., yearbooks), abundant records and information exist regarding parameter . For example, this information includes the average/mean amount of major indoor durable goods owned per 100 households, the values of these goods, and the average/mean living space per capita for the various types of buildings addressed. Based on the same rationale, we assigned a value to parameter by primarily reviewing the local Construction Statistical Yearbook, the local Price Yearbook and consultations with local building construction companies. The values of parameters ε and ε were based on the comprehensive analysis of relevant data in the available loss assessment reports of numerous past earthquakes nationwide [133][134][135], in relevant research books [136], and in available earthquake loss estimation project reports [128][129][130]137].
After assigning values to the variables in Equation (10) based on estimates or calculations following the above steps and after selecting various specific intensity scenarios, the potential direct monetary losses due to residential and public office building damage in the study area in 2009 can be displayed. Please see Appendix of this paper for an example calculation.

Field-Based Verifications
From previous sections, the reliability and accuracy of the integrated method developed in the study is primarily determined by the following: (1) the accuracies of the estimated building heights and footprint areas and (2) the reliabilities of the estimated correlations between the building height and number of floors and between the number of floors and structural type in the study area. To test the accuracy and reliability of our estimates, a careful field sampling investigation was conducted.
The first step of this field sampling investigation was to select an adequate set of sample buildings using a mixed sampling method. First, we delineated various sample sub-areas by cluster sampling. As explained earlier, the accuracies of estimated building heights and footprint areas depend on the accuracies of the SIPDs measurements and the roof outlines in the Hr-ORSIs. The primary factors that influence the clear identification and accurate measurement of these two parameters include the following: (1) various shading or interference effects in the images of a building's shadow or roof that stem from the building's surroundings (e.g., trees); (2) high gray-tone similarities between the images of a building's shadow or roof and its background (e.g., water bodies); and (3) a building's irregular or unique shape (relevant to measuring roof outlines only). Based on the significance of these three aspects regarding their combined effects on image interpretation, we divided the GeoEye-1 images spanning the entire study area into three clusters. (1) A low-effect cluster, namely a low-effect sub-area, was characterized by low building densities with very little or no effects from trees or other ground objects. Building-occupied land in the low-effect cluster comprised approximately 30% of the total building-occupied land in the study area; (2) A moderate-effect cluster (i.e., moderate-effect sub-area) was characterized by a moderate building density. A few or several of the buildings' shadows and/or roof outlines were shaded or interfered with, but the distinct roof edge points that were used for determining the SIPDs, such as the corner of a building's roof edge, were always clearly visible. The land occupied by buildings in this sub-area accounted for approximately 55% of the total land occupied by buildings in the study area; (3) A high-effect cluster (i.e., high-effect sub-area) of land occupied by buildings accounted for approximately 15% of the total land occupied by buildings in the study area. This sub-area always had a very high building density and significant shadowing and roof shading or interference effects. The rooftop edge points in this sub-area for determining the SIPDs frequently cannot be easily identified. In some cases, irregular-shaped buildings may be included. Next, we selected specific sample sites and selected a site size by using the stratified sampling concept, which involved two steps. (1) We determined the total size of the specific sample site(s) of each cluster (sub-area). In this case, the building-occupied lands of the low, moderate, and high-effect clusters (sub-areas) comprised approximately 30%, 55% and 15%, respectively, of the total building-occupied land in the study area. To ensure that the building sample is adequate for testing the overall accuracy of the estimated heights and footprint areas of all of the assessed buildings in the study area, we concluded that the total sizes of the specific sample site(s) with low, moderate and high effects should be 2 km 2 , 3.6 km 2 , and 1 km 2 , respectively; (2) Next, we selected the locations of these specific sites. The determining factor was the representativeness of the sites regarding the above-mentioned interference effects, gray-tone similarities, and the number of irregular-shaped buildings. Finally, we designated all of the individual residential and public office buildings in these specific sites as our final sample buildings for analysis, which numbered 954.
The second step was to perform verification through fieldwork. For this step, we used a handheld laser rangefinder to measure the true floor height, total building height and footprint area of each sample building. In addition, we used a handheld GPS to locate the exact spatial positions of the sample buildings and to record their structural types and numbers of floors. The accuracy or reliability of the integrated method developed in this study was evaluated using these data.

Accuracy of the Building Height Estimates
After comparing the true heights of the buildings measured in the field with the estimated heights of the same buildings obtained using the method we proposed in this study, we concluded that the overall accuracy of the building height estimation algorithm we developed was satisfactory (Figures 9-11; Table 5). Overall, a good linear relationship was observed between the measured/true heights obtained from the fieldwork and the estimated heights that were derived from the images (Figure 9). The mean measured height of the sample buildings was 19.72 m, and the mean estimated height of these same buildings was 19.39 m. Furthermore, the mean absolute error (MAE) between the estimated and measured heights was 0.55 m, and the mean absolute percentage error (MAPE) was 3.23% ( Table 5). The mean measured floor height of the sample buildings was 3.09 m, which closely matched our earlier estimated height. In terms of estimation error differences in the buildings with different actual floor numbers, Figure 10 shows that larger errors mainly occurred when estimating the heights of the buildings that were 4 through 6 and 11, 12, 17 and 18 floors in height. Regarding the differences between the measured and the estimated height of each individual sample building, only 27 buildings had absolute differences that were near or greater than 3 m, including only 2.83% of the 954 sample buildings (Table 5) and occurring at the moderate and high effect sites (Table 5; Figure 11). In this integrated method, the height of a building is an important reference by which a building's floor number is judged. The large differences between the buildings' mean height (19.72 m) and/or mean floor height (3.09 m) that were measured in the field and the estimating errors (e.g., MAE of 0.55 m) indicated that these estimation errors did not significantly influence the reliability of using these estimated heights to judge the number of floors in a building, except in a few extreme cases. Such an accuracy level is high enough to produce reliable estimates of the number of floors in a building by using the estimated building heights. Figure 12 provides several examples of measured vs. estimated height differences that were much larger than 3 m, which could yield an incorrect estimate of number of floors. The building with the largest absolute height estimation error was building No. 189 at the high effect site. Figure 12a shows that most of its shadow was blocked by a building in the front, which significantly distorted its shadow, made it difficult to identify a reliable SIPD and yielded a large height estimation error (15.75 m). Figure 12b shows the cases in which large portions of the buildings' shadows were projected onto surrounding trees, which rendered the shadows of the upper potions of the buildings nearly indiscernible. This type of interference can also produce serious estimation errors. However, these extreme instances of shading or interferences were rare across the study area and were primarily observed in high-effect sub-area ( Figure 11) in which the land occupied by buildings only accounted for approximately 15% of the total land occupied by buildings in the study area.    1) On all x-axes, the sample buildings at each site were sorted based on their height. Specifically, the buildings were arranged from low to high, with larger building (serial) numbers corresponding with larger actual/measured building heights; (2) Relative errors were obtained by subtracting the measured heights from their estimated counterparts. Thus, one can observe that a systematical under-estimation occurred in the building height estimations in this study, especially at the moderate and high effect sites. Note: MAE represents the mean absolute error; MAPE represents the mean absolute percentage error; r is the Pearson's product-moment correlation coefficient; and RMSE is the root mean square error.

Accuracy of the Building Footprint Area Estimates
The overall accuracies of the buildings' estimated footprint areas were satisfactory (Figures 13 and 14; Table 6). Overall, a good linear relationship was observed between the measured/true footprint areas obtained in the field and the estimated footprint areas derived from images ( Figure 13). The average measured/true footprint area of the sample buildings was 622.77 m 2 , and the average estimated footprint area was 623.55 m 2 . The MAE between these measured and estimated areas was 16.54 m 2 , and the MAPE was 2.99% (Table 6). Regarding the differences between the measured and the estimated areas of each individual sample building, the buildings that had relative differences of less than 2.99% accounted for 63.94% of the total sample buildings. Only 43 buildings had relative differences that exceeded 10%, which accounted for 4.51% of all of the sample buildings ( Table 6). The larger errors nearly all occurred at the moderate and high effect sites (Table 6; Figure 14).    (Figure 14), which mainly occurred because the roof of the building was in the shadow of a building in front of it (No. 13) (Figure 15a). Apart from image shading, interference, or the high gray-tone similarity effects, the irregular shape of a building's roof was another source of error when estimating building footprint areas. Figure 15b shows examples of this error for two rows of uniquely shaped buildings. It was not easy to precisely digitize these irregular roof outlines, which means that certain errors were inevitable. However, the occurrence of such extreme errors when estimating building footprint areas were rare, as were the extreme errors that occurred when estimating the building heights as discussed above. x y Figure 14. Relative differences between the measured and estimated footprint areas of the sample buildings. Note: (1) On all x-axes, the sample buildings at each effect site were arranged by footprint size, with larger (serial) numbers of buildings corresponding with larger building footprint areas; (2) Relative errors were obtained by subtracting the measured footprint areas from their estimated footprints. Note: MAE is the mean absolute error; MAPE is the mean absolute percentage error; r is the Pearson's product-moment correlation coefficient; and RMSE is the root mean square error. Figure 15. Images of buildings with large footprint area estimation errors.

Reliabilities of the Estimated Correlations between Building Height, Number of Floors and Structural Type
In this study, we established two sets of general rules using the proposed integrated method composed of a series of estimated or prescribed correlations regarding building height, the number of floors and structural type (Sections 4.1 and 4.2) to rapidly and generally determine the floor numbers and structure types of all residential and public office buildings in the study area. The data obtained in the field can also serve to test the reliabilities of these general rules or estimated/prescribed correlations.
In terms of this specific test, the 954 buildings used for field-based verification are actually a simple random sample.
The field data from these sample buildings showed that all of the single-floor sample buildings were between 3.00 and 3.90 m tall, all of the buildings with two floors were 6.4 to 6.8 m tall, and so forth. After comparing these actual relationships between the number of floors and the building height with our pre-set rules or estimated/prescribed correlations (e.g., buildings with an estimated height of approximately 3 to 4 m have one floor), it was found that these pre-set or pre-assumed rules agreed well with the corresponding actual relationships. The confusion matrix showed that the overall agreement between our pre-set rules/estimated correlations and the actual relationships was 91.61%, with a kappa coefficient of 0.89 ( Table 7).
The field data of the sample showed that all of the single-floor sample buildings were of S/B type; 95.16% of the sample buildings with two to six floors were B/C type, and 4.84% of the sample buildings with two to six floors were R/C type (Table 8). In addition, all of the sample buildings with seven or more floors were R/C type, except very few special cases of buildings taller than seven floors that were constructed before 2000, which were B/C type (Table 8). By comparing these facts with our estimated/pre-set correlations (i.e., the residential and public office buildings with one estimated floor are all S/B type structures, the buildings with two to six estimated floors are all B/C type, and with seven or more estimated floors are all R/C type), we found that the overall agreement between our general and rapid building structural type estimate rules and the corresponding facts was 95.81%, with a kappa coefficient of 0.88 (Table 8). This value was slightly higher than that regarding the rapid estimation of the buildings' floor numbers. The errors in this estimation primarily occurred when estimating the structure type of the building with six or fewer floors. This estimated bias occurred because some of the relatively lower buildings are also R/C type structures in the study area due to recent economy development. Table 7. Summary of the confusion matrix between the actual floor numbers of the sample buildings and those estimated by our pre-set rules.

Actual Earthquake Loss Estimation Project-Based Verification
In 2009, an actual earthquake loss estimation project that used the conventional data collection method was launched in the study area [130]. This project is one important reason for our selection of 2009 as the year to develop our integrated method. The target area of this project was the central part of Tangshan (i.e., its most urbanized area). This area overlapped with the area of the present study but was much smaller. This project not only included the building types that we addressed in our study, specifically the R/C, B/C, and S/B buildings, but also involved special and less prevalent structures, including steel buildings, single-floor plants, single-floor spacious buildings (e.g., cinemas), and temporary buildings. This project spanned more than three years and was completed in late 2012.
This verification method focused on comparing the constructed/floor areas of the R/C, B/C, and S/B buildings that were established by Li et al. with the corresponding measurements that were estimated in this study. We found that the relative difference between these two studies regarding the total floor area of all three types of buildings was small (i.e., −3.996%). The relative differences in the R/C and B/C buildings were larger, reaching −16.667% and 12.029%, respectively (Table 9). These two larger differences were observed because a small number of newly constructed buildings with six or fewer floors were assigned to the R/C category in the work of Li et al. but were generally classified as B/C structures in our study.
Given this particular verification, the accuracy and reliability of the integrated method proposed in this study for rapid and large scale estimations of seismic risk to the residential and public office buildings are considered satisfactory and can attain the accuracy levels provided by conventional data collection methods. Table 9. Comparisons between the determined building floor areas of Li et al. [130] and the floor areas that were estimated in the current study. * In China, particularly in its eastern half, many single-floor houses consist of two parts, the main house and its auxiliary structures. Among these houses, the main house is used as a living room, and the auxiliary structures are often used as kitchens, storerooms, and so forth. In general, the main rooms are taller and much larger than the auxiliary rooms. In actual pre-earthquake loss estimates in China based on conventional data collection, these auxiliary structures have often been neglected. This exclusion occurred in the research conducted by Li et al. [130]. This oversight is due to the much lower importance of these smaller auxiliary structures. More fundamentally and frequently, this omission is due to the difficulties that occur when gathering pertinent data using conventional methods. However, in our present study, these auxiliary structures were included and treated as S/B buildings. The total floor area of these auxiliary structures obtained using our integrated method was 2,126,325 m 2 . After evaluating these differences based on the floor area of the S/B buildings and on the total floor area of all of these three types of buildings, this value was excluded from our data when performing related comparisons to ensure that the findings of Li et al. could be compared with those of this study.

Simulating the Impacts of the Ms 7.8 Tangshan Earthquake of 1976
The Ms 7.8 Tangshan earthquake of 1976 claimed over 240,000 lives and caused over 0.82 billion USD (in constant 1990 valuations) in total property losses [110,113]. This earthquake is the deadliest earthquake in the history of modern China. In recent years since this earthquake, the following question has often been posed: What would occur if an earthquake of similar magnitude recurred in the 21st century in the same place, where rapid, or even accelerated, socioeconomic growth and fixed assets accumulation have lasted for many years? However, there is no specific answer to this question because ample and detailed data, applicable methods, or even a basic understanding of this extremely difficult problem are lacking.
The integrated method developed in this study can contribute to solving this complex problem to a certain degree. Therefore, a simulation of the impacts from the recurrence of this large earthquake was performed. This simulation provided a basic understanding of what would occur if a strong earthquake were to occur in the same location today. In addition, this simulation offered an opportunity to apply the proposed new method. For this simulation, we assumed that a similarly strong earthquake recurred in 2009 at the same location as the 1976 earthquake. The location of this simulation corresponded exactly with the area of the current study (i.e., the central part of the disaster area associated with the Ms 7.8 Tangshan earthquake of 1976) (Figure 16). The scenario intensities were assumed to be the same as those experienced in the study area in 1976 ( Figure 16). The affected objects of the simulation only consisted of R/C, B/C, and S/B buildings. Moreover, to fully use the results of this current study, we combined the data in Table 3 with similar data from Section 5.1. Based on this combination, we determined that all of the single-story residential and public office buildings addressed in the study area were S/B type. Of the buildings with two to six floors, we determined that 95.20% were B/C type and 4.80% were R/C type. Of the buildings with seven or more floors, we determined that 1.52% were B/C type and 98.48% were R/C type. We developed certain impact simulation scenarios after assigning specific values (e.g., the total floor areas and DPMs of each type of building) that were obtained earlier in this study to the corresponding variables in Equations (8)-(10) and after considering the actual intensities that occurred in 1976. The results indicated that if the 1976 Tangshan earthquake had recurred in 2009, the total damage to the residential and public office buildings in the study area would have been 8.00 times greater than in 1976 (Figure 17a), despite the significant structural improvements that were made to these types of buildings since 1976 (Figure 8d).
The total direct economic losses due to damaged buildings for these two usages would have been 28.73 times greater in 2009 than in 1976 (Figure 17a). The spatial distributions of these simulated impact scenarios and their contrasts with the spatial differences of the actual losses in 1976 are shown in Figure 17b,c. The most notable difference was the large increase in the loss risk potential in the periurban area and in several "satellite" towns. For example, the simulated building damage and direct economic loss in the Guoyuan area in 2009 were 23.57 and 95.00 times greater than their respective values associated with the 1976 event (Figure 17a).
This particular simulation only addressed the increases in the seismic disaster risk potential due to building stock growth. To our knowledge, building damage is a beginning component and only one component of a complete earthquake disaster scenario. Therefore, these simulated risk increases only represent a first chapter and a portion of the area's rapid seismic disaster risk accumulations that have resulted from rapid socioeconomic development. For several decades, GDP growth has been a top priority in the study area. This building-focused simulation is a strong warning to the booming economy and the crowded residents in this area regarding the current and future challenges that they may face.

Discussion
The previous sections of the paper described the development and application of the proposed integrated method. This section will discuss the general advantages, limitations, application prospects, and near-future extendible or scalable directions of the method.

General Advantages
The first advantage of this integrated method is its speed and labor savings compared with conventional LRC or DPM-based quantitative earthquake loss estimate methods. This advantage, particularly its speed, is highly critical for analyzing seismic loss risks in an area where socioeconomic growth is rapid. We based this study on the central Tangshan area, which contains densely packed buildings, to demonstrate the development and application of this method. The study area measured 456 km 2 , which greatly exceeded the areas considered in all other conventional earthquake loss estimation projects conducted in Mainland China that we reviewed [128][129][130][137][138][139][140]. After the Hr-ORSI and Br-LK sources were compiled, the potential seismic risks to the residential and public office buildings in this study area were calculated over approximately ten working days by two skilled individuals. In contrast with durations of several years (typically 2 to 3 years in mainland China) and the great deal of labor involved in conventional earthquake loss estimation projects, the speed and labor savings of this method stand out. This advantage is ensured by three aspects of the proposed integrated method. (1) The Hr-ORSIs sources are employed as a starting point to replace time and labor-intensive conventional methods of collecting building data and preparing building inventories; (2) Next, the Br-LK sources are reviewed to further accomplish this replacement. As is commonly known, large gaps exist between the information that Hr-ORSIs can directly provide and the specific parameters or data that specific seismic loss calculations for buildings require. Thus, a Br-LK review was proposed in this study to bridge these gaps. In other words, the Br-LK review was introduced to tightly link up remote sensing technologies with technical cores of mainstream quantitative seismic loss risk estimates. Our work, which employed the Br-LK review method to help determine the required parameters, is similar to the studies conducted during the previous decade that summarize baseline building DPMs for various regions or locations in China [71,141]. Among these studies, Yin and Yang [71] divided the entire Chinese mainland into three large regions and developed a set of baseline building DPMs for each of these regions. Hu et al. [141] attempted to sum the baseline building DPMs of much smaller and specific locations (e.g., various counties). The difference between these studies and our work is that we addressed, for example, the correlations between building heights and the number of floors and between the number of floors and the structural type in addition to developing location-specific DPMs. The current study indicates that this review is not a time-or labor-consuming task once the various Br-LK sources have been compiled; (3) The speed of the integrated method is also evident in our new algorithm for the large-scale estimation of building heights from single Hr-ORSIs. Although this algorithm involves manual work, it does not require substantial labor and time.
In addition to being rapid and saving labor, this integrated method provides an accuracy that is comparable with that of conventional quantitative earthquake loss estimation methods, as indicated by our relevant detailed verification presented earlier in this paper. This accuracy is primarily a function of the following two aspects of this integrated method: (1) The methods by which the heights and footprint areas of a large number of buildings of interest were estimated, including the new algorithm for extracting building heights and the old but realistic manual methods of obtaining building footprint areas that we adopted for the purpose of actual application; (2) The comprehensive and careful manner in which various Br-LKs were reviewed to obtain reliable correlations between the building height and the number of floors and between the number of floors and the structural types, reliable DPMs of the buildings, and reliable economic parameters for calculating monetary loss due to building damage specific to a certain area and a certain time context. Similarly to the location-specific baseline DPM summarization work in China that was conducted by the Chinese earthquake engineering experts mentioned above, this approach ensures that the reliabilities of these correlations and parameters are reasonable. Compared with the extreme difficulties of preparing and updating the building data and inventories of large areas with large numbers of buildings and rapid socioeconomic development using conventional methods, establishing and updating a small specialized database consisting of various Br-LKs of such an area is much more practical and easier. By using such a database and increasingly available and advanced Hr-ORSIs, the correlations, metrics, and parameters of the buildings in the area (e.g., correlations between the number of floors and structural types) that are directly or indirectly required in conventional seismic building loss risk calculations can be conveniently updated. Therefore, if a comprehensive Br-LK database of certain areas can be developed and maintained properly, the proposed integrated method may be applied more conveniently in these areas in a timely and accurate manner.
In addition to offering noteworthy speed, substantial labor savings, and good accuracy, this integrated method is also less costly than conventional seismic loss estimation methods, although Hr-ORSIs are still expensive. In this study, we employed GeoEye-1 images with panchromatic and multispectral resolutions of 0.41 m and 1.65 m, respectively, and with coverage spanning an area of 456 km 2 to demonstrate how this method is developed and applied. The cost of these images was approximately $9000, which accounted for nearly 90% of the total financial investment of this study. By comparing this current study with certain conventional earthquake loss estimation projects performed on the Chinese mainland [128][129][130]137,139,140], we noted that our study area was larger (in most cases, much larger) than the target areas considered in each of these conventional projects. However, the total financial investment of the integrated method addressed in this current study was less than one-tenth of the mean cost of these conventional method-based earthquake loss estimation projects.

Applicability to Other Locations
The proposed integrated method begins with estimating the building heights by using Hr-ORSIs. After considering the availability and costs of these images, we developed an algorithm that is based on single images to accomplish this task. This algorithm is theoretically sound if it is assumed that the targeted buildings are located on a level site (i.e., a site with flat terrain). Otherwise, the height estimates would be biased. For example, if the imaging points of in Figure 3 were projected onto a hillside or a low-lying area, the SIPDs would be somewhat shortened or lengthened. Therefore, this building height estimation method is only applicable in locations that are relatively level/flat, including plains, the low-lying parts of basins, the flat parts of plateaus, and even certain locally flat locations in mountainous and hilly regions.
The specific correlations between the numbers of floors, building heights and structural types and the location-specific DPMs of the residential and public office buildings established in the context of the study area can be applied to the same building groups in other areas of Tangshan. In addition, these correlations and DPMs are probably reliable in some other parts of North China and in some of its neighboring areas, including the cities of Beijing and Tianjin, the south-central region of Inner Mongolia, and the provinces of Hebei, Shanxi, Henan, and Shandong. This trend occurs because the longstanding dwelling and architectural traditions and modern construction practices in these areas (as in the Tangshan area) all essentially belong to the North China architectural style [102][103][104][105]. This consistency was also used by Chinese earthquake engineering experts to develop baseline DPMs for these locations [71,141]. The probable and general reliabilities and applicabilities of these correlations and DPMs across these mentioned areas are also supported by several seismic loss estimation studies that have used conventional building inventory preparation methods in recent years in the southern parts of the North China Plain, such as Nanyang in Henan province [128] and Handan in Hebei province [129]. For example, the data from the Nanyang area indicated that 92.6% of the buildings with seven or more floors are R/C type buildings and 96.4% of the buildings with two to six floors are B/C type buildings, which is consistent with the correlations in Tangshan. Tangshan and its surrounding areas underwent a large-scale, unified, systematic reconstruction campaign after the 1976 earthquake; thus, the buildings in and around Tangshan are more uniform than those in other parts of North China and its neighboring locations. However, this uniformity does not necessarily mean that the correlations and DPMs from this region are not applicable to these other locations. First, the high uniformities in Tangshan are mainly manifested as consistencies or similarities regarding certain aspects of land-use and land coverage (e.g., certain lands used for housing) rather than, for example, as a correlation between building height and the number of floors. These land use/coverage-focused consistencies allow us to process building images in bulk (i.e., in location-specific building groups). Second, since the completion of the post-disaster reconstruction after the 1976 earthquake (middle to late 1980s), the buildings in Tangshan have undergone significant changes, particularly during the previous decade, as have experienced the buildings in many other regions of the Chinese mainland. These changes have altered these land use/coverage-focused uniformities (i.e., not as uniform as they were before). Some information regarding building stock changes in the study area was presented earlier in this paper.
Except in North China and in some of its neighboring areas, the specific correlations regarding the building heights, the numbers of floors, and structural types and the DPMs of residential and public office buildings established in the context of the study area most likely do not conform to the correlations of the buildings with the same usage designations at different other places on the Chinese mainland, especially in villages and small towns and particularly regarding height/floor number-structure type correlations and DPMs. This difference occurs because the long-existing dwelling and architectural traditions and modern construction practices in these other areas belong to other large architectural styles (i.e., Northeast China, Northwest China, East China, and Southwest China style) [102][103][104][105]. For example, thicker walls are used in cold northeast areas of China [142], clay is often used to make walls for public houses in the vast dry areas in northwestern China (esp. rural areas) [102,143], and hollow bricks are frequently used to build residential houses in many areas of the Qinghai-Tibet Plateau due to either the lack of natural resources (i.e., clay) for locally producing regular brick or because of the high costs of transporting in the bricks [144].
Furthermore, the height/floor number-structure type correlations and the DPMs of the residential and public office buildings (esp. the former) in many other mainland areas of China, especially in the western half, are likely not as clear or regular as those shown for Tangshan. Such a likely complex is associated with the more common co-existence of multi-cultural dwellings and architectural conventions, less developed land and urban planning practices [105,145], and/or other actualities in these western places. In view of more insufficient exploration than that conducted on DPMs (e.g., Yin and Yang [71] have summarized several baseline building DPMs for different large regions of the Chinese mainland; and Hu et al. [141] attempted to develop baseline DPMs for smaller areas), examining the height/floor number-structure type correlations is probably an important first step for applying the proposed integrated method in this type of area. The Br-LK review method is probably an effective method for solving this issue. As demonstrated previously, Br-LK reviews can help address many aspects of such correlations of a place. For example, (1) in terms of the Chinese mainland, Br-LK reviews can help by judging whether the overall height/floor number-structure type correlations in the residential and public office buildings of a specific area of interest are complex or not. Generally, the areas with good land and urban planning practices, mono-dwelling and architectural cultures, and relatively developed economies present relatively clear correlations, and vice versa; (2) For a specific Chinese mainland location, Br-LK reviews can further help check buildings for residential and public usages, for which the building heights may present complex relationships with structure types. In general, apart from a few exceptions that usually result from very old buildings, the residential and public office buildings on the Chinese mainland with seven or more floors in areas with a seismic fortification intensity of VIII should be R/C type buildings according to the relevant seismic building codes. However, varied relationships in such areas often exist for buildings with six or fewer floors. These relatively short buildings, especially those with one to three floors, were mainly built by the public by "acting on their own", and usually lacked unified planning and management. Consequently, wealthy families (even wealthy institutions) can build houses (or offices) that meet higher standards than poor families (institutions). As presented in the verification portion of this paper (Section 5.2), such variations in buildings with six or fewer floors also exist in the study area, but are not significant. Thus, building height or the number of floors are generally critical indicators of the buildings' seismic structural type, but not always and not all. (3) Br-LK reviews can considerably complement such indicators or references, such as building uses, ages, and so forth, as discussed in Section 4.2 of this paper. Systematically compiling and reviewing the Br-LKs of places characterized by varied or complex floor number-structure type relationships, and subsequently, carefully developing certain specific multi-indicator-based methods to better examine these relationships may strongly promote broader applications of the proposed integrated method.
For this particular issue, observations from remote sensing images could be very helpful. First, as mentioned previously, Hr-ORSIs can directly provide rudimentary, general and wide-reaching impressions if the correlations in a specific area are complex or simple (for examples, see Figure 7). More importantly, through careful analysis and data-mining, remote sensing images can also provide many of these multi-indicators for such correlation analyses. Apart from heights, these indicators include many other features of the buildings themselves and of the building environments that can be derived from remote sensing images. Regarding the issue of this aspect, Geiß et al. [35] recently conducted an enlightening attempt to determine how buildings' structure types could be statistically correlated with the various remote sensing features of these buildings themselves and of these buildings' surrounding environments. If this valuable attempt can further combine other useful information that can be derived from Br-LK reviews, including those of the buildings themselves (e.g., building archives and urban planning reports with recorded usage, age, and ownership, seismic building codes, and other relevant regulations and policies) and those of the buildings' context (e.g., local land and urban planning, local economy and local dwelling conventions that are contained in various local documents), developing certain multi-indicator-based methods (as an alternative to conventional earthquake engineering method) to rapidly examine building structure types may become more promising. This scenario may occur because some of the Br-LK-derived indicators could be boundary or controlling indicators, for example, seismic building codes of China clearly regulate that the residential and public office buildings with seven or more floors in areas where the seismic fortification criterion is intensity VIII should be R/C type on the Chinese mainland.
Because of the various economic conditions that are prevalent in various locations, the 2009 parameters we developed based on the central portion of the Tangshan administrative region for calculating monetary losses due to building damage may not be applicable to all of North China and its neighboring locations and are probably much less applicable for other areas of Mainland Chinese, especially the western areas. Among these locations and areas, some are wealthier, more industrialized, and more developed than other areas, resulting in variations in building contents and building replacement and/or repair costs. However, these economy-relevant parameters are easier to obtain from Br-LK reviews than are to obtain those regarding the physical buildings themselves (e.g., the height/floor number-structure type correlations and DPMs). Therefore, this aspect of analysis should not serve as a critical obstacle when applying this integrated method across the Chinese mainland.
Although the proposed integrated method was developed in the context of the Chinese mainland, its applicability is not restricted to the Chinese mainland. Rather, the general ideas and overall guidelines that closely and organically combine Hr-ORSIs and various Br-LKs to address the technical cores for quantitatively and specifically estimating seismic loss risks to buildings are likely universally applicable. However, except the worldwide applicability of the new algorithm developed herein for estimating the heights of buildings, all of the other key elements that underlie this integrated method are highly likely to vary, even sharply, from one location to another. These varied elements include the correlations between building heights and numbers of floors and between numbers of floors and structural types, DPMs, economic metrics for calculating monetary losses due to a buildings' damage, and the specific forms and contents of Br-LKs that serve as the bridging sources for determining these elements. To apply this integrated method to an area outside of the Chinese mainland, Br-LKs must be accumulated for the area first, and then using Br-LK reviews to carefully examine the associated elements specific to that area. This general integrated approach can be employed when addressing earthquake loss risks outside of China, where active seismic activities occur in areas with rapid socioeconomic developments. The overall guidelines of the integrated idea could also be tested in earthquake-prone areas outside of China, where socioeconomic changes are not as rapid. However, colleting and/or updating detailed building data through conventional methods remain a significant challenge.

Applicability to a Range of Building Types
The proposed integrated method is applicable to the first four types of buildings in Yin's structural vulnerability classification system [68]. These buildings primarily consist of residential buildings and public service offices. Our algorithm for estimating building heights using single Hr-ORSIs is applicable to the fifth (i.e., the last) type of building (i.e., specially structured buildings for commercial, industrial, or other special use). However, the DPMs, correlations, and parameters involving building heights, numbers of floors, structural types, and the economic metrics for calculating the monetary losses that were established in this study are not applicable to this fifth type of building. (1) The correlations of the building heights, numbers of floors, and structural types are significantly different for the fifth type of building relative to buildings used for residential purposes or public service offices; (2) Particularly, the DPMs of these specially structured buildings are unique. Moreover, inadequate DPMs exist for many of these special structures and inadequate baseline DPMs exist that could be used to analyze such buildings in groups; (3) Furthermore, the contents of such buildings differ substantially from those built for residential and public service purposes. These three aspects are beyond the scope of this paper and many of them may be difficult to resolve based on currently available science and technology knowledge and actual seismic damage experiences. Thus, it is still very difficult to determine specific seismic loss risks rapidly, accurately, and at a large scale for this fifth type of building relative to the first four types of buildings. As a result, the applicability of this integrated method inevitably decreases to a certain degree in certain modern urban centers or historic places, such as certain locations in Beijing and Tianjin, where many are specially structured buildings.

Applicability to Various Scenarios
From a disaster cycle perspective, the integrated method proposed in this study is primarily intended to serve during normal, or non-disaster times and may be referred to whenever needed during such a period. For example, we could use the method to yearly estimate the loss risks to residential and public office buildings in each Yearly Key Seismically Risky Area on the Chinese mainland (At the end of each year in China, the China Earthquake Administration (CEA) identifies a certain number of specific locations with relatively high levels of seismic activity as Key Seismically Risky Areas for the following year. Usually, several tens of locations are identified with spatial coverage of each location that usually spans several counties.) [146]. This type of analysis would be very helpful for making yearly adjustments to existing hazard mitigation and disaster preparedness plans in the areas or for developing new plans where necessary. Because of its speed, this new method could also play an important role before an imminent earthquake. The Law of the People's Republic of China on Protecting Against and Mitigating Earthquake Disasters (2008 edition) has clearly stipulated that "After the release of earthquake forecast opinions, the people's governments of relevant provinces, autonomous regions and municipalities directly under the Central Government may, in accordance with the earthquake situation forecasted, declare that the relevant areas enter an emergency period of an imminent earthquake; the relevant local people's governments shall, in accordance with the earthquake emergency plan, organize the relevant departments to do a good job in the emergency prevention and preparations for earthquake relief." This scenario is very rare due to the extreme difficulties of short-term and imminent predictions of actual earthquakes and the difficulty of officially issuing an earthquake forecast to the public. However, such forecasts have been made several times for the Chinese mainland. Among these, the most famous and most successful forecast involved the Ms 7.3 earthquake in Haicheng, Liaoning province, in 1975. Approximately nine hours before the earthquake, a timely official forecast and warning was broadcast [147,148]. It was reported that this well-timed forecast and warning saved over 100,000 lives [149]. If a similar situation occurs again in the future, we are certain that this integrated method could play a role in providing more timely and specific preparations.
In addition to these two scenarios in a disaster cycle, the proposed integrated method can also be tested in the following two cases. (1) As in the particular simulation involving the Ms 7.8 Tangshan earthquake of 1976, this method can be used to simulate at least a portion of the modern effect scenarios of any historical earthquake if its ground motion intensities are available. This simulation helps form a conceivable and convincing reference for relevant areas to design more evidence-based seismic impact lessening strategies; (2) In addition, with the multi-temporal Hr-ORSIs and Br-LK, this method can serve as a tool for evaluating the year-to-year changes in seismic disaster risks and exploring their developmental mechanisms in different areas, thereby helping to adjust the corresponding disaster risk reduction strategies.

Applicability to Larger Areas
Theoretically, this integrated method is applicable to a level/flat location with any geographic size. However, this method likely cannot be applied to areas that are too large because Hr-ORSIs are currently expensive. Targeting areas with geographical sizes that are too large may be financially unrealistic. In this paper, the study area measured 456 km 2 , which exceeded (in most cases, greatly exceeded) the areas considered in each of the conventional earthquake loss estimation projects that we reviewed. The total expense for the method addressed in this paper was less than one-tenth of the mean cost of these projects conducted on the Chinese mainland, and nearly 90% of the total cost was for the GeoEye-1 images. In the economic context of the Chinese mainland, it may be financially feasible and acceptable to use GeoEye-1 images to provide coverage of areas measuring up to approximately 4560 km 2 (i.e., approximately 10 times of our 456 km 2 study area). Thus, from the perspectives of the current economic context of the Chinese mainland and the prices of Hr-ORSIs that are currently prevalent there, the proposed integrated method involving remote sensing images with meter or sub-meter resolutions may be applicable to flat portions of the Chinese mainland that cover up to several thousand square kilometers.
In the Chinese mainland, the application of this integrated method that employs remote sensing images with meter or sub-meter resolutions as underlying data to areas that are much larger than this limit is likely to be financially difficult. However, this integrated method may still play an important role in these cases. For example, this approach can be used to assess the risks to important but small sites in large target areas, thereby helping to optimize the allocation of risk reduction resources.

Future Directions
Obtaining timely, accurate and large-scale characterization of seismic loss risks and determining their changes in areas of rapid socioeconomic growth is a work in progress. Although many challenges await future solutions, we focus on the following four issues.

Possible Methods for Estimating Building Heights in Hilly Areas
The estimation of building heights using Hr-ORSIs is an important initial step in the proposed method. Building height is one key index for directly or indirectly obtaining data, including the number of floors, floor area, structural type, and DPM of a building. However, the algorithm we developed in this study for estimating building heights using single optical images alone is only reliable for buildings located on level/flat ground, such as plains, where the distortion of a building's image by the local topography is minimal. However, in mountainous or hilly areas, serious image distortions due to the local topography are highly likely (e.g., in images of buildings next to a steep hillside). To solve this problem, multi-source remote sensing data fusion techniques or stereo images could be introduced. However, these techniques may result in a greater monetary investment, which is a large problem when the area of interest is considerably large and contains a large number of packed buildings. Thus, the solution to this hilly area-focused issue may be to apply the proposed algorithm to buildings in the flat portions of the areas of interest. For buildings located in the portions with rugged topography, multi-source remote-sensing data-fusion techniques or stereo images should be applied.

Accumulating Br-LKs of Other Earthquake-Prone Areas of China
Except in North China and in some of its neighboring areas, correlations between parameters such as the heights, numbers of floors, and structural types of buildings; the DPMs; and the coefficients for calculating the monetary losses that developed in this study in Tangshan are unlikely to match those of other Chinese mainland regions, particularly those of rural and suburban locations. These differences from one region of China to another have evolved over many years due to either the diversities of local architectural styles and cultures or to the differences of local macro and micro economic conditions. However, the levels of seismic activity in certain regions, such as some regions of western China, can significantly exceed those in North China. Moreover, as in North China, nearly all of these regions have experienced and are experiencing rapid socioeconomic development. This situation implies that the work similar to this study that targets these places is even more urgent. Acquiring and accumulating such purpose and location-targeted sources of Br-LKs (1) to enable careful determination of required parameters; correlations, and coefficients and then (2) enabling timely and accurate large-scale assessments of seismic building loss risks and their changes in such regions will be one of our next efforts.

Examining the Prospect of Using Images with Lower Resolutions in the Proposed Method
The proposed integrated method employed remote-sensing images with meter or sub-meter resolutions (i.e., very high-resolution optical remote sensing images) as one source of input data. Currently, these very high resolution images are considerably expensive. Thus, it is worth considering the use of images with lower resolutions to decrease costs. The first priority of such a replacement should be that the images have a resolution that is exactly next to that of the highest resolution images, such as the Chinese-Brazil Earth Resources Satellite GF-1 (CBERS GF-1) images with a panchromatic resolution of 2 m and a multispectral resolution of 8 m [150], the SPOT-5 images recorded at a panchromatic resolution of 2.5 m [151], and the FORMOSAT-2 images with a panchromatic resolution of 2 m and a multispectral resolution of 8 m [152]. These images are the most promising for playing the same roles as those with meter or sub-meter resolutions. Particularly, these images are more easily available and cheaper (sometimes even free), thus possess the potential to allow the proposed integrated method to be more broadly applicable because covering a larger geographic area is less expensive and a more flexibility is allowed when using existing imagery resources. However, for assessing seismic loss risks to buildings, several uncertainties remain when using these images with slightly lower resolutions. The focus of these uncertainties is whether these images can ensure that the "revised" integrated method produces an assessment result with a similar level of accuracy. Specifically, focus is given to whether the buildings' heights and footprint areas can be extracted one by one from these images as accurately as expected when using the newly developed methods in this study (for assessing heights) or the already existing methods (for assessing footprint areas). If not, what could solve this problem? Thus, examining these pending uncertainties is an important future direction of our research.

Testing the Applicability of this Hybrid Method to a Highway System
Acquiring knowledge of building loss risks and their changes alone does not meet our goal of a comprehensively understanding earthquake loss risks and their rapidly changing scenarios in areas of rapid socioeconomic development. We hope that the hybrid approach and its general guidelines presented in this paper, which combine remote sensing imagery and local knowledge and employ a careful Br-LK review to bridge the gaps between remote sensing and technical cores of conventional seismic loss risk estimations, can address the many challenges posed by this goal. Next, we aim to test the applicability of this hybrid method to the seismic risks of highway systems. According to the basic logic of this integrated approach, the highway system is similar to buildings in several respects.
(1) First, roads and bridges are clearly visible in high-resolution remote sensing images; (2) As building-relevant local knowledge (Br-LK) was used to assess, for example, building structural types, highway system-relevant local knowledge (HSr-LK) (e.g., transportation yearbooks) could be used to assess information regarding highway quality and quantity that is critical for estimating potential earthquake-induced highway damage; (3) Significant and rapid changes have occurred to the highway system on the Chinese mainland during previous decades [153,154]. This integrated approach is likely to play a role in the large-scale rapid estimation of seismic loss risks to highway transportation systems. Detailed and specific confirmation in this regard is worth careful examination in the future.

Final Remarks
Rapid socioeconomic growth in an earthquake-prone area can cause rapid seismic loss risk changes. These changes challenge seismic disaster risk reduction strategies in the area, which must be updated frequently to remain realistic, practical, and effective over time. The top priority for overcoming these difficulties is to capture ongoing rapid risk changes in a timely, definitive, and accurate manner. However, methods for rapidly and accurately characterizing such seismic disaster risks over a large area have been insufficiently explored. This current study developed an integrated or hybridized method that combines remote sensing data and local knowledge to resolve this issue. The focus of the proposed integrated method is not solely on remote sensing technology or independently on mainstream technical cores of LRC or DPM-based seismic loss risk estimates. Instead, this method focuses on how the large gaps between these two sides may be bridged properly and quickly at large scales. This bridging topic is an issue that has been insufficiently explored relative to the fast and significant progress of the two individual sides on their own. Thus, this method may be important for promoting deep and actual practice-aimed remote sensing applications as remote sensing data becomes more readily available and as remote sensing technology becomes more advanced for use in pressing seismic loss risk reduction battlefields of earthquake active areas with rapid socioeconomic development.
By focusing on the seismic loss risks to buildings, the proposed integrated method started with employing either the new algorithm developed herein or the old but actual practice-favorable method to accurately and quickly extract the heights and footprint areas of a large number of residential and public office buildings from single high-resolution optical remote sensing images (Hr-ORSIs). Based on these two beginning parameters, this integrated method adopted careful and comprehensive building-relevant local knowledge (Br-LK) reviews as a bridge to help quickly and accurately determine or estimate technical cores or key data for performing specific mainstream quantitative seismic loss risk calculations to these buildings, including their floor areas, structural types, damage probability matrixes (DPMs), and economic parameters for calculating monetary losses due to their structural and indoor content damage. In doing so, the conventional building inventory preparation and vulnerability analysis processes were properly replaced, thus the timeliness, definitiveness, and accuracy were all considerably enhanced when estimating earthquake-induced building loss risks across large areas with rapid socioeconomic development.
The proposed integrated method was demonstrated in Tangshan, China, an area that has experienced the severe impacts of the Ms 7.8 earthquake in 1976 and significant expansion of the population and economy since that event. The total floor area of the residential and public office buildings in central Tangshan in 2009, which was estimated using this integrated method employing GeoEye-1 images with meter and sub-meter resolutions, was 3.99% lower than the corresponding number obtained by a conventional earthquake loss estimation project. Field-based verification indicated that the mean relative error of the method for estimating the floor areas of the addressed buildings was 2.99%. In addition to being capable of offering an accuracy that is comparable to that of conventional and mainstream quantitative earthquake loss estimate methods based on conventional building inventory preparation methods, this integrated method presents the following advantages relative to conventional methods: noteworthy speed, substantial labor savings, much lower cost, and significant serviceability when facing issues defined by time (non-disaster vs. impending disaster time, historical vs. current issue, status quo vs. situation changes). However, the timely, accurate, and all-round large-scale characterization of seismic loss risks and their changes in areas with rapid socioeconomic growth is a work that contains many challenges. The integrated method presented in this study that focused on seismic loss risks to buildings is only a beginning and a small portion of the entire approach. Even for estimating loss risks to buildings, several weaknesses exist in this integrated approach that challenge its more broad application and require to be explored in the future, such as the uncertainty when dealing with buildings in mountainous and hilly areas, the incapability of addressing specially structured buildings, the financial infeasibility of covering very large geographic areas, and the uncertainness involved with using remote sensing images with lower resolutions as alternatives.
As a particular application of the proposed integrated method, a simulation of the modern impacts of the 1976 Ms 7.8 Tangshan earthquake was performed. The results indicated that if this earthquake had recurred in 2009, the total damaged floor area of the residential and public office buildings and associated direct monetary losses in the study area could have been 8.00 and 28.73 times greater, respectively, than in 1976. This simulated scenario is a strong warning to the residents of crowded earthquake-prone areas regarding the increasing challenges that they may face and provides concrete and solid evidence of what may accompany a booming economy of an area with active seismic activities. In earthquake active areas with rapid socioeconomic development, these hidden/latent troubles must be dealt with properly and timely when considering economic growth.
This integrated method was developed based on the Chinese mainland and only focused on the seismic loss risks of buildings. However, the general ideas and overall guidelines that employ careful and comprehensive local knowledge reviews to bridge the large gaps between what remote sensing can directly provide and what mainstream quantitative seismic loss risk estimates specifically require may be universally applicable, either across boundaries or across different seismic impact-susceptible human objects. The general ideas and overall guidelines can be tested either outside of the Chinese mainland or on structures other than buildings. However, the specific risk analysis cores, the specific remote sensing applications, and the specific local knowledge defined by these cores and applications in these tests likely vary, even sharply, from one location or type of seismic-susceptible human object to another. Exactly what these diversities may be, particularly the place or object-specific local knowledge, what the specific gaps are between such risk analysis cores and the involved remote sensing applications, and in what ways and how much these gaps can be bridged by comprehensively and carefully reviewing relevant local knowledge, are really attractive and worthwhile careful explorations.

A2. Parameter Values
Based on the findings of this study, the values of all relevant parameters in the above equations are as follows.
(1) : Total floor area of the -type buildings. The specific numbers in 2009 are presented in the following table (Table A1). (2) | : DPMs of the buildings. For specific values, please see Table 4 in Section 4.3 of the main body of the paper.
(3) : Seismic intensity. As previously described in the appendix, we selected an intensity of VI to represent a moderate impact intensity and X to represent a heavy shock. The same intensity (VI or X) was maintained without any variations between the sites.
(4) : Mean monetary value of the building contents per unit floor area of the -type buildings. The specific values for the study area in 2009 are provided in Table A2.  Table A3. (6) ε Loss ratio of the building contents with the degree of building damage ( ). The specific values are presented in Table A4.  (7): ε : Loss ratio of the -type buildings with a damage degree of . The specific values are presented in Table A5.

A3.1. Building Damage Scenarios
If the study area suffered a seismic intensity of VI in 2009, only the S/B buildings would have sustained significant/substantial damage. The damage of many of these buildings could reach the Moderate degree of damage, and some could reach the Major degree of damage. However, the B/C and R/C buildings would only sustain a light impact (i.e., Slight degree of damage). The spatial distributions of the damages throughout the study area when subjected to the same impact intensity of VI indicated that the locations with larger building damage ratios would mainly be rural ( Figure A1a), where S/B buildings were dominant in 2009 (Table A1 and Figure 8a in the main body). However, if a seismic intensity of X affected the area in 2009, the residential and public office use buildings sustaining serious damage and destruction would have included all of the S/B, B/C and R/C buildings. In this scenario, few places in the study area would be free of heavy impacts ( Figure A1b).

A3.2. Direct Monetary Loss Scenarios
If a seismic intensity of VI was encountered in 2009, the direct monetary losses in the study area would mostly have resulted from the damage to the S/B buildings. In terms of spatial appearance, the places with larger loss ratios would be rural locations ( Figure A2a) where most of the buildings were S/B in 2009 (Table A1 and Figure 8a in the main body). However, if an intensity of X occurred, the direct monetary losses of the study area would have resulted from severe and heavy damages of all of the S/B, B/C and R/C buildings. In terms of spatial differences in this impact case, urban centers and some periurban regions would sustain the heaviest monetary losses within the study area ( Figure A2b). This significant feature is a combined product of the following factors: (1) Greater building densities occur and more good/high standard buildings exist in the urban and periurban locations (Table A1 and Figure 8a in the main body); (2) Good/high standard buildings require higher replacement and repair costs; and (3) The building contents in these locations are usually greater and more expensive relative to other locations in the study area.