The application of ichnofabrics towards bridging the dichotomy between siliciclastic and carbonate shelf facies: examples from the Upper Jurassic Fulmar Formation (UK) and the Jubaila Formation (Saudi Arabia)

An initial study of the ichnofabrics of the Upper Jurassic (Kimmeridgian) Jubaila Formation of Saudi Arabia shows that the ichnofabrics are closely matched to the relatively well-described ichnofabrics of the contemporary Fulmar Formation of the UK Continental Shelf (North Sea), in respect of the lower shoreface/offshore transition facies to offshore facies. The ichnology and ichnofabrics of the Lower Jubaila Formation show that deposition took place on an open-marine platform on the Arabian craton subject to periodic storm activity, but under a persisting equilibrium between sediment accumulation and subsidence. This is consistent with the moderately deep-marine foraminiferal assemblages and the presence of calcareous nannofossils. Cyclicity is absent, though storm beds may be grouped, in contrast with the genetic sequences present in the rift and halokinetic scenario of the North Sea. In contrast with the siliciclastic setting hardgrotinds (with Gastrochaenolites), more common firmground omission surfaces, and micritic mudstones with Chondrites and Zoophycos are notable features of the carbonate facies. In siliciclastic successions (parasequences) the latter ichnotaxa are generally regarded as having been deposited in rather deeper water, but in the carbonate Jubaila Formation are interpreted as being associated with local areas of lower turbulence. Likewise, the hardgrounds and firmgrounds, which have not been traced laterally, are tentatively regarded to be of local significance.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

Why the western Pacific subtropical high has extended westward since the late 1970s

The western Pacific subtropical high (WPSH) is closely related to Asian climate. Previous examination of changes in the WPSH found a westward extension since the late 1970s, which has contributed to the inter-decadal transition of East Asian climate. The reason for the westward extension is unknown, however. The present study suggests that this significant change of WPSH is partly due to the atmosphere's response to the observed Indian Ocean-western Pacific (IWP) warming. Coordinated by a European Union's Sixth Framework Programme, Understanding the Dynamics of the Coupled Climate System (DYNAMITE), five AGCMs were forced by identical idealized sea surface temperature patterns representative of the IWP warming and cooling. The results of these numerical experiments suggest that the negative heating in the central and eastern tropical Pacific and increased convective heating in the equatorial Indian Ocean/ Maritime Continent associated with IWP warming are in favor of the westward extension of WPSH. The SST changes in IWP influences the Walker circulation, with a subsequent reduction of convections in the tropical central and eastern Pacific, which then forces an ENSO/Gill-type response that modulates the WPSH. The monsoon diabatic heating mechanism proposed by Rodwell and Hoskins plays a secondary reinforcing role in the westward extension of WPSH. The low-level equatorial flank of WPSH is interpreted as a Kelvin response to monsoon condensational heating, while the intensified poleward flow along the western flank of WPSH is in accord with Sverdrup vorticity balance. The IWP warming has led to an expansion of the South Asian high in the upper troposphere, as seen in the reanalysis.

Atmosphere feedbacks during ENSO in a coupled GCM with a modified atmospheric convection scheme

The too diverse representation of ENSO in a coupled GCM limits one’s ability to describe future change of its properties. Several studies pointed to the key role of atmosphere feedbacks in contributing to this diversity. These feedbacks are analyzed here in two simulations of a coupled GCM that differ only by the parameterization of deep atmospheric convection and the associated clouds. Using the Kerry–Emanuel (KE) scheme in the L’Institut Pierre-Simon Laplace Coupled Model, version 4 (IPSL CM4; KE simulation), ENSO has about the right amplitude, whereas it is almost suppressed when using the Tiedke (TI) scheme. Quantifying both the dynamical Bjerknes feedback and the heat flux feedback in KE, TI, and the corresponding Atmospheric Model Intercomparison Project (AMIP) atmosphere-only simulations, it is shown that the suppression of ENSO in TI is due to a doubling of the damping via heat flux feedback. Because the Bjerknes positive feedback is weak in both simulations, the KE simulation exhibits the right ENSO amplitude owing to an error compensation between a too weak heat flux feedback and a too weak Bjerknes feedback. In TI, the heat flux feedback strength is closer to estimates from observations and reanalysis, leading to ENSO suppression. The shortwave heat flux and, to a lesser extent, the latent heat flux feedbacks are the dominant contributors to the change between TI and KE. The shortwave heat flux feedback differences are traced back to a modified distribution of the large-scale regimes of deep convection (negative feedback) and subsidence (positive feedback) in the east Pacific. These are further associated with the model systematic errors. It is argued that a systematic and detailed evaluation of atmosphere feedbacks during ENSO is a necessary step to fully understand its simulation in coupled GCMs.

The role of atmosphere feedbacks during ENSO in the CMIP3 models

