A decade of change in the Skidaway River estuary I. Hydrography and nutrients

The Skidaway River estuary is a tidally-dominated subtropical estuary in the southeastern USA surrounded by extensiveSpartina salt marshes. Weekly smapling at high and low tide began in 1986 for hydrography, nutrients, chlorophylla, particulate matter, and microbial and plankton biomass and composition; hydrographic and nutrient data during 1986–1996 are reported here. Salinity varied inversely with river discharge and exhibited variability at all time scales but with no long-term trend. Water temperature typically ranged over 25°C and was without apparent long-term frend. Seasonal cycles in concentrations of NO₃, NH₄, PO₄, Si(OH)₄, and DON were observed, with annual maxima generally occurring in late summer. Superimposed on seasonal cycles, all five nutrients exhibited steady increases in minimum, mean, and maximum concentrations; mean concentrations increased c. 50–150% during the decade. Nutrient concentrations were highly correlated with water temperature over the ten-year period, but weakly related to salinity and discharge. Nutrients were strongly correlated with one another, and the relative ratios among inorganic nutrients showed little long-term trend. Correlations among temperature and nutrient concentrations exhibited considerable inter-annual variability. Major spikes in organic and inorganic nutrient concentrations coincided with significant rainfall events; concentrations increased hyperbolically with rainfall. Although pristine compared to more heavily impacted waterways primarily outside the region, residential development and population density have been increasing rapidly during the past 15–20 years. Land use is apparently altering nutrient loading over the long-term (months-years), and superimposed on this are stochastic meteorological events that accelerate these changes over the short term (days-weeks).     

Numerical Homogenization of the Rigidity Tensor in Hooke's Law Using the Node-Based Finite Element Method

Combining a geological model with a geomechanical model, it generally turns out that the geomechanical model is built from units that are at least a 100 times larger in volume than the units of the geological model. To counter this mismatch in scales, the geological data model's heterogeneous fine-scale Young's moduli and Poisson's ratios have to be “upscaled” to one “equivalent homogeneous” coarse-scale rigidity. This coarse-scale rigidity relates the volume-averaged displacement, strain, stress, and energy to each other, in such a way that the equilibrium equation, Hooke's law, and the energy equation preserve their fine-scale form on the coarse scale. Under the simplifying assumption of spatial periodicity of the heterogeneous fine-scale rigidity, homogenization theory can be applied. However, even then the spatial variability is generally so complex that exact solutions cannot be found. Therefore, numerical approximation methods have to be applied. Here the node-based finite element method for the displacement as primary variable has been used. Three numerical examples showing the upper bound character of this finite element method are presented.     

Mathematical modeling of PCB bioaccumulation in Perna viridis

In the present work, we built a mathematical model of polychlorinated biphenyl (PCB) bioaccumulation in Perna viridis, namely, a one-compartment model with a time dependent incorporation rate R (μg g⁻¹ lipid per ppb water per day), with positive substrate cooperativity as the underlying physical mechanism. The temporal change of the PCB concentration Q (μg g⁻¹ lipid) in the soft tissues of the mussel depends on the competition of the input rate RW and the output rate kQ, where W is the concentration of PCB in water (ppb water) and k is the elimination rate (per day). From our experimental data, k=0.181±0.017 d⁻¹. The critical concentration in water W_c for positive substrate cooperativity was found to be 2.4 ppb. Below W_c, R is a constant. For a water concentration of 0.5 ppb Aroclor 1254, R=24.0±2.4 μg g⁻¹ lipid ppb⁻¹ d⁻¹. Above W_c, positive substrate cooperativity comes into effect and R becomes a function of time and dependent on the concentration Q in a form R=γQ/(Q+δ). This is the case for a water concentration of 5 ppb Aroclor 1254, where γ=15.1 μg g⁻¹ lipid ppb⁻¹ d¹ and δ≈200 μg g⁻¹ lipid. From this model, the uptake is exponentially increasing when the PCB concentration in the mussel is small compared to 200 μg g⁻¹ lipid, and hyperbolically increasing when the concentration is large compared to 200 μg g⁻¹ lipid, which are consistent with the experimental data. The model is useful for understanding the true processes taking place during the bioaccumulation and for risk assessment with higher confidence. Future experimental data which challenge the present model are anticipated and in fact desirable for improvement and perfection of the model.     

Determination of S-wave Site Response in Anchorage, Alaska in the 1–9 Hz Frequency Band

