Does variable channel morphology lead to dynamic salmon habitat?

Stream channel morphology forms the template upon which hydraulic aspects of aquatic habitat are created, yet spatial and temporal variability in habitat imposed by changing morphology is not well understood. This paper presents a conceptual model linking sediment supply patterns to spatial and temporal variability in channel form and aquatic habitat. To evaluate this model, change over time in three habitat variables is quantified using a 2D hydrodynamic modeling approach. A 45-year record of topographic data from Carnation Creek, a catchment in coastal British Columbia, is used for the flow modeling. Using the Nays2DH modeling platform, water depths and velocities are simulated in eight channel segments located at different positions relative to locations of historical colluvial input using seven flow levels ranging from 3% to 400% of mean annual discharge (0.02 to 3.31 m s ). Results indicate that habitat availability changes through time as a result of sediment supply-driven changes to channel morphology and wood loads, but patterns in habitat vary as a function of dominant channel segment morphology. Spatial and temporal variability in morphology also influences the relationship between habitat availability and river discharge, leading to non-stationary habitat-discharge rating curves. When habitat areas are predicted by applying these curves to daily flow series spanning annual dry seasons, over 50% of the variance in cumulative seasonal habitat area can be explained by year-to-year changes in channel morphology and wood loading, indicating that changing morphology is an important factor for driving temporal habitat variability. This variance is related to the morphological variability of a channel segment, which in turn is associated with the segment position relative to zones of colluvial input. Collectively, these results suggest that variability in habitat is impacted by channel morphology, and can be evaluated partly on the basis of a channel's sediment supply regime. © 2019 John Wiley & Sons, Ltd.     

Starved versus alluvial river bedforms: an experimental investigation

Laboratory experiments were conducted to investigate the formation of river bedforms under sediment supply-limited conditions, i.e. when a motionless substratum is bared by the dynamics of the mobile sediments. Three series of experiments were organized in a laboratory flume by fixing all the hydrodynamic and morphodynamic parameters but varying the thickness of the initial layer of mobile sediments which covers the rigid bottom of the flume. At the end of all the experiments, which lasted for the same amount of time, the formation of transverse sand dunes was observed. For decreasing , the rigid bottom of the flume was bared progressively earlier during the experiment and the measurements showed a clear tendency of the bedforms to lengthen, i.e. to increase their crest-to-crest distance. Moreover, under strong supply limitation, the two-dimensional transverse dunes turned into three-dimensional barchanoid forms and into isolated barchan dunes characterized by an abrupt reduction in bedform heights. A two-dimensional Fourier analysis of the bottom profile was performed, providing the amplitude of the main streamwise and spanwise harmonic components of the bottom morphology as a function of . © 2019 John Wiley & Sons, Ltd.     

Spatio-temporal characteristics of small-scale wave–current ripples on the Ameland ebb-tidal delta

Ebb-tidal deltas are highly dynamic environments affected by both waves and currents that approach the coast under various angles. Among other bedforms of various scales, these hydrodynamics create small-scale bedforms (ripples), which increase the bed roughness and will therefore affect hydrodynamics and sediment transport. In morphodynamic models, sediment transport predictions depend on the roughness height, but the accuracy of these predictors has not been tested for field conditions with strongly mixed (wave–current dominated) forcing. In this study, small-scale bedforms were observed in the field with a 3D Profiling Sonar at five locations on the Ameland ebb-tidal delta, the Netherlands. Hydrodynamic conditions ranged from wave dominated to current dominated, but were mixed most of the time. Small-scale ripples were found on all studied parts of the delta, superimposed on megaripples. Even though a large range of hydrodynamic conditions was encountered, the spatio-temporal variations in small-scale ripple dimensions were relatively small (height 0.015 m, length 0.11 m). Also, the ripples were always highly three-dimensional. These small dimensions are probably caused by the fact that the bed consists of relatively fine sediment. Five bedform height predictors were tested, but they all overestimated the ripple heights, partly because they were not created for small grain sizes. Furthermore, the predictors all have a strong dependence on wave- and current-related velocities, whereas the ripple heights measured here were only related to the near-bed orbital velocity. Therefore, ripple heights and lengths in wave–current-dominated, fine-grained coastal areas ( mm) may be best estimated by constant values rather than values dependent on the hydrodynamics. In the case of the Ameland ebb-tidal delta, these values were found to be m and m. ©2019 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Orthotropic anisotropy: on contributions of elasticity parameters to a difference in quasi-P-wave-squared velocities resulted from propagation in two orthogonal symmetry planes

