A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.     

Catastrophic,rainfall-induced debris flows in Andean villages of Tarapacá, Atacama Desert,northern Chile

In March 2012, during the rainy season in the Altiplano plateau, a >100-year return period rainfall event affected the deeply incised valleys of the Precordillera of the Tarapacá Region, northern Chile. This extreme event in a very arid region triggered a number of debris and mud flows that caused severe damage and destruction in several small villages along the Camiña and Tarapacá valleys. The highly vulnerable location of the villages on top of alluvial fans due to socioeconomic and cultural reasons is a key factor to explain the level of destruction in most villages. In this paper, this unusual, remarkable landslide event is described, and the hazard faced by these settlements for future rainfall episodes and possible mitigation measures are discussed.     

Late Quaternary vegetation and climate dynamics in central-eastern Brazil: insights from a ~35k cal a bp peat record in the Cerrado biome

The late Quaternary evolution of central-eastern Brazil has been under-researched. Questions remain as to the origin of the Cerrado, a highly endangered biome, and other types of vegetation, such as the Capões – small vegetation islands of semi-deciduous and mountain forests. We investigated the factors that influenced the expansion and contraction of the Cerrado and Capões during the late Quaternary (last ~35 ka), using a multi-proxy approach: stable isotopes (δ¹³C, δ¹⁵N), geochemistry, pollen and multivariate statistics derived from a peat core (Pinheiro mire, Serra do Espinhaço Meridional). Five major shifts in precipitation, temperature, vegetation and landscape stability occurred at different timescales. Our study revealed that changes in the South Atlantic Convergence Zone (SACZ) seem to have been coeval with these shifts: from the Late Glacial Maximum to mid-Holocene the SACZ remained near (~29.6 to ~16.5k cal a bp ) and over (~16.5 to ~6.1 k cal a bp ) the study area, providing humidity to the region. This challenges previous research which suggested that climate was drier for this time period. At present, the Capões are likely to be a remnant of a more humid climate; meanwhile, the Cerrado biome seems to have stablished in the late Holocene, after ~3.1 k cal a bp .     

Retrowedge-related Carboniferous units and coeval magmatism in the northwestern Neuquén province, Argentina

The studied Carboniferous units comprise metasedimentary (Guaraco Norte Formation), pyroclastic (Arroyo del Torreón Formation), and sedimentary (Huaraco Formation) rocks that crop out in the northwestern Neuquén province, Argentina. They form part of the basement of the Neuquén Basin and are mostly coeval with the Late Paleozoic accretionary prism complex of the Coastal Cordillera, south-central Chile. U–Pb SHRIMP dating of detrital zircon yielded a maximum depositional age of 374?Ma (Upper Devonian) for the Guaraco Norte Formation and 389?Ma for the Arroyo del Torreón Formation. Detrital magmatic zircon from the Guaraco Norte Formation are grouped into two main populations of Devonian and Ordovician (Famatinian) ages. In the Arroyo del Torreón Formation, zircon populations are also of Devonian and Ordovician (Famatinian), as well as of Late Neoproterozoic and Mesoproterozoic ages. In both units, there is a conspicuous population of Devonian magmatic zircon grains (from 406?±?4?Ma to 369?±?5?Ma), indicative of active magmatism at that time range. The ε_Hf values of this population range between ?2.84 and ?0.7, and the TDM-(Hf) are mostly Mesoproterozoic, suggesting that the primary sources of the Devonian magmatism contained small amounts of Mesoproterozoic recycled crustal components. The chemical composition of the Guaraco Norte Formation corresponds to recycled, mature polycyclic sediment of mature continental provenance, pointing to a passive margin with minor inputs from continental margin magmatic rocks. The chemical signature of the Huaraco Formation indicates that a magmatic arc was the main provenance for sediments of this unit, which is consistent with the occurrence of tuff—mostly in the Arroyo del Torreón Formation and very scarcely in the Huaraco Formation—with a volcanic-arc signature, jointly indicating the occurrence of a Carboniferous active arc magmatism during the deposition of the two units. The Guaraco Norte Formation is interpreted to represent passive margin deposits of mostly Lower Carboniferous age (younger than 374?Ma and older than 326?Ma) that precede the onset of the accretionary prism in Chile and extend into the earliest stage of the accretion, in a retrowedge position. The Arroyo del Torreón and Huaraco formations are considered to be retrowedge basin deposits to the early frontal accretionary prism (Eastern Series) of Chile. The presence of volcanism with arc signature in the units provides evidence of a Mississippian magmatic arc that can be correlated with limited exposures of the same age in the Frontal Cordillera (Argentina). The arc would have migrated to the West (Coastal Batholith) during Pennsylvanian–Permian times (coevally with the later basal accretionary prism/Western Series). The source of a conspicuous population of Devonian detrital zircon interpreted to be of magmatic origin in the studied units is discussed in various possible geotectonic scenarios, the preferred model being a magmatic arc developed in the Chilenia block, related to a west-dipping subduction beneath Chilenia before and shortly after its collision against Cuyania/Gondwana, at around 390?Ma and not linked to the independent, Devonian–Mississippian arc, developed to the south, in Patagonia.     

