Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

The structure of rotating stars from an extension of the Kippenhahn-Thomas method

A method is presented whereby the structure of rotating stars may be determined from an initial guess at the geometry of equipotential surfaces. The method may be considered an extension of the work of Kippenhahn and Thomas in that a uniformly continuous geometry is defined in terms of the appropriate spherical model with Roche characteristics at the surface of the configuration and sphericity at the centre. A simple Cowling model in uniform rotation is employed to illustrate the technique and for comparison purposes with previous work.     

The constancy of galactic cosmic rays as recorded by cosmogenic nuclides in iron meteorites

We measured the He, Ne, and Ar isotopic concentrations and the ¹⁰Be, ²⁶Al, ³⁶Cl, and ⁴¹Ca concentrations in 56 iron meteorites of groups IIIAB, IIAB, IVA, IC, IIA, IIB, and one ungrouped. From ⁴¹Ca and ³⁶Cl data, we calculated terrestrial ages indistinguishable from zero for six samples, indicating recent falls, up to 562 ± 86 ka. Three of the studied meteorites are falls. The data for the other 47 irons confirm that terrestrial ages for iron meteorites can be as long as a few hundred thousand years even in relatively humid conditions. The ³⁶Cl‐³⁶Ar cosmic ray exposure (CRE) ages range from 4.3 ± 0.4 Ma to 652 ± 99 Ma. By including literature data, we established a consistent and reliable CRE age database for 67 iron meteorites. The high quality of the CRE ages enables us to study structures in the CRE age histogram more reliably. At first sight, the CRE age histogram shows peaks at about 400 and 630 Ma. After correction for pairing, the updated CRE age histogram comprises 41 individual samples and shows no indications of temporal periodicity, especially not if one considers each iron meteorite group separately. Our study contradicts the hypothesis of periodic GCR intensity variations (Shaviv 2002, 2003), confirming other studies indicating that there are no periodic structures in the CRE age histogram (e.g., Rahmstorf et al. 2004; Jahnke 2005). The data contradict the hypothesis that periodic GCR intensity variations might have triggered periodic Earth climate changes. The ³⁶Cl‐³⁶Ar CRE ages are on average 40% lower than the ⁴¹K‐K CRE ages (e.g., Voshage 1967). This offset can either be due to an offset in the ⁴¹K‐K dating system or due to a significantly lower GCR intensity in the time interval 195–656 Ma compared to the recent past. A 40% lower GCR intensity, however, would have increased the Earth temperature by up to 2 °C, which seems unrealistic and leaves an ill‐defined ⁴¹K‐K CRE age system the most likely explanation. Finally, we present new ²⁶Al/²¹Ne and ¹⁰Be/²¹Ne production rate ratios of 0.32 ± 0.01 and 0.44 ± 0.03, respectively.     

Calibration of the kinematic method of studying solar wind disturbances on the basis of a one-dimensional MHD solution and a simulation study of the heliosphere disturbances between 22 november and 6 december 1977

An empirical kinematic method developed by Hakamada and Akasofu (1982) is calibrated on the basis of a one-dimensional MHD solution. The calibrated results are used to simulate the stream-stream interaction and the background corotating structure in a simple situation and also during 22 November–6 December 1977. The solar wind disturbances caused by solar activities during this period are then introduced into the above background stream in simulating the heliospheric disturbance event which was observed by an aligned set of spacecraft at distances between 0.6 and 1.6 a.u. The observations and the simulated results are satisfactory, and a little more refinement in the simulation could reconstruct reasonably well the data by filling the data gaps in the solar wind speed, the density and the IMF magnitude.     

Photoadaptation of carboxylating enzymes and photosynthesis during a spring bloom

Properties of the light saturation curve of photosynthesis and ribulose-1,5-bisphosphate carboxylase (RuBPC) activity are shown to change qualitatively in a natural population of marine phytoplankton during a spring bloom. Evidence is presented to show that these changes constitute photoadapative responses to increasing irradiance. As irradiance increased during the bloom, both the level of light-saturated photosynthesis (P_m) and the initial slope of the light saturation curve (α = photosynthetic efficiency) increased whether those parameters were normalized to chlorophyll a concentration (P_m^b, α^b) or to cell numbers (P_m^c, α^c). The magnitudes of these changes were such that I_k (= P_m/α, the photoadaptation parameter) did not change, but I_m, the light intensity at which photosynthesis becomes saturated, increased. RuBPC activity, both chlorophyll a (RuBPC^b) and cell number normalized (RuBPC^c), also increased during the bloom. We suggest that these adaptations were achieved by simultaneously increasing the number of photosynthetic units, proportionately decreasing the photosynthetic unit size, and increasing both the concentrations of the enzymes of the dark reactions and possibly also of photosynthetic electron transport components.We also observed diminished levels of photoinhibition in the high light adapted cells late in the bloom and have suggested that this was a consequence of the same suite of physiological changes.In situ carbon fixation per cell increased during the bloom whereas no change occurred in this parameter when normalized to chlorophyll a concentration. Although these photoadaptive responses thus permitted carbon to be fixed in situ more rapidly per cell, at a constant efficiency with respect to investment of energy in the photosynthetic apparatus, they did not result in a change in growth rate. Based on consideratios of the role of time scale in physiological adaptation, however, it is suggested that the observed alterations in photosynthesis with increasing irradiance might permit a cell to more rapidly fill an energy quota for division, possibly an advantage in a mixing environment in which energy is patchily distributed, both spatially and temporalyy.Phosphoenolpyruvate carboxylase activity when normalized to chlorophyll a (PEPC^b) did not change during the bloom while chlorophyll a normalized dark carbon fixation decreased sharply and was quantitatively small compared to PEPC^b. On this basis and considering that RuBPC^b increased during the bloom, it is suggested that, although PEPC may be involved in dark carbon fixation, its most important quantitative role is probably an indirect one in light dependent photosynthesis.We have also considered the relevance of laboratory results on photoadaptation to interpretations of field studies and have suggested that batch culture studies must be treated with caution but that turbidistat and semi-continuous methods provide reasonable simulations of natural conditions.     

