Dynamics of nonlinear error growth and season-dependent predictability of El Nino events in the Zebiak-Cane model

With the intermediate-complexity Zebiak-Cane model, we investigate the 'spring predictability barrier' (SPB) problem for El Nino events by tracing the evolution of conditional nonlinear optimal perturbation (CNOP), where CNOP is superimposed on the El Nino events and acts as the initial error with the biggest negative effect on the El Nino prediction. We show that the evolution of CNOP-type errors has obvious seasonal dependence and yields a significant SPB, with the most severe occurring in predictions made before the boreal spring in the growth phase of El Nino. The CNOP-type errors can be classified into two types: one possessing a sea-surface-temperature anomaly pattern with negative anomalies in the equatorial central-western Pacific, positive anomalies in the equatorial eastern Pacific, and a thermocline depth anomaly pattern with positive anomalies along the Equator, and another with patterns almost opposite to those of the former type. In predictions through the spring in the growth phase of El Nino, the initial error with the worst effect on the prediction tends to be the latter type of CNOP error, whereas in predictions through the spring in the decaying phase, the initial error with the biggest negative effect on the prediction is inclined to be the former type of CNOP error. Although the linear singular vector (LSV)-type errors also have patterns similar to the CNOP-type errors, they cover a more localized area than the CNOP-type errors and cause a much smaller prediction error, yielding a less significant SPB. Random errors in the initial conditions are also superimposed on El Nino events to investigate the SPB. We find that, whenever the predictions start, the random errors neither exhibit an obvious season-dependent evolution nor yield a large prediction error, and thus may not be responsible for the SPB phenomenon for El Nino events. These results suggest that the occurrence of the SPB is closely related to particular initial error patterns. The two kinds of CNOP-type error are most likely to cause a significant SPB. They have opposite signs and, consequently, opposite growth behaviours, a result which may demonstrate two dynamical mechanisms of error growth related to SPB: in one case, the errors grow in a manner similar to El Nino; in the other, the errors develop with a tendency opposite to El Nino. The two types of CNOP error may be most likely to provide the information regarding the 'sensitive area' of El Nino-Southern Oscillation (ENSO) predictions. If these types of initial error exist in realistic ENSO predictions and if a target method or a data assimilation approach can filter them, the ENSO forecast skill may be improved. Copyright (C) 2009 Royal Meteorological Society

Density-dependent effects on seston dynamics and rates of filtering and biodeposition of the suspension-cultured scallop Chlamys farreri in a eutrophic bay (northern China): An experimental study in semi-in situ flow-through systems

Effects of stocking density on seston dynamics and filtering and biodeposition by the suspension-cultured Zhikong scallop Chlamys farreri Jones et Preston in a eutrophic bay (Sishili Bay, northern China), were determined in a 3-month semi-field experiment with continuous flow-through seawater from the bay. Results showed that the presence of the scallops could strongly decrease seston and chlorophyll a concentrations in the water column. Moreover, in a limited water column, increasing scallop density could cause seston depletion due to scallop's filtering and biodeposition process, and impair scallop growth. Both filtration rate and biodeposition rate of C. farreri showed significant negative correlation with their density and positive relationship with seston concentration. Calculation predicts that the daily removal of suspended matter from water column by the scallops in Sishili Bay ecosystem can be as high as 45% of the total suspended matter; and the daily production of biodeposits by the scallops in early summer in farming zone may amount to 7.78 g m(-2), with daily C, N and P biodeposition rates of 3.06 x 10(-1), 3.86 x 10(-2) and 9.80 x 10(-3) g m(-2), respectively. The filtering and biodeposition by suspension-cultured scallops could substantially enhance the deposition of total suspended particulate material, suppress accumulation of particulate organic matter in water column, and increase the flux of C, N and P to benthos, strongly enhancing pelagic-benthic coupling. It was suggested that the filtering-biodeposition process by intensively suspension-cultured bivalve filter-feeders could exert strong top-down control on phytoplankton biomass and other suspended particulate material in coastal ecosystems. This study also indicated that commercially suspension-cultured bivalves may simultaneously and potentially aid in mitigating eutrophication pressures on coastal ecosystems subject to anthropogenic N and P loadings, serving as a eutrophic-environment bioremediator. The ecological services (e.g. filtering capacity, top-down control, and benthic-pelagic coupling) functioned by extractive bivalve aquaculture should be emphasized in coastal ecosystems. (c) 2005 Elsevier B.V. All rights reserved.