Relationships Between Community Structure of the Intertidal Macroinfauna and Sandy Beach Characteristics Along the Chilean Coast

Abstract. Eight sandy beaches were seasonally sampled along the coast of Chile, from ca. 21 to 42° S (about 3000 km) to study the relationship between community structure of the intertidal macroinfauna and beach characteristics. Sediment samples (0.1 m², 30 cm deep) were collected (July – September 1998 and December 1998 – January 1999) with plastic cylinders at 15 equally spaced levels along three replicated transects extending from above the drift line to the swash zone. The sediment was sieved through a 1 mm mesh and the organisms collected stored in 5 % formalin. To define beach types, Dean's parameter (Ω) was calculated from wave heights and periods, and fall velocity of sand particles from the swash zone. Crustaceans (mainly peracarids) were the most diverse group with 14 species, followed by polychaetes with 5 species. The talitrid amphipod Orchestoidea tuberculata , the cirolanid isopods Excirolana braziliensis and E. hirsuticauda and the anomuran decapod Emerita analoga were the most widely distributed and common species. Regression analyses between species richness, abundance and biomass of the whole macroinfauna versus sediment characteristics, beach face slopes and morphodynamic beach states showed no significant relationships. Thus, macroinfaunal community characteristics did not increase linearly from lower intermediate to higher intermediate or dissipative beach states as had been found before in Chile or in other coasts. A comparative analysis with data from sandy beaches of other world regions showed that the number of species inhabiting Chilean sandy beaches was generally lower, whereas total population abundances were generally higher compared with values reported elsewhere.     

Inelastic displacement ratios for evaluation of existing structures

Results of a detailed statistical study of constant relative strength inelastic displacement ratios to estimate maximum lateral inelastic displacement demands on existing structures from maximum lateral elastic displacement demands are presented. These ratios were computed for single‐degree‐of‐freedom systems with different levels of lateral strength normalized to the strength required to remain elastic when subjected to a relatively large ensemble of recorded earthquake ground motions. Three groups of soil conditions with shear wave velocities higher than 180m/s are considered. The influence of period of vibration, level of lateral yielding strength, site conditions, earthquake magnitude, distance to the source, and strain‐hardening ratio are evaluated and discussed. Mean inelastic displacement ratios and those associated with various percentiles are presented. A special emphasis is given to the dispersion of these ratios. It is concluded that distance to the source has a negligible influence on constant relative strength inelastic displacement ratios. However, for periods smaller than 1s earthquake magnitude and soil conditions have a moderate influence on these ratios. Strain hardening decreases maximum inelastic displacement at a fairly constant rate depending on the level of relative strength for periods of vibration longer than about 1.0s while it decreases maximum inelastic displacement non‐linearly as the period of vibration shortens and as the relative‐strength ratio increases for periods of vibration shorter than 1.0s. Finally, results from non‐linear regression analyses are presented that provide a simplified expression to be used to approximate mean inelastic displacement ratios during the evaluation of existing structures built on firm sites. Copyright © 2003 John Wiley & Sons, Ltd.     

Exercise in downscaling on sea surface temperature along Chinese coast

