Present-day oxidative subsidence of organic soils and mitigation in the Sacramento-San Joaquin Delta,California, USA

Subsidence of organic soils in the Sacramento-San Joaquin Delta threatens sustainability of the California (USA) water supply system and agriculture. Land-surface elevation data were collected to assess present-day subsidence rates and evaluate rice as a land use for subsidence mitigation. To depict Delta-wide present-day rates of subsidence, the previously developed SUBCALC model was refined and calibrated using recent data for CO₂ emissions and land-surface elevation changes measured at extensometers. Land-surface elevation change data were evaluated relative to indirect estimates of subsidence and accretion using carbon and nitrogen flux data for rice cultivation. Extensometer and leveling data demonstrate seasonal variations in land-surface elevations associated with groundwater-level fluctuations and inelastic subsidence rates of 0.5–0.8 cm yr^–1. Calibration of the SUBCALC model indicated accuracy of ±0.10 cm yr^–1 where depth to groundwater, soil organic matter content and temperature are known. Regional estimates of subsidence range from <0.3 to >1.8 cm yr^–1. The primary uncertainty is the distribution of soil organic matter content which results in spatial averaging in the mapping of subsidence rates. Analysis of leveling and extensometer data in rice fields resulted in an estimated accretion rate of 0.02–0.8 cm yr^–1. These values generally agreed with indirect estimates based on carbon fluxes and nitrogen mineralization, thus preliminarily demonstrating that rice will stop or greatly reduce subsidence. Areas below elevations of –2 m are candidate areas for implementation of mitigation measures such as rice because there is active subsidence occurring at rates greater than 0.4 cm yr^–1.     

Models of travel demand with endogenous preference change and heterogeneous agents

In the literatures of regional science, urban economics, and urban development planning, a working assumption is that individuals respond to incentives and regulations, given their preferences. Models for planning and policy analyses are used to consider what might occur if the incentives or regulations were different. In these models, however, preferences are usually assumed to be given and stable, and agents are usually assumed to be homogeneous. This paper focuses on the implications of making preferences in models of policy implementation endogenously determined and time varying heterogeneous agents. We consider first the recent literature on intertemporal choice and preference change, which cuts across many disciplines, and more briefly the literature on norm-regarding behavior. We then elaborate a simple model of transportation demand—from a static to a dynamic orientation, from fixed and exogenously given preferences of strictly self-regarding agents to endogenously determined and policy-induced preferences of heterogeneous agents—and illustrate its characteristics with simple numerical examples.     

The formation of porphyry copper deposits

Copper is a moderately incompatible chalcophile element. Its behavior is strongly controlled by sulfides. The speciation of sulfur is controlled by oxygen fugacity. Therefore, porphyry Cu deposits are usually oxidized (with oxygen fugacities > ΔFMQ +2) (Mungall 2002; Sun et al. 2015). The problem is that while most of the magmas at convergent margins are highly oxidized, porphyry Cu deposits are very rare, suggesting that high oxygen fugacity alone is not sufficient. Partial melting of mantle peridotite even at very high oxygen fugacities forms arc magmas with initial Cu contents too low to form porphyry Cu deposits directly (Lee et al. 2012; Wilkinson 2013). Here we show that partial melting of subducted young oceanic slabs at high oxygen fugacity (>ΔFMQ +2) may form magmas with initial Cu contents up to >500 ppm, favorable for porphyry mineralization. Pre-enrichment of Cu through sulfide saturation and accumulation is not necessarily beneficial to porphyry Cu mineralization. In contrast, re-melting of porphyritic hydrothermal sulfide associated with iron oxides may have major contributions to porphyry deposits. Thick overriding continental crust reduces the “leakage” of hydrothermal fluids, thereby promoting porphyry mineralization. Nevertheless, it is also more difficult for ore forming fluids to penetrate the thick continental crust to reach the depths of 2–4 km where porphyry deposits form.     

Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Lateral continuity of basement seismic reflections in 15 Ma ultrafast-spreading crust at ODP Site 1256

