Circumglobal wave train and the summer monsoon over northwestern India and Pakistan: the explicit role of the surface heat low

This study examines the influence of the mid-latitude circulation on the surface heat low (HL) and associated monsoon rainfall over northwestern India and Pakistan using the ERA40 data and high resolution (T106L31) climate model ECHAM5 simulation. Special emphasis is given to the surface HL which forms over Pakistan and adjoining areas of India, Iran and Afghanistan during the summer season. A heat low index (HLI) is defined to depict the surface HL. The HLI displays significant correlations with the upper level mid-latitude circulation over western central Asia and low level monsoon circulation over Arabian Sea and acts as a bridge connecting the mid-latitude wave train to the Indian summer monsoon. A time-lagged singular value decomposition analysis reveals that the eastward propagation of the mid-latitude circumglobal wave train (CGT) influences the surface pressure anomalies over the Indian domain. The largest low (negative) pressure anomalies over the western parts of the HL region (i.e., Iran and Afghanistan) occur in conjunction with the upper level anomalous high that develops over western-central Asia during the positive phase of the CGT. The composite analysis also reveals a significant increase in the low pressure anomalies over Iran and Afghanistan during the positive phase of CGT. The westward increasing low pressure anomalies with its north?Csouth orientation provokes enormous north?Csouth pressure gradient (lower pressure over land than over sea). This in turn enables the moist southerly flow from the Arabian Sea to penetrate farther northward over northwestern India and Pakistan. A monsoon trough like conditions develops over northwestern India and Pakistan where the moist southwesterly flow from the Arabian Sea and the Persian Gulf converge. The convergence in association with the orographic uplifting expedites convection and associated precipitation over northwestern India and Pakistan. The high resolution climate model ECHAM5 simulation also underlines the proposed findings and mechanism.     

Diagnosing GCM errors over West Africa using relaxation experiments. Part I: summer monsoon climatology and interannual variability

The CNRM atmospheric general circulation model Arpege-Climat is relaxed towards atmospheric reanalyses outside the 10°S?C32°N 30°W?C50°E domain in order to disentangle the regional versus large-scale sources of climatological biases and interannual variability of the West African monsoon (WAM). On the one hand, the main climatological features of the monsoon, including the spatial distribution of summer precipitation, are only weakly improved by the nudging, thereby suggesting the regional origin of the Arpege-Climat biases. On the other hand, the nudging technique is relatively efficient to control the interannual variability of the WAM dynamics, though the impact on rainfall variability is less clear. Additional sensitivity experiments focusing on the strong 1994 summer monsoon suggest that the weak sensitivity of the model biases is not an artifact of the nudging design, but the evidence that regional physical processes are the main limiting factors for a realistic simulation of monsoon circulation and precipitation in the Arpege-Climat model. Sensitivity experiments to soil moisture boundary conditions are also conducted and highlight the relevance of land?Catmosphere coupling for the amplification of precipitation biases. Nevertheless, the land surface hydrology is not the main explanation for the model errors that are rather due to deficiencies in the atmospheric physics. The intraseasonal timescale and the model internal variability are discussed in a companion paper.     

Processes controlling the surface temperature signature of the Madden�CJulian Oscillation in the thermocline ridge of the Indian Ocean