ＩＮＴＲＯＤＵＣＴＩＯＮＡｎｔａｒｃｔｉｃｋｒｉｌｌ(ＥｕｐｈａｕｓｉａｓｕｐｅｒｂａＤａｎａ)ｉｓａｋｅｙｏｒｇａｎｉｓｍｏｆｔｈｅＡｎｔａｒｃｔｉｃｍａｒｉｎｅｅｃｏｓｙｓｔｅｍａｎｄａｐｏｔｅｎｔｉａｌｆｉｓｈｅｒｙｒｅｓｏｕｒｃｅ.Ｉｔｈａｓｂｅｅｎｔｈｅｓｕｂｊｅｃｔｏｆｍｕｃｈｒｅｓｅａｒｃｈ.Ｄｅｓｐｉｔｅｔｈｅｌｏｎｇｈｉｓｔｏｒｙｏｆｒｅｓｅａｒｃｈｄａｔｉｎｇｂａｃｋｔｏｔｈｅ1930’ｓ,ｔｈｅｒｅａｒｅｓｔｉｌｌｑｕｅｓｔｉｏｎｓａｂｏｕｔｉｔｓｐｏｐｕｌａｔｉｏｎｓｔｒｕｃｔｕ…     

Modelling considerations in probabilistic performance‐based seismic evaluation: case study of the I‐880 viaduct

Limitations associated with deterministic methods to quantify demands and develop rational acceptance criteria have led to the emergence of probabilistic procedures in performance‐based seismic engineering. The Pacific Earthquake Engineering Research performance‐based methodology is one such approach. In this paper, the impact of certain modelling decisions made at different stages of the evaluation process on the performance assessment of a typical multi‐bent viaduct is examined. Modelling, in the context of this paper, covers hazard modelling, structural modelling and loss modelling. The specific application considered in this study is a section of an existing viaduct in California: the I‐880 interstate highway. Several simulation models of the viaduct are developed, a series of nonlinear time‐history analyses are carried out to predict demands, measures of damage are evaluated and the probability of closure of the viaduct is estimated using the specified hazard for the site. It is concluded that the methodology offers several advantages over existing deterministic performance‐based procedures. Results of the investigation indicate that the assessment methodology is particularly sensitive to the reliability of decisions made by bridge inspectors following a seismic event, and to the dispersion in the demand estimation, which in turn is influenced by several factors including soil–structure interaction effects and ground motion scaling procedures. Copyright © 2005 John Wiley & Sons, Ltd.     

A framework for creating and validating a non-linear spectrum-biomass model to estimate the secondary succession biomass in moist tropical forests

Uncertainties remain in the use of remote sensing technologies to provide validated model-derived estimates of the biomass of the secondary succession (SS) forests in the Amazon Basin. The objectives of this study were to develop a modeling framework for creating a valid spectrum-biomass model to estimate the SS biomass, to assess the utility of the framework and the accuracy and validity of the model, and to identify the model’s determinants. Data sources for this study include 1992–1993 vegetation inventory data and 1991 Landsat Thematic Mapper (TM) data on the Altamira region of Para, Brazil, and 1994–1995 vegetation inventory data and 1994 Landsat TM data on the nearby Bragantina region. The allometric approach was used to estimate the biomass of the sampled sites based on the vegetation inventory data. A framework for the spectrum-biomass regression model was developed based on the estimated biomass of the sampled sites and the Landsat data. The framework includes (1) the pooling of data from Bragantina and the use of ANCOVA to justify this approach; (2) image calibration; (3) biomass data age-adjustment, (4) selection of independent variables, (5) regression model development, and (6) model assessment and validation. The cubic regression model with TM Band5-related predictors was found to best fit the data as evidenced by an adjusted R-squared value of 0.865, mean square error (MSE) of the model, and the analysis of residuals. Residual analysis showed that the model might yield a better estimation on a lower biomass values than on higher biomass values. In addition, further analyses identified several determinants that can impact the accuracy of the spectrum-biomass model. ANCOVA confirmed that the relationship between the biomass and the spectrum is independent of the Altamira and Bragantina regions, and that it was appropriate to pool sampled data from both regions in the proposed model. The model development methodology and the model produced from this research will be of value to researchers using the spectrum-biomass modeling approaches to estimate the biomass and study the SS rates in moist tropical forests.