Impact of dynamic feedbacks between sedimentation,sea-level rise,and biomass production on near-surface marsh stratigraphy and carbon accumulation

Salt marshes accrete both organic and inorganic sediments. Here we present analytical and numerical models of salt marsh sedimentation that, in addition to capturing inorganic processes, explicitly account for above- and belowground organic processes including root growth and decay of organic carbon. The analytical model is used to examine the bias introduced by organic processes into proxy records of sedimentation, namely ¹³⁷Cs and ²¹⁰Pb. We find that accretion rates estimated using ²¹⁰Pb will be less than accretion rates estimated using the ¹³⁷Cs peak in steadily accreting marshes if (1) carbon decay is significant and (2) data for ²¹⁰Pb extend below the ¹³⁷Cs peak. The numerical model expands upon the analytical model by including belowground processes such as compaction and root growth, and by explicitly tracking the evolution of aboveground biomass and its effect on sedimentation rates. Using the numerical model we explore how marsh stratigraphy responds to sediment supply and the rate of sea-level rise. It is calibrated and tested using an extensive data set of both marsh stratigraphy and measurements of vegetation dynamics in a Spartina alterniflora marsh in South Carolina, USA. We find that carbon accumulation in marshes is nonlinearly related to both the supply of inorganic sediment and the rate of sea-level rise; carbon accumulation increases with sea-level rise until sea-level rise reaches a critical rate that drowns the marsh vegetation and halts carbon accumulation. The model predicts that changes in carbon storage resulting from changing sediment supply or sea-level rise are strongly dependent on the background sediment supply: if inorganic sediment supply is reduced in an already sediment poor marsh the storage of organic carbon will increase to a far greater extent than in a sediment-rich marsh, provided that the rate of sea-level rise does not exceed a threshold. These results imply that altering sediment supply to estuaries (e.g., by damming upstream rivers or altering littoral sediment transport) could lead to significant changes in the carbon budgets of coastal salt marshes.     

Application of a nonlinear frequency domain wave–current interaction model to shallow water recurrence effects in random waves

The effect of ambient currents on nearshore nonlinear wave–wave energy transfer in random waves is studied with the use of a nonlinear frequency domain wave–current interaction model. We focus on the phenomenon of wave recurrence as a classical nonlinear phenomenon whose characteristics are well established for systems truncated to small numbers of frequency modes. The model used for this study is first extended to enhance accuracy; comparisons of permanent form solutions to analytical forms confirm the model accuracy. Application of the model to a highly truncated system confirmed the model’s consistency with published results for both positive (following) and negative (adverse) currents. Propagation of random wave spectra over a flat bottom was performed with the model, with the intent of determining the prevalence of recurrence between the spectral peak and its harmonics. For spectra of moderate Ursell number, it was found that positive currents extended the length scale of recurrence relative to the case with no currents; conversely, negative currents reduced the recurrence lengths. However, beyond a propagation distance of ≈40 wavelengths of the spectral peak, recurrence becomes almost completely damped as the spectra becomes broad and the spectral energies equilibrate. For spectra of high Ursell number, in contrast, recurrence is almost immediately damped, suggesting that the nonlinearity is sufficient to allow immediate spectral broadening and equilibration and overwhelming any preferential interactions among the spectral peak and its harmonics, regardless of current magnitude or direction.     

Nearshore suspended sediment variations,central Surinam coast

Nearshore suspended sediment concentration along the muddy Surinam coast is highly variable; maximum values are many times greater than on other muddy coasts. Water samples taken at four field stations during various stages of the tide range in concentration from 15 to 3,700 mg/l near the surface and from 100 to 30,000 mg/l near the bottom. Highest overall concentrations and greatest variability in concentration occur in water over large banks of fluid mud (thixotropic gel) that extend 2–3 km offshore and 5–10 km along-shore. On both the intertidal and subaqueous portions of these mudbanks highest concentrations are found at low tide. Results provide evidence that an exchange between fluid mud and suspended sediment takes place during each tidal cycle.     