We investigate the dependence of quasi P-wave phase velocity propagating in orthotropic media on particular elasticity parameters. Specifically, due to mathematical facilitation, we consider the squared-velocity difference, , resulted from propagation in two mutually perpendicular symmetry planes. In the context of the effective medium theory, may be viewed as a parameter evaluating the influence of cracks – embedded in the background medium – parallel to one or both aforementioned planes. Our investigation is both theoretical and numerical. Based on Christoffel's equations, we propose two accurate approximations of . Due to them, we interpret the aforementioned squared-velocity difference as being twice more dependent on , than on . To describe the magnitude of the dependence, we consider the proportions between the partial derivatives of . Further, it occurs that is influenced by the ratio of vertically propagating quasi P-wave to vertically propagating quasi S-wave. Anomalously high might be caused by the low P/S ratio, which in turn can be an indicator of the presence of gas in natural fractures or aligned porosity. Also, we carry out numerical sensitivity study, according to which is approximately twice more dependent on than on , twice more sensitive to than to , and equally dependent on as on . The dependence on and can be neglected, especially for small phase angles. We verify the approximations and perform the sensitivity study, using eight examples of the elasticity tensors.     

Estimating the timescale of fluvial response to anthropogenic disturbance using two generations of dams on the South River,Massachusetts, USA

Centuries-long intensive land-use change in the north-eastern United States provides the opportunity to study the timescale of geomorphic response to anthropogenic disturbances. In this region, forest-clearing and agricultural practices following EuroAmerican settlement led to deposition of legacy sediment along valley bottoms, including behind mill dams. The South River in western Massachusetts experienced two generations of damming, beginning with mill dams up to 6-m high in the eighteenth–nineteenth century, and followed by construction of the Conway Electric Dam (CED), a 17-m-tall hydroelectric dam near the watershed outlet in 1906. We use the mercury (Hg) concentration in upstream deposits along the South River to constrain the magnitude, source, and timing of inputs to the CED impoundment. Based on cesium-137 (¹³⁷Cs) chronology and results from a sediment mixing model, remobilized legacy sediment comprised % of the sediment load in the South River prior to 1954; thereafter, from 1954 to 1980s, erosion from glacial deposits likely dominated (63 ± 14%), but with legacy sediments still a substantial source (37 ± 14%). We also use the CED reservoir deposits to estimate sediment yield through time, and find it decreased after 1952. These results are consistent with high rates of mobilization of legacy sediment as historic dams breached in the early twentieth century, and suggest rapid initial response to channel incision, followed by a long decay in the second half of the century, that is likely dependent on large flood events to access legacy sediment stored in banks. Identifying sources of sediment in a watershed and quantifying erosion rates can help to guide river restoration practices. Our findings suggest a short fluvial recovery time from the eighteenth–nineteenth century to perturbation during the first half of the twentieth century, with subsequent return to a dominant long-term signal from erosion of glacial deposits, with anthropogenic sediment persisting as a secondary source. © 2020 John Wiley & Sons, Ltd.     

Multi-parameter reflection waveform inversion for acoustic transversely isotropic media with a vertical symmetry axis

Full waveform inversion in transversely isotropic media with a vertical symmetry axis provides an opportunity to better match the data at the near and far offsets. However, multi-parameter full waveform inversion, in general, suffers from serious cycle-skipping and trade-off problems. Reflection waveform inversion can help us recover a background model by projecting the residuals of the reflected wavefield along the reflection wavepath. Thus, we extend reflection waveform inversion to acoustic transversely isotropic media with a vertical symmetry axis utilizing the proper parameterization for reduced parameter trade-off. From a radiation patterns analysis, an acoustic transversely isotropic media with a vertical symmetry axis is better described by a combination of the normal-moveout velocity and the anisotropic parameters η and δ for reflection waveform inversion applications. We design a three-stage inversion strategy to construct the optimal resulting model. In the first stage, we only invert for the background by matching the simulated reflected wavefield from the perturbations of and δ with the observed reflected wavefield. In the second stage, the background and η are optimized simultaneously and the far-offset reflected wavefield mainly contribute to their updates. We perform Born modelling to compute the reflected wavefield for the two stages of reflection waveform inversion. In the third stage, we perform full waveform inversion for the acoustic transversely isotropic media with a vertical symmetry axis to delineate the high-wavenumber structures. For this stage, the medium is described by a combination of the horizontal velocity , η and ε instead of , η and δ. The acoustic multi-parameter full waveform inversion utilizes the diving waves to improve the background as well as utilizes reflection for high-resolution information. Finally, we test our inversion algorithm on the modified Sigsbee 2A model (a salt free part) and a two-dimensional line from a three-dimensional ocean bottom cable dataset. The results demonstrate that the proposed reflection waveform inversion approach can recover the background model for acoustic transversely isotropic media with a vertical symmetry axis starting from an isotropic model. This recovered background model can mitigate the cycle skipping of full waveform inversion and help the inversion recover higher resolution structures.     

