Recursive computation of finite difference of associated Legendre functions     

Toshio Fukushima 《Journal of Geodesy》2012,86(9):745-754

The existing methods to compute the definite integral of associated Legendre function (ALF) with respect to the argument suffer from a loss of significant figures independently of the latitude. This is caused by the subtraction of similar quantities in the additional term of their recurrence formulas, especially the finite difference of their values between two endpoints of the integration interval. In order to resolve the problem, we develop a recursive algorithm to compute their finite difference. Also, we modify the algorithm to evaluate their definite integrals assuming that their values at one endpoint are known. We numerically confirm a significant increase in computing precision of the integral by the new method. When the interval is one arc minute, for example, the gain amounts to 2–4 digits for the degree of harmonics in the range 2 ≤ n ≤ 2,048. This improvement in precision is achieved at a negligible increase in CPU time, say less than 5%.  相似文献   

Equivalent forcing depth in tropical oceans     

Takeshi Doi  Tomoki Tozuka  Toshio Yamagata 《Dynamics of Atmospheres and Oceans》2010,50(3):415-423

The forcing efficiency for the first and the second baroclinic modes by the wind stress in tropical oceans has been discussed by calculating equivalent forcing depth from annual mean, seasonal, and pentadal density profiles of the observational data. In the annual mean field, the first mode is forced preferentially in the western Pacific and the Indian Ocean, whereas the second mode is more strongly excited in the Atlantic and the eastern Pacific. This difference is mostly due to the pycnocline depth; the second mode is more dominantly forced where the pycnocline depth is shallower. We also revealed large seasonal variations of the second mode's equivalent forcing depth in the western Indian Ocean. The first mode is more dominantly forced during boreal spring and fall in the western Indian Ocean, while the second mode becomes more dominantly forced during boreal summer and winter. Those are due to seasonal variations of both the zonal wind and the pycnocline depth. Moreover, we show that the excitation of the second mode in the western Pacific increases after the late 1970s, which is associated with the decreasing trend of the zonal pycnocline gradient. Revealing the variation of the equivalent forcing depth will be useful for understanding the oceanic response to winds in tropical oceans and the improvement in the predictability of air-sea coupled climate variability in the tropics.  相似文献   

How are seasonal prediction skills related to models’ performance on mean state and annual cycle?     

June-Yi Lee  Bin Wang  I.-S. Kang  J. Shukla  A. Kumar  J.-S. Kug  J. K. E. Schemm  J.-J. Luo  T. Yamagata  X. Fu  O. Alves  B. Stern  T. Rosati  C.-K. Park 《Climate Dynamics》2010,35(2-3):267-283

Given observed initial conditions, how well do coupled atmosphere–ocean models predict precipitation climatology with 1-month lead forecast? And how do the models’ biases in climatology in turn affect prediction of seasonal anomalies? We address these questions based on analysis of 1-month lead retrospective predictions for 21 years of 1981–2001 made by 13 state-of-the-art coupled climate models and their multi-model ensemble (MME). The evaluation of the precipitation climatology is based on a newly designed metrics that consists of the annual mean, the solstitial mode and equinoctial asymmetric mode of the annual cycle, and the rainy season characteristics. We find that the 1-month lead seasonal prediction made by the 13-model ensemble has skills that are much higher than those in individual model ensemble predictions and approached to those in the ERA-40 and NCEP-2 reanalysis in terms of both the precipitation climatology and seasonal anomalies. We also demonstrate that the skill for individual coupled models in predicting seasonal precipitation anomalies is positively correlated with its performances on prediction of the annual mean and annual cycle of precipitation. In addition, the seasonal prediction skill for the tropical SST anomalies, which are the major predictability source of monsoon precipitation in the current coupled models, is closely link to the models’ ability in simulating the SST mean state. Correction of the inherent bias in the mean state is critical for improving the long-lead seasonal prediction. Most individual coupled models reproduce realistically the long-term annual mean precipitation and the first annual cycle (solstitial mode), but they have difficulty in capturing the second annual (equinoctial asymmetric) mode faithfully, especially over the Indian Ocean (IO) and Western North Pacific (WNP) where the seasonal cycle in SST has significant biases. The coupled models replicate the monsoon rain domains very well except in the East Asian subtropical monsoon and the tropical WNP summer monsoon regions. The models also capture the gross features of the seasonal march of the rainy season including onset and withdraw of the Asian–Australian monsoon system over four major sub-domains, but striking deficiencies in the coupled model predictions are observed over the South China Sea and WNP region, where considerable biases exist in both the amplitude and phase of the annual cycle and the summer precipitation amount and its interannual variability are underestimated.  相似文献   

