Evidence of dispersion relations for the nonlinear response of the Lorenz 63 system

Along the lines of the nonlinear response theory developed by Ruelle, in a previous paper we have proved under rather general conditions that Kramers-Kronig dispersion relations and sum rules apply for a class of susceptibilities describing at any order of perturbation the response of Axiom A non equilibrium steady state systems to weak monochromatic forcings. We present here the first evidence of the validity of these integral relations for the linear and the second harmonic response for the perturbed Lorenz 63 system, by showing that numerical simulations agree up to high degree of accuracy with the theoretical predictions. Some new theoretical results, showing how to derive asymptotic behaviors and how to obtain recursively harmonic generation susceptibilities for general observables, are also presented. Our findings confirm the conceptual validity of the nonlinear response theory, suggest that the theory can be extended for more general non equilibrium steady state systems, and shed new light on the applicability of very general tools, based only upon the principle of causality, for diagnosing the behavior of perturbed chaotic systems and reconstructing their output signals, in situations where the fluctuation-dissipation relation is not of great help.

Precise calculations of the existence of multiple AMOC equilibria in coupled climate models Part I: equilibrium states

This study examines criteria for the existence of two stable states of the Atlantic Meridional Overturning Circulation (AMOC) using a combination of theory and simulations from a numerical coupled atmosphere–ocean climate model. By formulating a simple collection of state parameters and their relationships, the authors reconstruct the North Atlantic Deep Water (NADW) OFF state behavior under a varying external salt-flux forcing. This part (Part I) of the paper examines the steady-state solution, which gives insight into the mechanisms that sustain the NADW OFF state in this coupled model; Part II deals with the transient behavior predicted by the evolution equation. The nonlinear behavior of the Antarctic Intermediate Water (AAIW) reverse cell is critical to the OFF state. Higher Atlantic salinity leads both to a reduced AAIW reverse cell and to a greater vertical salinity gradient in the South Atlantic. The former tends to reduce Atlantic salt export to the Southern Ocean, while the latter tends to increases it. These competing effects produce a nonlinear response of Atlantic salinity and salt export to salt forcing, and the existence of maxima in these quantities. Thus the authors obtain a natural and accurate analytical saddle-node condition for the maximal surface salt flux for which a NADW OFF state exists. By contrast, the bistability indicator proposed by De Vries and Weber does not generally work in this model. It is applicable only when the effect of the AAIW reverse cell on the Atlantic salt budget is weak.

The influence of ozone precursor emissions from four world regions on tropospheric composition and radiative climate forcing

Ozone (O3) precursor emissions influence regional and global climate and air quality through changes in tropospheric O3 and oxidants, which also influence methane (CH4) and sulfate aerosols (SO42−). We examine changes in the tropospheric composition of O3, CH4, SO42− and global net radiative forcing (RF) for 20% reductions in global CH4 burden and in anthropogenic O3 precursor emissions (NOx, NMVOC, and CO) from four regions (East Asia, Europe and Northern Africa, North America, and South Asia) using the Task Force on Hemispheric Transport of Air Pollution Source-Receptor global chemical transport model (CTM) simulations, assessing uncertainty (mean ± 1 standard deviation) across multiple CTMs. We evaluate steady state O3 responses, including long-term feedbacks via CH4. With a radiative transfer model that includes greenhouse gases and the aerosol direct effect, we find that regional NOx reductions produce global, annually averaged positive net RFs (0.2 ± 0.6 to 1.7 ± 2 mWm−2/Tg N yr−1), with some variation among models. Negative net RFs result from reductions in global CH4 (−162.6 ± 2 mWm−2 for a change from 1760 to 1408 ppbv CH4) and regional NMVOC (−0.4 ± 0.2 to −0.7 ± 0.2 mWm−2/Tg C yr−1) and CO emissions (−0.13 ± 0.02 to −0.15 ± 0.02 mWm−2/Tg CO yr−1). Including the effect of O3 on CO2 uptake by vegetation likely makes these net RFs more negative by −1.9 to −5.2 mWm−2/Tg N yr−1, −0.2 to −0.7 mWm−2/Tg C yr−1, and −0.02 to −0.05 mWm−2/Tg CO yr−1. Net RF impacts reflect the distribution of concentration changes, where RF is affected locally by changes in SO42−, regionally to hemispherically by O3, and globally by CH4. Global annual average SO42− responses to oxidant changes range from 0.4 ± 2.6 to −1.9 ± 1.3 Gg for NOx reductions, 0.1 ± 1.2 to −0.9 ± 0.8 Gg for NMVOC reductions, and −0.09 ± 0.5 to −0.9 ± 0.8 Gg for CO reductions, suggesting additional research is needed. The 100-year global warming potentials (GWP100) are calculated for the global CH4 reduction (20.9 ± 3.7 without stratospheric O3 or water vapor, 24.2 ± 4.2 including those components), and for the regional NOx, NMVOC, and CO reductions (−18.7 ± 25.9 to −1.9 ± 8.7 for NOx, 4.8 ± 1.7 to 8.3 ± 1.9 for NMVOC, and 1.5 ± 0.4 to 1.7 ± 0.5 for CO). Variation in GWP100 for NOx, NMVOC, and CO suggests that regionally specific GWPs may be necessary and could support the inclusion of O3 precursors in future policies that address air quality and climate change simultaneously. Both global net RF and GWP100 are more sensitive to NOx and NMVOC reductions from South Asia than the other three regions.