During boreal winter, there is a prominent maximum of intraseasonal sea-surface temperature (SST) variability associated with the Madden?CJulian Oscillation (MJO) along a Thermocline Ridge located in the southwestern Indian Ocean (5°S?C10°S, 60°E?C90°E; TRIO region). There is an ongoing debate about the relative importance of air-sea heat fluxes and oceanic processes in driving this intraseasonal SST variability. Furthermore, various studies have suggested that interannual variability of the oceanic structure in the TRIO region could modulate the amplitude of the MJO-driven SST response. In this study, we use observations and ocean general circulation model (OGCM) experiments to quantify these two effects over the 1997?C2006 period. Observational analysis indicates that Ekman pumping does not contribute significantly (on average) to intraseasonal SST variability. It is, however, difficult to quantify the relative contribution of net heat fluxes and entrainment to SST intraseasonal variability from observations alone. We therefore use a suite of OGCM experiments to isolate the impacts of each process. During 1997?C2006, wind stress contributed on average only about 20% of the intraseasonal SST variability (averaged over the TRIO region), while heat fluxes contributed about 70%, with forcing by shortwave radiation (75%) dominating the other flux components (25%). This estimate is consistent with an independent air-sea flux product, which indicates that shortwave radiation contributes 68% of intraseasonal heat flux variability. The time scale of the heat-flux perturbation, in addition to its amplitude, is also important in controlling the intraseasonal SST signature, with longer periods favouring a larger response. There are also strong year-to-year variations in the respective role of heat fluxes and wind stress. Of the five strong cooling events identified in both observations and the model (two in 1999 and one in 2000, 2001 and 2002), intraseasonal-wind stress dominates the SST signature during 2001 and contributes significantly during 2000. Interannual variations of the subsurface thermal structure associated with the Indian Ocean Dipole or El Ni?o/La Ni?a events modulate the MJO-driven SST signature only moderately (by up to 30%), mainly by changing the temperature of water entrained into the mixed layer. The primary factor that controls year-to-year changes in the amplitude of TRIO, intraseasonal SST anomalies is hence the characteristics of intraseasonal surface flux perturbations, rather than changes in the underlying oceanic state.     

Tropical influence on boreal summer mid-latitude stationary waves

While organized tropical convection is a well-known source of extratropical planetary waves, state-of-the-art climate models still show serious deficiencies in simulating accurately the atmospheric response to tropical sea surface temperature (SST) anomalies and the associated teleconnections. In the present study, the remote influence of the tropical atmospheric circulation is evaluated in ensembles of global boreal summer simulations in which the Arpege-Climat atmospheric General Circulation Model (GCM) is nudged towards 6-h reanalyses. The nudging is applied either in the whole tropical band or in a regional summer monsoon domain. Sensitivity tests to the experimental design are first conducted using prescribed climatological SST. They show that the tropical relaxation does not improve the zonal mean extratropical climatology but does lead to a significantly improved representation of the mid-latitude stationary waves in both hemispheres. Low-pass filtering of the relaxation fields has no major effect on the model response, suggesting that high-frequency tropical variability is not responsible for extratropical biases. Dividing the nudging strength by a factor 10 only decreases the magnitude of the response. Model errors in each monsoon domain contribute to deficiencies in the model??s mid-latitude climatology, although an exaggerated large-scale subsidence in the central equatorial Pacific appears as the main source of errors for the representation of stationary waves in the Arpege-Climat model. Case studies are then conducted using either climatological or observed SST. The focus is first on summer 2003 characterized by a strong and persistent anticyclonic anomaly over western Europe. This pattern is more realistic in nudging experiments than in simulations only driven by observed SST, especially when the nudging domain is centred over Central America. Other case studies also show a significant tropical forcing of the summer mid-latitude stationary waves and suggest a weak influence of prescribed observed SST in the northern extratropics. Results therefore indicate that improving the tropical divergent circulation and its response to tropical SST anomalies remains a key issue for increasing the skill of extratropical seasonal predictions, not only in the winter hemisphere but also in the boreal summer hemisphere where the prediction of heatwave and drought likelihood is expected to become an important challenge with increasing concentrations of greenhouse gases.     

Double-difference Relocation of the Aftershocks of the Tecom��n, Colima, Mexico Earthquake of 22 January 2003