Several studies using ocean-atmosphere GCMs suggest that the atmospheric component plays a dominant role in the modelled ENSO. To help elucidate these findings, the two main atmosphere feedbacks relevant to ENSO, the Bjerknes positive feedback (µ) and the heat flux negative feedback (), are analysed here in 12 coupled GCMs. We find that the models generally underestimate both feedbacks, leading to an error compensation. The strength of is inversely related to the ENSO amplitude in the models and the latent heat and shortwave flux components of this feedback dominate. Furthermore, the shortwave component could help explain the model diversity in both overall and ENSO amplitude.

Coupled simulations of the Indian monsoon intraseasonal oscillation using a fine-resolution mixed-layer ocean model

The intraseasonal variability (ISV) of the Indian summer monsoon is dominated by a 30–50 day oscillation between “active” and “break” events of enhanced and reduced rainfall over the subcontinent, respectively. These organized convective events form in the equatorial Indian Ocean and propagate north to India. Atmosphere–ocean coupled processes are thought to play a key role the intensity and propagation of these events. A high-resolution, coupled atmosphere–mixed-layer-oceanmodel is assembled: HadKPP. HadKPP comprises the Hadley Centre Atmospheric Model (HadAM3) and the K Profile Parameterization (KPP) mixed-layer ocean model. Following studies that upper-ocean vertical resolution and sub-diurnal coupling frequencies improve the simulation of ISV in SSTs, KPP is run at 1 m vertical resolution near the surface; the atmosphere and ocean are coupled every three hours. HadKPP accurately simulates the 30–50 day ISV in rainfall and SSTs over India and the Bay of Bengal, respectively, but suffers from low ISV on the equator. This is due to the HadAM3 convection scheme producing limited ISV in surface fluxes. HadKPP demonstrates little of the observed northward propagation of intraseasonal events, producing instead a standing oscillation. The lack of equatorial ISV in convection in HadAM3 constrains the ability of KPP to produce equatorial SST anomalies, which further weakens the ISV of convection. It is concluded that while atmosphere–ocean interactions are undoubtedly essential to an accurate simulation of ISV, they are not a panacea for model deficiencies. In regions where the atmospheric forcing is adequate, such as the Bay of Bengal, KPP produces SST anomalies that are comparable to the Tropical Rainfall Measuring Mission Microwave Imager (TMI) SST analyses in both their magnitude and their timing with respect to rainfall anomalies over India. HadKPP also displays a much-improved phase relationship between rainfall and SSTs over a HadAM3 ensemble forced by observed SSTs, when both are compared to observations. Coupling to mixed-layer models such as KPP has the potential to improve operational predictions of ISV, particularly when the persistence time of SST anomalies is shorter than the forecast lead time.

The basic ingredients of the North Atlantic storm track. Part I: land-sea contrast and orography

Understanding and predicting changes in storm tracks over longer time scales is a challenging problem, particularly in the North Atlantic. This is due in part to the complex range of forcings (land–sea contrast, orography, sea surface temperatures, etc.) that combine to produce the structure of the storm track. The impact of land–sea contrast and midlatitude orography on the North Atlantic storm track is investigated through a hierarchy of GCM simulations using idealized and “semirealistic” boundary conditions in a high-resolution version of the Hadley Centre atmosphere model (HadAM3). This framework captures the large-scale essence of features such as the North and South American continents, Eurasia, and the Rocky Mountains, enabling the results to be applied more directly to realistic modeling situations than was possible with previous idealized studies. The physical processes by which the forcing mechanisms impact the large-scale flow and the midlatitude storm tracks are discussed. The characteristics of the North American continent are found to be very important in generating the structure of the North Atlantic storm track. In particular, the southwest–northeast tilt in the upper tropospheric jet produced by southward deflection of the westerly flow incident on the Rocky Mountains leads to enhanced storm development along an axis close to that of the continent’s eastern coastline. The approximately triangular shape of North America also enables a cold pool of air to develop in the northeast, intensifying the surface temperature contrast across the eastern coastline, consistent with further enhancements of baroclinicity and storm growth along the same axis.

Very large inverse problems in atmosphere and ocean modelling

For the very large nonlinear dynamical systems that arise in a wide range of physical, biological and environmental problems, the data needed to initialize a numerical forecasting model are seldom available. To generate accurate estimates of the expected states of the system, both current and future, the technique of ‘data assimilation’ is used to combine the numerical model predictions with observations of the system measured over time. Assimilation of data is an inverse problem that for very large-scale systems is generally ill-posed. In four-dimensional variational assimilation schemes, the dynamical model equations provide constraints that act to spread information into data sparse regions, enabling the state of the system to be reconstructed accurately. The mechanism for this is not well understood. Singular value decomposition techniques are applied here to the observability matrix of the system in order to analyse the critical features in this process. Simplified models are used to demonstrate how information is propagated from observed regions into unobserved areas. The impact of the size of the observational noise and the temporal position of the observations is examined. The best signal-to-noise ratio needed to extract the most information from the observations is estimated using Tikhonov regularization theory. Copyright © 2005 John Wiley & Sons, Ltd.