Modeling electrocortical activity through improved local approximations of integral neural field equations

Neural field models of firing rate activity typically take the form of integral equations with space-dependent axonal delays. Under natural assumptions on the synaptic connectivity we show how one can derive an equivalent partial differential equation (PDE) model that properly treats the axonal delay terms of the integral formulation. Our analysis avoids the so-called long-wavelength approximation that has previously been used to formulate PDE models for neural activity in two spatial dimensions. Direct numerical simulations of this PDE model show instabilities of the homogeneous steady state that are in full agreement with a Turing instability analysis of the original integral model. We discuss the benefits of such a local model and its usefulness in modeling electrocortical activity. In particular, we are able to treat “patchy” connections, whereby a homogeneous and isotropic system is modulated in a spatially periodic fashion. In this case the emergence of a “lattice-directed” traveling wave predicted by a linear instability analysis is confirmed by the numerical simulation of an appropriate set of coupled PDEs.

A model of the hemodynamic response and oxygen delivery to brain

A recent nonlinear system by Friston et al. (2000. NeuroImage 12: 466–477) links the changes in BOLD response to changes in neural activity. The system consists of five subsystems, linking: (1) neural activity to flow changes; (2) flow changes to oxygen delivery to tissue; (3) flow changes to changes in blood volume and venous outflow; (4) changes in flow, volume, and oxygen extraction fraction to deoxyhemoglobin changes; and finally (5) volume and deoxyhemoglobin changes to the BOLD response. Friston et al. exploit, in subsystem 2, a model by Buxton and Frank coupling flow changes to changes in oxygen metabolism which assumes tissue oxygen concentration to be close to zero. We describe below a model of the coupling between flow and oxygen delivery which takes into account the modulatory effect of changes in tissue oxygen concentration. The major development has been to extend the original Buxton and Frank model for oxygen transport to a full dynamic capillary model making the model applicable to both transient and steady state conditions. Furthermore our modification enables us to determine the time series of CMRO2 changes under different conditions, including CO2 challenges. We compare the differences in the performance of the “Friston system” using the original model of Buxton and Frank and that of our model. We also compare the data predicted by our model (with appropriate parameters) to data from a series of OIS studies. The qualitative differences in the behaviour of the models are exposed by different experimental simulations and by comparison with the results of OIS data from brief and extended stimulation protocols and from experiments using hypercapnia.

Impact of oceanic processes on the carbon cycle during the last termination

During the last termination (from ~18 000 years ago to ~9000 years ago), the climate significantly warmed and the ice sheets melted. Simultaneously, atmospheric CO2 increased from ~190 ppm to ~260 ppm. Although this CO2 rise plays an important role in the deglacial warming, the reasons for its evolution are difficult to explain. Only box models have been used to run transient simulations of this carbon cycle transition, but by forcing the model with data constrained scenarios of the evolution of temperature, sea level, sea ice, NADW formation, Southern Ocean vertical mixing and biological carbon pump. More complex models (including GCMs) have investigated some of these mechanisms but they have only been used to try and explain LGM versus present day steady-state climates. In this study we use a coupled climate-carbon model of intermediate complexity to explore the role of three oceanic processes in transient simulations: the sinking of brines, stratification-dependent diffusion and iron fertilization. Carbonate compensation is accounted for in these simulations. We show that neither iron fertilization nor the sinking of brines alone can account for the evolution of CO2, and that only the combination of the sinking of brines and interactive diffusion can simultaneously simulate the increase in deep Southern Ocean δ13C. The scenario that agrees best with the data takes into account all mechanisms and favours a rapid cessation of the sinking of brines around 18 000 years ago, when the Antarctic ice sheet extent was at its maximum. In this scenario, we make the hypothesis that sea ice formation was then shifted to the open ocean where the salty water is quickly mixed with fresher water, which prevents deep sinking of salty water and therefore breaks down the deep stratification and releases carbon from the abyss. Based on this scenario, it is possible to simulate both the amplitude and timing of the long-term CO2 increase during the last termination in agreement with ice core data. The atmospheric δ13C appears to be highly sensitive to changes in the terrestrial biosphere, underlining the need to better constrain the vegetation evolution during the termination.