On 22 January 2003, the M _w?=?7.6 Tecomán earthquake struck offshore of the state of Colima, Mexico, near the diffuse triple junction between the Cocos, Rivera, and North American plates. Three-hundred and fifty aftershocks of the Tecomán earthquake with magnitudes between 2.6 and 5.8, each recorded by at least 7 stations, are relocated using the double difference method. Initial locations are determined using P and S readings from the Red Sismológica Telemétrica del Estado de Colima (RESCO) and a 1-D velocity model. Because only eight RESCO stations were operating immediately following the Tecomán earthquake, uncertainties in the initial locations and depths are fairly large, with average uncertainties of 8.0?km in depth and 1.4?km in the north?Csouth and east?Cwest directions. Events occurring between 24 January and 31 January were located using not only RESCO phase readings but also additional P and S readings from 11 temporary stations. Average uncertainties decrease to 0.8?km in depth, 0.3?km in the east?Cwest direction, and 0.7?km in the north?Csouth direction for events occurring while the temporary stations were deployed. While some preliminary studies of the early aftershocks suggested that they were dominated by shallow events above the plate interface, our results place the majority of aftershocks along the plate interface, for a slab dipping between approximately 20° and 30°. This is consistent with the slab positions inferred from geodetic studies. We do see some upper plate aftershocks that may correspond to forearc fault zones, and faults inland in the upper plate, particularly among events occurring more than 3?months after the mainshock.     

Volcanic Markers of the Post-Subduction Evolution of Baja California and Sonora, Mexico: Slab Tearing Versus Lithospheric Rupture of the Gulf of California

The study of the geochemical compositions and K-Ar or Ar-Ar ages of ca. 350 Neogene and Quaternary lavas from Baja California, the Gulf of California and Sonora allows us to discuss the nature of their mantle or crustal sources, the conditions of their melting and the tectonic regime prevailing during their genesis and emplacement. Nine petrographic/geochemical groups are distinguished: ??regular?? calc-alkaline lavas; adakites; magnesian andesites and related basalts and basaltic andesites; niobium-enriched basalts; alkali basalts and trachybasalts; oceanic (MORB-type) basalts; tholeiitic/transitional basalts and basaltic andesites; peralkaline rhyolites (comendites); and icelandites. We show that the spatial and temporal distribution of these lava types provides constraints on their sources and the geodynamic setting controlling their partial melting. Three successive stages are distinguished. Between 23 and 13 Ma, calc-alkaline lavas linked to the subduction of the Pacific-Farallon plate formed the Comondú and central coast of the Sonora volcanic arc. In the extensional domain of western Sonora, lithospheric mantle-derived tholeiitic to transitional basalts and basaltic andesites were emplaced within the southern extension of the Basin and Range province. The end of the Farallon subduction was marked by the emplacement of much more complex Middle to Late Miocene volcanic associations, between 13 and 7 Ma. Calc-alkaline activity became sporadic and was replaced by unusual post-subduction magma types including adakites, niobium-enriched basalts, magnesian andesites, comendites and icelandites. The spatial and temporal distribution of these lavas is consistent with the development of a slab tear, evolving into a 200-km-wide slab window sub-parallel to the trench, and extending from the Pacific coast of Baja California to coastal Sonora. Tholeiitic, transitional and alkali basalts of subslab origin ascended through this window, and adakites derived from the partial melting of its upper lip, relatively close to the trench. Calc-alkaline lavas, magnesian andesites and niobium-enriched basalts formed from hydrous melting of the supraslab mantle triggered by the uprise of hot Pacific asthenosphere through the window. During the Plio-Quaternary, the ??no-slab?? regime following the sinking of the old part of the Farallon plate within the deep mantle allowed the emplacement of alkali and tholeiitic/transitional basalts of deep asthenospheric origin in Baja California and Sonora. The lithospheric rupture connected with the opening of the Gulf of California generated a high thermal regime associated to asthenospheric uprise and emplaced Quaternary depleted MORB-type tholeiites. This thermal regime also induced partial melting of the thinned lithospheric mantle of the Gulf area, generating calc-alkaline lavas as well as adakites derived from slivers of oceanic crust incorporated within this mantle.     

Overview of Recent Coastal Tectonic Deformation in the Mexican Subduction Zone