—?Site response was estimated at 19 sites in the Anchorage basin in south-central Alaska, using 15 local earthquakes recorded with good signal-to-noise ratio by a temporary weak motion network. The receiver-function-type horizontal-to-vertical spectral ratios (HVSR) were computed at 1–9 Hz frequency band and the resulting HVSR contour maps at 1, 5 and 9 Hz are presented here. The spatial site response distribution shows considerable variation from the foothills of the Chugach Mountains in the east to the western part of Anchorage. The site response increases by a factor of 3 and 3.5 at 1 and 5 Hz, respectively, from the area of older glacial deposits in the eastern part of the city to the area occupied by the Bootlegger Cove formation, particularly in the section adjoining Knik Arm. At 9 Hz, the variation of HVSR from the east to the west is smaller, approximately by a factor of 2. Moreover, the trend of the HVSR variation at 1 and 5 Hz shows good correlation with that of the soil class obtained from surface measurements of S-wave velocity in the 0–30?m depth range and available results on ground failure susceptibility of Anchorage.     

Seismic evaluation of existing nuclear facility using ambient vibration test to characterize dynamic behavior of the structure and microtremor measurements to characterize the soil: a case study

The methods used for a building seismic hazard evaluation are presented with the associated results. The goals of the study are (1) to check the soil nature and the existence or not of a possible site effect around the installation and (2) to characterize the dynamic behavior of the building using ambient vibration records.

The results of the soil study with the Nakamura method are very difficult to interpret because they are not stable in space and time. The spectral ratios method has been used with regional earthquake records. The results of the application of this method allowed us to conclude that the installation was free of site effect.

The ambient vibration measurements on the building brought the conclusion to determine the first and second modes of the structure. These results have been used to calibrate numerical model. The modal shapes in plan (high roof) and in elevation (main column) have been evaluated. The damping of the building has been computed using ambient vibration records.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

The Influence of 3-D Structure on the Propagation of Seismic Waves Away from Earthquakes

-- The seismic records from significant earthquakes are profoundly affected by 3-D variations in crustal structure both in the source zone itself and in propagation to some distance. Even in structurally complex zones such as Japan and Mexico relatively coherent arrivals are found associated with different classes of propagation paths. The presence of strong lateral variations can disrupt the arrivals, and impose significant variations in propagation characteristics for different directions from the source as illustrated by observations for the 1995 Kobe and 2000 Tottori-ken Seibu earthquakes in western Japan. Such effects can be modelled in 3 dimensions using a hybrid scheme with a pseudospectral representation for horizontal coordinates and finite differences in depth. This arrangement improves parallel implementation by minimising communication costs. For a realistic 3-D model for the structure in western Japan the 3-D simulations to frequencies close to 1 Hz provide a good representation of the observations from subduction zones events such as the 1946 Nankai earthquake and the 2000 Tottori-ken Seibu earthquake. The model can therefore be used to investigate the pattern of ground motion expected for future events e.g., in current seismic gaps.     

Studying Thermal Eutrophication in Cooling Reservoirs of Nuclear Power Plants

The processes of thermal eutrophication of the cooling reservoirs of the Smolensk and Kursk nuclear power plants (NPP) were studied. Data characterizing changes in the concentrations of inorganic forms of nitrogen and phosphates in water during its heating in the NPP cooling system were obtained. Different aspects of the effect of heated water on the production–destruction processes are discussed. The conditions under which the intensity of thermal eutrophication can substantially increase were determined.     

Optimized Seismic Threshold Monitoring – Part 2: Teleseismic Processing

v--vThis second paper (Part 2) pertaining to optimized site-specific threshold monitoring addresses the application of the method to regions covered by a teleseismic or a combined regional-teleseismic network. In the first paper (Part 1) we developed the method for the general case, and demonstrated its application to an area well-covered by a regional network (the Novaya Zemlya nuclear test site). In the present paper, we apply the method to the Indian and Pakistani nuclear test sites, and show results during the periods of nuclear testing by these two countries in May 1998. Since the coverage by regional stations in these areas is poor, an optimized approach requires the use of selected, high-quality stations at teleseismic distances.¶To optimize the threshold monitoring of these test sites, we use as calibration events either one of the nuclear explosions or a nearby earthquake. From analysis of the calibration events we derive values for array beamforming steering delays, filter bands, short-term averages (STA) lengths, phase travel times (P waves), and amplitude-magnitude relationships for each station. By applying these parameters, we obtain a monitoring capability of both test sites ranging from m_b 2.8-3.0 using teleseismic stations only. When including the nearby Nilore station to monitor the Indian tests, we show that the threshold can be reduced by about 0.4 magnitude units. In particular, we demonstrate that the Indian tests on 13 May, 1998, which were not detected by any known seismic station, must have corresponded to a magnitude (m_b) of less than 2.4.¶We also discuss the effect of a nearby aftershock sequence on the monitoring capability for the Pakistani test sites. Such an aftershock sequence occurred in fact on the day of the last Pakistani test (30 May, 1998), following a large (m_b 5.5) earthquake in Afghanistan located about 1100 km from the test site. We show that the threshold monitoring technique has sufficient resolution to suppress the signals from these interfering aftershocks without significantly affecting the true peak of the nuclear explosion on the threshold trace.