K-T transition in Deccan traps of central India marks major marine seaway across India.

Deccan intertrappean sediments in central India are generally considered as terrestrial deposits of Maastrichtian age, but the Cretaceous–Tertiary (K–T) position is still unknown. Here we report the discovery of the K–T transition, a marine incursion and environmental changes preserved within the intertrappean sediments at Jhilmili, Chhindwara District, Madhya Pradesh. Integrative biostratigraphic, sedimentologic, mineralogic and chemostratigraphic analyses reveal the basal Danian in the intertrappean sediments between lower and upper trap basalts that regionally correspond to C29r and the C29R/C29N transition, respectively. Intertrappean deposition occurred in predominantly terrestrial semi-humid to arid environments. But a short aquatic interval of fresh water ponds and lakes followed by shallow coastal marine conditions with brackish marine ostracods and early Danian zone P1a planktic foraminifera mark this interval very close to the K–T boundary. This marine incursion marks the existence of a nearby seaway, probably extending inland from the west through the Narmada and Tapti rift valleys. The Jhilmili results thus identify the K–T boundary near the end of the main phase of Deccan eruptions and indicate that a major seaway extended at least 800 km across India.

A closer look at chaotic advection in the stratosphere: part I: geometric structure

The relevance of chaotic advection to stratospheric mixing and transport is addressed in the context of (i) a numerical model of forced shallow-water flow on the sphere, and (ii) a middle-atmosphere general circulation model. It is argued that chaotic advection applies to both these models if there is suitable large-scale spatial structure in the velocity field and if the velocity field is temporally quasi-regular. This spatial structure is manifested in the form of “cat’s eyes” in the surf zone, such as are commonly seen in numerical simulations of Rossby wave critical layers; by analogy with the heteroclinic structure of a temporally aperiodic chaotic system the cat’s eyes may be thought of as an “organizing structure” for mixing and transport in the surf zone. When this organizing structure exists, Eulerian and Lagrangian autocorrelations of the velocity derivatives indicate that velocity derivatives decorrelate more rapidly along particle trajectories than at fixed spatial locations (i.e., the velocity field is temporally quasi-regular). This phenomenon is referred to as Lagrangian random strain.

A closer look at chaotic advection in the stratosphere: part II: statistical diagnostics

Statistical diagnostics of mixing and transport are computed for a numerical model of forced shallow-water flow on the sphere and a middle-atmosphere general circulation model. In particular, particle dispersion statistics, transport fluxes, Liapunov exponents (probability density functions and ensemble averages), and tracer concentration statistics are considered. It is shown that the behavior of the diagnostics is in accord with that of kinematic chaotic advection models so long as stochasticity is sufficiently weak. Comparisons with random-strain theory are made.

Lateral boundary contributions to wave-activity invariants and nonlinear stability theorems for balanced dynamics

The slow advective-timescale dynamics of the atmosphere and oceans is referred to as balanced dynamics. An extensive body of theory for disturbances to basic flows exists for the quasi-geostrophic (QG) model of balanced dynamics, based on wave-activity invariants and nonlinear stability theorems associated with exact symmetry-based conservation laws. In attempting to extend this theory to the semi-geostrophic (SG) model of balanced dynamics, Kushner & Shepherd discovered lateral boundary contributions to the SG wave-activity invariants which are not present in the QG theory, and which affect the stability theorems. However, because of technical difficulties associated with the SG model, the analysis of Kushner & Shepherd was not fully nonlinear. This paper examines the issue of lateral boundary contributions to wave-activity invariants for balanced dynamics in the context of Salmon's nearly geostrophic model of rotating shallow-water flow. Salmon's model has certain similarities with the SG model, but also has important differences that allow the present analysis to be carried to finite amplitude. In the process, the way in which constraints produce boundary contributions to wave-activity invariants, and additional conditions in the associated stability theorems, is clarified. It is shown that Salmon's model possesses two kinds of stability theorems: an analogue of Ripa's small-amplitude stability theorem for shallow-water flow, and a finite-amplitude analogue of Kushner & Shepherd's SG stability theorem in which the ‘subsonic’ condition of Ripa's theorem is replaced by a condition that the flow be cyclonic along lateral boundaries. As with the SG theorem, this last condition has a simple physical interpretation involving the coastal Kelvin waves that exist in both models. Salmon's model has recently emerged as an important prototype for constrained Hamiltonian balanced models. The extent to which the present analysis applies to this general class of models is discussed.

Observed and simulated precipitation responses in wet and dry regions 1850–2100

Global warming is expected to enhance fluxes of fresh water between the surface and atmosphere, causing wet regions to become wetter and dry regions drier, with serious implications for water resource management. Defining the wet and dry regions as the upper 30% and lower 70% of the precipitation totals across the tropics (30° S–30° N) each month we combine observations and climate model simulations to understand changes in the wet and dry regions over the period 1850–2100. Observed decreases in precipitation over dry tropical land (1950–2010) are also simulated by coupled atmosphere–ocean climate models (−0.3%/decade) with trends projected to continue into the 21st century. Discrepancies between observations and simulations over wet land regions since 1950 exist, relating to decadal fluctuations in El Niño southern oscillation, the timing of which is not represented by the coupled simulations. When atmosphere-only simulations are instead driven by observed sea surface temperature they are able to adequately represent this variability over land. Global distributions of precipitation trends are dominated by spatial changes in atmospheric circulation. However, the tendency for already wet regions to become wetter (precipitation increases with warming by 3% K−1 over wet tropical oceans) and the driest regions drier (precipitation decreases of −2% K−1 over dry tropical land regions) emerges over the 21st century in response to the substantial surface warming.