Holocene and Pleistocene tectonic deformation of the coast in the Mexico subudction margin is recorded by geomorphic and stratigraphic markers. We document the spatial and temporal variability of active deformation on the coastal Mexican subduction margin. Pleistocene uplift rates are estimated using wave-cut platforms at ca. 0.7?C0.9?m/ka on the Jalisco block coast, Rivera-North America tectonic plate boundary. We examine reported measurements from marine notches and shoreline angle elevations in conjunction with their radiocarbon ages that indicate surface uplift rates increasing during the Holocene up to ca. 3?±?0.5?m/ka. In contrast, steady rates of uplift (ca. 0.5?C1.0?m/ka) in the Pleistocene and Holocene characterize the Michoacan coastal sector, south of El Gordo graben and north of the Orozco Fracture Zone (OFZ), incorporated within the Cocos-North America plate boundary. Significantly higher rates of surface uplift (ca. 7?m/ka) across the OFZ subduction may reflect the roughness of subducting plate. Absence of preserved marine terraces on the coastal sector across El Gordo graben likely reflects slow uplift or coastal subsidence. Stratigraphic markers and their radiocarbon ages show late Holocene (ca. last 6?ka bp) coastal subsidence on the Guerrero gap sector in agreement with a landscape barren of marine terraces and with archeological evidence of coastal subsidence. Temporal and spatial variability in recent deformation rates on the Mexican Pacific coast may be due to differences in tectonic regimes and to localized processes related to subduction, such as crustal faults, subduction erosion and underplating of subducted materials under the southern Mexico continental margin.     

Scenario-Based Seismic Risk Analysis: An Engineering Approach to the Development of Source and Site-Specific Ground Motion Time Histories in Areas of Low Seismicity

Modern engineering design methods require ground motion time histories as input for non-linear dynamic structural analysis. Non-linear dynamic methods of analysis are increasingly applied in the context of probabilistic risk assessments and for cost-effective design of critical infrastructures. In current engineering practice artificial time histories matching deterministic design spectra or probabilistic uniform hazard spectra are most frequently used for engineering analysis. The intermediate step of generation of response spectra can lead to a biased estimate of the potential damage from earthquakes because of insufficient consideration of the true energy content and strong motion duration of earthquakes. Thus, assessment of seismic risk may seem unrealistic. An engineering approach to the development of three-component ground motion time histories has been established which enables consideration of the typical characteristics of seismic sources, regional ground motion attenuation, and the main geotechnical characteristics of the target site. Therefore, the approach is suitable for use in scenario-based risk analysis a larger number of time histories are required for representation of the seismic hazard. Near-field effects are implemented in the stochastic source model using engineering approximations. The approach is suggested for use in areas of low seismicity where ground motion records of larger earthquakes are not available. Uncertainty analysis indicates that ground motions generated by individual earthquakes are well constrained and that the usual lognormal model is not the best choice for predicting the upper tail of the distribution of the ground motions.     

A Downscaling Approach Toward High-Resolution Surface Mass Balance Over Antarctica

The Antarctic ice sheet surface mass balance shows high spatial variability over the coastal area. As state-of-the-art climate models usually require coarse resolutions to keep computational costs to a moderate level, they miss some local features that can be captured by field measurements. The downscaling approach adopted here consists of using a cascade of atmospheric models from large scale to meso-?? scale. A regional climate model (Modèle Atmosphérique Régional) forced by meteorological reanalyses provides a diagnostic physically-based rain- and snowfall downscaling model with meteorological fields at the regional scale. Although the parameterizations invoked by the downscaling model are fairly simple, the knowledge of small-scale topography significantly improves the representation of spatial variability of precipitation and therefore that of the surface mass balance. Model evaluation is carried out with the help of shallow firn cores and snow height measurements provided by automatic weather stations. Although downscaling of blowing snow still needs to be implemented in the model, the net accumulation gradient across Law Dome summit is shown to be induced mostly by orographic effects on precipitation.