Relaxing the Tropics to an ‘observed’ state: analysis using a simple baroclinic model

The technique of relaxation of the tropical atmosphere towards an analysis in a month-season forecast model has previously been successfully exploited in a number of contexts. Here it is shown that when tropical relaxation is used to investigate the possible origin of the observed anomalies in June–July 2007, a simple dynamical model is able to reproduce the observed component of the pattern of anomalies given by an ensemble of ECMWF forecast runs. Following this result, the simple model is used for a range of experiments on time-scales of relaxation, variables and regions relaxed based on a control model run with equatorial heating in a zonal flow. A theory based on scale analysis for the large-scale tropics is used to interpret the results. Typical relationships between scales are determined from the basic equations, and for a specified diabatic heating a chain of deductions for determining the dependent variables is derived. Different critical time-scales are found for tropical relaxation of different dependent variables to be effective. Vorticity has the longest critical time-scale, typically 1.2 days. For temperature and divergence, the time-scales are 10 hours and 3 hours, respectively. However not all the tropical fields, in particular the vertical motion, are reproduced correctly by the model unless divergence is heavily damped. To obtain the correct extra-tropical fields, it is crucial to have the correct rotational flow in the subtropics to initiate the Rossby wave propagation from there. It is sufficient to relax vorticity or temperature on a time-scale comparable or less than their critical time-scales to obtain this. However if the divergent advection of vorticity is important in the Rossby Wave Source then strong relaxation of divergence is required to accurately represent the tropical forcing of Rossby waves.

Importance of oceanic resolution and mean state on the extra-tropical response to El Niño in a matrix of coupled models

The extra-tropical response to El Niño in configurations of a coupled model with increased horizontal resolution in the oceanic component is shown to be more realistic than in configurations with a low resolution oceanic component. This general conclusion is independent of the atmospheric resolution. Resolving small-scale processes in the ocean produces a more realistic oceanic mean state, with a reduced cold tongue bias, which in turn allows the atmospheric model component to be forced more realistically. A realistic atmospheric basic state is critical in order to represent Rossby wave propagation in response to El Niño, and hence the extra-tropical response to El Niño. Through the use of high and low resolution configurations of the forced atmospheric-only model component we show that, in isolation, atmospheric resolution does not significantly affect the simulation of the extra-tropical response to El Niño. It is demonstrated, through perturbations to the SST forcing of the atmospheric model component, that biases in the climatological SST field typical of coupled model configurations with low oceanic resolution can account for the erroneous atmospheric basic state seen in these coupled model configurations. These results highlight the importance of resolving small-scale oceanic processes in producing a realistic large-scale mean climate in coupled models, and suggest that it might may be possible to “squeeze out” valuable extra performance from coupled models through increases to oceanic resolution alone.

Trefftz discontinuous Galerkin methods for acoustic scattering on locally refined meshes

We extend the a priori error analysis of Trefftz-discontinuous Galerkin methods for time-harmonic wave propagation problems developed in previous papers to acoustic scattering problems and locally refined meshes. To this aim, we prove refined regularity and stability results with explicit dependence of the stability constant on the wave number for non convex domains with non connected boundaries. Moreover, we devise a new choice of numerical flux parameters for which we can prove L2-error estimates in the case of locally refined meshes near the scatterer. This is the setting needed to develop a complete hp-convergence analysis.

Uncertainties and assessments of chemistry-climate models of the stratosphere

In recent years a number of chemistry-climate models have been developed with an emphasis on the stratosphere. Such models cover a wide range of time scales of integration and vary considerably in complexity. The results of specific diagnostics are here analysed to examine the differences amongst individual models and observations, to assess the consistency of model predictions, with a particular focus on polar ozone. For example, many models indicate a significant cold bias in high latitudes, the “cold pole problem”, particularly in the southern hemisphere during winter and spring. This is related to wave propagation from the troposphere which can be improved by improving model horizontal resolution and with the use of non-orographic gravity wave drag. As a result of the widely differing modelled polar temperatures, different amounts of polar stratospheric clouds are simulated which in turn result in varying ozone values in the models. The results are also compared to determine the possible future behaviour of ozone, with an emphasis on the polar regions and mid-latitudes. All models predict eventual ozone recovery, but give a range of results concerning its timing and extent. Differences in the simulation of gravity waves and planetary waves as well as model resolution are likely major sources of uncertainty for this issue. In the Antarctic, the ozone hole has probably reached almost its deepest although the vertical and horizontal extent of depletion may increase slightly further over the next few years. According to the model results, Antarctic ozone recovery could begin any year within the range 2001 to 2008. The limited number of models which have been integrated sufficiently far indicate that full recovery of ozone to 1980 levels may not occur in the Antarctic until about the year 2050. For the Arctic, most models indicate that small ozone losses may continue for a few more years and that recovery could begin any year within the range 2004 to 2019. The start of ozone recovery in the Arctic is therefore expected to appear later than in the Antarctic.