The lateral and vertical growth of laterite weathering profiles,Hawaiian Islands,USA

The Hawaiian Islands permit investigation of tropical chemical weathering rates and processes on a single rock type, basalt. Chronosequences are investigated as a function of rainfall due to the varying age of each island, including Kauai (~4 Ma), Oahu (~2 Ma), and Hawaii's Kohala Peninsula (~0.3 to 0.17 Ma). Understanding tropical critical zone (CZ) development is vital given the large populations in developing countries that rely on it. HVSR (horizontal-to-vertical spectral ratio) seismic soundings on Kauai indicate that ~60% of the variability in laterite thickness is due to gradients in precipitation, with errors in erosion corrections and variability in the original permeability structure of the volcanic sequence playing important roles. Basalts have higher horizontal than vertical hydraulic conductivity (K_h > K_v) , and local variability in likely drives much of the remaining differences in laterite thickness. HVSR is well suited for estimating laterite thickness as it has been shown to reliably detect the base of the weathering profile, is rapid (20 min/sounding), highly portable, and occupies a very small footprint. Comparison of Kauai and Oahu weathering profiles suggests that the Oahu laterites are fully or nearly fully formed, despite being half the age of Kauai. By contrast, the young laterites on Kohala (~170 to ~300 ka) exhibit greatly contrasting thicknesses, where coastal laterites are thick and interior laterites are thin, suggesting that early weathering on shield volcanoes produces wedge-shaped laterites near the coast. With time, the thick (coastal) end of the wedge propagates upslope such that a fully developed, constant-thickness laterite carapace can form in ~2 Ma or less. The development of thickened coastal laterites on young substrates depends on greater water–rock ratios as vertically infiltrating water upslope is diverted laterally. This view of laterite development is very different compared to endmember models of continental weathering and CZ development. © 2020 John Wiley & Sons, Ltd.     

Experiment and Simulation of Non-Reactive Solute Transport in Porous Media

To more accurately predict the migration behavior of pollutants in porous media, we conduct laboratory scale experiments and model simulation. Aniline (AN) is used in one-dimensional soil column experiments designed under various media and hydrodynamic conditions. The advection-dispersion equation (ADE) and the continuous-time random walk (CTRW) were used to simulate the breakthrough curves (BTCs) of the solute transport. The results show that the media and hydrodynamic conditions are two important factors affecting solute transport and are related to the degree of non-Fickian transport. The simulation results show that CTRW can more effectively describe the non-Fickian phenomenon in the solute transport process than ADE. The sensitive parameter in the CTRW simulation process is , which can reflect the degree of non-Fickian diffusion in the solute transport. Understanding the relationship of with velocity and media particle size is conducive to improving the reactive solute transport model. The results of this study provide a theoretical basis for better prediction of pollutant transport in groundwater.     

Separation of multi-mode surface waves by supervised machine learning methods

Logistic regression, neural networks and support vector machines are tested for their effectiveness in isolating surface waves in seismic shot records. To distinguish surface waves from other arrivals, we train the algorithms on three distinguishing features of surface-wave dispersion curves in the domain: spectrum coherency of the trace's magnitude spectrum, local dip and the frequency range for a fixed wavenumber k in the spectrum. Numerical tests on synthetic data show that the kernel-based support vector machines algorithm gives the highest accuracy in predicting the surface-wave window in the domain compared to neural networks and logistic regression. This window is also used to automatically pick the fundamental dispersion curve. The other two methods correctly pick the low-frequency part of the dispersion curve but fail at higher frequencies where there is interference with higher-order modes.     