Developing a reliable strategy to infer the effective soil hydraulic properties from field evaporation experiments for agro-hydrological models

The Richards equation has been widely used for simulating soil water movement. However, the take-up of agro-hydrological models using the basic theory of soil water flow for optimizing irrigation, fertilizer and pesticide practices is still low. This is partly due to the difficulties in obtaining accurate values for soil hydraulic properties at a field scale. Here, we use an inverse technique to deduce the effective soil hydraulic properties, based on measuring the changes in the distribution of soil water with depth in a fallow field over a long period, subject to natural rainfall and evaporation using a robust micro Genetic Algorithm. A new optimized function was constructed from the soil water contents at different depths, and the soil water at field capacity. The deduced soil water retention curve was approximately parallel but higher than that derived from published pedo-tranfer functions for a given soil pressure head. The water contents calculated from the deduced soil hydraulic properties were in good agreement with the measured values. The reliability of the deduced soil hydraulic properties was tested in reproducing data measured from an independent experiment on the same soil cropped with leek. The calculation of root water uptake took account for both soil water potential and root density distribution. Results show that the predictions of soil water contents at various depths agree fairly well with the measurements, indicating that the inverse analysis is an effective and reliable approach to estimate soil hydraulic properties, and thus permits the simulation of soil water dynamics in both cropped and fallow soils in the field accurately. (C) 2009 Elsevier B.V. All rights reserved.

Sea ice–ocean coupling using a rescaled vertical coordinate z*

Realistic representation of sea ice in ocean models involves the use of a non-linear free-surface, a real freshwater flux and observance of requisite conservation laws. We show here that these properties can be achieved in practice through use of a rescaled vertical coordinate ‘‘z*” in z-coordinate models that allows one to follow undulations in the free-surface under sea ice loading. In particular, the adoption of "z*" avoids the difficult issue of vanishing levels under thick ice. Details of the implementation within MITgcm are provided. A high resolution global ocean sea ice simulation illustrates the robustness of the z* formulation and reveals a source of oceanic variability associated with sea ice dynamics and ice-loading effects. The use of the z* coordinate allows one to achieve perfect conservation of fresh water, heat and salt, as shown in extended integration of coupled ocean sea ice atmospheric model.

Extratropical cyclones and the projected decline of winter Mediterranean precipitation in the CMIP5 models

The Mediterranean region has been identified as a climate change "hot-spot" due to a projected reduction in precipitation and fresh water availability which has potentially large socio-economic impacts. To increase confidence in these projections, it is important to physically understand how this precipitation reduction occurs. This study quantifies the impact on winter Mediterranean precipitation due to changes in extratropical cyclones in 17 CMIP5 climate models. In each model, the extratropical cyclones are objectively tracked and a simple approach is applied to identify the precipitation associated to each cyclone. This allows us to decompose the Mediterranean precipitation reduction into a contribution due to changes in the number of cyclones and a contribution due to changes in the amount of precipitation generated by each cyclone. The results show that the projected Mediterranean precipitation reduction in winter is strongly related to a decrease in the number of Mediterranean cyclones. However, the contribution from changes in the amount of precipitation generated by each cyclone are also locally important: in the East Mediterranean they amplify the precipitation trend due to the reduction in the number of cyclones, while in the North Mediterranean they compensate for it. Some of the processes that determine the opposing cyclone precipitation intensity responses in the North and East Mediterranean regions are investigated by exploring the CMIP5 inter-model spread.

Impact of cyclones on sea ice in the Fram Strait: numerical simulations and observations

The sea ice export from the Arctic is of global importance due to its fresh water which influences the oceanic stratification and, thus, the global thermohaline circulation. This study deals with the effect of cyclones on sea ice and sea ice transport in particular on the basis of observations from two field experiments FRAMZY 1999 and FRAMZY 2002 in April 1999 and March 2002 as well as on the basis of simulations with a numerical sea ice model. The simulations realised by a dynamic-thermodynamic sea ice model are forced with 6-hourly atmospheric ECMWF- analyses (European Centre for Medium-Range Weather Forecasts) and 6-hourly oceanic data of a MPI-OM-simulation (Max-Planck-Institute Ocean Model). Comparing the observed and simulated variability of the sea ice drift and of the position of the ice edge shows that the chosen configuration of the model is appropriate for the performed studies. The seven observed cyclones change the position of the ice edge up to 100 km and cause an extensive decrease of sea ice coverage by 2 % up to more than 10 %. The decrease is only simulated by the model if the ocean current is strongly divergent in the centre of the cyclone. The impact is remarkable of the ocean current on divergence and shear deformation of the ice drift. As shown by sensitivity studies the ocean current at a depth of 6 m – the sea ice model is forced with – is mainly responsible for the ascertained differences between simulation and observation. The simulated sea ice transport shows a strong variability on a time scale from hours to days. Local minima occur in the time series of the ice transport during periods with Fram Strait cyclones. These minima are not caused by the local effect of the cyclone’s wind field, but mainly by the large-scale pattern of surface pressure. A displacement of the areas of strongest cyclone activity in the Nordic Seas would considerably influence the ice transport.