Calculation of beach change under interacting cross-shore and longshore processes

This paper presents a mathematical approach and numerical model that simulates beach and dune change in response to cross-shore processes of dune growth by wind and dune erosion by storms, and by gradients in longshore sand transport that will alter shoreline position. Sub-aerial transport processes are represented, whereas sub-aqueous transport is neglected. The system is tightly coupled morphologically, with the berm playing a central role. For example, the potential for sand to be transported to the dune by wind depends on berm width, and sand lost in erosion of the dune during storms can widen the berm. Morphologic equilibrium considerations are introduced to improve reliability of predictions and stability of the non-linear model. An analytical solution is given under simplification to illustrate properties of the model. Sensitivity tests with the numerical solution of the coupled equations demonstrate model performance, with one test exploring beach and dune response to potential increase in storm-wave height with global warming. Finally, the numerical model is applied to examine the consequences of groin shortening at Westhampton Beach, Long Island, New York, as an alternative for providing a sand supply to the down-drift beach. Results indicate that the sand will be released over several decades as the shoreline and dune move landward in adjustment to the new equilibrium condition with the shortened groins.     

Chromite oxidation by manganese oxides in subseafloor basalts and the presence of putative fossilized microorganisms

Chromite is a mineral with low solubility and is thus resistant to dissolution. The exception is when manganese oxides are available, since they are the only known naturally occurring oxidants for chromite. In the presence of Mn(IV) oxides, Cr(III) will oxidise to Cr(VI), which is more soluble than Cr(III), and thus easier to be removed. Here we report of chromite phenocrysts that are replaced by rhodochrosite (Mn(II) carbonate) in subseafloor basalts from the Koko Seamount, Pacific Ocean, that were drilled and collected during the Ocean Drilling Program (ODP) Leg 197. The mineral succession chromite-rhodochrosite-saponite in the phenocrysts is interpreted as the result of chromite oxidation by manganese oxides. Putative fossilized microorganisms are abundant in the rhodochrosite and we suggest that the oxidation of chromite has been mediated by microbial activity. It has previously been shown in soils and in laboratory experiments that chromium oxidation is indirectly mediated by microbial formation of manganese oxides. Here we suggest a similar process in subseafloor basalts.     

A Numerical Framework for Modeling Folds in Structural Geology

A numerical framework for modeling folds in structural geology is presented. This framework is based on a novel and recently published Hamilton–Jacobi formulation by which a continuum of layer boundaries of a fold is modeled as a propagating front. All the fold classes from the classical literature (parallel folds, similar folds, and other fold types with convergent and divergent dip isogons) are modeled in two and three dimensions as continua defined on a finite difference grid. The propagating front describing the fold geometry is governed by a static Hamilton–Jacobi equation, which is discretized by upwind finite differences and a dynamic stencil construction. This forms the basis of numerical solution by finite difference solvers such as fast marching and fast sweeping methods. A new robust and accurate scheme for initialization of finite difference solvers for the static Hamilton–Jacobi equation is also derived. The framework has been integrated in simulation software, and a numerical example is presented based on seismic data collected from the Karama Block in the North Makassar Strait outside Sulawesi.     

Finite element modeling of hydraulic fracturing in 3D

A procedure based on the finite element method is suggested for modeling of 3D hydraulic fracturing in the subsurface. The proposed formulation partitions the stress field into the initial stress state and an additional stress state caused by pressure buildup. The additional stress is obtained as a solution of the Biot equations for coupled fluid flow and deformations in the rock. The fluid flow in the fracture is represented on a regular finite element grid by means of “fracture” porosity, which is the volume fraction of the fracture. The use of the fracture porosity allows for a uniform finite element formulation for the fracture and the rock, both with respect to fluid pressure and displacement. It is demonstrated how the fracture aperture is obtained from the displacement field. The model has a fracture criterion by means of a strain limit in each element. It is shown how this criterion scales with the element size. Fracturing becomes an intermittent process, and each event is followed by a pressure drop. A procedure is suggested for the computation of the pressure drop. Two examples of hydraulic fracturing are given, when the pressure buildup is from fluid injection by a well. One case is of a homogeneous rock, and the other case is an inhomogeneous rock. The fracture geometry, well pressure, new fracture area, and elastic energy released in each event are computed. The fracture geometry is three orthogonal fracture planes in the homogeneous case, and it is a branched fracture in the inhomogeneous case.     

Effect of sub‐core scale heterogeneities on acoustic and electrical properties of a reservoir rock: a CO2 flooding experiment of brine saturated sandstone in a computed tomography scanner