A Molecular Dynamics Simulation: the Effect of Finite Size on the Thermal Conductivity in a Single Crystal Silicon

Non-equilibrium molecular dynamics (NEMD) simulations are performed to calculate thermal conductivity. The environment-dependent interatomic potential (EDIP) potential on crystal silicon is adopted as a model system. The issues are related to nonlinear response, local thermal equilibrium and statistical averaging. The simulation results by non-equilibrium molecular dynamics show that the calculated thermal conductivity decreases almost linearly as the film thickness reduced at the nanometre scale. The effect of size on the thermal conductivity is also obtained by a theoretic analysis of the kinetic theory and formulas of the heat capacity. The analysis reveals that the contributions of phonon mean free path (MFP) and phonon number in a finite cell to thermal conductivity are very important.

Atomistic Origin Of Rate-Dependent Serrated Plastic Flow In Metallic Glasses

Nanoindentation simulations on a binary metallic glass were performed under various strain rates by using molecular dynamics. The rate-dependent serrated plastic flow was clearly observed, and the spatiotemporal behavior of its underlying irreversible atomic rearrangement was probed. Our findings clearly validate that the serration is a temporally inhomogeneous characteristic of such rearrangements and not directly dependent on the resultant shear-banding spatiality. The unique spatiotemporal distribution of shear banding during nanoindentation is highlighted in terms of the potential energy landscape (PEL) theory.

Size-Dependent Interface Phonon Transmission And Thermal Conductivity Of Nanolaminates

An analytical model for size-dependent interface phonon transmission and thermal conductivity of nanolaminates is derived based on the improved acoustic mismatch theory and the Lindemann melting theory by considering the size effect of phonon velocity and the interface lattice mismatch effect. The model suggests that the interface phonon transmission is dominant for the cross-plane thermal conductivity of nanolaminates and superlattices, and the intrinsic variety of size effect of thermal conductivity for different systems is proposed based on the competition mechanism of size effect of phonon transport between two materials constituting the interfaces. The model's prediction for thermal conductivity of nanolaminates agrees with the experimental results. (C) 2008 American Institute of Physics.

Molecular dynamics simulation of crack-tip processes in copper

The crack tip processes in copper under mode II loading have been simulated by a molecular dynamics method. The nucleation, emission, dislocation free zone (DFZ) and pile-up of the dislocations are analyzed by using a suitable atom lattice configuration and Finnis & Sinclair potential. The simulated results show that the dislocation emitted always exhibits a dissociated fashion. The stress intensity factor for dislocation nucleation, DFZ and dissociated width of partial dislocations are strongly dependent on the loading rate. The stress distributions are in agreement with the elasticity solution before the dislocation emission, but are not in agreement after the emission. The dislocation can move at subsonic wave speed (less than the shear wave speed) or at transonic speed (greater than the shear wave speed but less than the longitudinal wave speed), but at the longitudinal wave speed the atom lattice breaks down.

Atomistic Simulation On Size-Dependent Yield Strength And Defects Evolution Of Metal Nanowires

A simple derivation based on continuum mechanics is given, which shows the surface stress is critical for yield strength at ultra-small scales. Molecular dynamics (MD) simulations with modified embedded atom method (MEAM) are employed to investigate the mechanical behaviors of single-crystalline metal nanowires under tensile loading. The calculated yield strengths increasing with the decrease of the cross-sectional area of the nanowires are in accordance with the theoretical prediction. Reorientation induced by stacking faults is observed at the nanowire edge. In addition. the mechanism of yielding is discussed in details based on the snapshots of defects evolution. The nanowires in different crystallographic orientations behave differently in stretching deformation. This study on the plastic properties of metal nanowires will be helpful to further understanding of the mechanical properties of nanomaterials. (C) 2009 Elsevier B.V. All rights reserved.