Low and high frequency Madden–Julian oscillations in austral summer: interannual variations     

Takeshi Izumo  Sébastien Masson  Jérome Vialard  Clément de Boyer Montegut  Swadhin K. Behera  Gurvan Madec  Keiko Takahashi  Toshio Yamagata 《Climate Dynamics》2010,35(4):669-683

The Madden–Julian oscillation (MJO) is the main component of intraseasonal variability of the tropical convection, with clear climatic impacts at an almost-global scale. Based on satellite observations, it is shown that there are two types of austral-summer MJO events (broadly defined as 30–120 days convective variability with eastward propagation of about 5 m/s). Equatorial MJO events have a period of 30–50 days and tend to be symmetric about the equator, whereas MJO events centered near 8°S tend to have a longer period of 55–100 days. The lower-frequency variability is associated with a strong upper-ocean response, having a clear signature in both sea surface temperature and its diurnal cycle. These two MJO types have different interannual variations, and are modulated by the Indian Ocean Dipole (IOD). Following a negative IOD event, the lower-frequency southern MJO variability increases, while the higher-frequency equatorial MJO strongly diminishes. We propose two possible explanations for this change in properties of the MJO. One possibility is that changes in the background atmospheric circulation after an IOD favour the development of the low-frequency MJO. The other possibility is that the shallower thermocline ridge and mixed layer depth, by enhancing SST intraseasonal variability and thus ocean–atmosphere coupling in the southwest Indian Ocean (the breeding ground of southern MJO onset), favour the lower-frequency southern MJO variability.  相似文献   

Climate variability in the southern Indian Ocean as revealed by self-organizing maps     

Yushi Morioka  Tomoki Tozuka  Toshio Yamagata 《Climate Dynamics》2010,35(6):1059-1072

Using a non-linear statistical analysis called “self-organizing maps”, the interannual sea surface temperature (SST) variations in the southern Indian Ocean are investigated. The SST anomalies during austral summer from 1951 to 2006 are classified into nine types with differences in the position of positive and negative SST anomaly poles. To investigate the evolution of these SST anomaly poles, heat budget analysis of mixed-layer using outputs from an ocean general circulation model is conducted. The warming of the mixed-layer by the climatological shortwave radiation is enhanced (suppressed) as a result of negative (positive) mixed-layer thickness anomaly over the positive (negative) SST anomaly pole. This contribution from shortwave radiation is most dominant in the growth of SST anomalies. In contrast to the results reported so far, the contribution from latent heat flux anomaly is not so important. The discrepancy in the analysis is explained by the modulation in the contribution from the climatological heat flux by the interannual mixed-layer depth anomaly that was neglected in the past studies.  相似文献   

Generation and termination of Indian Ocean dipole events in 2003, 2006 and 2007     

Suryachandra A. Rao  Jing-Jia Luo  Swadhin K. Behera  Toshio Yamagata 《Climate Dynamics》2009,33(6):751-767

