Response to Comment on "Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply".

Tolstoy reports the existence of a characteristic 100 thousand year (ky) period in the bathymetry of fast-spreading seafloor but does not argue that sea level change is a first-order control on seafloor morphology worldwide. Upon evaluating the overlap between tectonic and Milankovitch periodicities across spreading rates, we reemphasize that fast-spreading ridges are the best potential recorders of a sea level signature in seafloor bathymetry.

Response to Comment on "Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply".

Huybers et al present new bathymetric spectra from an intermediate-spreading ridge as evidence for a primary contribution of sea level cycles to the morphology of the seafloor. Although we acknowledge the possibility that sea level-modulated magmatic constructions may be superimposed on a first-order tectonic fabric, we emphasize the difficulty of deciphering these different contributions in the frequency domain alone.

Sensitivity of seafloor bathymetry to climate-driven fluctuations in mid-ocean ridge magma supply.

Recent studies have proposed that the bathymetric fabric of the seafloor formed at mid-ocean ridges records rapid (23,000 to 100,000 years) fluctuations in ridge magma supply caused by sealevel changes that modulate melt production in the underlying mantle. Using quantitative models of faulting and magma emplacement, we demonstrate that, in fact, seafloor-shaping processes act as a low-pass filter on variations in magma supply, strongly damping fluctuations shorter than about 100,000 years. We show that the systematic decrease in dominant seafloor wavelengths with increasing spreading rate is best explained by a model of fault growth and abandonment under a steady magma input. This provides a robust framework for deciphering the footprint of mantle melting in the fabric of abyssal hills, the most common topographic feature on Earth.

Towards time-dependent current-density-functional theory in the non-linear regime.

Time-Dependent Density-Functional Theory (TDDFT) is a well-established theoretical approach to describe and understand irradiation processes in clusters and molecules. However, within the so-called adiabatic local density approximation (ALDA) to the exchange-correlation (xc) potential, TDDFT can show insufficiencies, particularly in violently dynamical processes. This is because within ALDA the xc potential is instantaneous and is a local functional of the density, which means that this approximation neglects memory effects and long-range effects. A way to go beyond ALDA is to use Time-Dependent Current-Density-Functional Theory (TDCDFT), in which the basic quantity is the current density rather than the density as in TDDFT. This has been shown to offer an adequate account of dissipation in the linear domain when the Vignale-Kohn (VK) functional is used. Here, we go beyond the linear regime and we explore this formulation in the time domain. In this case, the equations become very involved making the computation out of reach; we hence propose an approximation to the VK functional which allows us to calculate the dynamics in real time and at the same time to keep most of the physics described by the VK functional. We apply this formulation to the calculation of the time-dependent dipole moment of Ca, Mg and Na2. Our results show trends similar to what was previously observed in model systems or within linear response. In the non-linear domain, our results show that relaxation times do not decrease with increasing deposited excitation energy, which sets some limitations to the practical use of TDCDFT in such a domain of excitations.

Modeling workplace bullying using catastrophe theory.

Workplace bullying is defined as negative behaviors directed at organizational members or their work context that occur regularly and repeatedly over a period of time. Employees' perceptions of psychosocial safety climate, workplace bullying victimization, and workplace bullying perpetration were assessed within a sample of nearly 5,000 workers. Linear and nonlinear approaches were applied in order to model both continuous and sudden changes in workplace bullying. More specifically, the present study examines whether a nonlinear dynamical systems model (i.e., a cusp catastrophe model) is superior to the linear combination of variables for predicting the effect of psychosocial safety climate and workplace bullying victimization on workplace bullying perpetration. According to the AICc, and BIC indices, the linear regression model fits the data better than the cusp catastrophe model. The study concludes that some phenomena, especially unhealthy behaviors at work (like workplace bullying), may be better studied using linear approaches as opposed to nonlinear dynamical systems models. This can be explained through the healthy variability hypothesis, which argues that positive organizational behavior is likely to present nonlinear behavior, while a decrease in such variability may indicate the occurrence of negative behaviors at work.

On transition rates in surface hopping.