Investigating young water fractions in a small Mediterranean mountain catchment: Both precipitation forcing and sampling frequency matter

The proportion of water younger than 2–3 months (young water fraction, F_yw) has become increasingly investigated in catchment hydrology. F_yw is typically estimated by comparing seasonal tracer cycles in precipitation and streamflow, through water sampling. However, some open research questions remain, such as: (i) whether part of the summer precipitation should be discarded because the high evapotranspiration demand, (ii) how well F_yw serves as a metric to compare catchments, and (iii) how sampling frequency affects F_yw estimates. To address these questions, we investigated F_yw in soil-, ground- and stream waters for the small Mediterranean Can Vila catchment. Rainfall was sampled at 5-mm intervals. Mobile soil water and groundwater were sampled fortnightly. Stream water was sampled depending on flow at variable time intervals (30 min to 1 week). Over 58 months, this sampling provided 1,529 δ¹⁸O determinations. Isotopic analyses results led us to include summer precipitation in the input signal. We found the highest F_yw in mobile soil waters (34%), while this was almost zero for groundwater except during wet periods. For stream waters, F_yw depended on the discharge variations, so that the flow-weighted young water fraction () was 22.6%, whereas the time-weighted F_yw was just 6.2%. Both and its discharge sensitivity (S_d) varied when different 12-month sampling periods were investigated. The young water fraction that would be obtained from a virtual thorough sampling () was estimated from the S_d and the observed stream flow. This showed an underestimation of by 25% for the frequent dynamic sampling and 66% for weekly sampling, due to missing high flows. Our results confirm that F_yw and its discharge sensitivity are metrics very sensitive to meteorological forcing during the analysed period. Thus, comparisons between catchments need long-term mean annual values and their variability. Our findings also support the dependence of F_yw estimates on the sampling rate and show the advantages of flow-weighted analysis. Finally, catchment water turnover investigations should be accompanied by the analysis of flow duration curves.     

Post-fire management effects on sediment (dis)connectivity in Mediterranean forest ecosystems: Channel and catchment response

Forest fires and post-fire practices influence sediment connectivity (SC). In this study, we use the ‘aggregated index of connectivity’ (AIC) to assess SC in five Mediterranean catchments (198–1090 ha) affected by a wildfire in 2012 in south-eastern Spain. Two temporal scenarios were considered, immediately after the fire and before post-fire management, and 2 years after the fire including all practices (hillslope barriers, check-dams, afforestation, salvage logging and skid trails). One LiDAR (light detection and ranging)-derived digital elevation model (DEM, 2 m × 2 m resolution) was generated, per scenario. The five catchment outlets were established as the computation target (AIC_OUT), and structural and functional SC were calculated. Index outputs were normalized to make the results of the non-nested catchments comparable (AIC_N-OUT). The output analysis includes the SC distribution along the catchments and at local scale (929 sub-catchments, 677 in the burned area), the hillslope and channel measures' effect on SC, and a sedimentological analysis using observed area-specific sediment yield (SSY) at 10 new (built after post-fire practices) concrete check-dams located in the catchments (SSY = 1.94 Mg ha⁻¹ yr⁻¹; σ = 1.22). The catchments with more circular shapes and steeper slopes were those with higher AIC_N-OUT. The structural SC maps – removing the rainfall erosivity influence – allowed evaluating the actual role played by the post-fire practices that reduced SC ( $\overline{x}$= − 1.19%; σ = 0.41); while functional SC was linked to the actual change of SC ( $\overline{x}$= + 5.32%; σ = 0.62). Hillslope treatments resulted in significant changes on AIC_N-OUT at sub-catchment scale with certain disconnectivity. A good and positive correlation was found between the SSY and the changes of AIC_N-OUT. However, the coarse DEM resolution explained the lack of effect of the rock check-dams – located on the secondary channels – on AIC_N-OUT. AIC_N-OUT proved to be a useful tool for decision making in post-fire restoration, but an optimal input data is still necessary to refine calculations.     

The role of bathymetry and directional wave conditions on observed crescentic bar dynamics

