Uncertainty analysis and probabilistic segmentation of electrical resistivity images: the 2D inverse problem

In this paper, we present the uncertainty analysis of the 2D electrical tomography inverse problem using model reduction and performing the sampling via an explorative member of the Particle Swarm Optimization family, called the Regressive‐Regressive Particle Swarm Optimization. The procedure begins with a local inversion to find a good resistivity model located in the nonlinear equivalence region of the set of plausible solutions. The dimension of this geophysical model is then reduced using spectral decomposition, and the uncertainty space is explored via Particle Swarm Optimization. Using this approach, we show that it is possible to sample the uncertainty space of the electrical tomography inverse problem. We illustrate this methodology with the application to a synthetic and a real dataset coming from a karstic geological set‐up. By computing the uncertainty of the inverse solution, it is possible to perform the segmentation of the resistivity images issued from inversion. This segmentation is based on the set of equivalent models that have been sampled, and makes it possible to answer geophysical questions in a probabilistic way, performing risk analysis.     

Objectives of fisheries management: case studies from the UK, France, Spain and Denmark

The main objectives of fisheries management are generally similar throughout the world. These are often stated in policy documents such as the Common Fisheries Policy and the Magnuson–Stevens Fishery Conservation and Management Act. However, at the local level often the key objectives of management are more detailed, characterised by both the overriding management structure and the status and type of fishery concerned. In this paper, we consider case study fisheries from the UK, France, Spain and Denmark to compare some of the various types of fisheries and fisheries management systems that exist in the European Union. From this, we define the key objectives for each management system.     

Numerical estimation of carbonate rock properties using multiscale images

Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples.     

Current-influenced depositional provinces, continental margin off Cape Hatteras, identified by petrologic method

Analyses of DSRV “Alvin” core samples on the Cape Hatteras margin indicate major textural and compositional changes at depths of about 1000 and well below 2500 m. The distribution patterns of petrologic parameters correlate well with water mass flow and suspended-sediment plumes measured on this margin by other workers. Our study also shows: (a) vigorous erosion and sediment transport at depths of less than 400 m resulting from the NE-trending Gulf Stream flow; (b) deposition, largely planktonic-rich sediment released from the Gulf Stream, on the upper- to mid-slope, to depths of about 800–1200 m; (c) winnowing, resuspension and deposition induced by periodically intensified slope currents on the mid-slope to uppermost rise, between about 1000 and 2500 m; and (d) prevailing deposition on the upper rise proper (below 2500 m), from transport by the SW-trending Western Boundary Undercurrent. Sediments moved by bottom currents have altered the composition and distribution patterns of material transported downslope by offshelf spillover; this mixing of gravity-emplaced and bottom-current-transported sediment obscures depositional boundaries. Moreover, reworking of the seafloor by benthic organisms alters physical properties and changes erodability of surficial sediments by bottom currents. Measurement of current flow above the seafloor and direct observation of the bottom are insufficient to delineate surficial sediment boundaries. Detailed petrologic analyses are needed to recognize the long-term signature of processes and define depositional provinces.     

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

Semi‐active tuned mass damper building systems: Application

Seismic performance attributes of multi‐story passive and semi‐active tuned mass damper (PTMD and SATMD) building systems are investigated for 12‐story moment resisting frames modeled as ‘10+2’ stories and ‘8+4’ stories. Segmented upper portion of the stories are isolated as a tuned mass, and a passive viscous damper or semi‐active resetable device is adopted as energy dissipation strategy. The semi‐active approach uses feedback control to alter or manipulate the reaction forces, effectively re‐tuning the system depending on the structural response. Optimum tuned mass damper control parameters and appropriate matching SATMD configurations are adopted from a companion study on a simplified two‐degree‐of‐freedom system. Statistical performance metrics are presented for 30 probabilistically scaled earthquake records from the SAC project. Time history analyses are used to compute response reduction factors across a wide range of seismic hazard intensities. Results show that large SATMD systems can effectively manage seismic response for multi‐degree‐of freedom systems across a broad range of ground motions in comparison to passive solutions. Specific results include the identification of differences in the mechanisms by which SATMD and PTMD systems remove energy, based on the differences in the devices used. Additionally, variability is seen to be tighter for the SATMD systems across the suites of ground motions used, indicating a more robust control system. While the overall efficacy of the concept is shown the major issues, such as isolation layer displacement, are discussed in detail not available in simplified spectral analyses, providing further insight into the dynamics of these issues for these systems. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydrological response of an afforested catchment in a Mediterranean humid mountain area: a comparative study with a natural forest

The effects of afforestation on water resources are still controversial. The aim of this paper is to (i) analyse the hydrological response of an afforested area in the Central Pyrenees and (ii) compare the hydrological response of an afforested area with the response observed in a natural undisturbed forest. The Araguás catchment was cultivated until the 1950s, and then afforested with pines in an effort to control the active degradation processes. The hydrological response was variable and complex, because the discharge was generated by a combination of distinct runoff processes. The hydrological response showed that (i) afforestation produced moderate peak discharges, stormflows and recession limbs, and long rising limbs; (ii) no one single variable was able to explain the hydrological response: rainfall volume and intensity did not explain the hydrological response and antecedent rainfall and initial discharge (indicating antecedent moisture conditions) did play an important role; (iii) seasonal differences were observed suggesting different runoff generating processes; and (iv) the effect of forest cover on peak discharges became less important as the size of the hydrological event increased. The stormflow coefficient showed a clear seasonal pattern with an alternation between a wet period, when the catchment was hydrologically responsive, and a dry summer period when the catchment rarely responded to any event, and two transitional periods (wetting and drying). Compared with a natural forest, the afforested area recorded greater flows and peak discharges, faster response times and shorter recession limbs. Afforestation reduces the water yield and the number of floods compared with non‐vegetated areas and abandoned lands. Copyright © 2016 John Wiley & Sons, Ltd.     

