The Huashan piedmont fault, forming a part of the southern margin of the Weihe graben, is one of the important normal faults that control the subsidence of the intracontinental rift. Developing on the footwall of the fault, the Huashan block has experienced rapid cooling during the Cenozoic, especially since the early-middle Miocene. Mountain exhumation causes and transports a great amount of sediments to the adjacent hanging wall, setting a typical case of mountain-basin coupling system. Studies on active tectonics, historical and paleo earthquakes and field investigations reveal that the middle section(Huaxian-Huayin)of the fault is much more active than the west(Lantian-Huaxian)and east(Huayin-Lingbao)sections.
    We extracted channel profiles of rivers that originate from the main water divide of the northern flank of the Huashan Mountain. Based on the method of slope-area analysis and the integral approach, we identified knickpoints, calculated channel concavity and steepness indices, and constructed paleo river profiles. Of most rivers, the concavities are within a relatively narrow range of 0.3~0.6, with no obvious correlation with tectonics. However, channel steepness and knickpoint distribution vary spatially. In the east section, rivers are under steady-state with smooth, concave-up channels and lower steepness((104±30)m^0.9). In the other two sections, rivers are mainly under transient state with slope-break knickpoints. For the channel segments below knickpoints, steepness indices are much higher in the middle section((230±92)m^0.9)than in the west((152±53)m^0.9). Thus, the variance of fault activity can be reflected by channel steepness pattern. Above the knickpoints, channel steepness indices are much lower(middle(103±23)m^0.9, west(60±14)m^0.9). What's more, we found a statistically significant power-law scaling between knickpoint retreat distance and catchment drainage area. Thus, we attributed these knickpoints to be the results of recent rapid uplift of the Huashan block. The relief of paleo channels(middle(1000±153)m, west(751±170)m)accounts for~60%~80% of the relief of modern rivers(middle(1323±249)m, west(1057±231)m), which means that ~20%~40% of modern channel relief was caused by the episode of the rapid uplift. Assuming a balance between the rates of rock uplift and downstream river incision, a power-law function between uplift rates and channel steepness can be derived. According to the fault throw rates of the middle section 1.5~3mm/a(since late Pleistocene), we constrained slope exponent n~0.5 and channel erodibility K~1.5×10^-4m^0.55/a. Combining the knickpoint age formula, we estimated that the rapid mountain uplift/fault throw began at ~(0.55±0.25)Ma BP. Therefore, the middle of the Huashan piedmont fault is more active than the west and east sections. The fast fault throw of the west and middle sections since the middle Pleistocene has caused rapid mountain uplift and high topographic relief.

LA-ICP-MS(laser ablation-inductively coupled-mass spectrometry)has been recently used for rapid, accurate and precise U-Pb geochronology on zircon grains. In this paper, we adopted an Agilent 7900 quadrupole ICP-MS coupled with a Resolution M50-LR 193nm excimer laser system to establish integrated measurement procedures. Before analysis, the system is tuned to achieve sensitivities better than 30000 cps/s for ²³⁸ U with a 40μm spot size, at~3.5J/cm². Detailed parameters for laser system and ICP-MS are presented here. Then, we analyzed five reference zircons(91500, GJ-1, Plesovice, FCT, Penglai)with a wide range in age from~1064 to~4.4Ma. Two standard zircons, 91500 and GJ-1, are employed as external reference standards. Generally, second zircon standard is analyzed in an effort to ensure accuracy and evaluate reproducibility. A typical analysis sequence includes one international glass standard(NIST610), two external reference standards, five grains of unknown zircon with every eight ablations. Laser induced time-dependent elemental fractionation is corrected using the intercept method, whereas instrument drift, mass bias and elemental fractional caused by ionization differences are corrected by external reference standard 91500 or GJ-1. Compared with 91500 and GJ-1, common Pb content of Plesovice, FCT, Penglai can't be ignored. Thus, we did common Pb correction for the above three standard zircons. The performance of the established procedure was assessed by analyzing zircon range in age from~1 064 to~4Ma. The results show that the ages of these five references are consistent with the ages of published studies with accuracy for three international references(91500, GJ-1, Plesovice)better than 3% and two young secondary references(FCT, Penglai)lower than 7% at the 2 sigma level, which indicates that our analytical procedure is reliable. For individual laser analysis, the uncertainties are mainly from three sources:Measurement error of isotope ratio, error of correction factors for instrument drift and element fractionation, and error of recommended age of external references. Compared to three international references, there are three extra uncertainties for young reference zircons, including:1)little radioactive isotopes closing to blank level increase the measurement error of isotope ratio; 2) effect of common lead becomes more significant;3) the nonhomogeneous samples couldn't match references well. Therefore, accuracy and precision of measurement depend on absolute age, content of common lead and matching degree between references and samples. In summary, the accuracy and precision obtained using the technique presented in this study are similar to those of other LA-ICP-MS laboratories.

