A coupling of mixed and continuous Galerkin finite element methods for poroelasticity I: the continuous in time case

In this paper, we formulate a finite element procedure for approximating the coupled fluid and mechanics in Biot’s consolidation model of poroelasticity. Here, we approximate the pressure by a mixed finite element method and the displacements by a Galerkin method. Theoretical convergence error estimates are derived in a continuous in-time setting for a strictly positive constrained specific storage coefficient. Of particular interest is the case when the lowest-order Raviart–Thomas approximating space or cell-centered finite differences are used in the mixed formulation, and continuous piecewise linear approximations are used for displacements. This approach appears to be the one most frequently applied to existing reservoir engineering simulators.     

Stabilization of Expansive Clays Using Granulated Blast Furnace Slag (GBFS) and GBFS-Cement

Expansive clays undergo swelling when subjected to water. This can cause damage, especially to light weight structures, water conveyance canals, lined reservoirs, highways, and airport runways unless appropriate measures are taken. In this study, granulated blast furnace slag (GBFS) and GBFS-cement (GBFSC) were utilized to overcome or to limit the expansion of an artificially prepared expansive soil sample (sample A). GBFS and GBFSC were added to sample A in proportions of 5–25% by weight. The effects of these stabilizers on grain size distribution, Atterberg limits, swelling percentage and rate of swell of soil samples were determined. GBFS and GBFSC were shown to successfully decreasing the total amount of swell while increasing the rate of swell.     

Temporal variability in summertime bottom hypoxia in shallow areas of Mobile Bay,Alabama

This paper addresses temporal variability in bottom hypoxia in broad shallow areas of Mobile Bay, Alabama. Time-series data collected in the summer of 2004 from one station (mean depth of 4 m) exhibit bottom dissolved oxygen (DO) variations associated with various time scales of hours to days. Despite a large velocity shear, stratification was strong enough to suppress vertical mixing most of the time. Bottom DO was closely related to the vertical salinity gradient (ΔS). Hypoxia seldom occurred when ΔS (over 2.5 m) was <2 psu and occurred almost all the time when ΔS was >8 psu in the absence of extreme events like hurricanes. Oxygen balance between vertical mixing and total oxygen demand was considered for bottom water from which oxygen demand and diffusive oxygen flux were estimated. The estimated decay rates at 20°C ranging between 0.175–0.322 d⁻¹ and the corresponding oxygen consumption as large as 7.4 g O₂ m⁻² d⁻¹ fall at the upper limit of previously reported ranges. The diffusive oxygen flux and the corresponding vertical diffusivity estimated for well mixed conditions range between 8.6–9.5 g O₂ m⁻² d⁻¹ and 2.6–2.9 m² d⁻¹, respectively. Mobile Bay hypoxia is likely to be associated with a large oxygen demand, supported by both water column and sediment oxygen demands, so that oxygen supply from surface water during destratification events would be quickly exhausted to return to hypoxic conditions within a few hours to days after destratification events are terminated.     

Probabilistic seismic hazard assessment for South-Eastern France

The accurate evaluation and appropriate treatment of uncertainties is of primary importance in modern probabilistic seismic hazard assessment (PSHA). One of the objectives of the SIGMA project was to establish a framework to improve knowledge and data on two target regions characterized by low-to-moderate seismic activity. In this paper, for South-Eastern France, we present the final PSHA performed within the SIGMA project. A new earthquake catalogue for France covering instrumental and historical periods was used for the calculation of the magnitude-frequency distributions. The hazard model incorporates area sources, smoothed seismicity and a 3D faults model. A set of recently developed ground motion prediction equations (GMPEs) from global and regional data, evaluated as adequately representing the ground motion characteristics in the region, was used to calculate the hazard. The magnitude-frequency distributions, maximum magnitude, faults slip rate and style-of-faulting are considered as additional source of epistemic uncertainties. The hazard results for generic rock condition (Vs30 = 800 m/s) are displayed for 20 sites in terms of uniform hazard spectra at two return periods (475 years and 10,000 years). The contributions of the epistemic uncertainties in the ground motion characterizations and in the seismic source characterization to the total hazard uncertainties are analyzed. Finally, we compare the results with existing models developed at national scale in the framework of the first generation of models supporting the Eurocode 8 enforcement, (MEDD 2002 and AFPS06) and at the European scale (within the SHARE project), highlighting significant discrepancies at short return periods.     

Interpretation of Gravity Anomalies with the Normalized Full Gradient (NFG) Method and an Example

The Normalized Full Gradient (NFG) method which was put forward about 50 years ago has been used for downward continuation of gravity potential data, especially in the former Union of Soviet Socialist Republics. This method nullifies perturbations due to the passage of mass depth during downward continuation. The method depends on the downwards analytical continuation of normalized full gradient values of gravity data. Analytical continuation discriminates certain structural anomalies which cannot be distinguished in the observed gravity field. This method has been used in various petroleum and tectonic studies. The Trapeze method was used for the determination of Fourier coefficients during the application of this method. No other techniques for calculating these coefficients have been used. However, the Filon method was used for the determination of Fourier coefficients during the application of the NFG method in this work. This method, rather than the Trapeze method, should be preferred for indicating abnormal mass resources at the lower harmonics. In this study, the NFG method using the Filon method has been applied the first time to theoretical models of gravity profiles as example field at the Hasankale-Horasan petroleum exploration province where successful results were achieved. Hydrocarbon presence was shown on the NFG sections by the application of NFG downward continuation operations on theoretical models. Important signs of hydrocarbon structure on the NFG section for field and model data at low harmonics are obtained more effectively using this method.     