Shape effects on the yield stress and deformation of silicon nanowires: A molecular dynamics simulation

The tension and compression of single-crystalline silicon nanowires (SiNWs) with different cross-sectional shapes are studied systematically using molecular dynamics simulation. The shape effects on the yield stresses are characterized. For the same surface to volume ratio, the circular cross-sectional SiNWs are stronger than the square cross-sectional ones under tensile loading, but reverse happens in compressive loading. With the atoms colored by least-squares atomic local shear strain, the deformation processes reveal that the failure modes of incipient yielding are dependent on the loading directions. The SiNWs under tensile loading slip in {111} surfaces, while the compressive loading leads the SiNWs to slip in the {110} surfaces. The present results are expected to contribute to the design of the silicon devices in nanosystems.

Multi-scale Calculation of Settling Speed of Coarse Particles by Accelerated Stokesian Dynamics without Adjustable Parameter

The calculation of settling speed of coarse particles is firstly addressed, with accelerated Stokesian dynamics without adjustable parameters, in which far field force acting on the particle instead of particle velocity is chosen as dependent variables to consider inter-particle hydrodynamic interactions. The sedimentation of a simple cubic array of spherical particles is simulated and compared to the results available to verify and validate the numerical code and computational scheme. The improvedmethod keeps the same computational cost of the order O(N log N) as usual accelerated Stokesian dynamics does. Then, more realistic random suspension sedimentation is investigated with the help ofMont Carlo method. The computational results agree well with experimental fitting. Finally, the sedimentation of finer cohesive particle, which is often observed in estuary environment, is presented as a further application in coastal engineering.

The dynamics of hippocampal sensory gating during the development of morphine dependence and withdrawal in rats

The effects of morphine on hippocampal sensory gating (N40) during the development of morphine dependence and withdrawal were investigated in the double click auditory evoked potential (EP) suppression paradigm. Rats were made dependent upon morphine hydrochloride by a series of injections (every 12h) over 6 days, followed by withdrawal after stopping morphine administration. Hippocampal gating was examined during the development of dependence and withdrawal. Moreover, the DA antagonist haloperidol was used to assess the contribution of dopamine to hippocampal gating induced by morphine. Our results showed that the morphine-treated rats exhibited significantly disrupted hippocampal gating during the development of morphine dependence and this disrupted gating was partially reversed by haloperidol pretreatment. In contrast, there was significantly enhanced hippocampal gating at the fifth and sixth days of withdrawal. The dynamics of hippocampal gating during the development of morphine dependence and withdrawal suggests the interaction between the hippocampus and opioids. (c) 2005 Elsevier Ireland Ltd. All rights reserved.

Prediction and elucidation of the population dynamics of Microcystis spp. in Lake Dianchi (China) by means of artificial neural networks

Lake Dianchi is a shallow and turbid lake, located in Southwest China. Since 1985, Lake Dianchi has experienced severe cyanabacterial blooms (dominated by Microcystis spp.). In extreme cases, the algal cell densities have exceeded three billion cells per liter. To predict and elucidate the population dynamics ofMicrocystis spp. in Lake Dianchi, a neural network based model was developed. The correlation coefficient (R 2) between the predicted algal concentrations by the model and the observed values was 0.911. Sensitivity analysis was performed to clarify the algal dynamics to the changes of environmental factors. The results of a sensitivity analysis of the neural network model suggested that small increases in pH could cause significantly reduced algal abundance. Further investigations on raw data showed that the response of Microcystis spp. concentration to pH increase was dependent on algal biomass and pH level. When Microcystis spp. population and pH were moderate or low, the response of Microcystis spp. population would be more likely to be positive in Lake Dianchi; contrarily, Microcystis spp. population in Lake Dianchi would be more likely to show negative response to pH increase when Microcystis spp. population and pH were high. The paper concluded that the extremely high concentration of algal population and high pH could explain the distinctive response of Microcystis spp. population to +1 SD (standard deviation) pH increase in Lake Dianchi. And the paper also elucidated the algal dynamics to changes of other environmental factors. One SD increase of water temperature (WT) had strongest positive relationship with Microcystis spp. biomass. Chemical oxygen demand (COD) and total phosphorus (TP) had strong positive effect on Microcystis spp. abundance while total nitrogen (TN), biological oxygen demand in five days (BOD5), and dissolved oxygen had only weak relationship with Microcystis spp. concentration. And transparency (Tr) had moderate positive relationship with Microcystis spp. concentration.