Non-destructive estimation of potato leaf chlorophyll from canopy hyperspectral reflectance using the inverted PROSAIL model

Optimizing nitrogen (N) fertilization in crop production by in-season measurements of crop N status may improve fertilizer N use efficiency. Hyperspectral measurements may be used to assess crop N status by estimating leaf chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content. Trials were conducted with two potato cultivars under different N fertility rates (0–300 kg N ha⁻¹). Canopy reflectance, leaf area index (LAI) and leaf chlorophyll and N contents were measured. The PROSAIL model was able to predict leaf chlorophyll content with reasonable accuracy later in the growing season. The low estimation accuracy earlier in the growing season could be due to model sensitivity to non-homogenous canopy architecture and soil background interference before full canopy closure. Canopy chlorophyll content (leaf chlorophyll content × LAI) was predicted less accurately than leaf chlrophyll content due to the low estimation accuracy of LAI for values higher than 4.5.     

Correction to: Projection of Iron Ore Production

Natural Resources Research - The original version of this article unfortunately contained an error in Equation 3.     

Geomorphic constraints on the active tectonics of southern Taiwan

Taiwan is a region of rapid active tectonics, yet the study of the tectonic processes that shape the interior of the island is difficult due to the high rates of erosion and dense vegetation. We use digital topography to look for indications of active deformation preserved in the local geomorphology. In particular, anomalies in the regional pattern of drainage are used to infer zones of enhanced tectonic activity. The apparent anticlockwise rotation of major river systems in plan view indicates the presence of a diffuse zone of left-lateral shear running down the southeastern side of Taiwan. Asymmetries in the catchments of individual drainage basins show the influence of varying rates of uplift across southern Taiwan, with the most rapid uplift close to Taitung at the indentation point of the Luzon arc with the Chinese continental margin. Our interpretations, though based predominantly on remote-sensing observations, are consistent with the available field evidence. This study demonstrates the usefulness of drainage basins as tectonic markers in the quantification of regional strain and uplift, which may have wider applicability in other deforming parts of the world.     

Reduced complexity strategies for modelling urban floodplain inundation

Significant advances in flood inundation modelling have been made in the last decade through the use of a new generation of 2D hydraulic numerical models. These offer the potential to predict the local pattern and timing of flood depth and velocity, enabling informed flood risk zoning and improved emergency planning. With the availability of high resolution DEMs derived from airborne lidar, these models can theoretically now be routinely parameterized to represent considerable topographic complexity, even in urban areas where the potential exists to represent flows at the scale of individual buildings. Currently, however, computational constraints on conventional finite element and volume codes typically require model discretization at scales well below those achievable with lidar and are thus unable to make optimal use of this emerging data stream.In this paper we review two strategies that attempt to address this mismatch between model and data resolution in an effort to improve urban flood forecasts. The first of these strives for a solution by simplifying the mathematical formulation of the numerical model by using a computationally efficient 2D raster storage cell approach coupled to a 1D channel model. This parsimonious model structure enables simulations over large model domains offering the opportunity to employ a topographic discretization strategy which explicitly represents the built environment. The second approach seeks to further reduce the computational overhead of this raster method by employing a subgrid parameterization to represent the effect of buildings and micro-relief on flow pathways and floodplain storage. This multi-scale methodology enables highly efficient model applications at coarse spatial resolutions while retaining information about the complex geometry of the built environment.These two strategies are evaluated through numerical experiments designed to reconstruct a flood in the small town of Linton in southern England, which occurred in response to a 1 in 250 year rainfall event in October 2001. Results from both approaches are encouraging, with the spatial pattern of inundation and flood wave propagation matching observations well. Both show significant advantages over a coarse resolution model without subgrid parameterisation, particularly in terms of their ability to reproduce both hydrograph and inundation depth measurements simultaneously, without need for recalibration. The subgrid parameterization is shown to achieve this without contributing significant computational complexity and reduces model run-times by an order of magnitude.     