Impact of atmospheric forcing on Antarctic continental shelf water masses

The Antarctic continental shelf seas feature a bimodal distribution of water mass temperature, with the Amundsen and Bellingshausen Seas flooded by Circumpolar Deep Water that is several degrees Celsius warmer than the cold shelf waters prevalent in the Weddell and Ross Seas. This bimodal distribution could be caused by differences in atmospheric forcing, ocean dynamics, ocean and ice feedbacks, or some combination of these factors. In this study, a highly simplified coupled sea ice–mixed layer model is developed to investigate the physical processes controlling this situation. Under regional atmospheric forcings and parameter choices the 10-yr simulations demonstrate a complete destratification of the Weddell Sea water column in winter, forming cold, relatively saline shelf waters, while the Amundsen Sea winter mixed layer remains shallower, allowing a layer of deep warm water to persist. Applying the Weddell atmospheric forcing to the Amundsen Sea model destratifies the water column after two years, and applying the Amundsen forcing to the Weddell Sea model results in a shallower steady-state winter mixed layer that no longer destratifies the water column. This suggests that the regional difference in atmospheric forcings alone is sufficient to account for the bimodal distribution in Antarctic shelf-sea temperatures. The model prediction of mixed layer depth is most sensitive to the air temperature forcing, but a switch in all forcings is required to prevent destratification of the Weddell Sea water column.

Coalescence of particles by differential sedimentation

We consider a three dimensional system consisting of a large number of small spherical particles, distributed in a range of sizes and heights (with uniform distribution in the horizontal direction). Particles move vertically at a size-dependent terminal velocity. They are either allowed to merge whenever they cross or there is a size ratio criterion enforced to account for collision efficiency. Such a system may be described, in mean field approximation, by the Smoluchowski kinetic equation with a differential sedimentation kernel. We obtain self-similar steady-state and time-dependent solutions to the kinetic equation, using methods borrowed from weak turbulence theory. Analytical results are compared with direct numerical simulations (DNS) of moving and merging particles, and a good agreement is found.

Intercomparison of general circulation models for hot extrasolar planets

We compare five general circulation models (GCMs) which have been recently used to study hot extrasolar planet atmospheres (BOB, CAM, IGCM, MITgcm, and PEQMOD), under three test cases useful for assessing model convergence and accuracy. Such a broad, detailed intercomparison has not been performed thus far for extrasolar planets study. The models considered all solve the traditional primitive equations, but employ di↵erent numerical algorithms or grids (e.g., pseudospectral and finite volume, with the latter separately in longitude-latitude and ‘cubed-sphere’ grids). The test cases are chosen to cleanly address specific aspects of the behaviors typically reported in hot extrasolar planet simulations: 1) steady-state, 2) nonlinearly evolving baroclinic wave, and 3) response to fast timescale thermal relaxation. When initialized with a steady jet, all models maintain the steadiness, as they should—except MITgcm in cubed-sphere grid. A very good agreement is obtained for a baroclinic wave evolving from an initial instability in pseudospectral models (only). However, exact numerical convergence is still not achieved across the pseudospectral models: amplitudes and phases are observably di↵erent. When subject to a typical ‘hot-Jupiter’-like forcing, all five models show quantitatively di↵erent behavior—although qualitatively similar, time-variable, quadrupole-dominated flows are produced. Hence, as have been advocated in several past studies, specific quantitative predictions (such as the location of large vortices and hot regions) by GCMs should be viewed with caution. Overall, in the tests considered here, pseudospectral models in pressure coordinate (PEBOB and PEQMOD) perform the best and MITgcm in cubed-sphere grid performs the worst.

Cropping systems, carbon sequestration and erosion in Brazil, a review