The Ocean Drilling Program (ODP) initiated drilling at Site 1256D in the Guatemala Basin, about 1,000 km off the East Pacific Rise to penetrate plutonic rocks, anticipated to be relatively shallow in this region, formed at an ultra-fast spreading rate. IODP Expedition E312 successfully drilled into gabbros at ~1,150 m in basement. Multi-channel seismic traces show weak laterally coherent sub-basement reflections at borehole depths. Synthetic reflectivity seismograms were computed using a Ricker wavelet and impedance profiles from borehole sonic logs. These seismograms show significant sub-basement amplitude peaks. A zero-offset vertical seismic profile, shot on E312, was processed to investigate the authenticity of these reflections and their relationship to borehole geology. A dual scheme of the median filtering and F–K dip filtering was used. Tests with synthetic seismograms indicate the approach is effective at reasonable SNR levels. Downgoing energy is clearly identified but negligible upgoing energy is visible over random noise. These results indicate that lava flows and igneous contacts in upper ocean crust have significant topography on lateral scales less than the Fresnel Zone (~300 m) due to igneous and tectonic processes.     

Mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradation

The mass-front velocities of granular flows results from the joint action of particle size gradations and the underlying surfaces.However,because of the complexity of friction during flow movement,details such as the slope-toe impedance effects and momentum-transfer mechanisms have not been completely explained by theoretical analyses,numerical simulations,or field investigations.To study the mass-front velocity of dry granular flows influenced by the angle of the slope to the runout plane and particle size gradations we conducted model experiments that recorded the motion of rapid and long-runout rockslides or avalanches.Flume tests were conducted using slope angles of 25°,35°,45°,and 55° and three particle size gradations.The resulting mass-front motions consisted of three stages:acceleration,velocity maintenance,and deceleration.The existing methods of velocity prediction could not explain the slowing effect of the slope toe or the momentum-transfer steady velocity stage.When the slope angle increased from 25° to 55°,the mass-front velocities dropped significantly to between 44.4% and59.6% of the peak velocities and energy lossesincreased from 69.1% to 83.7% of the initial,respectively.The velocity maintenance stages occurred after the slope-toe and mass-front velocity fluctuations.During this stage,travel distances increased as the angles increased,but the average velocity was greatest at 45°.At a slope angle of 45°,as the median particle size increased,energy loss around the slope toe decreased,the efficiency of momentum transfer increased,and the distance of the velocity maintenance stage increased.We presented an improved average velocity formula for granular flow and a geometrical model of the energy along the flow line.     

An Analytical Model For Threshold Crossing Rates Of Concentration Fluctuations In Dispersing Plumes

A simple analytical model is developed for the meanupcrossing rate of plume concentration fluctuations assuming that thisprocess can be well approximated by a lognormal process. The resultingexpression requires only the specification of the in-plume fluctuationintensity and in-plume Taylor micro-time scale and, hence, does notexplicitly involve the joint probability density function of theconcentration and its derivative. The analytical model provides agood fit to some field measurements of the mean upcrossing rate ina dispersing plume.     

Study on uncertainty of inversion and resolution limit in the presence of noise

This study analyzed the uncertainty of inversion and the resolution limit in the presence of noise by means of statistical experiments. The exhaustive method is adopted to obtain the global optimal solution in each experiment. We found that even with small level of noise, solutions fluctuate in a large range for the thin bed. The distribution of solutions in the presence of noise is closely related to the spread of the cost function in the absence of noise. As a result, the area of a certain neighborhood around the true solution on the spread of the cost function in the absence of noise is used to evaluate the uncertainty of inversion and the resolution limit in the presence of noise. In the case that the SNR (signal-to-noise ratio) is 5 in this study, solutions focus around the true solution with a very small uncertainty only when the bed thickness is greater than the reciprocal of the double predominant frequency of the convoluting wavelet.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.