Storm rainfall conditions for floods and debris flows from recently burned areas in southwestern Colorado and southern California

Debris flows generated during rain storms on recently burned areas have destroyed lives and property throughout the Western U.S. Field evidence indicate that unlike landslide-triggered debris flows, these events have no identifiable initiation source and can occur with little or no antecedent moisture. Using rain gage and response data from five fires in Colorado and southern California, we document the rainfall conditions that have triggered post-fire debris flows and develop empirical rainfall intensity–duration thresholds for the occurrence of debris flows and floods following wildfires in these settings. This information can provide guidance for warning systems and planning for emergency response in similar settings.Debris flows were produced from 25 recently burned basins in Colorado in response to 13 short-duration, high-intensity convective storms. Debris flows were triggered after as little as six to 10 min of storm rainfall. About 80% of the storms that generated debris flows lasted less than 3 h, with most of the rain falling in less than 1 h. The storms triggering debris flows ranged in average intensity between 1.0 and 32.0 mm/h, and had recurrence intervals of two years or less. Threshold rainfall conditions for floods and debris flows sufficiently large to pose threats to life and property from recently burned areas in south-central, and southwestern, Colorado are defined by: I = 6.5D^− 0.7 and I = 9.5D^− 0.7, respectively, where I = rainfall intensity (in mm/h) and D = duration (in hours).Debris flows were generated from 68 recently burned areas in southern California in response to long-duration frontal storms. The flows occurred after as little as two hours, and up to 16 h, of low-intensity (2–10 mm/h) rainfall. The storms lasted between 5.5 and 33 h, with average intensities between 1.3 and 20.4 mm/h, and had recurrence intervals of two years or less. Threshold rainfall conditions for life- and property-threatening floods and debris flows during the first winter season following fires in Ventura County, and in the San Bernardino, San Gabriel and San Jacinto Mountains of southern California are defined by I = 12.5D^−0.4, and I = 7.2D^−0.4, respectively. A threshold defined for flood and debris-flow conditions following a year of vegetative recovery and sediment removal for the San Bernardino, San Gabriel and San Jacinto Mountains of I = 14.0D^−0.5 is approximately 25 mm/h higher than that developed for the first year following fires.The thresholds defined here are significantly lower than most identified for unburned settings, perhaps because of the difference between extremely rapid, runoff-dominated processes acting in burned areas and longer-term, infiltration-dominated processes on unburned hillslopes.     

Decadal‐Scale Changes of the Ödenwinkelkees,Central Austria,Suggest Increasing Control of Topography and Evolution Towards Steady State

Small mountain glaciers have short mass balance response times to climate change and are consequently very important for short‐term contributions to sea level. However, a distinct research and knowledge gap exists between (1) wider regional studies that produce overview patterns and trends in glacier changes, and (2) in situ local scale studies that emphasise spatial heterogeneity and complexity in glacier responses to climate. This study of a small glacier in central Austria presents a spatiotemporally detailed analysis of changes in glacier geometry and changes in glaciological behaviour. It integrates geomorphological surveys, historical maps, aerial photographs, airborne LiDAR data, ground‐based differential global positioning surveys and Ground Penetrating Radar surveys to produce three‐dimensional glacier geometry at 13 time increments spanning from 1850 to 2013. Glacier length, area and volume parameters all generally showed reductions with time. The glacier equilibrium line altitude increased by 90 m between 1850 and 2008. Calculations of the mean bed shear stress rapidly approaching less than 100 kPA, of the volume–area ratio fast approaching 1.458, and comparison of the geometric reconstructions with a 1D theoretical model could together be interpreted to suggest evolution of the glacier geometry towards steady state. If the present linear trend in declining ice volume continues, then the Ödenwinkelkees will disappear by the year 2040, but we conceptualise that non‐linear effects of bed overdeepenings on ice dynamics, of supraglacial debris cover on the surface energy balance, and of local topographically driven controls, namely wind‐redistributed snow deposition, avalanching and solar shading, will become proportionally more important factors in the glacier net balance.