Coherence loss and recovery of an electron spin coupled inhomogeneously to a one-dimensional interacting spin bath: An adaptive time-dependent density-matrix renormalization group study

Coherence evolution and echo effect of an electron spin, which is coupled inhomogeneously to an interacting one-dimensional finite spin bath via hyperfine-type interaction, are studied using the adaptive time-dependent density-matrix renormalization group method. It is found that the interplay of the coupling inhomogeneity and the transverse intrabath interactions results in two qualitatively different coherence evolutions, namely, a coherence-preserving evolution characterized by periodic oscillation and a complete decoherence evolution. Correspondingly, the echo effects induced by an electron-spin flip at time tau exhibit stable recoherence pulse sequence for the periodic evolution and a single peak at root 2 tau for the decoherence evolution, respectively. With the diagonal intrabath interaction included, the specific feature of the periodic regime is kept, while the root 2 tau-type echo effect in the decoherence regime is significantly affected. To render the experimental verifications possible, the Hahn echo envelope as a function of tau is calculated, which eliminates the inhomogeneous broadening effect and serves for the identification of the different status of the dynamic coherence evolution, periodic versus decoherence.

High-temperature electron-hole liquid and dynamics of Fermi excitons in a In0.65Al0.35As/Al0.4Ga0.6As quantum dot array

State-filling effects of the exciton in a In0.65Al0.35As/Al0.4Ga0.6As quantum dot array are observed by quantum dot array photolumineseence at a sample temperature of 77 K. The exciton emission at low excitation density is dominated by the radiative recombination of the states in the s shell and at high excitation density the emission mainly results from the radiative recombination of the exciton state in the p shell. The spectral interval between the states in the s and p shells is about 30-40 mcV. The time resolved photoluminescence shows that the decay time of exciton states in the p shell is longer than that of exciton states in the s shell, and the emission intensity of the exciton state in the p shell is superlinearly dependent on excitation density. Furthermore, electron-hole liquid in the quantum dot array is observed at 77 K, which is a much higher temperature than that in bulk. The emission peak of the. recombination, of electron-hole liquid has an about 200 meV redshift from the exciton fluorescence. Two excitation density-dependent emission peaks at 1.56 and 1.59 eV are observed, respectively, which result from quantum confinement effects in QDs. The emission intensity of electron-hole liquid is directly proportional to the cubic of excitation densities and its decay time decreases significantly at the high excitation density.

Depth dependent elastic strain in ZnO epilayer: combined Rutherford backscattering/channeling and X-ray diffraction

A ZnO layer was grown by metalorganic chemical vapor deposition (MOCVD) on a sapphire (0 0 0 1) substrate. The perpendicular and parallel elastic strain of the ZnO epilayer, e(perpendicular to) = 0.19%, e(parallel to) = -0.29%, respectively, were derived by using the combination of Rutherford backscattering (RBS)/channeling and X-ray diffraction (XRD). The ratio vertical bar e(parallel to)/ e(perpendicular to)vertical bar = 1.5 indicates that ZnO layer is much stiffer in the a-axis direction than in the c-axis direction. By using RBS/C, the depth dependent elastic strain was deduced. The strain is higher at the depth close to the interface and decreases towards the surface. The negative tetragonal distortion was explained by considering the lattice mismatch and thermal mismatch in ZnO thin film. (c) 2004 Elsevier B.V. All rights reserved.