Nearshore sandbars are important features in the surf zone of many beaches because they strongly influence the mean circulation and evolving morphology. Due to variations in wave conditions, sandbars can experience cross-shore migration and vary in shape from alongshore uniform (shore-parallel) to alongshore rhythmic (crescentic). Sandbar dynamics have been studied extensively, but existing observational studies usually do not quantify the processes leading to crescentic bar formation and straightening. This study analyses the dynamics of crescentic bar events at the fetch-limited beach of Castelldefels (northwestern Mediterranean Sea, Spain) using 7.5 years of hourly time-exposure video images and detailed wave conditions. The results show that, despite the generally calm wave conditions, the sandbars were very dynamic in the cross-shore and longshore directions. They often migrated rapidly offshore during storms (up to 70 m in one day) and more slowly onshore during post-storm conditions. Crescentic bars were often present at the study site (48% of the time), but only when the sandbar was at least 10 m from the shoreline. They displayed a large variability in wavelengths (100–700 m), alongshore migration speeds (0–50 m/day) and cross-shore amplitudes (5–20 m). Wavelengths increased for larger bar–shoreline distances and the alongshore migration speeds were strongly correlated with the alongshore component of the radiation stresses. Crescentic patterns typically developed during low–medium energetic waves with limited obliquity ( $\theta \lesssim 20$° at 10 m depth), while bar straightening occurred during medium–high energetic waves with strong oblique angles of incidence ( $\theta \gtrsim 15$°). Overall, this study provides further proof for the important role of wave direction in crescentic bar dynamics and highlights the strong dependence of crescentic bar development on the initial bathymetric configuration.     

Deterministic ground motion simulations with shallow crust nonlinearity at Garner Valley in Southern California

We present deterministic ground motion simulations that account for the cyclic multiaxial response of sediments in the shallow crust. We use the Garner Valley in Southern California as a test case. The multiaxial constitutive model is based on the bounding surface plasticity theory in terms of total stress and is implemented in a high‐performance computing finite‐element parallel code. A major advantage of this model is the small number of free parameters that need to be calibrated given a shear modulus reduction curve and the ultimate soil strength. This, in turn, makes the model suitable for regional‐scale simulations, where geotechnical data in the shallow crust are scarce. In this paper, we first describe a series of numerical experiments designed to verify the model implementation. This is followed by a series of idealized large‐scale simulations in a 35 26 4.5 km domain that encompasses the Garner Valley downhole array site, which is an instrumented and well‐characterized site in Southern California. Material properties were extracted from the Southern California Earthquake Center Community velocity model, CVM‐S4.26, considering its optional geotechnical layer, while the modulus reduction curves and soil strength were selected empirically to constrain the nonlinear soil model parameters. Our nonlinear simulations suggest that peak ground displacements within the valley increase relative to the linear case, while peak ground accelerations can increase or decrease, depending on the frequency content of the excitation. The comparisons of our simulations against hybrid three‐dimensional–one‐dimensional site response analyses suggest the inadequacy of the latter to capture the complexity of fully three‐dimensional simulations.     

Floodplain downed wood volumes: a comparison across three biomes

Downed large wood (LW) in floodplains provides habitat and nutrients for diverse organisms, influences hydraulics and sedimentation during overbank flows, and affects channel form and lateral migration. Very few studies, however, have quantified LW volumes in floodplains that are unaltered by human disturbance. We compare LW volumes in relatively unaltered floodplains of semiarid boreal lowland, subtropical lowland, and semiarid temperate mountain rivers in the United States. Average volumes of downed LW are 42.3 m³ ha^?1, 50.4 m³ ha^?1, and 116.3 m³ ha^?1 in the semiarid boreal, subtropical, and semiarid temperate sites, respectively. Observed patterns support the hypothesis that the largest downed LW volumes occur in the semiarid temperate mountain sites, which is likely linked to a combination of moderate‐to‐high net primary productivity, temperature‐limited decomposition rates, and resulting slow wood turnover time. Floodplain LW volumes differ among vegetation types within the semiarid boreal and semiarid temperate mountain regions, reflecting differences in species composition. Lateral channel migration and flooding influence vegetation communities in the semiarid boreal sites, which in turn influences floodplain LW loads. Other forms of disturbance such as fires, insect infestations, and blowdowns can increase LW volumes in the semiarid boreal and semiarid temperate mountain sites, where rates of wood decay are relatively slow compared with the subtropical lowland sites. Although sediment is the largest floodplain carbon reservoir, floodplain LW stores substantial amounts of organic carbon and can influence floodplain sediment storage. In our study sites, floodplain LW volumes are lower than those in adjacent channels, but are higher than those in upland (i.e. non‐floodplain) forests. Given the important ecological and physical effects of floodplain LW, efforts to add LW to river corridors as part of restoration activities, and the need to quantify carbon stocks within river corridors, we urge others to quantify floodplain and instream LW volumes in diverse environments. Copyright © 2016 John Wiley & Sons, Ltd.     