Deep-Sea Sediment Compression Curves: Some Controlling Factors, Spurious Overconsolidation, Predictions, and Geophysical Reproduction

Oedometer tests have been carried out on 70 undisturbed surficial clays (at approximately 250 mm below the mudline), mostly collected by free-fall corers from sites widely scattered throughout the deep-sea North Atlantic. Acoustic measurements were also made, initially on contiguous samples and ultimately on the same sample using a geophysically instrumented oedometer which also collected electrical resistivity data. Apart from those quiescent areas below the carbonate compensation depth, such as north of the West Indies where very fine clays exist, most of the samples are silty clays whose geotechnical-geophysical properties are dependent on the type of clay minerals present (and their ability to take in moisture), the sand-size fraction, and the quantity of carbonate present. Thus the pure clays have high compressibilities which decrease on the addition of coarse particles, while the converse is true for the acoustic parameters, these increasing with the sand fraction. Using the notion of the intrinsic compression line for all samples, and comparison to it of the measured compression curves, it is clear that, contrary to some previously held ideas, most deep-sea clays are normally consolidated; the addition of carbonate has the effect of creating an open, stronger sediment skeleton. Interestingly, where information is available, the variation with depth of a sample's acoustic velocity follows the void ratio pressure relationship of the compression curve. This allows the construction of an in-situ sediment compression curve using the in-situ geophysical observations.     

Variability of extreme precipitation in coastal river basins of the southern mexican Pacific region

Extreme wet and dry years (± 1 standard deviation, respectively), as well as the top 95 percentile (P95) of daily precipitation events, derived from tropical cyclone (TC) and nontropical cyclone (NTC) rainfall, were analyzed in coastal river basins in Southern Oaxaca, Mexico (Río Verde, Río Tehuantepec, and the Southern Coast). The study is based on daily precipitation records from 47 quality-controlled stations for the 1961 to 1990 period and TC data for the Eastern Tropical Pacific (EPAC). The aim of this study was to evaluate extreme (dry and wet) trends in the annual contribution of daily P95 precipitation events and to determine the relationship of summer precipitation with El Niño Southern Oscillation (ENSO) and the Pacifical Decadal Oscillation (PDO). A regionalization based on a rotated principal component analysis (PCA) was used to produce four precipitation regions in the coastal river basins. A significant negative correlation (significance at the 95% level) was only found with ONI in rainfall Region 3, nearest to the Gulf of Tehuantepec. Wet years, mainly linked to TC-derived P95 precipitation events, were associated with SST anomalies (≥?0.6°C) similar to weak La Niña and Neutral cool conditions, while dry years were associated with SST positive anomalies similar to Neutral warm conditions (≤?0.5°C). The largest contribution of extreme P95 precipitation derived from TCs to the annual precipitation was observed in Region 3. A significant upward trend in the contribution of TC-derived precipitation to the annual precipitation was found only in Region 1, low Río Verde.     

Multi‐Criteria Analysis for Optimization of Sodium Chromate Production from Chromic Waste

The proposal of application of multi‐criteria analysis and technological quality method for optimization of technological processes on the case study of sodium chromate production variants from waste was presented. The research on the use of chromic waste for the production of sodium chromate made it possible to determine the optimal process parameters. Technological quality method was used to define the value of criterions for multi‐criteria analysis while multi‐criteria analysis allowed determining the impact of criterions weights for ranking variants of technological process. Both methods in a complex way evaluate the degree of modernization of technologies, describing them qualitatively within the area of environmental, technical, and the economic effects of the proposed technological variants. Based on the proposed evaluating criteria both methods yield comparable results. The most advantageous alternatives provided realization of the sodium chromate production process with the use of in‐process recycling of chromic mud, off‐site recycling of chromic tannery waste and on‐site recycling of waste from the old chromic heaps.     

Using constrained inversion of gravity and magnetic field to produce a 3D litho‐prediction model

Geologically constrained inversion of gravity and magnetic field data of the Victoria property (located in Sudbury, Canada) was undertaken in order to update the present three‐dimensional geological model. The initial and reference model was constructed based on geological information from over 950 drillholes to constrain the inversion. In addition, downhole density and magnetic susceptibility measured in six holes were statistically analysed to derive lower and upper bounds on the physical properties attributed to the lithological units in the reference model. Constrained inversion of the ground gravity and the airborne magnetic data collected at the Victoria property were performed using GRAV3D and MAG3D, respectively. A neural network was trained to predict lithological units from the physical properties measured in six holes. Then, the trained network was applied on the three‐dimensional distribution of physical properties derived from the inversion models to produce a three‐dimensional litho‐prediction model. Some of the features evident in the lithological model are remnants of the constraints, where the data did not demand a significant change in the model from the initial constraining model (e.g., the thin pair of diabase dykes). However, some important changes away from the initial model are evident; for example, a larger body was predicted for quartz diorite, which may be related to the prospective offset dykes; a new zone was predicted as sulfide, which may represent potential mineralisation; and a geophysical subcategory of metabasalt was identified with high magnetic susceptibility and high density. The litho‐prediction model agrees with the geological expectation for the three‐dimensional structure at Victoria and is consistent with the geophysical data, which results in a more holistic understanding of the subsurface lithology.