Rare-earth elements and isotope ratios at DSDP Leg 82 Holes

We report 48 analyses of rare-earth elements (REE) and 15 143Nd/144Nd and 87Sr/86Sr analyses for basalts from the eight holes drilled during Leg 82. Discrete and distinct REE patterns and 143Nd/144Nd ratios characterize the eight holes, with little variation observed downhole except in Holes 561 and 558, thus suggesting dominantly long-term temporal and large-scale spatial variations in the mantle source of these basalts beneath the Mid-Atlantic Ridge over the last 35 Ma of its spreading activity. There is a good inverse correlation between 143Nd/144Nd and (La/Sm)EF with one exception in Hole 558 (approximately 35 Ma), the latter suggesting a recent (35 Ma) light REE depletion event, perhaps caused by dynamic or fractional melting. Short-term temporal and small-scale spatial mantle source variability is also evident in Hole 561 (approximately 18 Ma), which has rapid fluctuations in REE patterns and 143Nd/144Nd ratios (suggesting rapid transfer of magma from the time of melting) and is evidence contrary to the presence of a well-mixed magma chamber at this particular site and time. The mantle source variations noted can be interpreted within two extreme models. The first model invokes a convecting mantle depleted in large ion lithophile elements (LILE) and containing lumps (or veins) of LILE-enriched material of various shapes and sizes, passively and randomly distributed throughout. A second more restrictive model considers the interaction of fixed mantle plumes and the LILE-depleted asthenosphere flowing towards a migrating Mid- Atlantic Ridge (MAR) axis. With the exception of Hole 558 and the uncertainties of reconstructions of absolute plate movements in the region, the observed variations can be explained by two hot spots; the nearly ridge-centered Azores hot spot (plume) and another hot spot located beneath the African plate that may be affecting the source of basalts currently erupting at the MAR axis at 35°N and which, in the past, would have produced the New England chain of seamounts on the North American plate and (later) the Atlantis-Great Meteor chain on the African plate. Basalts erupted south of the Hayes Fracture Zone have not been affected by either of these two hot spots over the last 35 Ma and appear to have been continuously derived from the LILE-depleted source. Subaxial flow downridge from the Azores plume appears to have started 9 Ma, on the basis of the southward converging V-shaped time-transgressive ridges branching from the Pico and Corves Island, or not earlier than 16 Ma, on the basis of the geochemical results. Variations within Hole 558 remains unexplained by the latter model, unless we hypothesize a third hot spot.

Initiation of clement surface conditions on the earliest Earth

In the beginning the surface of the Earth was extremely hot, because the Earth as we know it is the product of a collision between two planets, a collision that also created the Moon. Most of the heat within the very young Earth was lost quickly to space while the surface was still quite hot. As it cooled, the Earth's surface passed monotonically through every temperature regime between silicate vapor to liquid water and perhaps even to ice, eventually reaching an equilibrium with sunlight. Inevitably the surface passed through a time when the temperature was around 100°C at which modern thermophile organisms live. How long this warm epoch lasted depends on how long a thick greenhouse atmosphere can be maintained by heat flow from the Earth's interior, either directly as a supplement to insolation, or indirectly through its influence on the nascent carbonate cycle. In both cases, the duration of the warm epoch would have been controlled by processes within the Earth's interior where buffering by surface conditions played little part. A potentially evolutionarily significant warm period of between 105 and 107 years seems likely, which nonetheless was brief compared to the vast expanse of geological time.

Theoretical and numerical analysis of large-scale heat transfer problems with temperature-dependent pore-fluid densities

Purpose - In many scientific and engineering fields, large-scale heat transfer problems with temperature-dependent pore-fluid densities are commonly encountered. For example, heat transfer from the mantle into the upper crust of the Earth is a typical problem of them. The main purpose of this paper is to develop and present a new combined methodology to solve large-scale heat transfer problems with temperature-dependent pore-fluid densities in the lithosphere and crust scales. Design/methodology/approach - The theoretical approach is used to determine the thickness and the related thermal boundary conditions of the continental crust on the lithospheric scale, so that some important information can be provided accurately for establishing a numerical model of the crustal scale. The numerical approach is then used to simulate the detailed structures and complicated geometries of the continental crust on the crustal scale. The main advantage in using the proposed combination method of the theoretical and numerical approaches is that if the thermal distribution in the crust is of the primary interest, the use of a reasonable numerical model on the crustal scale can result in a significant reduction in computer efforts. Findings - From the ore body formation and mineralization points of view, the present analytical and numerical solutions have demonstrated that the conductive-and-advective lithosphere with variable pore-fluid density is the most favorite lithosphere because it may result in the thinnest lithosphere so that the temperature at the near surface of the crust can be hot enough to generate the shallow ore deposits there. The upward throughflow (i.e. mantle mass flux) can have a significant effect on the thermal structure within the lithosphere. In addition, the emplacement of hot materials from the mantle may further reduce the thickness of the lithosphere. Originality/value - The present analytical solutions can be used to: validate numerical methods for solving large-scale heat transfer problems; provide correct thermal boundary conditions for numerically solving ore body formation and mineralization problems on the crustal scale; and investigate the fundamental issues related to thermal distributions within the lithosphere. The proposed finite element analysis can be effectively used to consider the geometrical and material complexities of large-scale heat transfer problems with temperature-dependent fluid densities.

