Regularized sparse‐grid geometric sampling for uncertainty analysis in non‐linear inverse problems

This paper introduces an efficiency improvement to the sparse‐grid geometric sampling methodology for assessing uncertainty in non‐linear geophysical inverse problems. Traditional sparse‐grid geometric sampling works by sampling in a reduced‐dimension parameter space bounded by a feasible polytope, e.g., a generalization of a polygon to dimension above two. The feasible polytope is approximated by a hypercube. When the polytope is very irregular, the hypercube can be a poor approximation leading to computational inefficiency in sampling. We show how the polytope can be regularized using a rotation and scaling based on principal component analysis. This simple regularization helps to increase the efficiency of the sampling and by extension the computational complexity of the uncertainty solution. We demonstrate this on two synthetic 1D examples related to controlled‐source electromagnetic and amplitude versus offset inversion. The results show an improvement of about 50% in the performance of the proposed methodology when compared with the traditional one. However, as the amplitude versus offset example shows, the differences in the efficiency of the proposed methodology are very likely to be dependent on the shape and complexity of the original polytope. However, it is necessary to pursue further investigations on the regularization of the original polytope in order to fully understand when a simple regularization step based on rotation and scaling is enough.     

Denudation and retreat of the Serra do Mar escarpment in southern Brazil derived from in situ‐produced 10Be concentration in river sediment

The Serra do Mar escarpment, located along the southeastern coast of Brazil, is a high‐elevation passive margin escarpment. This escarpment evolved from the denudation of granites, migmatites and gneisses. The granites outcrop in the form of a ridge along the escarpment crest, due to its differential erosion (‘sugarloaf’ hills) from the surrounding lithologies. Several studies suggest that the passive margin escarpments are actively retreating toward the interior of the continent. However, no prior study has calculated the long‐term denudation rates of Serra do Mar to test this hypothesis. In this study, we measured the in situ‐produced ¹⁰Be concentration in fluvial sediments to quantify the catchment‐wide long‐term denudation rates of the Serra do Mar escarpment in southern Brazil. We sampled the fluvial sediments from ten watersheds that drain both sides of the escarpment. The average long‐term denudation rate of the oceanic side is between 2.1‐ and 2.6‐fold higher than the rate of the continental side: 26.04 ± 1.88 mm ka^‐1 (integrating over between 15.8 ka^‐1 and 46.6 ka^‐1) and 11.10 ± 0.37 mm ka^‐1 (integrating over between 52.9 ka^‐1 and 85.4 ka^‐1), respectively. These rates indicate that the coastal base level is controlling the escarpment retreat toward the continental high lands, which is consistent with observations made at other high‐elevation passive margins around the globe. The results also demonstrate the differential erosion along the Serra do Mar escarpment in southern Brazil during the Quaternary, where drainages over granites had lower average denudation rates in comparison with those over migmatites and gneisses. Moreover, the results demonstrate that the ocean‐facing catchments have been eroded more intensely than those facing the continent. The results also reveal that drainage over the granites decreases the average denudation rates of the ocean‐facing catchments and the ‘sugarloaf’ hills therefore are natural barriers that slowly retreat once they are exhumed. Copyright © 2013 John Wiley & Sons, Ltd.     

Forest Fire Alert System: a Geo Web GIS prioritization model considering land susceptibility and hotspots – a case study in the Carajás National Forest,Brazilian Amazon

To increase the monitoring potential of forest fires, an alert classification methodology using satellite-mapped hotspots has been established to help forest managers in the prioritization of which hotspot to be verified in the field, thus potentially improving the distribution of fire-fighting resources. A computer application was developed based on web-distributed geographical information technology whose main function is to interact automatically generated satellite hotspots and risk areas indicated in fire-susceptibility maps and classify them into five alert levels. The location of the hotspots is available continuously every 4 h, and a susceptibility map is produced daily through map algebra algorithm, which uses static (topography, vegetation and land use) and dynamic (weather) variables. Every process runs through automated geoprocessing routines. The methodology was tested during the dry period of 2007 in the Carajás National Forest, in the Brazilian Amazon, within an area of 400,000 ha. It is a critical area constantly threatened by fires caused by invasions and deforestation owing to intense agribusiness advances and mining activities in its surroundings. This situation results in observations of many hotspots inside the study area for the same day and almost the same time period, in places of extreme opposites, demanding complex rapid analysis and hindering the decision of the displacement of fire-fighting teams. Further, a major mining company operates within the National Forest area, maintaining actions of protection as part of its environmental mining license. Results are presented under three aspects: (i) the credibility of the daily susceptibility map (algorithm), which showed strong correlation between areas of greatest risks and the confirmed forest fires; (ii) the reliability of hotspots (alert levels), confirming 71% of fires; (iii) accuracy in the decision of which hotspot to be checked, which revealed the same number of verifications at different alert levels, 82% confirmed alert 5 hotspots (maximum) and only 50% from alert 1 (minimum), resulting in faster fire-fighting actions, minimizing burned areas and, in some cases, allowing fire control before its spreading. Therefore, the methodology demonstrated that GIS routines are able to determine the relationship between a reality-based, interpreted susceptibility map of the area and satellite-generated hotspots, highlighting the ones of highest hazard level through the alert classification, becoming an important tool to help decisions from the fire-control center, especially for high-risk regions. The methodology may be extrapolated to other forested areas.     