Soils represent a large carbon pool, approximately 1500 Gt, which is equivalent to almost three times the quantity stored in terrestrial biomass and twice the amount stored in the atmosphere. Any modification of land use or land management can induce variations in soil carbon stocks, even in agricultural systems that are perceived to be in a steady state. Tillage practices often induce soil aerobic conditions that are favourable to microbial activity and may lead to a degradation of soil structure. As a result, mineralisation of soil organic matter increases in the long term. The adoption of no-tillage systems and the maintenance of a permanent vegetation cover using Direct seeding Mulch-based Cropping system or DMC, may increase carbon levels in the topsoil. In Brazil, no-tillage practices (mainly DMC), were introduced approximately 30 years ago in the south in the Parana state, primarily as a means of reducing erosion. Subsequently, research has begun to study the management of the crop waste products and their effects on soil fertility, either in terms of phosphorus management, as a means of controlling soil acidity, or determining how manures can be applied in a more localised manner. The spread of no-till in Brazil has involved a large amount of extension work. The area under no-tillage is still increasing in the centre and north of the country and currently occupies ca. 20 million hectares, covering a diversity of environmental conditions, cropping systems and management practices. Most studies of Brazilian soils give rates of carbon storage in the top 40 cm of the soil of 0.4 to 1.7 t C ha(-1) per year, with the highest rates in the Cerrado region. However, caution must be taken when analysing DMC systems in terms of carbon sequestration. Comparisons should include changes in trace gas fluxes and should not be limited to a consideration of carbon storage in the soil alone if the full implications for global warming are to be assessed.

The influence of Fe-rich coatings on the dissolution of anorthite at pH 2.6

Surface coatings are very common on mineral grains in soils but most laboratory dissolution experiments are carried out on pristine, uncoated mineral grains. An experiment designed to unambiguously isolate the effect of surface coatings on mineral dissolution from any influence of solution saturation state is reported. Two aliquots of 53 to 63 mum anorthite feldspar powder were used. One was dissolved in pH 2.6 HCl, the other in pH 2.6 FeCl3 solution, both for similar to6000 h in flow-through reactors. An amorphous Fe-rich, Al-, Ca- and Si-free orange precipitate coated the anorthite dissolved in the FeCl3 solution. BET surface area of the anorthite increased from 0.16 to 1.65 m(2) g(-1) in the HCl experiment and to 3.89 m(2) g(-1) in the FeCl3 experiment. The increase in surface area in the HCl experiment was due to the formation of etch pits on the anorthite grain surface whilst the additional increase in the FeCl3 experiment was due to the micro- and meso-porous nature of the orange precipitate. This precipitate did not inhibit or slow the dissolution of the anorthite. Steady state dissolution rates for the anorthite dissolved in the HCl and FeCl3 were similar to2.5 and 3.2 X 10(-10) mol(feldspar) m(-2) s(-1) respectively. These rates are not significantly different after the cumulative uncertainty of 17% in their value due to uncertainty in the inputs parameters used in their calculation is taken into account. Results from this experiment support previous theoretical and inference-based conclusions that porous coatings should not inhibit mineral dissolution. Copyright (C) 2003 Elsevier Ltd.

Does reactive surface area depend on grain size? Results from pH 3, 25 circle C far-from-equilibrium flow-through dissolution experiments on anorthite and biotite

Laboratory determined mineral weathering rates need to be normalised to allow their extrapolation to natural systems. The principle normalisation terms used in the literature are mass, and geometric- and BET specific surface area (SSA). The purpose of this study was to determine how dissolution rates normalised to these terms vary with grain size. Different size fractions of anorthite and biotite ranging from 180-150 to 20-10 mu m were dissolved in pH 3, HCl at 25 degrees C in flow through reactors under far from equilibrium conditions. Steady state dissolution rates after 5376 h (anorthite) and 4992 h (biotite) were calculated from Si concentrations and were normalised to initial- and final- mass and geometric-, geometric edge- (biotite), and BET SSA. For anorthite, rates normalised to initial- and final-BET SSA ranged from 0.33 to 2.77 X 10(-10) mol(feldspar) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 5.74 to 8.88 X 10(-10) mol(feldspar) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.11 to 1.65 mol(feldspar) g(-1) s(-1). For biotite, rates normalised to initial- and final-BET SSA ranged from 1.02 to 2.03 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric SSA ranged from 3.26 to 16.21 X 10(-12) mol(biotite) m(-2) s(-1), rates normalised to initial- and final-geometric edge SSA ranged from 59.46 to 111.32 x 10(-12) mol(biotite) m(-2) s(-1) and rates normalised to initial- and final-mass ranged from 0.81 to 6.93 X 10(-12) mol(biotite) g(-1) s(-1). For all normalising terms rates varied significantly (p <= 0.05) with grain size. The normalising terms which gave least variation in dissolution rate between grain sizes for anorthite were initial BET SSA and initial- and final-geometric SSA. This is consistent with: (1) dissolution being dominated by the slower dissolving but area dominant non-etched surfaces of the grains and, (2) the walls of etch pits and other dissolution features being relatively unreactive. These steady state normalised dissolution rates are likely to be constant with time. Normalisation to final BET SSA did not give constant ratios across grain size due to a non-uniform distribution of dissolution features. After dissolution coarser grains had a greater density of dissolution features with BET-measurable but unreactive wall surface area than the finer grains. The normalising term which gave the least variation in dissolution rates between grain sizes for biotite was initial BET SSA. Initial- and final-geometric edge SSA and final BET SSA gave the next least varied rates. The basal surfaces dissolved sufficiently rapidly to influence bulk dissolution rate and prevent geometric edge SSA normalised dissolution rates showing the least variation. Simple modelling indicated that biotite grain edges dissolved 71-132 times faster than basal surfaces. In this experiment, initial BET SSA best integrated the different areas and reactivities of the edge and basal surfaces of biotite. Steady state dissolution rates are likely to vary with time as dissolution alters the ratio of edge to basal surface area. Therefore they would be more properly termed pseudo-steady state rates, only appearing constant because the time period over which they were measured (1512 h) was less than the time period over wich they would change significantly. (c) 2006 Elsevier Inc. All rights reserved.