A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity

Models of the air-sea transfer velocity of gases may be either empirical or mechanistic. Extrapolations of empirical models to an unmeasured gas or to another water temperature can be erroneous if the basis of that extrapolation is flawed. This issue is readily demonstrated for the most well-known empirical gas transfer velocity models where the influence of bubble-mediated transfer, which can vary between gases, is not explicitly accounted for. Mechanistic models are hindered by an incomplete knowledge of the mechanisms of air-sea gas transfer. We describe a hybrid model that incorporates a simple mechanistic view—strictly enforcing a distinction between direct and bubble-mediated transfer—but also uses parameterizations based on data from eddy flux measurements of dimethyl sulphide (DMS) to calibrate the model together with dual tracer results to evaluate the model. This model underpins simple algorithms that can be easily applied within schemes to calculate local, regional, or global air-sea fluxes of gases.

A reconciliation of empirical and mechanistic models of the air-sea gas transfer velocity

Models of the air-sea transfer velocity of gases may be either empirical or mechanistic. Extrapolations of empirical models to an unmeasured gas or to another water temperature can be erroneous if the basis of that extrapolation is flawed. This issue is readily demonstrated for the most well-known empirical gas transfer velocity models where the influence of bubble-mediated transfer, which can vary between gases, is not explicitly accounted for. Mechanistic models are hindered by an incomplete knowledge of the mechanisms of air-sea gas transfer. We describe a hybrid model that incorporates a simple mechanistic view—strictly enforcing a distinction between direct and bubble-mediated transfer—but also uses parameterizations based on data from eddy flux measurements of dimethyl sulphide (DMS) to calibrate the model together with dual tracer results to evaluate the model. This model underpins simple algorithms that can be easily applied within schemes to calculate local, regional, or global air-sea fluxes of gases.

Air-sea transfer of gas phase controlled compounds

Gases in the atmosphere/ocean have solubility that spans several orders of magnitude. Resistance in the molecular sublayer on the waterside limits the air-sea exchange of sparingly soluble gases such as SF6 and CO2. In contrast, both aerodynamic and molecular diffusive resistances on the airside limit the exchange of highly soluble gases (as well as heat). Here we present direct measurements of air-sea methanol and acetone transfer from two open cruises: the Atlantic Meridional Transect in 2012 and the High Wind Gas Exchange Study in 2013. The transfer of the highly soluble methanol is essentially completely airside controlled, while the less soluble acetone is subject to both airside and waterside resistances. Both compounds were measured concurrently using a proton-transfer-reaction mass spectrometer, with their fluxes quantified by the eddy covariance method. Up to a wind speed of 15 m s-1, observed air-sea transfer velocities of these two gases are largely consistent with the expected near linear wind speed dependence. Measured acetone transfer velocity is ~30% lower than that of methanol, which is primarily due to the lower solubility of acetone. From this difference we estimate the "zero bubble" waterside transfer velocity, which agrees fairly well with interfacial gas transfer velocities predicted by the COARE model. At wind speeds above 15 m s-1, the transfer velocities of both compounds are lower than expected in the mean. Air-sea transfer of sensible heat (also airside controlled) also appears to be reduced at wind speeds over 20 m s-1. During these conditions, large waves and abundant whitecaps generate large amounts of sea spray, which is predicted to alter heat transfer and could also affect the air-sea exchange of soluble trace gases. We make an order of magnitude estimate for the impacts of sea spray on air-sea methanol transfer.

Air-sea transfer of gas phase controlled compounds

Gases in the atmosphere/ocean have solubility that spans several orders of magnitude. Resistance in the molecular sublayer on the waterside limits the air-sea exchange of sparingly soluble gases such as SF6 and CO2. In contrast, both aerodynamic and molecular diffusive resistances on the airside limit the exchange of highly soluble gases (as well as heat). Here we present direct measurements of air-sea methanol and acetone transfer from two open cruises: the Atlantic Meridional Transect in 2012 and the High Wind Gas Exchange Study in 2013. The transfer of the highly soluble methanol is essentially completely airside controlled, while the less soluble acetone is subject to both airside and waterside resistances. Both compounds were measured concurrently using a proton-transfer-reaction mass spectrometer, with their fluxes quantified by the eddy covariance method. Up to a wind speed of 15 m s-1, observed air-sea transfer velocities of these two gases are largely consistent with the expected near linear wind speed dependence. Measured acetone transfer velocity is ~30% lower than that of methanol, which is primarily due to the lower solubility of acetone. From this difference we estimate the "zero bubble" waterside transfer velocity, which agrees fairly well with interfacial gas transfer velocities predicted by the COARE model. At wind speeds above 15 m s-1, the transfer velocities of both compounds are lower than expected in the mean. Air-sea transfer of sensible heat (also airside controlled) also appears to be reduced at wind speeds over 20 m s-1. During these conditions, large waves and abundant whitecaps generate large amounts of sea spray, which is predicted to alter heat transfer and could also affect the air-sea exchange of soluble trace gases. We make an order of magnitude estimate for the impacts of sea spray on air-sea methanol transfer.