Evolution of Indian Ocean Dipole (IOD) events in 2003, 2006 and 2007 is investigated using observational and re-analysis data products. Efforts are made to understand various processes involved in three phases of IOD events; activation, maturation and termination. Three different triggers are found to activate the IOD events. In preceding months leading to the IOD evolution, the thermocline in southeastern Indian Ocean shoals by reflection of near equatorial upwelling Rossby waves at the East African coast into anomalous upwelling equatorial Kelvin waves. Strengthening (weakening) of northern (southern) portion of ITCZ in March/April and May/June of IOD years, leads to strengthening of alongshore winds along Sumatra/Java coasts. With the combined shallow thermocline and increased latent heat flux due to enhanced wind speeds, the SST in the southeastern Indian Ocean cools in following months. On intraseasonal time scales convection-suppressing phase of Madden-Julian oscillation (MJO) propagates from west to east in May/June of IOD year, and easterlies associated with this phase of MJO causes further shoaling of thermocline in southeastern Indian Ocean, through anomalous upwelling Kelvin wave. All these three mechanisms appear to be involved in initiating IOD event in 2006. On the other hand, except the strengthening/weakening of ITCZ, all other mechanisms are involved in activation of 2003 IOD event. Activation of 2007 IOD event was due to propagation of convection-suppressing MJO in May/June and strengthening of mean winds along Sumatra/Java coast from March to June through changes in convection. The IOD events matured into full-fledged events in the following months after activation, by surface heat fluxes, vertical and horizontal advection of cool waters supported by local along-shore upwelling favorable winds and remote equatorial easterly wind anomalies through excitation of upwelling Kelvin waves. Propagating MJO signals in the tropical Indian Ocean brings significant changes in evolution of IOD events on MJO time scales. Termination of 2003 and 2007 IOD events is achieved by strong convection-enhancing MJOs propagating from west to east in the tropical Indian Ocean which deepen the thermocline in the southeastern equatorial Indian Ocean. IOD event in 2006 was terminated by seasonal reversal of monsoon winds along Sumatra/Java coasts which stops the local coastal upwelling.  相似文献   

An index for tropical temperate troughs over southern Africa     

Satyaban B. Ratna  Swadhin Behera  J. Venkata Ratnam  Keiko Takahashi  Toshio Yamagata 《Climate Dynamics》2013,41(2):421-441

Strong cases of the tropical temperate troughs (TTT) that are responsible for the most of the summer rainfall over subtropical southern Africa are analyzed. An index for identifying the TTT is introduced for the first time using anomalies of outgoing longwave radiation (OLR) and the wind. The TTT is associated with a ridge-trough-ridge wave-like structure in the lower troposphere over southern Africa and the adjoining Indian Ocean. Therefore, the index considers physical processes that occur over southern Africa, adjoining the Atlantic and Indian Oceans to depict the variability of the TTT events. Unusually strong TTT events are identified when the standard deviations of the TTT indices defined by the OLR and wind anomalies in the selected regions are above 1.5 and 0.5 respectively. After applying this criterion and filtering out consecutive events, 55 TTT events are identified during the study period of December–January–February seasons from 1980–1981 to 2009–2010. From the composite analyses of those 55 events, it is found that the TTTs evolve with suppressed (enhanced) convection over the southwest Indian Ocean adjacent to Madagascar (southern Africa). The suppressed convection is, in turn, found to be associated with the enhanced convection around Sumatra in the southeast tropical Indian Ocean. This may explain why more TTT events occur in La Niña years as compared to El Niño years. Time evolution of the canonical TTT event shows that it starts 3 days prior to the mature phase of the event, suggesting possible predictability. After reaching a matured state, the system moves east toward the Indian Ocean and decays within the subsequent couple of days. In addition, the intertropical convergence zone (ITCZ) structure changes over Southern Africa/Madagascar during the TTT event and remains similar to climatology over other regions. The results indicate that the continental part of the ITCZ intensifies prior to the TTT event and then spreads southward following the mid-latitude influence during and after the event.  相似文献   

Clashing strategic cultures and climate policy     

Megan Ceronsky  Cameron Hepburn  Michael Obersteiner  Yoshiki Yamagata 《Climate Policy》2013,13(4):347-357

Abstract

Kagan (2002) argues that the different responses of Europeans and Americans to major strategic and international challenges is not simply due to differences in the current administrations, but rather results from (i) a power gap and (ii) differing ideologies. This article applies Kagan's theory to climate policy, employing terrorism policy as a point of comparison. We argue that the power gap between Europe and America is unable to explain the differences in climate policy. In contrast, the ideology gap may indeed have some explanatory value. Furthermore, we argue that one additional feature is critical—the costs and benefits imposed by climate change and terrorism prevention, and the process by which such costs and benefits are evaluated, differ between America and Europe.  相似文献   

Poleward propagation of boreal summer intraseasonal oscillations in a coupled model: role of internal processes     

R. S. Ajayamohan  H. Annamalai  Jing-Jia Luo  Jan Hafner  Toshio Yamagata 《Climate Dynamics》2011,37(5-6):851-867