A computational workflow for rupture‐to‐structural‐response simulation and its application to Istanbul

Scenario‐based earthquake simulations at regional scales hold the promise in advancing the state‐of‐the‐art in seismic risk assessment studies. In this study, a computational workflow is presented that combines (i) a broadband Green's function‐based fault‐rupture and ground motion simulation—herein carried out using the “UCSB (University of California at Santa Barbara) method”, (ii) a three‐dimensional physics‐based regional‐scale wave propagation simulation that is resolved at  Hz, and (iii) a local soil‐foundation‐structure finite element analysis model. These models are interfaced with each other using the domain reduction method. The innermost local model—implemented in ABAQUS—is additionally enveloped with perfectly matched layer boundaries that absorb outbound waves scattered by the structures contained within it. The intermediate wave propagation simulation is carried out using Hercules , which is an explicit time‐stepping finite element code that is developed and licensed by the CMU‐QUAKE group. The devised workflow is applied to a  km region on the European side of Istanbul, which was modeled using detailed soil stratigraphy data and realistic fault rupture properties, which are available from prior microzonation surveys and earthquake scenario studies. The innermost local model comprises a chevron‐braced steel frame building supported by a shallow foundation slab, which, in turn, rests atop a three‐dimensional soil domain. To demonstrate the utility of the workflow, results obtained using various simplified soil‐structure interaction analysis techniques are compared with those from the detailed direct model. While the aforementioned demonstration has a limited scope, the devised workflow can be used in a multitude of ways, for example, to examine the effects of shallow‐layer soil nonlinearities and surface topography, to devise site‐ and structure‐specific seismic fragilities, and for calibrating regional loss models, to name a few.     

Wood process domains and wood loads on floodplains

Downed large wood on floodplains creates similar geomorphic and ecological effects as wood in the active channel, but has been the subject of fewer geomorphic studies. I propose floodplain large wood process domains that are distinguished based on recruitment source at the reach to river-length scale. Wood recruited to the floodplain can be autochthonous (individual or mass recruitment from floodplain forest), fluvially transported, or transported from adjacent hillslopes via mass movements that come down the valley side slopes or down the main channel. Fluvially transported wood can be further distinguished as being deposited: within the channel and subsequently accreted to the floodplain; marginal to the channel; on the floodplain during overbank flow; or on tributary fans. The mechanism of wood recruitment to a floodplain influences the spatial distribution of the wood across the floodplain and the proportion of wood pieces within jams, which in turn influences geomorphic and ecological effects of the floodplain wood. Using published studies of floodplain wood load for unmanaged river corridors, I hypothesize that the climate-controlled balance between forest primary productivity and decay rates of downed wood is the first-order control on floodplain large wood loads. Disturbance regime and wood recruitment mechanism are second-order controls on wood load and primary controls on the spatial distribution of large wood. Understanding of floodplain large wood can be applied to quantifying the effect of large wood on river corridors; river restoration; paleoenvironmental inferences; and estimation of organic carbon stock in river corridors. © 2019 John Wiley & Sons, Ltd.     

How geomorphic context governs the influence of wildfire on floodplain organic carbon in fire-prone environments of the western United States