Rational Function Interpolation of Electromagnetic Transfer Functions of High-Speed Interconnect Systems from Discrete Time-Domain and Frequency-Domain Data

We present new methodologies to generate rational function approximations of broadband electromagnetic responses of linear and passive networks of high-speed interconnects, and to construct SPICE-compatible, equivalent circuit representations of the generated rational functions. These new methodologies are driven by the desire to improve the computational efficiency of the rational function fitting process, and to ensure enhanced accuracy of the generated rational function interpolation and its equivalent circuit representation. Toward this goal, we propose two new methodologies for rational function approximation of high-speed interconnect network responses. The first one relies on the use of both time-domain and frequency-domain data, obtained either through measurement or numerical simulation, to generate a rational function representation that extrapolates the input, early-time transient response data to late-time response while at the same time providing a means to both interpolate and extrapolate the used frequency-domain data. The aforementioned hybrid methodology can be considered as a generalization of the frequency-domain rational function fitting utilizing frequency-domain response data only, and the time-domain rational function fitting utilizing transient response data only. In this context, a guideline is proposed for estimating the order of the rational function approximation from transient data. The availability of such an estimate expedites the time-domain rational function fitting process. The second approach relies on the extraction of the delay associated with causal electromagnetic responses of interconnect systems to provide for a more stable rational function process utilizing a lower-order rational function interpolation. A distinctive feature of the proposed methodology is its utilization of scattering parameters. For both methodologies, the approach of fitting the electromagnetic network matrix one element at a time is applied. It is shown that, with regard to the computational cost of the rational function fitting process, such an element-by-element rational function fitting is more advantageous than full matrix fitting for systems with a large number of ports. Despite the disadvantage that different sets of poles are used in the rational function of different elements in the network matrix, such an approach provides for improved accuracy in the fitting of network matrices of systems characterized by both strongly coupled and weakly coupled ports. Finally, in order to provide a means for enforcing passivity in the adopted element-by-element rational function fitting approach, the methodology for passivity enforcement via quadratic programming is modified appropriately for this purpose and demonstrated in the context of element-by-element rational function fitting of the admittance matrix of an electromagnetic multiport.

The spectral albedo of sea ice and salt crusts on the tropical ocean of Snowball Earth: II. Optical modeling

During the Snowball Earth events of the Neoproterozoic, tropical regions of the ocean could have developed a precipitated salt lag deposit left behind by sublimating sea ice. The major salt would have been hydrohalite, NaCl•2H2O. The crystals in such a deposit can be small and highly scattering, resulting in an allwave albedo similar to that of snow. The snow-free sea ice from which such a crust could develop has a lower albedo, around 0.5, so the development of a crust would substantially increase the albedo of tropical regions on Snowball Earth. Hydrohalite crystals are much less absorptive than ice in the near- infrared part of the solar spectrum, so their presence at the surface would increase the overall albedo as well as altering its spectral distribution. In this paper, we use laboratory measurements of the spectral albedo of a hydrohalite lag deposit, in combination with a radiative transfer model, to infer the inherent optical properties of hydrohalite as functions of wavelength. Using this result, we model mixtures of hydrohalite and ice representing both artificially created surfaces in the laboratory and surfaces relevant to Snowball Earth. The model is tested against sequences of laboratory measurements taken during the formation and the dissolution of a lag deposit of hydrohalite. We present a parameterization for the broadband albedo of cold, sublimating sea ice as it forms and evolves a hydrohalite crust, for use in climate models of Snowball Earth.

Micropolar flow in a porous channel with high mass transfer

Two dimensional flow of a micropolar fluid in a porous channel is investigated. The flow is driven by suction or injection at the channel walls, and the micropolar model due to Eringen is used to describe the working fluid. An extension of Berman's similarity transform is used to reduce the governing equations to a set of non-linear coupled ordinary differential equations. The latter are solved for large mass transfer via a perturbation analysis where the inverse of the cross-flow Reynolds number is used as the perturbing parameter. Complementary numerical solutions for strong injection are also obtained using a quasilinearisation scheme, and good agreement is observed between the solutions obtained from the perturbation analysis and the computations.