Development of convection along the SPCZ within a Madden-Julian oscillation

A subtropical Rossby-wave propagation mechanism is proposed to account for the poleward and eastward progression of intraseasonal convective anomalies along the South Pacific convergence zone (SPCZ) that is observed in a significant proportion of Madden–Julian oscillations (MJOs). Large-scale convection, associated with an MJO, is assumed to be already established over the Indonesian region. The latent heating associated with this convection forces an equatorial Rossby-wave response with an upper-tropospheric anticyclone centred over, or slightly to the west of, the convection. Large potential-vorticity (PV) gradients, associated with the subtropical jet and the tropopause, lie just poleward of the anticyclone, and large magnitude PV air is advected equatorwards on the eastern side of the anticyclone. This ‘high’ PV air, or upper-tropospheric trough, is far enough off the equator that it has associated strong horizontal temperature gradients, and it induces deep ascent on its eastern side, at a latitude of about 15–30°. If this deep ascent is over a region susceptible to deep convection, such as the SPCZ, then convection may be forced or triggered. Hence convection develops along the SPCZ as a forced response to convection over Indonesia. The response mechanism is essentially one of subtropical Rossby-wave propagation. This hypothesis is based on a case study of a particularly strong MJO in early 1988, and is tested by idealized modelling studies. The mechanism may also be relevant to the existence of the mean SPCZ, as a forced response to mean Indonesian convection.

Atmospheric response in summer linked to recent Arctic sea ice loss

Since 2007 a large decline in Arctic sea ice has been observed. The large-scale atmospheric circulation response to this decline is investigated in ERA-Interim reanalyses and HadGEM3 climate model experiments. In winter, post-2007 observed circulation anomalies over the Arctic, North Atlantic and Eurasia are small compared to interannual variability. In summer, the post-2007 observed circulation is dominated by an anticyclonic anomaly over Greenland which has a large signal-to-noise ratio. Climate model experiments driven by observed SST and sea ice anomalies are able to capture the summertime pattern of observed circulation anomalies, although the magnitude is a third of that observed. The experiments suggest high SSTs and reduced sea ice in the Labrador Sea lead to positive temperature anomalies in the lower troposphere which weaken the westerlies over North America through thermal wind balance. The experiments also capture cyclonic anomalies over Northwest Europe, which are consistent with downstream Rossby wave propagation

Efficient incorporation of channel cross-section geometry uncertainty into regional and global scale flood inundation models

This paper investigates the challenge of representing structural differences in river channel cross-section geometry for regional to global scale river hydraulic models and the effect this can have on simulations of wave dynamics. Classically, channel geometry is defined using data, yet at larger scales the necessary information and model structures do not exist to take this approach. We therefore propose a fundamentally different approach where the structural uncertainty in channel geometry is represented using a simple parameterization, which could then be estimated through calibration or data assimilation. This paper first outlines the development of a computationally efficient numerical scheme to represent generalised channel shapes using a single parameter, which is then validated using a simple straight channel test case and shown to predict wetted perimeter to within 2% for the channels tested. An application to the River Severn, UK is also presented, along with an analysis of model sensitivity to channel shape, depth and friction. The channel shape parameter was shown to improve model simulations of river level, particularly for more physically plausible channel roughness and depth parameter ranges. Calibrating channel Manning’s coefficient in a rectangular channel provided similar water level simulation accuracy in terms of Nash-Sutcliffe efficiency to a model where friction and shape or depth were calibrated. However, the calibrated Manning coefficient in the rectangular channel model was ~2/3 greater than the likely physically realistic value for this reach and this erroneously slowed wave propagation times through the reach by several hours. Therefore, for large scale models applied in data sparse areas, calibrating channel depth and/or shape may be preferable to assuming a rectangular geometry and calibrating friction alone.

The combined influences of westerly phase of the quasi biennial oscillation and 11-year solar maximum conditions on the Northern Hemisphere extratropical winter circulation