The effect of sub‐core scale heterogeneity on fluid distribution pattern, and the electrical and acoustic properties of a typical reservoir rock was studied by performing drainage and imbibition flooding tests with CO₂ and brine in a laboratory. Moderately layered Rothbach sandstone was used as a test specimen. Two core samples were drilled; one perpendicular and the other parallel to the layering to allow injection of fluids along and normal to the bedding plane. During the test 3D images of fluid distribution and saturation levels were mapped by an industrial X‐ray CT‐scanner together with simultaneous measurement of electrical resistivity, ultrasonic velocities as well as amplitudes. The results showed how the layering and the flooding direction influenced the fluid distribution pattern and the saturation level of the fluids. For a given fluid saturation level, the measured changes in the acoustic and electrical parameters were affected by both the fluid distribution pattern and the layering orientation relative to the measurement direction. The P‐wave amplitude and the electrical resistivity were more sensitive to small changes in the fluid distribution patterns than the P‐wave velocity. The change in amplitude was the most affected by the orientation of the layering and the resulting fluid distribution patterns. In some instances the change due to the fluid distribution pattern was higher than the variation caused by the change in CO₂ saturation. As a result the Gassmann relation based on ‘uniform' or ‘patchy' saturation pattern was not suitable to predict the P‐wave velocity variation. Overall, the results demonstrate the importance of core‐imaging to improve our understanding of fluid distribution patterns and the associated effects on measured rock‐physics properties.     

Changing patterns of residential and workplace segregation in the Stockholm metropolitan area

Immigrant–native segregation is present in the spaces in which individuals from different ethnic/racial groups practice their everyday lives; interact with others and develop their ethnic, social and spatial networks. The overwhelming majority of academic research on immigrant segregation has focused on the residential domain, thus largely overlooking other arenas of daily interaction. The present study contributes to the emerging literature on immigrant residential and workplace segregation by examining changes in patterns of residential and workplace segregation over time. We draw our data from the Stockholm metropolitan region, Sweden’s main port of entry for immigrants. The results suggest a close association between residential and workplace segregation. Immigrant groups that are more segregated at home are also more segregated in workplace neighborhoods. More importantly, we found that a changing segregation level in one domain tends to involve a similar trend in the other domain.     

Recolonisation by Macrobenthos Mobilises Organic Phosphorus from Reoxidised Baltic Sea Sediments

In recent decades, eutrophication has increased the extent of hypoxic and anoxic conditions in many coastal marine environments. In such conditions, the nutrient flux across the sediment?Cwater interface is a key process controlling the biogeochemical dynamics, and thereby the level and character of biological production. In some areas, management attempts to drive the ecosystem towards phosphorus (P) limitation, which calls for reliable knowledge on the mechanisms controlling P-cycling. We report a well-controlled laboratory experiment on benthic fluxes of P, when shifting from a state of hypoxic and azoic sediments to oxic and zoic bottom conditions. Adding any of three types of macrobenthic fauna (mysid shrimp, pontoporeid amphipod and tellinid clam) to oxygenated aquarium sections resulted in benthic P fluxes that differed consistently from the azoic control sections. All species caused liberation of dissolved organically bound P (DOP) from the sediment, in contrast to the azoic systems. The shrimp and the amphipod also resuspended the sediment, which resulted in a release of P bound to particles (>0.45???m). Dissolved inorganic phosphate (DIP) was released during hypoxic conditions, but was taken up after oxygenation, irrespective of the presence or absence of bottom fauna. In the presence of fauna, the uptake of DIP roughly equalled the release of DOP, suggesting that the benthic efflux of DOP following oxygenation and bottom fauna (re)colonisation might be considerable. This is an hitherto overlooked animal-controlled nutrient flux, which is missing from coastal marine P budgets.     

Accuracy and efficiency of stencils for the eikonal equation in earth modelling

Motivated by the needs for creating fast and accurate models of complex geological scenarios, accuracy and efficiency of three stencils for the isotropic eikonal equation on rectangular grids are evaluated using a fast marching implementation. The stencils are derived by direct modelling of the wave front, resulting in new and valuable insight in terms of improved upwind and causality conditions. After introducing a method for generalising first-order upwind stencils to higher order, a new second-order diagonal stencil is presented. Similarly to the multistencil fast marching approach, the diagonal stencil makes use of nodes in the diagonal directions, whereas the traditional Godunov stencil uses solely edge-connected neighbours. The diagonal stencil uses nodes close to each other, reaching upwind, to get a more accurate estimate of the angle of incidence of the arriving wave front. Although the stencils are evaluated in a fast marching setting, they can be adapted to other efficient eikonal solvers. All first- and second-order stencils are evaluated in a range of tests. The first test case models a folded structure from the Zagros fold belt in Iran. The other test cases are constructed to investigate specific properties of the examined stencils. The numerical investigation considers convergence rates and CPU times for non-constant and constant speed first-arrival computations. In conclusion, the diagonal stencil is the most efficient and accurate of the three alternatives.