The topography and geomorphology of active orogens result from the interaction of tectonics and climate. In most orogens, a fluvial channel is most sensitive to the coupling between tectonics, lithology, and climate. Meanwhile, the related signals have been recorded by both the drainage geometry and channel longitudinal profile. Thus, how to extract tectonic information from fluvial channels has been a focused issue in geologic and geomorphologic studies.
The well known stream-power river incision model bridges the gap between tectonic uplift, river incision and channel profile change, making it possible to retrieve rock uplift pattern from river profiles. In this model, the river incision rate depends on the rock erodibility, contributing drainage area and river gradient. The steady-state form of the river incision model predicts a power-law scaling between the drainage area and channel gradient. Via a linear regression to the log-transformed slope-area data, the slope and intercept are channel concavity and steepness indices, respectively. The concavity relates to lithology, climatic setting and incision process while the channel steepness can be used to map the spatial pattern of rock uplift. For its simple calculation process, the slope-area analysis has been widely used in the study of tectonic geomorphology during past decades.
However, to calculate river slope, the coarse channel elevation data must be smoothed, re-sampled, and differentiated without any reasonable smooth window or rigid mathematical fundamentals. One may lose important information and derive stream-power parameters with high uncertainties. In this paper, we introduce the integral approach, a procedure that has been widely used in the latest four years and demonstrated to be a better method for river profile analysis than the traditional slope-area analysis. Via the integration to the steady-state form of the stream-power river incision equation, the river longitudinal profile can be converted into a straight line of which the independent variable is the integral quantity χ with the unit of distance and the dependent variable is the relative channel elevation. We can calculate the linear correlation coefficient between elevation and χ based on a series of concavity values and find the best linear fit to be the reasonable channel concavity index. The slope of the linear fit to the χ value and elevation is simply related to the ratio of the uplift rate to the erodibility.
Without calculating channel slope, the integral approach makes up for the drawback of the slope-area analysis. Meanwhile, via the integral approach, a steady-state river profile can be expressed as a continuous function, which can provide theoretical principle for some geomorphic parameters (e.g., slope-length index, hypsometric integral). In addition, we can determine the drainage network migration direction using this method. Therefore, the integral approach can be used as a better method for tectonogeomorphic research.

The M_S 8.0 Wenchuan earthquake and its 2216 aftershocks were relocated using the double difference algorithm.The horizontal and vertical errors of the 2061 relocated hypocenters are approximately 1～2km and 2～3km,respectively.The epicenter of main shock is approximately 31.00°N,103.38°E,the focal depth is about 13km and the seismogenic structure is the central fault of Longmenshan Fault zone.The total length of spatial distribution of aftershocks along the strike of the fault is about 330km and the predominance distribution of focal depth is 3～20km,which shows obviously the characteristic of segmented activity.The seismicity of the southern part mainly concentrates on the central fault of Longmenshan Fault zone,and some earthquakes occurred on the range-front and range-back faults;the dip of the three faults seems to become gentler gradually from west to east,forming imbricate ruptures.The central fault and Pingwu-Qingchuan Fault of the Longmenshan Fault zone are involved in the seismogenic process,and the seismic rupture is both of thrust napping and right-lateral strike-slipping.