The combined influences of the westerly phase of the quasi-biennial oscillation (QBO-W) and solar maximum (Smax) conditions on the Northern Hemisphere extratropical winter circulation are investigated using reanalysis data and Center for Climate System Research/National Institute for Environmental Studies chemistry climate model (CCM) simulations. The composite analysis for the reanalysis data indicates strengthened polar vortex in December followed by weakened polar vortex in February–March for QBO-W during Smax (QBO-W/Smax) conditions. This relationship need not be specific to QBO-W/Smax conditions but may just require strengthened vortex in December, which is more likely under QBO-W/Smax. Both the reanalysis data and CCM simulations suggest that dynamical processes of planetary wave propagation and meridional circulation related to QBO-W around polar vortex in December are similar in character to those related to Smax; furthermore, both processes may work in concert to maintain stronger vortex during QBO-W/Smax. In the reanalysis data, the strengthened polar vortex in December is associated with the development of north–south dipole tropospheric anomaly in the Atlantic sector similar to the North Atlantic oscillation (NAO) during December–January. The structure of the north–south dipole anomaly has zonal wavenumber 1 (WN1) component, where the longitude of anomalous ridge overlaps with that of climatological ridge in the North Atlantic in January. This implies amplification of the WN1 wave and results in the enhancement of the upward WN1 propagation from troposphere into stratosphere in January, leading to the weakened polar vortex in February–March. Although WN2 waves do not play a direct role in forcing the stratospheric vortex evolution, their tropospheric response to QBO-W/Smax conditions appears to be related to the maintenance of the NAO-like anomaly in the high-latitude troposphere in January. These results may provide a possible explanation for the mechanisms underlying the seasonal evolution of wintertime polar vortex anomalies during QBO-W/Smax conditions and the role of troposphere in this evolution.

ENSO related variation of equatorial MRG and Rossby waves and forcing from higher latitudes

The contrasting behaviour of westward-moving mixed Rossby-gravity (WMRG) and the first Rossby (R1) waves in El Niño (EN) and La Niña (LN) seasons is documented with a focus on the Northern Hemisphere winter. The eastward-moving variance in the upper troposphere is dominated by WMRG and R1 structures that appear to be Doppler-shifted by the flow and are referred to as WMRG-E and R1-E. In the East Pacific and Atlantic the years with stronger equatorial westerly winds have the stronger WMRG and WMRG- E. In the East Pacific, R1 is also a maximum in LN. However, R1-E exhibits an eastward-shift between LN and EN. The changes with ENSO phase provide a test-bed for the understanding of these waves. In the East Pacific and Atlantic, the stronger WMRG-E and WMRG with stronger westerlies are in accord with the dispersion relation with simple Doppler-shifting by the zonal flow. The possible existence of free waves can also explain stronger R1 in EN in the Eastern Hemisphere. 1-D free wave propagation theory based on wave activity conservation is also important for R1. However, this theory is unable to explain the amplitude maxima for other waves observed in the strong equatorial westerly regions in the Western Hemisphere, and certainly not their ENSO-related variation. The forcing of equatorial waves by higher latitude wave activity and its variation with ENSO phase is therefore examined. Propagation of extratropical eastward-moving Rossby wave activity through the westerly ducts into the equatorial region where it triggers WMRG-E is favoured in the stronger westerlies, in LN in the East Pacific and EN in the Atlantic. It is also found that WMRG is forced by Southern Hemisphere westward-moving wavetrains arching into the equatorial region where they are reflected. The most significant mechanism for both R1 and R1-E appear to be lateral forcing by subtropical wavetrains.

Climatological features of cutoff low systems in the Southern Hemisphere

Cutoff lows (COLs) pressure systems climatology for the Southern Hemisphere (SH), between 10 degrees S and 50 degrees S, using the National Center for Environmental Prediction-National Center for Atmospheric Research (NCEP-NCAR) and the ERA-40 European Centre for Medium Range Weather Forecast (ECMWF) reanalyses are analyzed for the period 1979-1999. COLs were identified at three pressure levels (200, 300, and 500 hPa) using an objective method that considers the main physical characteristics of the conceptual model of COLs. Independently of the pressure level analyzed, the climatology from the ERA-40 reanalysis has more COLs systems than the NCEP-NCAR. However, both reanalyses present a large frequency of COLs at 300 hPa, followed by 500 and 200 hPa. The seasonality of COLs differs at each pressure level, but it is similar between the reanalyses. COLs are more frequent during summer, autumn, and winter at 200, 300, and 500 hPa, respectively. At these levels, they tend to occur around the continents, preferentially from southeastern Australia to New Zealand, the south of South America, and the south of Africa. To study the COLs at 200 and 300 hPa from a regional perspective, the SH was divided in three regions: Australia-New Zealand (60 E-130 W), South America (130 degrees W-20 degrees W), and southern Africa (20 degrees W-60 degrees E). The common COLs features in these sectors for both reanalyses are a short lifetime (similar to 80.0% and similar to 70.0% of COLs at 200 and 300 hPa, respectively, persisting for up to 3 days), mobility (similar to 70.0% and similar to 50% of COLs at 200 and 300 hPa, respectively, traveling distances of up to 1200 km), and an eastward propagation.

The role of the South Indian and Pacific oceans in South American monsoon variability