Flood Damage Reduction Program (FDRP) in Québec: Case Study of the Chaudiére River

This study presents the evolution of agreements between the governments of Canadaand Québec on flood damage reduction. In Québec, the implementation of a regulation about building in floodplains came about in 1983–1984. Today, this regulation takes the form of a policy called ``Policy of shores, littoral and floodplain protection'. Municipalities must adopt rules that concur with the principles of this policy.The Chaudiére River basin was selected for analysis of urban developmentduring the period following the application of building rules in flood-risk areas.Despite the ban on building in the strong current zone (0–20 year return periodflood zone), many buildings, essentially residential, have been erected in thiszone. These new constructions generally account for a low percentage of thetotal property value in the 0–100 year flooding area, but are legal since theyare connected to a water and sewage network that existed prior to the officialfloodplain regulation.Flood damage along the Chaudiére River will tend to increase for two reasons.Firstly, while respecting the policy mentioned above, sites are still available infloodplains for future development. Secondly, no structural flood protection workshave been erected in view of the fact they are only marginally profitable from anbenefit-cost point of view.     

Authigenic carbonates related to active seepage of methane-rich hot brines at the Cheops mud volcano,Menes caldera (Nile deep-sea fan,eastern Mediterranean Sea)

On the passive margin of the Nile deep-sea fan, the active Cheops mud volcano (MV; ca. 1,500 m diameter, ~20–30 m above seafloor, 3,010–3,020 m water depth) comprises a crater lake with hot (up to ca. 42 °C) methane-rich muddy brines in places overflowing down the MV flanks. During the Medeco2 cruise in fall 2007, ROV dives enabled detailed sampling of the brine fluid, bottom lake sediments at ca. 450 m lake depth, sub-surface sediments from the MV flanks, and carbonate crusts at the MV foot. Based on mineralogical, elemental and stable isotope analyses, this study aims at exploring the origin of the brine fluid and the key biogeochemical processes controlling the formation of these deep-sea authigenic carbonates. In addition to their patchy occurrence in crusts outcropping at the seafloor, authigenic carbonates occur as small concretions disseminated within sub-seafloor sediments, as well as in the bottom sediments and muddy brine of the crater lake. Aragonite and Mg-calcite dominate in the carbonate crusts and in sub-seafloor concretions at the MV foot, whereas Mg-calcite, dolomite and ankerite dominate in the muddy brine lake and in sub-seafloor concretions near the crater rim. The carbonate crusts and sub-seafloor concretions at the MV foot precipitated in isotopic equilibrium with bottom seawater temperature; their low δ¹³C values (–42.6 to –24.5‰) indicate that anaerobic oxidation of methane was the main driver of carbonate precipitation. By contrast, carbonates from the muddy lake brine, bottom lake concretions and crater rim concretions display much higher δ¹³C (up to –5.2‰) and low δ¹⁸O values (down to –2.8‰); this is consistent with their formation in warm fluids of deep origin characterized by ¹³C-rich CO₂ and, as confirmed by independent evidence, slightly higher heavy rare earth element signatures, the main driver of carbonate precipitation being methanogenesis. Moreover, the benthic activity within the seafloor sediment enhances aerobic oxidation of methane and of sulphide that promotes carbonate dissolution and gypsum precipitation. These findings imply that the coupling of carbon and sulphur microbial reactions represents the major link for the transfer of elements and for carbon isotope fractionation between fluids and authigenic minerals. A new challenge awaiting future studies in cold seep environments is to expand this work to oxidized and reduced sulphur authigenic minerals.     

Geophysical constraints on deep weathering and water storage potential in the Southern Sierra Critical Zone Observatory

The conversion of bedrock to regolith marks the inception of critical zone processes, but the factors that regulate it remain poorly understood. Although the thickness and degree of weathering of regolith are widely thought to be important regulators of the development of regolith and its water‐storage potential, the functional relationships between regolith properties and the processes that generate it remain poorly documented. This is due in part to the fact that regolith is difficult to characterize by direct observations over the broad scales needed for process‐based understanding of the critical zone. Here we use seismic refraction and resistivity imaging techniques to estimate variations in regolith thickness and porosity across a forested slope and swampy meadow in the Southern Sierra Critical Zone Observatory (SSCZO). Inferred seismic velocities and electrical resistivities image a weathering zone ranging in thickness from 10 to 35 m (average = 23 m) along one intensively studied transect. The inferred weathering zone consists of roughly equal thicknesses of saprolite (P‐velocity < 2 km s^?1) and moderately weathered bedrock (P‐velocity = 2–4 km s^?1). A minimum‐porosity model assuming dry pore space shows porosities as high as 50% near the surface, decreasing to near zero at the base of weathered rock. Physical properties of saprolite samples from hand augering and push cores are consistent with our rock physics model when variations in pore saturation are taken into account. Our results indicate that saprolite is a crucial reservoir of water, potentially storing an average of 3 m³ m^?2 of water along a forested slope in the headwaters of the SSCZO. When coupled with published erosion rates from cosmogenic nuclides, our geophysical estimates of weathering zone thickness imply regolith residence times on the order of 10⁵ years. Thus, soils at the surface today may integrate weathering over glacial–interglacial fluctuations in climate. Copyright © 2013 John Wiley & Sons, Ltd.