Potassium leaching in undisturbed soil cores following surface applications of gypsum

Displacement studies on leaching of potassium (K+) were conducted under unsaturated steady state flow conditions in nine undisturbed soil columns (15.5 cm in diameter and 25 cm long). Pulses of K+ applied to columns of undisturbed soil were leached with distilled water or calcium chloride (CaCl2) at a rate of 18 mm h(-1). The movement of K+ in gypsum treated soil leached with distilled water was at a similar rate to that of the untreated soil leached with 15 mM CaCl2. The Ca2+ concentrations in the leachates were about 15 mM, the expected values for the dissolution of the gypsum. When applied K+ was displaced with the distilled water, K+ was retained in the top 10-12.5 cm depth of soil. In the undisturbed soil cores there is possibility of preferential flow and lack of K+ sorption. The application of gypsum and CaCl2 in the reclamation of sodic soils would be expected to leach K+ from soils. It can also be concluded that the use of sources of water for irrigation which have a high Ca2+ concentration can also lead to leaching of K+ from soil. Average effluent concentration of K+ during leaching period was 30.2 and 28.6 mg l(-1) for the gypsum and CaCl2 treated soils, respectively. These concentrations are greater than the recommended guideline of the World Health Organisation (12 mg K+ l(-1)).

Uptake kinetics of metals by the earthworm Eisenia fetida exposed to field-contaminated soils

It is well known that earthworms can accumulate metals. However, most accumulation studies focus on Cd-, Cu-, Pb- or Zn-amended soils, additionally few studies consider accumulation kinetics. Here we model the accumulation kinetics of 18 elements by Eisenia fetida, exposed to 8 metal-contaminated and 2 uncontaminated soils. Tissue metal concentration was determined after 3, 7,14, 21, 28 and 42 days. Metal elimination rate was important in determining time to reach steady-state tissue metal concentration. Uptake flux to elimination rate ratios showed less variation and lower values for essential than for non-essential metals. In theory kinetic rate constants are dependent only on species and metal. Therefore it should be possible to predict steady-state tissue metal concentrations on the basis of very few measurements using the rate constants. However, our experiments show that it is difficult to extrapolate the accumulation kinetic constants derived using one soil to another. (C) 2009 Elsevier Ltd. All rights reserved.

The impact of fallowing and green manuring on soil conditions and the growth of sugarcane

There are currently concerns within some sugar industries that long-term monoculture has led to soil degradation and consequent yield decline. An investigation was conducted in Swaziland to assess the effects of fallowing and green manuring practices, over a seven-month period, on sugarcane yields and the physical properties of a poorly draining clay soil. In the subsequent first sugarcane crop after planting, yields were improved from 129 t ha(-1) under continuous sugarcane to 141-144 t ha(-1) after fallowing and green manuring, but there were no significant responses in the first and second ratoon crops. Also, in the first crop after planting, root length index increased from 3.5 km m(-2) under continuous sugarcane to 5.2-6.8 km m(-2) after fallowing, and improved rooting was still evident in the first ratoon crop where there had been soil drying during the fallow period. Soil bulk density, total porosity and water-holding capacity were not affected by the fallowing practices. However, air-filled porosity increased from 11% under continuous sugarcane to 16% after fallowing, and steady state ponded infiltration rates were increased from 0.61 mm h(-1) to 1.34 mm h(-1), but these improvements were no longer evident after a year back under sugarcane. Levels of soil organic matter were reduced in all cases, probably as a result of the tillage operations involved. In the plant crop, root length was well correlated with air-filled porosity, indicating the importance of improving belowground air supply for crop production on poorly draining clay soils.