We draw on published studies of floodplain organic carbon storage, wildfire-related effects on floodplains in temperate and high latitudes, and case studies to propose a conceptual model of the effects of wildfire on floodplain organic carbon storage in relation to climate and valley geometry. Soil organic carbon typically constitutes the largest carbon stock in floodplains in fire-prone regions, although downed wood can contain significant organic carbon. We focus on the influence of wildfire on soil organic carbon and downed wood as opposed to standing vegetation to emphasize the geomorphic influences resulting from wildfire on floodplain organic carbon stocks. The net effect of wildfire varies depending on site-specific characteristics including climate and valley geometry. Wildfire is likely to reduce carbon stock in steep, confined valley segments because increased water and sediment yields following fire create net floodplain erosion. The net effect of fire in partly confined valleys depends on site-specific interactions among floodplain aggradation and erosion, and, in high-latitude regions, permafrost degradation. In unconfined valleys in temperate latitudes, wildfire is likely to slightly increase floodplain organic carbon stock as a result of floodplain aggradation and wood deposition. In unconfined valleys in high latitudes underlain by permafrost, wildfire is likely in the short-term to significantly decrease floodplain organic carbon via permafrost degradation and reduce organic-layer thickness. Permafrost degradation reduces floodplain erosional resistance, leading to enhanced stream bank erosion and greater carbon fluxes into channels. The implications of warming climate and increased wildfires for floodplain organic carbon stock thus vary. Increasing wildfire extent, frequency, and severity may result in significant redistribution of organic carbon from floodplains to the atmosphere via combustion in all environments examined here, as well as redistribution from upper to lower portions of watersheds in the temperate zone and from floodplains to the oceans via riverine transport in the high-latitudes. © 2019 John Wiley & Sons, Ltd.     

Pre-stack seismic inversion based on -norm regularized logarithmic absolute misfit function

Ill-posedness is one of the most common and intractable issues that arise when solving geophysical inverse problems. Ill-posedness could be induced by various factors such as noise, band-limited intrinsic property of seismic data and inappropriate forward operators. Regularization has been proven to be an effective method widely accepted for mitigating the adverse effects of ill-posedness. Aiming to improve the stability and fidelity of the pre-stack seismic inversion process, we implement the inversion in a Bayesian framework, with a logarithmic absolute criterion taken as a likelihood function, and an -norm metric as a priori constraint. Here, we exploit the linear approximation as the forward operator, and optimize the regularized misfit function by the alternating direction method of multipliers. Applications of the method to synthetic and real data sets yielded improved inversion results in terms of accuracy and resolution, and demonstrated the robustness of the method to noise.     

Random noise attenuation via the randomized canonical polyadic decomposition

Tensor algebra provides a robust framework for multi-dimensional seismic data processing. A low-rank tensor can represent a noise-free seismic data volume. Additive random noise will increase the rank of the tensor. Hence, tensor rank-reduction techniques can be used to filter random noise. Our filtering method adopts the Candecomp/Parafac decomposition to approximates a N-dimensional seismic data volume via the superposition of rank-one tensors. Similar to the singular value decomposition for matrices, a low-rank Candecomp/Parafac decomposition can capture the signal and exclude random noise in situations where a low-rank tensor can represent the ideal noise-free seismic volume. The alternating least squares method is adopted to compute the Candecomp/Parafac decomposition with a provided target rank. This method involves solving a series of highly over-determined linear least-squares subproblems. To improve the efficiency of the alternating least squares algorithm, we uniformly randomly sample equations of the linear least-squares subproblems to reduce the size of the problem significantly. The computational overhead is further reduced by avoiding unfolding and folding large dense tensors. We investigate the applicability of the randomized Candecomp/Parafac decomposition for incoherent noise attenuation via experiments conducted on a synthetic dataset and field data seismic volumes. We also compare the proposed algorithm (randomized Candecomp/Parafac decomposition) against multi-dimensional singular spectrum analysis and classical prediction filtering. We conclude the proposed approach can achieve slightly better denoising performance in terms of signal-to-noise ratio enhancement than traditional methods, but with a less computational cost.     

Multi-objective fitting of concentration-discharge relationships

Concentration–discharge ($C–Q$) relationships are widely used to assess the link between hydrological and biogeochemical processes at the catchment scale. $C–Q$ relationships are mainly calibrated using mono-objective methods to represent, either concentrations or discharge-weighted concentrations (i.e., load). Based on its wide use in hydrological modelling, we test a multi-objective calibration for the $C–Q$ relationship parameters, using both concentration and load, and compare it to a mono-objective calibration applied on either concentrations or load. This work is carried out on a high-frequency dataset (ORACLE-Orgeval Observatory, France). Our findings show that the multi-objective calibration yield a better representation of $C–Q$ relationships parameters during the high and low-flow events. The multi-objective calibration can be used for all forms of $C–Q$ relationships and avoids issues of under-representation of dilution processes characterized by high-discharge, low-concentration periods.