Experimental investigation on the seismic performance of masonry buildings using shaking table testing

Masonry buildings worldwide exhibited severe damage and collapse in recent strong earthquake events. It is known that their brittle behavior, which is mainly due to the combination of low tensile strength, large mass and insufficient connection between structural elements, is the main limitation for their structural implementation in residential buildings. A new construction system for masonry buildings using concrete blocks units and trussed reinforcement is presented here and its seismic behavior is validated through shaking table tests. Dynamic tests of two geometrically identical two-story reduced scale (1:2) models have been carried out, considering artificial accelerograms compatible with the elastic response spectrum defined by the Eurocode 8. The first model was reinforced with the new proposed system while the second model was built with unreinforced masonry. The experimental analysis encompasses local and global parameters such as cracking patterns, failure mechanisms, and in-plane and out-of-plane behavior in terms of displacements and lateral drifts from where the global dynamic behavior of the two buildings is analyzed comparatively. Finally, behavior factors for the design recommendations in case of unreinforced masonry are also evaluated.     

Numerical relationships between magnetic parameters measured in Quaternary sediments and global paleoclimatic proxies

The complexity of most geological and geophysical problems prompts sometimes the use of non linear mathematical methods to handle them. An adaptive neuro fuzzy inference system (ANFIS) that combines fuzzy logic with neural networks, is applied here to study a paleoclimate section from the Quaternary sedimentary fill of the Lake Mucubají (western Venezuela). The purpose of this work is to find a set of numerical relationships that could predict the possible connections between oxygen isotope (δ¹⁸O) values from two different locations in the northern hemisphere (Ammersee in southern Germany and an ice core from the Greenland Ice Core Project — GRIP) and rock-magnetic parameters measured in Mucubají samples (i.e. mass-specific magnetic susceptibility — χ, magnetic remanence S-ratio, mass-specific saturation isothermal remanent magnetization — SIRM and anhysteretic remanent magnetization — ARM). The best inferences in terms of coefficient of determionation R² and the Root Mean-Square Error (RMSE) are obtained using those magnetic data as input that include information about magnetite grain size distributions, e.g., SIRM and ARM in FIS structures [1χ, 4ARM] and [4ARM, 1SIRM]. A comparison between Ammersee and GRIP actual data, as well as their corresponding inferences for the FIS structure [4ARM, 1SIRM], reveals a reasonable good inference of global trends for both records, overlooking the regional and/or local paleoclimate forcings that might have affected Ammersee. A better correlation between global isotope paleoclimate records and magnetic proxies, is perhaps prevented by the role played by local and regional paleoclimate and tectonism in Mucubají. We also argue that the ratio of ARM over SIRM appears to be related in a complex way to the onset and to the end of the Younger Dryas. Our novel approach to the assessment of a specific paleoclimate case study shows the potential of the ANFIS technique in solving problems where traditional univariate and multivariate linear regression methods could prove inadequate.     

Disentangling the role of the Pacific and Atlantic Oceans during the Amazonian droughts in 2015

Theoretical and Applied Climatology - El Niño events lead to worldwide changes in the inter-annual variability of precipitation, and more in particular to South America droughts. In the year...     

Rockfall risk assessment to persons travelling in vehicles along a road: the case study of the Amalfi coastal road (southern Italy)

The paper deals with the assessment of rockfall risk to persons travelling in vehicles along the SS163 road, an important transportation corridor supporting a high vehicle traffic within the well-known tourist area of the Amalfi Coast (southern Italy). To this aim, the Rockfall Hazard Rating System (RHRS) and quantitative risk assessment (QRA) procedures, in this latter case for three rockfall risk scenarios, are applied. With reference to a large portion (33.820 out of a total of 50.365 km) of the SS163 road, the obtained QRA results highlight that, although the estimated individual risk to life satisfies the adopted tolerable risk criterion, the computed societal risk cannot be tolerated. Starting from this result, site-specific QRA analyses—carried out with reference to some road sections chosen on the basis of the RHRS results—allow the detection of the SS163 portions where the individual risk to life exceeds the tolerable risk threshold and, then, the recourse to mitigation measures could reveal necessary. In this regard, RHRS and QRA methods can be considered complementary tools in prioritizing the road sections where construction funds can be profitably spent in order to mitigate the rockfall risk with reference to both direct consequences (life loss) and indirect ones (traffic delay and diversions).     