Trajectory surface hopping (TSH) is one of the most widely used quantum-classical algorithms for nonadiabatic molecular dynamics. Despite its empirical effectiveness and popularity, a rigorous derivation of TSH as the classical limit of a combined quantum electron-nuclear dynamics is still missing. In this work, we aim to elucidate the theoretical basis for the widely used hopping rules. Naturally, we concentrate thereby on the formal aspects of the TSH. Using a Gaussian wave packet limit, we derive the transition rates governing the hopping process at a simple avoided level crossing. In this derivation, which gives insight into the physics underlying the hopping process, some essential features of the standard TSH algorithm are retrieved, namely (i) non-zero electronic transition rate ("hopping probability") at avoided crossings; (ii) rescaling of the nuclear velocities to conserve total energy; (iii) electronic transition rates linear in the nonadiabatic coupling vectors. The well-known Landau-Zener model is then used for illustration.

Central role of detachment faults in accretion of slow-spreading oceanic lithosphere.

The formation of oceanic detachment faults is well established from inactive, corrugated fault planes exposed on sea floor formed along ridges spreading at less than 80 km Myr(-1) (refs 1-4). These faults can accommodate extension for up to 1-3 Myr (ref. 5), and are associated with one of the two contrasting modes of accretion operating along the northern Mid-Atlantic Ridge. The first mode is asymmetrical accretion involving an active detachment fault along one ridge flank. The second mode is the well-known symmetrical accretion, dominated by magmatic processes with subsidiary high-angle faulting and the formation of abyssal hills on both flanks. Here we present an examination of approximately 2,500 km of the Mid-Atlantic Ridge between 12.5 and 35 degrees N, which reveals asymmetrical accretion along almost half of the ridge. Hydrothermal activity identified so far in the study region is closely associated with asymmetrical accretion, which also shows high levels of near-continuous hydroacoustically and teleseismically recorded seismicity. Increased seismicity is probably generated along detachment faults that accommodate a sizeable proportion of the total plate separation. In contrast, symmetrical segments have lower levels of seismicity, which occurs primarily at segment ends. Basalts erupted along asymmetrical segments have compositions that are consistent with crystallization at higher pressures than basalts from symmetrical segments, and with lower extents of partial melting of the mantle. Both seismic evidence and geochemical evidence indicate that the axial lithosphere is thicker and colder at asymmetrical sections of the ridge, either because associated hydrothermal circulation efficiently penetrates to greater depths or because the rising mantle is cooler. We suggest that much of the variability in sea-floor morphology, seismicity and basalt chemistry found along slow-spreading ridges can be thus attributed to the frequent involvement of detachment faults in oceanic lithospheric accretion.

A study of dental staining among competitive swimmers.

OBJECTIVE: To estimate the prevalence of dental stains (DS) in competitive swimmers and quantify the risk of these stains compared with sportsmen in a non-swimmers group in Castellón, Spain. METHODS: Cross-sectional and case-control designs. Between July 1996 and March 1997, 404 subjects, (171 enrolled in two clubs of competitive swimming and 233 sportsmen from two schools), were examined in order to detect and classify yellowish-brown or dark-brown stains on the facial surface of the eight incisors. Participation rates were 88.6% for swimmers, and 95.7% for sportsmen. Mean of participants' age was 12 years, range 7-22 years. Castellón has three public competition swimming pools, two of which are indoors. Two of the pools used chloride products, and the third bromine for the disinfection of water. The recommended hygiene regulations were adhered to. RESULTS: Prevalence of DS was 60.2% in swimmers and 12.9% in sportsmen (P= 0.0001). Risk factors for DS included: use of competition swimming pools, age, gender, years of competition, daily consumption of coffee, red wine, and iron supplement during the last year. Professional dental cleanliness was a protective factor. In a logistic regression analysis, the use of competition swimming pools maintained a high risk of DS, odds ratio (OR)=9.28; 95% confidence interval (CI) 5.21-16.5, adjusted by the other variables. Amongst swimmers, more than 6 h of training a week increased the risk of these stains (OR=3.51; 95% CI 1.35-9.10). CONCLUSION: The study indicated a high risk of DS in competitive swimmers.