The study compares the simulated poleward migration characteristics of boreal summer intraseasonal oscillations (BSISO) in a suite of coupled ocean?Catmospheric model sensitivity integrations. The sensitivity experiments are designed in such a manner to allow full coupling in specific ocean basins but forced by temporally varying monthly climatological sea surface temperature (SST) adopted from the fully coupled model control runs (ES10). While the local air?Csea interaction is suppressed in the tropical Indian Ocean and allowed in the other oceans in the ESdI run, it is suppressed in the tropical Pacific and allowed in the other oceans in the ESdP run. Our diagnostics show that the basic mean state in precipitation and easterly vertical shear as well as the BSISO properties remain unchanged due to either inclusion or exclusion of local air?Csea interaction. In the presence of realistic easterly vertical shear, the continuous emanation of Rossby waves from the equatorial convection is trapped over the monsoon region that enables the poleward propagation of BSISO anomalies in all the model sensitivity experiments. To explore the internal processes that maintain the tropospheric moisture anomalies ahead of BSISO precipitation anomalies, moisture and moist static energy budgets are performed. In all model experiments, advection of anomalous moisture by climatological winds anchors the moisture anomalies that in turn promote the northward migration of BSISO precipitation. While the results indicate the need for realistic simulation of all aspects of the basic state, our model results need to be taken with caution because in the ECHAM family of coupled models the internal variance at intraseasonal timescales is indeed very high, and therefore local air?Csea interactions may not play a pivotal role.  相似文献   

Characteristic Features of Tin–iron–copper Mineralization in the Anle-Huanggangliang Mining Area, Inner Mongolia, China     

Daizo ISHIYAMA  Ryoji SATO  Toshio MIZUTA  Yohei ISHIKAWA  Jingbin WANG 《Resource Geology》2001,51(4):377-392

Abstract: The Anle Sn‐Cu and Huanggangliang Fe‐Sn deposits have been exploited in the Linxi district, which is located 165 km northwest of Chifeng City in northern China. In this study the formation mechanisms of the tin deposits in the Anle and Huanggangliang mining area were investigated to understand the mechanisms of tin mineralization in northern China. The veins of the Anle deposit are divided into cassiterite–quartz–chlorite veins, chalcopyrite‐bearing quartz veins, cassi–terite–chalcopyrite–bearing quartz veins and sphalerite‐quartz veins. The sequence of mineralization is tin mineralization (stage I), copper mineralization (stage II), and lead‐zinc mineralization (stage III). The Huanggangliang tin deposit consists of magnetite skarn orebodies and many cassiterite‐bearing feldspar–fluorite veins and veinlets cutting the magnetite orebodies. The fluid inclusions in quartz and fluorite in ores from the Anle and Huanggangliang tin deposits are divided into two‐phase fluid inclusions, vapor‐rich fluid inclusions and poly‐phase fluid inclusions. The final homogenization temperatures of fluid inclusions of quartz in the ores of the Anle deposit and fluorite of tin‐bearing feldspar veins in the Huanggangliang tin deposit range from 195 to 425C and from 215 to 450C, respectively. The fluids responsible for the Anle and Huanggangliang tin deposits were of very high temperature and NaCl‐rich ones containing K, Ca, Al, Si, Ti, Fe and Cl in addition to ore metals such as Sn and Cu. The temperature and chemical composition of fluid in fluid inclusions of igneous rocks in the mining area are very similar to those of fluid in fluid inclusions in the ores of these deposits. The fluid for these ore deposits had a close relation with the fluid coexisting with melt of Late Jurassic granitic rocks in this mining area. Salinities of fluid inclusions from these ore deposits and granitic rocks in the mining area were estimated to range from 35 to 50 wt % NaCl equivalent. Based on arsenopy‐rite geothermometry and fluid inclusion studies, a fluid containing 40 wt% NaCl (eq.) could be formed by phase separation of fluid having 6 wt% NaCl (eq.) at a temperature of 420 to 500C and a pressure of 0.3 to 0.4 kb. The temperatures and pressures presented above indicate an NaCl‐rich magmatic fluid derived from granitic melt that had intruded into a shallow level of crust caused the Sn–Fe–Cu mineralization of the mining area. The geological relationship between these ore deposits and granitic bodies around the ore deposits, and the similarity of fluids forming these ore deposits and coexisting with granitic melt, suggest that these ore deposits were formed by the activity of fluid derived from granitic melt in Late Jurassic age.  相似文献