Is the cutting of oil contaminated marshes an efficient clean-up technique in a subtropical estuary?

Cutting and removal of oil-impacted marsh plants are still used worldwide as a clean-up and recovery technique. To experimentally test the efficacy of cutting and removing marsh plants under subtropical conditions, we simulated an oil spill (Bunker MF-180) in Spartina alterniflora marshes and compared the responses of plant height, biomass, density of culms and number of flowering plants in high and low energy areas in Paranaguá Bay (S Brazil) for about 9 months. Cutting and removal were inefficient in promoting or accelerating the recovery of the impacted areas. Cut or uncut impacted marshes fully recovered within 6 months, both in low and high energy areas. Plant cutting should be practiced only when there is an effective risk of contamination of groundwater near urban areas, when obvious aesthetical issues are involved in areas of touristic interest or when there are real short-term conservation risks to threatened species.     

Geology and geochronology of the Betara region in south-southeastern Brazil: Evidence for possible Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extension events

The results of geological mapping, chemical analysis and radiometric dating of metabasic rocks of Betara Formation, and mapping and dating of those present in the Betara basement nucleus together with mylonitic granodiorite and syenogranite are reported here. U–Pb analysis of bulk zircon fractions from the metabasic rocks of the basement nucleus yielded a Statherian age of 1790 ± 22 Ma, while the metabasic rocks from the upper part of the Betara Formation yielded a Calymmian age between 1500 and 1450 Ma. This age is a minimum for the deposition of the Betara Formation. The older metabasic rocks are associated with post-tectonic, possibly anorogenic syenogranite, while the younger ones are gabbro or very porphyritic ankaramite whose REE patterns are consistent with crystallization from an N-MORB parent magma. The observations and data point to the probable events associated with extensional processes of the end of Paleoproterozoic and early Mesoproterozoic. Similar registers of Statherian (1.80–1.75 Ga) and Calymmian (1.50–1.45 Ga) extensional events are recorded in other parts of the South American and African continents. The Neoproterozoic witnessed the formation and junction of the tectonic slices which formed the Apiaí domain during the assemblage of western Gondwana.     

Collision of the Ganges–Brahmaputra Delta with the Burma Arc: Implications for earthquake hazard

We take a fresh look at the topography, structure and seismicity of the Ganges–Brahmaputra Delta (GBD)–Burma Arc collision zone in order to reevaluate the nature of the accretionary prism and its seismic potential. The GBD, the world's largest delta, has been built from sediments eroded from the Himalayan collision. These sediments prograded the continental margin of the Indian subcontinent by  400 km, forming a huge sediment pile that is now entering the Burma Arc subduction zone. Subduction of oceanic lithosphere with > 20 km sediment thickness is fueling the growth of an active accretionary prism exposed on land. The prism starts at an apex south of the GBD shelf edge at  18°N and widens northwards to form a broad triangle that may be up to 300 km wide at its northern limit. The front of the prism is blind, buried by the GBD sediments. Thus, the deformation front extends 100 km west of the surface fold belt beneath the Comilla Tract, which is uplifted by 3–4 m relative to the delta. This accretionary prism has the lowest surface slope of any active subduction zone. The gradient of the prism is only  0.1°, rising to  0.5° in the forearc region to the east. This low slope is consistent with the high level of overpressure found in the subsurface, and indicates a very weak detachment. Since its onset, the collision of the GBD and Burma Arc has expanded westward at  2 cm/yr, and propagated southwards at  5 cm/yr. Seismic hazard in the GBD is largely unknown. Intermediate-size earthquakes are associated with surface ruptures and fold growth in the external part of the prism. However, the possibility of large subduction ruptures has not been accounted for, and may be higher than generally believed. Although sediment-clogged systems are thought to not be able to sustain the stresses and strain-weakening behavior required for great earthquakes, some of the largest known earthquakes have occurred in heavily-sedimented subduction zones. A large earthquake in 1762 ruptured  250 km of the southern part of the GBD, suggesting large earthquakes are possible there. A large, but poorly documented earthquake in 1548 damaged population centers at the northern and southern ends of the onshore prism, and is the only known candidate for a rupture of the plate boundary along the subaerial part of the GBD–Burma Arc collision zone.