From paleo to policy: partitioning the historical point and nonpoint phosphorus loads to the St. Croix River,Minnesota-Wisconsin,USA

Paleolimnological studies show that phosphorus (P) loads to the federally protected St. Croix River, a tributary of the Upper Mississippi River, have increased about threefold over the last century. Ongoing management efforts to protect and restore the river hinge on the question of whether the increased nutrient load results from point-source discharges or nonpoint runoff from agricultural intensification and urban expansion. Here we determine the historical contribution of point source phosphorus (P) loads to the St. Croix watershed from 1900–2000 A.D. Historical point source loads were estimated based on discharge volumes, demographics, industrial sources, wastewater technologies, and facility discharge records, where available. Sewering in the basin began in 1905, and since that time, there have been as many as 169 permitted point source dischargers basinwide, including municipal, industrial, and agricultural facilities. Early wastewater management typically discharged untreated sewage; technological advances had secondary treatment in place at most facilities by the 1960s–1970s and much of the municipal population was served by tertiary treatment by the 1990s. Peak nutrient discharges from point sources occurred in the 1960s–1970s. Detergent phosphorus bans instituted in the late 1970s for Minnesota and Wisconsin, greater use of land and groundwater effluent disposal, and improvements in treatment technology brought about decreases in P loads in the 1980s and 1990s. Point-source discharges were compared to historical total phosphorus loads estimated in a whole-basin phosphorus mass balance to calculate the historical contribution of point sources, anthropogenic nonpoint sources, and natural or background sources. We estimated 1990s point source loads at 48 t P yr⁻¹, which represents about 10% of the total phosphorus load (459 t P yr⁻¹, flow-corrected to 412 t P yr⁻¹) to the basin. Without further controls on nutrient inputs to the St. Croix River, annual flow-corrected P loads are projected to increase to 498 t P yr⁻¹ by the 2020s with point source phosphorus loading contributions at 65 t P yr⁻¹ or 13% of the total load. However, if we exclude background P loads to the St. Croix (166 t P yr⁻¹), recent nutrient loads are primarily from anthropogenic nonpoint sources. Point sources also contribute over 19% of the current and future phosphorus load that can be attributed to human activities in the watershed. Interstate and federal efforts to decrease P loading to the St. Croix River by 20% will need to target both point and nonpoint sources. This is one of eight papers dedicated to the “Recent Environmental History of the Upper Mississippi River” published in this special issue of the Journal of Paleolimnology. D.R. Engstrom served as guest editor of the special issue.     

Identification of forcing mechanisms of convective initiation over mountains through high-resolution numerical simulations

正Convection and its ensuing severe weather,such as heavy rainfall,hail,tornado,and high wind,have significant impacts on our society and economy(e.g.,Cao et al.,2004;Fritsch and Carbone,2004;Verbout et al.,2006;Ashley and Black,2008;Cao,2008;Cao and Ma,2009;Zhang et al.,2014).Due to its localized and transient nature,the initiation of convection or convective initiation remains one of the least     

Parameterization of Sheared Convective Entrainment in the First-Order Jump Model: Evaluation Through Large-Eddy Simulation

In this note, two different approaches are used to estimate the entrainment-flux to surface-flux ratio for a sheared convective boundary layer (CBL); both are derived under the framework of the first-order jump model (FOM). That suggested by Sun and Wang (SW approach) has the advantage that there is no empirical constant included, though the dynamics are described in an implicit manner. The second, which was proposed by Kim et al. and Pino et al. (KP approach), explicitly characterizes the dynamics of the sheared entrainment, but uncertainties are induced through the empirical constants. Their performances in parameterizing the CBL growth rate are compared and discussed, and a new value of the parameter A ₃ in the KP approach is suggested. Large-eddy simulation (LES) data are employed to test both approaches: simulations are conducted for the CBL growing under varying conditions of surface roughness, free-atmospheric stratification, and wind shear, and data used when the turbulence is in steady state. The predicted entrainment rates in each case are tested against the LES data. The results show that the SW approach describes the evolution of the sheared CBL quite well, and the KP approach also reproduces the growth of the CBL reasonably, so long as the value of A ₃ is modified to 0.6.     

Effect of perturbations on the non linear stability of triangular points in the restricted three-body problem with variable mass

The effect of small perturbations ε and ε ^′ in the Coriolis and the centrifugal forces, respectively on the nonlinear stability of the triangular points in the restricted three-body problem with variable mass has been studied. It is found that, in the nonlinear sense, the triangular points are stable for all mass ratios in the range of linear stability except for three mass ratios, which depend upon ε, ε ^′ and β, the constant due to the variation in mass governed by Jeans’ law.