This work has investigated the impact of three different low-frequency sea surface temperature (SST) variability modes located in the Indian and the Pacific Oceans on the interannual variability of the South American Monsoon System (SAMS) using observed and numerical data. Rotated Empirical Orthogonal Function (REOF) analysis and numerical simulations with a General Circulation Model (GCM) were used. One of the three SST variability modes is located close to southeastern Africa. According to the composites, warmer waters over this region are associated with enhanced austral summer precipitation over the sub-tropics. The GCM is able to reproduce this anomalous precipitation pattern, simulating a wave train emanating from the Indian Ocean towards South America (SA). A second SST variability mode was located in the western Pacific Ocean. REOF analysis indicates that warmer waters are associated with drought conditions over the South Atlantic Convergence Zone (SACZ) and enhanced precipitation over the sub-tropics. The GCM indicates that the warmer waters over Indonesia generate drought conditions over tropical SA through a Pacific South America-like (PSA) wave pattern emanating from the western Pacific. Finally, the third SST variability mode is located over the southwestern South Pacific. The composites indicate that warmer waters are associated with enhanced precipitation over the SACZ and drought conditions over the sub-tropics. There is a PSA-like wave train emanating from Indonesia towards SA, and another crossing the Southern Hemisphere in the extra-tropics, probably associated with transient activity. The GCM is able to reproduce the anomalous precipitation pattern, although it is weaker than observed. The PSA-like pattern is simulated, but the model fails in reproducing the extra-tropical wave activity.

Moisture transport and intraseasonal variability in the South America monsoon system

This paper examines moisture transport on intraseasonal timescales over the continent and over the South Atlantic convergence zone (SACZ) during the South America (SA) summer monsoon. Combined Empirical Orthogonal Function analysis (EOFc) of Global Precipitation Climatology Project pentad precipitation, specific humidity, air temperature, zonal and meridional winds at 850 hPa (NCEP/NCAR reanalysis) are performed to identify the large-scale variability of the South America monsoon system and the SACZ. The first EOFc was used as a large-scale index for the South American monsoon (LISAM), whereas the second EOFc characterized the SACZ. LISAM (SACZ) index showed spectral variance on 30-90 (15-20) days and were both band filtered (10-100 days). Intraseasonal wet anomalies were defined when LISAM and SACZ anomalies were above the 75th percentile of their respective distribution. LISAM and SACZ wet events were examined independently of each other and when they occur simultaneously. LISAM wet events were observed with the amplification of wave activity in the Northern Hemisphere and the enhancement of northwesterly cross-equatorial moisture transport over tropical continental SA. Enhanced SACZ was observed with moisture transport from the extratropics of the Southern Hemisphere. Simultaneous LISAM and SACZ wet events are associated with cross-equatorial moisture transport along with moisture transport from Subtropical Southwestern Atlantic.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Sediment deformation in Miocene and post-Miocene strata, Northeastern Brazil: Evidence for paleoseismicity in a passive margin

Stratigraphic intervals characterized by varied and complex styles of soft-sediment deformation structures are well preserved in Miocene and Late Pleistocene to Holocene deposits of a sedimentary basin located in Northeastern Brazil. The Miocene strata, represented by the Barreiras Formation, record only brittle structures, including numerous faults and fractures with straight and high angle-dipping planes that are often filled with sands derived from overlying beds. Folds consisting of broad anticlines and synclines are also present in this unit. The late Pleistocene to Holocene deposits, named Post-Barreiras Sediments, contain an indurated sandy package with a large variety of ductile and brittle deformation structures (i.e., massive sandstones with isolated sand fragments and breccias, undulatory strata, sand dykes and diapirs, sinks and bowls, pebbly pockets, plunged sediment mixtures, fitted sand masses, cone-shaped cracks, fault grading and sedimentary enclaves). These features, confined to sharp-based stratigraphic horizons that progressively grade downward into undisturbed deposits, are related to seismic shocks of high surface-wave magnitude (i.e., Ms>5 or 6). Amalgamated seismites suggest that previously formed seismites were affected by subsequent seismic-wave propagation. Seismic waves caused by activity along one, or most likely, several tectonic structures would have propagated throughout the depositional environment, producing laterally extensive seismites. The close proximity to earthquake epicenters would have promoted pervasive re-sedimentation due to pore overpressure, resulting high volumes of massive sandstones and breccia. The similarity between deposits with correlatable strata from many other areas along the Brazilian coast allows raise the hypothesis that the seismic episodes might have affected sedimentation patterns in a large (i.e., extension of several hundreds of kilometers) geographic area. Thus, the modern seismicity recorded along Northeastern Brazil was recurrent during the Quaternary and, perhaps, also in the Pliocene. The estimated high magnitude of the seismic events and the great regional extent of the affected area demonstrate that the Brazilian coast experienced tectonic stress through the last geological episodes of its evolution, which would have favored sediment accumulation and penecontemporaneous re-sedimentation. This geological context is unexpected in a passive margin, inducing to revisit the debate on how active is a passive margin. (C) 2010 Elsevier B.V. All rights reserved.