large scale dynamics of the mesosphere and lower thermosphere an analysis using the extended canadian middle atmosphere model the extended canadian middle atmosphere model is used to investigate the large scale dynamics of the mesosphere and lower thermosphere mlt it is shown that the 4 day wave is substantially amplified in southern polar winter in the presence of instabilities arising from strong vertical shears in the mlt zonal mean zonal winds brought about by parameterized nonorographic gravity wave drag a weaker 4 day wave in northern polar winter is attributed to the weaker wind shears that result from weaker parameterized wave drag the 2 day wave also exhibits a strong dependence on zonal wind shears in agreement with previous modeling studies in the equatorial upper mesosphere the migrating diurnal tide provides most of the resolved westward wave forcing which varies semiannually in conjunction with the tide itself resolved forcing by eastward traveling disturbances is dominated by smaller scales nonmigrating tides and other planetary scale waves play only a minor role in the zonal mean zonal momentum budget in the tropics at these heights resolved waves are shown to play a significant role in the zonal mean meridional momentum budget in the mlt impacting significantly on gradient wind balance balance fails at low latitudes as a result of a strong reynolds stress associated with the migrating diurnal tide an effect which is most pronounced at equinox when the tide is strongest resolved and parameterized waves account for most of the imbalance at higher latitudes in summer this results in the gradient wind underestimating the actual eastward wind reversal by up to 40

inversion of complex body waves we have developed a method that inverts seismic body waves to determine the mechanism and rupture pattern of earthquakes the rupture pattern is represented as a sequence of subevents distributed on the fault plane this method is an extension of our earlier method in which the subevent mechanisms were fixed in the new method the subevent mechanisms are determined from the data and are allowed to vary during the sequence when subevent mechanisms are allowed to vary however the inversion often becomes unstable because of the complex trade offs between the mechanism the timing and the location of the subevents many different subevent sequences can explain the same data equally well and it is important to determine the range of allowable solutions some constraints must be imposed on the solution to stabilize the inversion we have developed a procedure to explore the range of allowable solutions and appropriate constraints in this procedure a network of grid points is constructed on the i plane where and i are respectively the onset time and the distance from the epicenter of a subevent the best fit subevent is determined at all grid points then the correlation is computed between the synthetic waveform for each subevent and the observed waveform the correlation as a function of and i and the best fit mechanisms computed at each i grid point depict the character of allowable solutions and facilitate a decision on the appropriate constraints to be imposed on the solution the method is illustrated using the data for the 1976 guatemala earthquake

a comparative study of mesospheric solar tides observed at adelaide and kyoto observations of winds in the upper mesosphere and lower thermosphere 80 110 km measured at the geographically conjugate locations of adelaide 35 s 138 e and kyoto 35 n 136 e in the years 1983 1985 are used to study hemispheric differences in the behavior of the solar diurnal and semidiurnal tides comparisons of the tidal amplitude and phase parameters show that there are distinct hemispheric asymmetries the diurnal tide is usually much larger in amplitude at adelaide and it exhibits a height structure which indicates a relatively larger contribution from the propagating 1 1 mode at this location comparisons with recent theoretical work suggest that while some of these differences may be due to hemispheric differences in forcing and or to interactions between the tide and the background winds in the middle atmosphere it is also possible that they indicate weaker turbulent induced damping of the 1 1 mode in the southern hemisphere the characteristics of the semidiurnal tide agree well at the two locations it is found that the zonal wind components are about 6 hours out of phase and the meridional components are in phase in all seasons this behavior taken together with the inferred vertical wavelengths which are greater than 100 km in summer and 45 80 km in winter indicates that there are strong contributions from the antisymmetric 2 3 mode probably forced by wind induced mode coupling from the symmetric 2 2 mode the phases near 90 km exhibit a bimodal or two state behavior with the phases being almost constant in one state near the june solstice and then making a rapid transition of between 4 and 6 hours to another phase stable state near the december solstice the duration of the state near the june solstice 6 8 months is longer than the duration of the state at the december solstice although the agreement between recent models and our observations is relatively good better agreement may be expected when more realistic wind distributions especially for the southern hemisphere are used in the models

observations of the two day wave in windii data during january 1993 analysis of variations in the zonal and meridional wind observed by the wind imaging interferometer windii on the upper atmosphere research satellite uars from january 19 to january 31 1993 reveal a strong quasi two day wave qtdw signature these observations were taken during the time of peak amplitude for this qtdw event wu et al 1993 and extend from 41 n to 71 s daytime and nighttime results are presented for the altitude ranges 90 to 150 km and 90 to 102 km respectively the qtdw is found to occur with significant amplitude well into the winter hemisphere and lower thermosphere its latitudinal structure below 103 km is found to be similar to that calculated by hagan et al 1993 using the background wind profile of portniagin and solov yeva 1992 the peak amplitude occurs at 94 km with a magnitude of about 50 m s for the meridional wind and 30 m s for the zonal wind mean profiles from the summer midlatitudes show a secondary maximum in amplitude near 115 km and a node at 130 km above 130 km the wave appears to be evanescent

momentum flux spectrum of convectively forced internal gravity waves and its application to gravity wave drag parameterization part ii impacts in a gcm waccm abstract impacts of a spectral parameterization of gravity wave drag gwd induced by cumulus convection gwdc in the ncar whole atmosphere community climate model waccm1b are investigated in the spectral gwdc parameterization reference wave momentum flux spectrum is launched at cloud top and analytically calculated based on the physical properties of convection and the large scale flow the cloud top wave momentum flux is strong mainly in the tropics and midlatitude storm track regions and exhibits anisotropy and spatiotemporal variability the anisotropy and variability are determined by the distributions and variations of convective activities the moving speed of convection and horizontal wind and stability in convection regions zonal mean zonal gwdc has a maximum of 13 27 37 50 m s 1 day 1 in the mesosphere in january july impacts of gwdc on zonal wind appear mainly in the low to midlatitudes of the upper stratosphere and mesosphere in these regions biases of zonal wind with respect to

source investigation of a small event using empirical green s functions and simulated annealing we propose a two step inversion of three component seismograms that 1 recovers the far field source time function at each station and 2 estimates the distribution of co seismic slip on the fault plane for small earthquakes magnitude 3 to 4 the empirical green s function egf method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path site and instrumental effects from the main event signal the deconvolution between the two earthquakes is an unstable procedure we have therefore developed a simulated annealing technique to recover a stable and positive source time function stf in the time domain at each station with an estimation of uncertainties given a good azimuthal coverage we can obtain information on the directivity effect as well as on the rupture process we propose an inversion method by simulated annealing using the stf to recover the distribution of slip on the fault plane with a constant rupture velocity model this method permits estimation of physical quantities on the fault plane as well as possible identification of the real fault plane we apply this two step procedure for an event of magnitude 3 recorded in the gulf of corinth in august 1991 a nearby event of magnitude 2 provides us with empirical green s functions for each station we estimate an active fault area of 0 02 to 0 15 km2 and deduce a stress drop value of 1 to 30 bar and an average slip of 0 1 to 1 6 cm the selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half graben

analysis of the reliability and resolution of the earthquake source history inferred from waveforms taking the chi chi earthquake as an example summary in the inversion of observed waveforms to obtain the earthquake source faulting process some uncertain factors may influence the reliability of the results several numerical tests have been analysed in this paper in order to understand the effect of some of these uncertain parameters such as the number and width of multiple time windows that define the source time function the assumed strike and dip of the fault the velocity structure and noise also the influence of the sampling rate the time window of the data and the distribution of observational stations is studied our research indicates the near source waveform data can adequately recover the source rupture process of the shallow part of the fault however they poorly constrain the slip distribution on the deep part of the fault using near source and teleseismic waveform data jointly can provide a more complete view of the rupture process of the whole fault the result using near fault data is very sensitive to the assumed strike and dip parameters of the fault model a 2 deviation of the presumed strike and a 5 deviation of the presumed dip from that of the true fault will lead to a distorted inversion result the influence of the uncertainty of the seismic velocity structure for the chi chi earthquake might be ignored if we can take the difference between ma et al s model and wu et al s model as the measure of the media uncertainty

correlation between the pacific north american pattern and the eastward propagation of sea surface temperature anomalies in the north pacific lag correlation and singular value decomposition analyses reveal that the pacific north american pattern is correlated with eastward propagating sea surface temperature sst anomalies in the north pacific the eastward movement of sst anomalies is also clear in one point lag correlation maps and longitude time plot of low pass filtered time series the timescale for the sst anomalies to travel across the north pacific is 1 year it is speculated that the propagation is related to the air sea interactions in the north pacific these findings provide observational support for the atmospheric bridge theory proposed by lau and nath 1994

gravity wave interactions with the diurnal propagating tide observations of the diurnal propagating solar tide in the mesosphere and lower thermosphere 50 150 km indicate substantial annual variations in tidal amplitudes observations and modeling efforts suggest that interactions with internal gravity waves can result in significant modifications of the diurnal tide herein a two dimensional 2 d linear steady state tidal model is used in conjunction with hybrid lindzen matsuno gravity wave parameterization to examine the role of gravity wave interactions in the seasonal variability of the diurnal tide in the mlt the analysis indicates that annual variations of the background atmosphere result in significant annual variations in the eddy diffusion and the diurnal harmonic of the momentum flux divergence wave stress produced by breaking gravity waves the seasonal variability in these two mechanisms can explain the observed semi annual variability in the diurnal tide amplitude minima at solstice and maxima at equinox the relative roles of wave stress and diffusion are examined and it is shown that both play important roles in tidal modification by damping its amplitude the effects of varying inputs to the gravity wave parameterization are explored

rupture process of the 1983 japan sea akita oki earthquake using a waveform inversion method the rupture process of the 1983 japan sea earthquake mjma 7 7 is determined by applying waveform inversion of displacement type strong motion records based on the bayesian method the synthesis of the seismic motions from the main shock is made using the records from the two aftershocks mjma 7 1 and mjma 6 1 as empirical green s functions the main shock fault consists of two subfaults the aftershock m 6 1 occurring south of the main shock fault is used for the southern subfault empirical green s function and the aftershock m 7 1 occurring north of the main shock fault is used for the northern subfault in this synthesis it is not necessary to calculate propagation path and local site effects therefore the information about the source process can be extracted from the comparison between observed seismic motions and the synthetic ones without detailed information on the velocity structure along the whole path from the source to each station the inversion for the rupture process is attempted for both a line source model and an area source model in both models the main shock fault surface is divided into several elements whose sizes are determined from the scaling relations between the main shock and its aftershocks the slip displacement and the rupture starting time on each element are estimated as model parameters in the inversion a good convergence is obtained after about 10 iteration steps for the line source case and after 5 iteration steps for the area source case the results are summarized 1 the average rupture velocity is approximately 2 5 km sec for the southern subfault and 2 0 km sec for the northern subfault and 2 there is a large slip displacement near the rupture initiation point and the north edge of the fault thus there is considerable heterogeneity in rupture propagation and slip distribution over the fault plane

time gcm results for the quasi two day wave to understand how the quasi two day wave qtdw interacts with the mean flow and migrating tides we perturbed the lower boundary of the national center for atmospheric research thermosphere ionosphere mesosphere electrodynamics general circulation model time gcm to excite the qtdw results from our initial model runs computed for a perpetual january 15th show wave amplitudes in wind and temperature that are in qualitative agreement with previous model results and observations the wave period and amplitude change over the duration of the model run and it is expected that this results from changes in the zonal mean zonal winds evidence supporting nonlinear interactions between the qtdw and the migrating tides is presented

the seasonal variation of the propagating diurnal tide in the mesosphere and lower thermosphere part i the role of gravity waves and planetary waves abstract the seasonal variation of the propagating diurnal tide in the mesosphere and lower thermosphere is examined using results from a 2 yr simulation of the extended version of the canadian middle atmosphere model cmam the cmam is shown to be able to reproduce not only the observed semiannual amplitude variation of the tide in the lower thermosphere but also more subtle features such as amplitude maxima that are stronger in march april than in september october a 4 to 6 h shift in phase between winter and summer in the northern hemisphere and a weaker seasonal variation of phase in the southern hemisphere part i of this two part series of papers investigates the importance of two of the mechanisms that have been proposed to explain the observed variation of tidal amplitude namely 1 interactions with small scale gravity waves and 2 interactions with planetary scale waves like the quasi 2 day wave analysis of the tidal momentum and thermodynamic budgets shows that the direct effects of param

extended canadian middle atmosphere model zonal mean climatology and physical parameterizations 1 this paper describes the energetics and zonal mean state of the upward extension of the canadian middle atmosphere model which extends from the ground to 210 km the model includes realistic parameterizations of the major physical processes from the ground up to the lower thermosphere and exhibits a broad spectrum of geophysical variability the rationale for the extended model is to examine the nature of the physical and dynamical processes in the mesosphere lower thermosphere mlt region without the artificial effects of an imposed sponge layer which can modify the circulation in an unrealistic manner the zonal mean distributions of temperature and zonal wind are found to be in reasonable agreement with observations in most parts of the model domain below 150 km analysis of the global average energy and momentum budgets reveals a balance between solar extreme ultraviolet heating and molecular diffusion and a thermally direct viscous meridional circulation above 130 km with the viscosity coming from molecular diffusion and ion drag below 70 km radiative equilibrium prevails in the global mean in the mlt region between 70 and 120 km many processes contribute to the global energy budget at solstice there is a thermally indirect meridional circulation driven mainly by parameterized nonorographic gravity wave drag this circulation provides a net global cooling of up to 25 k d 1

scale matters a multi resolution study of the determinants of patterned ground activity in subarctic finland abstract numerical studies of earth surface processes in relation to their environment are one of the central topics in physical geography however collinearity between explanatory variables and spatial autocorrelation can hinder the detection of key environmental correlates underlying response explanatory variables relationships identified by traditional regression methods moreover conclusions about the potential importance of environmental variables have generally made on analysis conducted only at one spatial scale resolution in this study a variation partitioning method provided a framework to obtain new insights into the relative roles of different factors determining patterned ground activity at multiple spatial resolutions the variation in the distribution of the sorted and non sorted patterned ground was decomposed into independent and joint effects of relief soil and spatial variables geographical location based on a multi resolution system of 1 ha 25 ha and 100 ha cells covering in total 100 km 2 of a subarctic landscape in northern finland the independent effects of relief and soil variables captured the largest fraction of the variation in the non sorted patterned ground distribution while relief had a major contribution for sorted patterned ground activity the independent effect of spatial autocorrelation on sorted patterned ground was higher than that on non sorted patterned ground however a considerable amount of variation in the distribution of both patterned ground types was accounted for by the joint effects of explanatory variables and may thus be causally related to two or all three groups of variables our analyses produced often contradicting results at different resolutions consequently this has substantial implications for the study of geomorphological systems since the choice of resolution can have a major effect on the inferences of analyses our results draw attention to the roles of resolution and spatial autocorrelation in the study of geomorphological systems

timescales for dust variability in the greenland ice core project grip ice core in the last 100 000 years the calcium representing dust concentration record of the last 100 000 years from the greenland ice core project grip ice core shows a huge dynamic range factor 100 the relationship between dust concentrations and temperature represented by the oxygen isotope ratio is not a simple one as has often been assumed a rapid alternation factor of 5 10 between low concentrations during the dansgaard oeschger interstadial periods and high levels in colder periods is superimposed on a long term trend encompassing a further factor of 5 10 within climate periods there is only a very weak relationship between ca concentration and temperature previous authors biscaye et al 1997 have suggested that the most likely source for the increased dust is eastern asia for the first time we consider each possible cause of both rapid and slow increases from source to deposition we suggest that to account for the size and rapidity of the fast changes significantly higher wind speed in the source area is required although changes in atmospheric residence time could also play a role for the slower long term variability changes in transport speed or possibly route are probably also involved changes in the size of the source area could give some change on longer time periods the probable importance of changes in source area wind speed almost simultaneous with greenland temperature changes confirms that climatic parameters in high and low latitudes were strongly coupled through the atmosphere during glacial climatic changes this adds to evidence that the atmospheric circulation system underwent almost instantaneous large scale changes during the last glacial period

modelling of active lineaments for predicting a possible earthquake scenario around dehradun garhwal himalaya india abstract the doon valley which consists of the doon gravels and siwalik rocks is surrounded by the main boundary thrust mbt in the north the main frontal thrust mft mohand thrust and bhimgoda thrust in the south the ganga tear fault in the east and the yamuna tear fault in the west lineaments identified from the aerial photographs show that a few of them trend parallel to the himalayan orogen and the others are transverse to it the effect of the neotectonic activity is clearly reflected in the doon valley implying that the doon valley is tectonically unstable today and that the possibility of earthquakes in this region cannot be ruled out the ruptures along identified active lineaments are modelled three dimensionally peak accelerations are obtained by using a methodology based on the semi empirical method of irikura which was used by midorikawa the efficacy of this modelling technique for the himalayan region is established by simulating peak ground accelerations for the uttarkashi earthquake of 20 october 1991 and comparing it with recorded one a seismic study of the region shows that it is prone to moderate as well as great earthquakes the technique of midorikawa has been applied for the preparation of zonation maps in this seismically active region two different zonation maps showing respectively the effects of moderate as well as great earthquakes in this region are presented in this paper the zonation map for moderate earthquakes shows that the region can experience peak ground accelerations around 310 cm s2 this value is close to the peak ground acceleration of 304 cm s2 recorded for the uttarkashi earthquake ms 7 1 at uttarkashi station which is about 60 km north of dehradun the prepared map shows that among all identified active lineaments rupture along four lineaments can cause major damage in the region among these four three lineaments are transverse to the major tectonic zones viz the mft and mbt in the region and one extends along the mft lying south of dehradun the zonation map for great earthquakes m 8 has been prepared by modelling a rupture in the central gap region of the himalaya seven different possibilities of rupture nucleation within the rupture plane are modelled and a zonation map containing three different zones is prepared the prepared zonation map for this worst possible scenario shows that cities like tehri uttarkashi srinagar deoprayag pauri and dehradun fall within zone 1 and can experience ground accelerations bigger than 600 cm s2 the region like hardwar and rishikesh fall within zone 2 and can experience ground accelerations around 300 cm s2 the zonation map shows that tehri can experience accelerations up to 0 84g and this value is comparable with that obtained by other researchers this study confirms that for very long periods the destruction in the region will be dominantly affected by ruptures having longer lengths and least effected by shorter ruptures therefore special consideration about the longer ruptures should be taken while considering effect of earthquake hazards in this region

the northern hemisphere tropospheric response to changes in the gravity wave drag scheme in a perpetual january gcm the sensitivity of a general circulation model climatology to changes in the gravity wave drag scheme is investigated with particular emphasis being placed upon the northern hemisphere wintertime tropospheric circulation it is found that by weakening the upper level drag the ferrel cell circulation intensifies to the north resulting in a greater than 10 k cooling in the region of the arctic tropopause the poleward side of the northern hemisphere sub tropical jet intensifies in the region of 60 70 n and a poleward shift is observed in the eddy meridional momentum flux a large part of the response is zonally asymmetric with the aleutian and icelandic lows becoming shallower with weakened upper level gravity wave drag by comparing the results of a gravity wave drag scheme having a wave breaking vertical profile with those of a simple uniform drag profile scheme it appears that the tropospheric response to weakened upper level drag depends very little upon the choice of vertical drag profile in the troposphere

rupture history and seismotectonics of the 1991 uttarkashi himalaya earthquake abstract the 19 october 1991 uttarkashi india earthquake occurred in the main thrust zone of the himalaya with a moment magnitude of 6 8 this event is characteristic of the present day motion on the thrust fault system we examine this earthquake using different sets of data in order to understand better the faulting process of a major earthquake in the himalayan region firstly the modeling of the teleseismic records indicates that the mechanism is similar to the published cmt and indicates a shallow between 10 and 15 km depth low angle thrust event in the vicinity of the source the earthquake was recorded by a network of accelerometers run by the university of roorkee six three component accelerometers were triggered within a radius of 60 km two of them were very close to the surface projection of the fault forward modeling of those records shows that the rupture propagated toward the west this forward modeling gives us the possibility to confirm the epicenter location and to evaluate the timing of the accelerograms the accelerogram records are inverted to obtain the distribution of slip on the fault plane the results show a complex rupture process the slip maxima 1 5 m occurred 10 km west and 15 km southwest of the hypocenter the slip source function obtained with near field data is similar to the function obtained from teleseismic records and shows a low moment release at the beginning of the rupture and a maximum rate of moment release 4 seconds after the relation between the slip distribution obtained by inversion isoseimals mapped faults and the aftershocks location is then discussed and we finally propose a seismotectonic interpretation of this earthquake the uttarkashi earthquake probably occurred along the detachment surface which coincides with the upper surface of the subducting indian lithosphere this detachment surface is gently dipping under the lesser himalaya and south of the vaikrita thrust the vaikrita thrust marks the line separating the very shallow dipping detachment along which big earthquakes like the uttarkashi earthquake could occur from the steeper dipping aseismic basement thrust this observation is important for correctly estimating the seismic hazard in the uttarkashi region

sensitivity of mesospheric mean flow planetary waves and tides to strength of gravity wave drag a global circulation model which extends from the surface to 125 km is used to study how the strength of gravity wave drag affects the dynamics of the mesosphere the strength of gravity wave drag has a strong influence on the zonal mean state of the mesosphere in particular the magnitude and variability of the summer mesopause temperature and the shear on the top of the mesospheric jets this change in the zonal mean state strongly affects the susceptibility of the mesosphere to baroclinic and or barotropic instability and hence the formation of the 2 day wave the 2 day wave in turn interacts nonlinearly with the diurnal tide producing secondary waves and a reduction in amplitude of the diurnal tide previous studies with quasi linear mechanistic tidal models have captured some semiannual variation in tidal amplitude through direct interactions between tides and gravity waves our fully nonlinear global circulation model results support an alternative explanation in terms of interactions between planetary waves and tides

the rupture process and asperity distribution of three great earthquakes from long period diffracted p waves moment tensor source mechanisms for the largest twenty nine shallow earthquakes from 1981 m s 6 5 have been determined using long period rayleigh wave spectra and p wave first motions recorded by the international deployment of accelerograms ida global digital seismic network gdsn and world wide standardised seismic network wwssn these are compared with the centroid moment tensor mechanisms presented by dziewonski and woodhouse which were obtained in a contrasting manner in most cases the source is well represented by a single double couple mechanisms with good agreement of scalar moment and orientation between the two studies apparent mismatches for other cases may be interpreted in terms of source complexity we emphasise the objectivity of comparing moment tensor elements apparent discrepancies in the moment and orientation of double couple fault plane solutions are usually due to differences in the mxz and myz moment tensor elements which are poorly constrained by long period data

shallow dynamic overshoot and energetic deep rupture in the 2011 mw 9 0 tohoku oki earthquake strong spatial variation of rupture characteristics in the moment magnitude mw 9 0 tohoku oki megathrust earthquake controlled both the strength of shaking and the size of the tsunami that followed finite source imaging reveals that the rupture consisted of a small initial phase deep rupture for up to 40 seconds extensive shallow rupture at 60 to 70 seconds and continuing deep rupture lasting more than 100 seconds a combination of a shallow dipping fault and a compliant hanging wall may have enabled large shallow slip near the trench normal faulting aftershocks in the area of high slip suggest dynamic overshoot on the fault despite prodigious total slip shallower parts of the rupture weakly radiated at high frequencies whereas deeper parts of the rupture radiated strongly at high frequencies

comparison of seismic waveform inversion results for the rupture history of a finite fault application to the 1986 north palm springs california earthquake the july 8 1986 north palm springs earthquake is used as a basis for comparison of several different approaches to the solution for the rupture history of a finite fault the inversion of different waveform data is considered both teleseismic p waveforms and local strong ground motion records linear parametrizations for slip amplitude are compared with nonlinear parametrizations for both slip amplitude and rupture time inversions using both synthetic and empirical green s functions are considered in general accurate green s functions are more readily calculable for the teleseismic problem where simple ray theory and flat layered velocity structures are usually sufficient however uncertainties in the variation in t with frequency most limit the resolution of teleseismic inversions a set of empirical green s functions that are well recorded at teleseismic distances could avoid the uncertainties in attenuation in the inversion of strong motion data the accurate calculation of propagation path effects other than attenuation effects is the limiting factor in the resolution of source parameters the assumption of a laterally homogeneous velocity structure is usually not a good one and the use of empirical green s functions is desirable considering the parametrization of the problem any degree of fault rupture complexity can be described in terms of a linear parametrization for slip amplitudes however a nonlinear parametrization for rupture times and slip amplitudes can have a distinct advantage over a simple linear one by limiting the number of unknown parameters regardless of the choice of data or the type of parametrization the model or solution will be affected by the choice of minimization norm and the type of stabilization used

dynamic rupture analysis inversion for the source process of the 1990 izu oshima japan earthquake m 6 5 a waveform inversion has been applied to strong motion data using a dynamic shear crack model we studied the 1990 izu oshima earthquake mjma 6 5 which has vertical strike slip faulting with unilateral rupture propagation the inversion has two steps a waveform inversion and a crack inversion that are applied iteratively a waveform inversion is used to determine the distribution of rupture starting times and slip dislocations using the slip functions calculated by the initial crack model or by previous crack inversion a crack inversion is used to calculate dynamic crack propagation that explains the results of the above inversion in this step we use the estimated rupture times as a locking fracture criterion the maximum shear stress attained before a fault segment breaks gives a lower bound estimate of the peak shear strength at each fault segment then the dynamic stress drop distribution is estimated from the slip distribution obtained from waveform inversion assuming a dynamic crack model from the results we determine the rise time distribution and the distribution of a dimensionless stress ratio s defined as strength excess stress drop our analysis gives the following picture of the rupture process of the 1990 izu oshima earthquake 1 an asperity type faulting having large slip and high stress drop was detected in the region around the initiation point of rupture 2 south of the asperity zone barrier type faulting characterized by incoherent propagating rupture small slip long rise time and high strength excess was detected this zone corresponds to the intersection of the fault with the 1978 earthquake mjma 7 0

dynamic rupture processes on a dipping fault and estimates of stress drop and strength excess from the results of waveform inversion summary the dynamic rupture processes on a dipping fault in a horizontally layered medium are investigated on 3 d spontaneous shear crack models the wave equations for the 3 d space are solved numerically by a finite difference method under appropriate boundary conditions the displacement components right above and right below the dipping fault plane have been obtained by solving the boundary conditions across the fault it was found that the final fault slip on the hanging wall side is appreciably larger than that on the foot wall side and that their amount is quite sensitive to the fault depth and the elastic heterogeneities of the medium particularly of the existence of low velocity surface layers on the basis of the above model the dynamic rupture process of the 1961 kita mino earthquake in central japan has been investigated with constraints of the fault parameters derived from the results of waveform inversion analysis for this purpose a locking fracture criterion was introduced for rupture propagation from which a lower bound of the maximum shear stress before the rupture and hence of the strength excess has been estimated the dynamic stress drop has also been evaluated by repeated iterations to minimize the difference between the dynamic and kinematic fault slips the results revealed quite heterogeneous distribution of strength excess indicating large values near the sw section of the fault at shallow depths and in the central bottom section in contrast to small values at mid depths in the ne section around the rupture nucleus zone the distribution of stress drop was also found to be heterogeneous with the existence of high stress drop zones in the ne shallow and bottom sections and a low stress drop zone at most of mid depths except in the sw section these results suggest that the dynamic rupture of this earthquake initiates at a small nucleus zone with low strength excess and small stress drop then broke the shallower and deeper fault sections with moderate strength excess and high stress drop and finally ruptured barrier zones with high strength excess these zones may be interpreted as weak and strong asperities and barriers respectively

rising slope of moment rate functions the 1989 earthquakes off east coast of honshu a series of large earthquakes occurred off east coast of honshu japan during the period of october to november 1989 mechanisms were studied for the foreshock m 6 5 the mainshock m 7 1 and the two aftershocks m 5 2 and 5 3 using the sts i broadband records at station inu about 700km away from the epicenters the lowpass filtered seismograms of these earthquakes show mutually very similar waveforms in good agreement with those calculated for a known mechanism by an extended reflectivity method the mainshock is 10 times as large as the foreshock and 2000 times as large as the aftershocks in terms of seismic moment assuming that the mainshock or foreshock is a collection of some elementary sources at a focal point and that aftershock seismograms represent the responses to such elementary sources the broadband seismogram of the mainshock or foreshock for the first 60s was deconvolved with that of an aftershock to obtain the moment rate function the moment rate functions of the mainshock and foreshock are both triangle like with durations of 14 and 9s respectively their rising slopes are mutually very different a greater slope is associated with the greater shock this and the similar results for other pairs of earthquakes suggest that the rising slope of moment rate function is causaly related to the eventual size of seismic rupture

analysis and interpretation of the heterogeneous rupture process application of the empirical green s function method and nonlinear inversion technique to large earthquakes abstract heterogeneous rupture processes were revealed using an empirical green s function egf method and a nonlinear inversion technique three events the 1980 east izu peninsula earthquake the 1983 japan sea earthquake and the 1987 east chiba earthquake are studied kernel events used as empirical green s functions were selected among the foreshocks and aftershocks the dimension and orientation of the fault plane were evaluated from the aftershock distribution and main shock fault plane solution a moment ratio between the main shock and the kernel event was estimated in order to identify the number of subfaults strong motion waveforms recorded at near distances were inverted to estimate moment release rupture time and rise time at each subfault with the stress drop calculated from the moment release and rise time distributions of these faulting parameters can be related to the aftershock distribution and geological structures general features are that the moment release was large where the aftershock activity was low that the irregular rupture propagations detected here were related to geological structures such as a ridge junction seamount chains and other tectonic plate boundaries and that stress drop distributions show some relation to the pattern of aftershocks an exception is the small moment release at the bending corner of the 1983 japan sea earthquake where the aftershock activity was low with this method it is possible to extract not only the detailed source history but also information on geophysical structures

observed temperature two day wave and its relatives near the stratopause the two day wave is observed in the upper atmosphere research satellitemicrowave limb sounder temperature data around 40 58 km between december 1991 and september 1994 the two day wave temperature signature is most significant after each solstice when the derived easterly winds near the stratopause extend across the equator to at least 208 latitude in the winter hemisphere and the zonal mean winds near the equator are inertially unstable with observed inertial instability disturbances the observed two day wave consists of a 2 0 day period zonal wavenumber 3 and a 1 8 day period zonal wavenumber 4 component named 3 2 0 and 4 1 8 respectively the 3 2 0 component is dominant during two of the three available austral summers but its amplitude is much weaker than the 4 1 8 component during the two observed boreal summers during the austral summers correspondence between amplification of the two day wave temperature signatures regions of reversed potential vorticity gradient due to meridional curvature of the zonal mean flow and the critical lines for the 3 2 0 and 4 1 8 modes suggest barotropic instability as a source of both wave components momentum redistribution by observed inertial instability appears to barotropically destabilize the equatorward flank of the easterly jet where the wave components subsequently grow during the boreal summers the 4 1 8 component appears to be excited by instability that is associated with vertical shear and curvature of the flow seated above the observational domain the boreal 3 2 0 mode appears unrelated to the zonal flow instability within the observational domain and may reflect a normal mode like response during these periods

intraseasonal and interannual variability of the quasi 2 day wave in the northern hemisphere summer mesosphere this study uses global synoptic meteorological fields from a high altitude data assimilation system to investigate the quasi 2 day wave q2dw and migrating diurnal tide during the northern hemisphere nh summers of 2007 2009 by applying a two dimensional fast fourier transform to meridional wind and temperature fields we identify q2dw source regions and diagnose propagation of q2dw activity into the upper mesosphere and lower thermosphere we find that the q2dw in nh summer is composed primarily of westward propagating zonal wave number 3 and wave number 4 components that originate within baroclinically unstable regions along the equatorward flank of the summer midlatitude easterly jet the amplitude of the wave number 3 q2dw tends to peak in july while the amplitude of the wave number 4 q2dw tends to peak in late june and again in early august the seasonal mean q2dw amplitudes are largest in 2009 when the amplitude of the migrating diurnal tide in the upper mesosphere near 30 n was relatively weak however there is no evidence of rapid amplification of the q2dw via nonlinear interaction with the diurnal tide instead variations of q2dw amplitudes during nh summer appear to be linked to variations in the strength and location of the mesospheric easterly jet from one summer to the next with a stronger jet producing larger q2dw amplitudes linear instability model calculations based on the assimilated wind fields indicate that the fastest growing modes are zonal wave numbers 3 and 4 with periods near 2 days that originate in the vicinity of the easterly jet

a three dimensional mountain wave imaged in satellite radiance throughout the stratosphere evidence of the effects of directional wind shear swath scanned thermal radiances from the advanced microwave sounding units amsu a on the noaa and aqua satellites are used to image the horizontal temperature structure of a long wavelength mountain wave that formed over southern scandinavia on 14 january 2003 data from all six stratospheric channels show this wave propagating through the full depth of the stratosphere in channels 9 11 altitudes 20 90 hpa the imaged wave has a phase structure broadly consistent with a stationary wave radiated by the southeastward flow over and above the quasi elliptical terrain of southern norway channel 12 radiances at 10 hpa however show a remarkable abrupt change in imaged wave structure the horizontal wavelength contracts phase lines rotate anticlockwise by 30 40 and peak activity migrates to the south to lie over denmark and northern germany similar structure persists in channels 13 and 14 at 2 5 hpa with progessively increasing radiance amplitudes these features are stable over the 10 hours of amsu a measurements from five separate overpasses and are validated against independent radiances acquired by the atmospheric infrared sounder airs and retrieved airs amsu a temperature profiles from the aqua overpass this change at the channel 11 12 interface coincides with an onset of anticlockwise rotation backing and intensification of background stratospheric winds with height fourier ray and spatial ray modelling incorporating these directionally sheared winds and simplified orographic forcing reproduce the salient features of the observations but only after wave induced temperature perturbations have been converted to channel radiances using a forward model the differential visibility of various components of this three dimensional mountain wave to the amsu a channel weighting functions has a first order impact on the observations at heights above 10 hpa once that is factored in the combined observations and modelling provide direct experimental support for the shutts model s predictions of how backing wind vectors affect the vertical evolution of three dimensional mountain waves implications of these observations for orographic gravity wave drag parametrization are briefly discussed copyright 2007 royal meteorological society

coordinated investigation of summer time mid latitude descending e layer es perturbations using na lidar ionosonde and meteor wind radar observations over logan utah 41 7 n 111 8 w 1 it is well known that there is a strong correlation between the formation of a descending sporadic e layer es and the occurrence of large upper atmospheric zonal wind shears most likely driven by solar thermal tides and or gravity waves we present new results of es perturbation events captured between 13 and 17 july 2011 ut days 194 198 as part of a coordinated campaign using a wind temperature na lidar at utah state university 41 7on 111 8 w and a canadian advanced digital ionosonde cadi scientific instrumentation ltd saskatoon saskatchewan canada and skiymet meteor wind radar both located at nearby bear lake observatory 41 9 n 111 4 w during this period the cadi detected strong descending es on 2 days 195 and 197 when large modulations of the top side mesospheric na layer occurred in synchronism with strong oscillations in the ionosonde e region echoes a weakening in the descending e layer echoes was observed on the other 2 days 196 and 198 coincident with a large reduction in the zonal diurnal and semidiurnal amplitudes above 95 km both tidal components were found to have comparable contributions to the total zonal wind shear that was critical for es formation and its downward propagation further investigation indicates that the weakening tidal amplitudes and the occurrence of the es events were also influenced by a strong quasi two day period modulation suggesting significant quasi two day wave qtdw interactions with the tides indeed a nonlinear wave wave interaction induced 16 hour period child wave was also detected with amplitude comparable to that of the prevailing tides these interaction processes and their associated effects are consistent with earlier thermosphere ionosphere mesosphere electrodynamics general circulation model studies of nonlinear interactions between the migrating tidal waves and the qtdw and were probably responsible for the observed damping of the tidal amplitudes resulting in the disruption of the es

mapping of the high frequency source radiation for the loma prieta earthquake california waveform inversion has been successfully applied to map the spatial distribution of earthquake fault slip from strong motion records but its application has been restricted to a relatively lower frequency range of the data 5 hz generation is incoherent between neighboring patches on the fault we formulated an analytical model for calculating the envelope time history of the squared high frequency displacement seismogram this envelope time history is found through the line integration of the squared energy radiation intensity over the isochron using a recursive stochastic inversion scheme we determined both the lower frequency slip distributions from the low pass filtered strong motion displacement data and the higher frequency energy radiation intensity from the envelope of the squared high frequency displacement data on an extended fault for the 1989 loma prieta earthquake respectively our resulting lower frequency slip model shows that there are two main slip patches and little slip directly updip of the hypocenter one major slip zone located 10 km northwest of the hypocenter with predominant dip slip and the other located 12 km southeast of the hypocenter with predominant oblique slip our preliminary results of the high frequency energy radiation intensity inversion indicate three large high frequency sources and few smaller ones located on the outer periphery of the large slip zones this suggests that the high frequency energy sources are located along or near the boundaries of large slip zones which is consistent with the theoretical consideration that high frequencies are primarily generated from the rupture stopping areas or places with large slip variation

simultaneous inversion of waveform and geodetic data for the rupture process of the 1984 naganoken seibu japan earthquake summary a simultaneous inversion method of strong motion seismograms and geodetic data for inferring earthquake rupture process was developed we assumed a source model composed of many point sources of which slip vectors and rupture times are unknown parameters in addition to observed data we introduced smoothing constraints to determine many parameters without instability since the leastsquares solution depends on the relative weights of the smoothing constraints to a data misfit function we determined their values by minimizing akaike s bayesian information criterion abic by using abic we can objectively choose a unique model this inversion method was applied to the naganoken seibu japan earthquake of 1984 dividing a 15 x 10 km fault into 150 subfaults we obtained the minimum abic model with a seismic moment of 2 7 x los dyne cm this model has a slip peak in the eastern and shallow area of the assumed fault the high slip area is just under the areas where thrown out boulders were found by umeda et al 1987 and iio yoshioka 1988 comparing the aftershock distribution with the estimated slip distribution we found that the high slip area had few aftershocks this suggests that aftershocks are caused by the stress generated in the marginal region around the main shock slip the inversion analysis also revealed that there was a slight change in slip angle during the earthquake the high slip area has some amount of dip slip component while the slips in the other areas are almost purely right lateral

deep learning of aftershock patterns following large earthquakes aftershocks are a response to changes in stress generated by large earthquakes and represent the most common observations of the triggering of earthquakes the maximum magnitude of aftershocks and their temporal decay are well described by empirical laws such as bath s law1 and omori s law2 but explaining and forecasting the spatial distribution of aftershocks is more difficult coulomb failure stress change3 is perhaps the most widely used criterion to explain the spatial distributions of aftershocks4 8 but its applicability has been disputed9 11 here we use a deep learning approach to identify a static stress based criterion that forecasts aftershock locations without prior assumptions about fault orientation we show that a neural network trained on more than 131 000 mainshock aftershock pairs can predict the locations of aftershocks in an independent test dataset of more than 30 000 mainshock aftershock pairs more accurately area under curve of 0 849 than can classic coulomb failure stress change area under curve of 0 583 we find that the learned aftershock pattern is physically interpretable the maximum change in shear stress the von mises yield criterion a scaled version of the second invariant of the deviatoric stress change tensor and the sum of the absolute values of the independent components of the stress change tensor each explain more than 98 per cent of the variance in the neural network prediction this machine learning driven insight provides improved forecasts of aftershock locations and identifies physical quantities that may control earthquake triggering during the most active part of the seismic cycle neural networks trained on data from about 130 000 aftershocks from around 100 large earthquakes improve predictions of the spatial distribution of aftershocks and suggest physical quantities that may control earthquake triggering

diurnal variations of temperature and winds inferred from timed and uars measurements 1 owing to their sampling properties in local solar time both the timed and uars spacecrafts provide the opportunity to quantitatively study the mesospheric thermal tides the information though is limited because it takes the timed instruments 2 months to cover the full range of 24 hours almost twice as long as uars we present here the results of an analysis in which the seasonal and global scale variations of amplitude and phase for the migrating diurnal tide are retrieved from atmospheric temperature and wind measurements on the timed saber tidi and uars mls hrdi satellites comparisons are made with results from the numerical spectral model nsm which is characterized in part by the parameterization of small scale gravity waves at 60 km and above the amplitude and phase variations of the diurnal tide from the nsm agree generally well with those derived from the saber temperatures in the saber temperature amplitudes at around 60 km there are evident asymmetries in latitude with respect to the equator which are not apparent in the mls data near 55 km or in the nsm the asymmetry in the saber temperatures disappears above 60 km for the meridional winds of the migrating diurnal tide derived from tidi the phase and relative amplitude variation with latitude and season are in qualitative agreement with the hrdi measurements and the nsm the amplitudes from tidi winds are systematically smaller but this may in part be due to interannual variations however the results for the zonal winds from tidi more often do not agree well with those from hrdi or the nsm

the character of polar tidal signatures in the extended canadian middle atmosphere model the characteristics of the diurnal semidiurnal and terdiurnal tides zonal wave numbers 5 to 5 in temperature and zonal wind in the polar mesosphere and lower thermosphere region as simulated by the extended canadian middle atmosphere model are examined the most significant diurnal semidiurnal and terdiurnal tides in the polar regions are ds0 dw1 and de1 sw3 sw2 sw1 ss0 se1 and se2 and tw3 ts0 and tw1 respectively and their latitudinal structures seasonal variations and hemispheric asymmetries noted of these components ds0 tw1 ts0 and tw3 exhibit a seasonally symmetric variation with both hemispheres strengthening simultaneously on the other hand dw1 strengthens asymmetrically so that when one hemisphere is strong the other is weak the remainder show no seasonal tendency but vacillate on shorter than seasonal time scales in a symmetric or antisymmetric manner at different times of the year global scale correlations of the amplitudes of the migrating tides dw1 sw2 and the stationary planetary wave 1 and the assumed child nonmigrating tides are also examined the results indicate that the correlations are highly time scale dependent and the significant correlations seen with the original time series are mainly due to longer term variations 18 days there are no consistent global correlations associated with the short term variations 18 days among these waves

observations of the quasi 2 day wave from the high resolution doppler imager on uars a strong westward traveling oscillation with a period of 2 days and zonal wave number 3 is observed in the mesospheric and lower thermospheric winds from the high resolution doppler imager on the upper atmosphere research satellite the important events happen in january july and september october of which the occurrence in january is the strongest with an amplitude over 60ms 1 detailed analyses for the periods of january 1992 and january 1993 reveal a cause and effect relationship in the wave developing process at 95km the global structures of the wave amplitude and phase are also presented

observations of the 2 day wave in nmc stratospheric analyses abstract observational characteristics of the 2 day wave a westward propagating zonal wave 3 oscillation in the summer subtropical upper stratosphere and mesosphere are studied based on five years of national meteorological center nmc operational stratospheric analyses these data show episodic occurrence of the 2 day wave in the upper stratosphere in january centered near 20 s and july august centered near 20 n these episodes are strongly correlated with observed reversals of the zonal mean potential vorticity gradient near the core of the summer easterly jet consistent with previous suggestions that the 2 day wave is generated by an in situ instability of this jet on the other hand the horizontal and vertical structure of the waves is very similar to that calculated by salby for a global normal mode rossby wave the combination of normal mode structure and instability signature suggests that the 2 day wave is a near resonant mode forced by dynamical instability

short term variability in the migrating diurnal tide caused by interactions with the quasi 2 day wave 1 the migrating diurnal tide is one of the dominant dynamical features in the low latitudes of the earth s mesosphere and lower thermosphere mlt region representing the atmospheric response to the largest component of solar forcing ground based observations of the tide have resolved short term variations attributed to nonlinear interactions between the tide and planetary waves that are also in the region using the ncar thermosphere ionosphere mesosphere electrodynamics general circulation model time gcm we simulate a quasi 2 day wave qtdw event under late january conditions in this case sideband sum and difference child waves are resolved indicating that a nonlinear interaction is occurring between the qtdw and the tide the migrating diurnal tide in the mlt displays local amplitude decreases of 20 40 as well as a shortening of vertical wavelength by roughly 4 km examining the physical mechanisms driving the interaction nonlinear advection is found to result in amplification of the tide in some regions and damping in others manifesting as increased smoothing of the tidal structure when the qtdw is present in the mlt additionally the qtdw also enhances the easterly summer mean wind jet that can also account for changes in tidal amplitude and vertical wavelength we find that qtdw induced background atmosphere changes in time gcm can drive tidal variability at levels greater than nonlinear advection a possibility not previously considered

source scaling relations of interface subduction earthquakes for strong ground motion and tsunami simulation the recording on high resolution broadband seismic networks of several great interface subduction earthquakes during the last decade provide an excellent opportunity to extend source scaling relations to very large magnitudes and to place constraints on the potential range of source parameters for these events at present there is a wide range of uncertainty in the median rupture areas predicted for any given seismic moment by current relationships between magnitude and rupture area for subduction interface earthquakes our goal is to develop an updated set of earthquake source scaling relations that will reduce this current large degree of epistemic uncertainty and improve the accuracy of seismic hazard analysis and the prediction of the strong motion characteristics and tsunamis of future subduction earthquakes to achieve this goal we compiled a database including slip models of interface earthquakes that occurred worldwide with moment magnitudes ranging from m 6 75 9 1 we characterized the seismic sources based on well established criteria to estimate the asperity areas as well as the average slip on the faults and we used these parameters to compute an updated set of magnitude scaling relations of the various characteristics of the fault additionally we followed an alternative approach to quantifying slip models for use in developing characteristic slip models of future earthquakes this involved analyzing the 2d fourier transforms of the slip functions in the compiled database and deriving a wavenumber spectral model of the slip distribution

alleviation of a systematic westerly bias in general circulation and numerical weather prediction models through an orographic gravity wave drag parametrization systematic westerly biases in the northern hemisphere wintertime flow of the meteorological office 15 layer operational model and 11 layer general circulation model are described evidence that the failure to parametrize subgrid scale orographic gravity wave drag may account for such biases is presented this evidence is taken from aircraft studies surface pressure drag measurements and studies of the zonally averaged momentum budget a parametrization scheme is described in which the surface stress is proportional to the near surface wind speed and static stability and to the variance of subgrid scale orography the stress is absorbed in the vertical by considering the influence of such gravity wave activity on static stability and vertical wind shear a richardson number dependent wave breaking formulation is devised and the vertical stress profile determined by a saturation hypothesis whereby the breaking waves are maintained at marginal stability it is shown that wave breaking preferentially occurs in the boundary layer and in the lower stratosphere results from a simple zonally symmetric model show how the adjustment to thermal wind balance with a wave drag in the stratosphere warms polar regions by adiabatic descent and decelerates the mean westerlies in the troposphere the influence of the parametrization scheme on integrations of the 11 layer model is described and found to be generally beneficial in a discussion of the reasons why this problem has only recently emerged it is suggested that the satisfactory northern hemisphere winter circulations of previous coarser general circulation models were due to a compensation implied by underestimating both the surface drag and the horizontal flux of momentum by explicitly resolved large scale eddies

high frequency radiation process during earthquake faulting envelope inversion of acceleration seismograms from the 1993 hokkaido nansei oki japan earthquake we investigate the process of high frequency 1 to 10 hz radiation on the fault plane of the 1993 hokkaido nansei oki japan earthquake mw 7 5 from the envelope inversion of strong motion acceleration seismograms for the analysis empirical green s functions are used because theoretical approach is not available for such high frequencies the source is modeled with two fault planes with different strike angles the rupture process of this earthquake is very complex in terms of high frequency wave generation the rupture which started on the northern fault plane had a delay of about 10 sec or propagated very slowly between the northern and southern fault planes high frequency wave radiation is large at the northern and southern edges of the source region deceleration of rupture is also observed there this is interpreted to be associated with stopping of rupture another high frequency wave radiation area is found at the center of the northern fault plane where discontinuity in the depth distribution of aftershocks suggests an existence of a barrier the areas of high and low frequency wave radiation are not correlated this is considered to result from the complexity of rupture process we cannot distinguish between westward and eastward dip of the southern fault plane because of one sided station distribution

seasonal variations of the diurnal tide induced by gravity wave filtering wind measurements with the uars spacecraft show that between 90 and 100 km altitude the fundamental diurnal tide is modulated semi annually with peak amplitudes of almost 100 m s near equinox more than a factor of two larger than during solstice applying hines doppler spread parameterization dsp for small scale gravity waves gw we present results from a numerical spectral model nsm to describe the seasonal variations in the mean zonal circulation interacting self consistently with the diurnal tide with the dsp the gw momentum source increases the amplitude of the tide at altitudes above 70 km in contrast to the rayleigh friction parameterization that causes damping since the tide absorbs gw momentum it reduces the wave generated semi annual oscillations sao in the mean zonal winds at equatorial latitudes outside the tropics gw momentum is absorbed by the mean zonal winds that peak near the summer and winter solstices as a result the gw reach the upper mesosphere with more momentum during equinox periods and thus cause a seasonal modulation of the tide that is predominately semi annual

rupture nucleations in the 1995 hyogo ken nanbu earthquake and its large aftershocks summary local seismograms of the hyogo ken nanbu earthquake m 7 2 and its largest aftershocks show that a slow rising initial phase precedes the main p phase high dynamic range short period seismograms of the foreshocks mainshock and aftershocks were recorded by the japan meteorological agency jma network and the net hyogo using these seismograms we investigate the initial rupture process of the mainshock and large aftershocks using an empirical green s function method the deconvolved moment rate function of the mainshock shows that it begins with a small pre event mw 4 0 then the moment releases gradually and after about 0 6 s it increases rapidly the first 0 7 s of the deconvolved moment rate functions are inverted to obtain the spatio temporal images of the initial rupture process the results show that at the hypocentral area less than 1 km from the hypocentre slip velocity rises slowly with a low peak in the surrounding area 1 km east and downward from the hypocentre slip velocity rises rapidly with a high peak it is also found that during the initial rupture process the slip velocity pulse grows continuously these observations are in accord with a theoretical model of the transition process from nucleation to high speed dynamic rupture in which slip velocity rises slowly at the nucleation zone we also analyse the initiation and the main dynamic rupture of large aftershocks spatio temporal images of the largest aftershocks m 5 1 and m 4 8 show that there is a clear initial rupture process with a low moment release and after about 0 2 s slip velocity increases rapidly we compare the initial rupture process of the mainshock with those of the large aftershocks m 4 0 5 1 the duration 0 6 s of the initial phase of the mainshock is longer than those 0 02 0 3 s of the large aftershocks furthermore the slip velocity of the initial rupture process of the mainshock is several times larger than those of the two largest aftershocks these results suggest that the mainshock starts in the region where the stress drop and fracture energy is higher

simulation of large scale dynamics in the mesosphere and lower thermosphere with the doppler spread parameterization of gravity waves 2 eddy mixing and the diurnal tide turbulence is thought to be produced by gravity waves that become unstable in the upper atmosphere owing to a relatively short vertical wavelength of the solar diurnal tide turbulent viscous type mixing is expected to be particularly important for tidal dissipation in the mesosphere and lower thermosphere mlt within the doppler spread parameterization dsp of gravity waves the vertical eddy exchange coefficients may be calculated from the gw energy deposition rate the dsp as implemented in the spectral mesosphere lower thermosphere model is shown to provide realistic simulations of both the diurnal tide and turbulent mixing in the mlt in particular the vertical structure and seasonal variations of eddy mixing agree well with in situ observations and independent estimates from global thermal balance the seasonal variations of turbulence are found to be primarily responsible for a strong semiannual variation of the diurnal tide in agreement with recent satellite observations

the 12 july 1993 hokkaido nansei oki japan earthquake effects of source complexity on surface wave radiation abstract we present an analysis of the data dependency of surface wave inversions for source parameters of the 12 july 1993 hokkaido nansei oki japan earthquake this dependency is demonstrated by differences between source mechanisms obtained from different data sets particularly from rayleigh waves or love waves although inversions of various data sets yield unanimously thrust fault mechanisms for the inversion of rayleigh waves the mechanism is characterized by a shallow dipping west plunging nodal plane which is consistent with p wave first motion data and the distribution of aftershocks close to the hypocenter however for the inversion of both rayleigh and love waves the mechanism is characterized by a shallow dipping east plunging nodal plane several numerical experiments were carried out to test a hypothesis that the data dependency is attributed to the effects of source complexity on the surface wave radiation with the source complexity being characterized by tempo spatial changes in focal mechanism and moment release during the rupture process the experiments show that source parameters obtained from inversions of surface waves with periods as long as a few hundred seconds can be largely affected by using various models of source complexity proposed previously for the hokkaido earthquake and the effects on love waves tend to be systematically different from those on rayleigh waves our preferred source model is characterized by significant tempo spatial changes in focal mechanism the experiments based on this model provide a good explanation for the data dependency although our analysis cannot completely reject the alternative explanation that attributes the dependency to the manifestation of errors in modeling the surface wave propagation and or inversion instability the model is consistent with the data of body waves and long period surface waves for a very broad frequency band implying the existence of weak zones with various fault geometries along the eastern margin of the japan sea

source models of two large intraslab earthquakes from broadband strong ground motions we construct source models of two large intraslab earthquakes the 1993 kushiro oki and 1994 hokkaido toho oki earthquakes along the southern kurile hokkaido arc based on strong motion records the source models well explain observed broadband strong ground motions we assume a distribution of rectangular shaped asperities on a fault plane and apply the empirical green9s function method to the synthesis of strong ground motions from the asperities the asperity parameters area and stress drop are estimated by matching the synthetic waveforms displacement velocity and acceleration time histories and spectra to the observed ones our source models have the following characteristics compared with those of inland and plate boundary earthquakes having a comparable seismic moment 1 the total asperity area is much smaller and 2 the stress drop is extremely high about 200 400 mpa these indicate that the large intraslab earthquakes radiate seismic energies from a small area in a short time this rupture process results in extremely strong radiation of short period seismic waves finally we confirm a correlation between outer entire fault size seismic moment and short period level of the acceleration source spectrum and inner fault parameters areas and stress drops of asperities for the two intraslab earthquakes

the quasi two day wave observed in the equatorial middle atmosphere an oscillation with a period near 2 days is found in radar wind measurements made in the upper mesosphere and lower thermosphere at christmas island 2 n 157 w in the central pacific the oscillation is particularly strong in the meridional wind component and seems to be present almost continuously ir the 80 100 km height region although some intermittency is evident wave amplitudes are especially large about one month after the solstices in july august and in january february at times coincident with large 2 day wave events in the summer mesosphere at extratropical latitudes power spectral and complex demodulation techniques are used to study the amplitude and frequency variations of the oscillation as a function of height and time downward phase propagation is found consistent with upward energy propagation and with a vertical wavelength of about 70 km an unusual feature is a shift in wave period from near 50 hours in july august to 48 hours in january february when the largest amplitudes up to 45 m s 1 are observed a 16 hour oscillation is also found at times when the 2 day wave amplitudes are largest the observations are consistent with the 2 day wave being a manifestation of the 3 0 rossby gravity normal mode another weaker oscillation is also found in the meridional wind field at a period near 44 hours this oscillation may be due to the 2 0 normal mode

the summer wind field between 80 and 93 km observed by the mst radar at poker flat alaska 65 n abstract results of an analysis of the summertime wind field between 80 93 km over poker flat alaska are presented the data were obtained using the large but as yet incomplete mst radar at poker flat and cover the periods 17 june 13 july 1979 and 22 june 22 july 1980 zonal and meridional mean wind profiles and tidal components during these periods are examined and are compared with other observations at comparable latitudes as well as with current theoretical profiles while our results agree reasonably closely with other observations and theory some discrepancies exist and are discussed in addition to the expected 8 12 and 24 h components of atmospheric motions we find strong evidence for a 16 h component particularly during the 1979 data period in addition we show examples of the average power spectrum of wind fluctuations from 3 min to 8 days these spectra not only show the expected tidal peaks but show also that the spectral energy density is continuous between the shortest 8 h t

kīholo bay hawai i earthquake sequence of 2006 relationship of the main shock slip with locations and source parameters of aftershocks 1 we study the source process of the k holo bay earthquake mw 6 7 which occurred beneath the northwest part of the island of hawai i on 15 october 2006 and static stress drops of small earthquakes that occurred in 2006 and 2007 around the main shock including aftershocks we relocate the aftershocks to determine the fault plane from the two nodal planes the relocated aftershocks define an e w trending plane that dips to the south in good agreement with one of the nodal planes given by the global centroid moment tensor solution waveform inversion is performed with multiple time windows to investigate the rupture speed and the slip distribution of the main shock waveforms of an aftershock with mw 5 2 are used to calculate empirical green s functions our results indicate that the rupture propagated unilaterally to the west with a rupture speed greater than 3 0 km s 63 of the shear wave velocity this westward rupture is consistent with the fact that aftershocks are distributed predominantly to the west of the main shock epicenter most aftershocks are located on the edge of patches with a large slip or asperities and some also occur inside the patches we also estimate static stress drops of 39 earthquakes 2 5 ml 3 5 that occurred in 2006 and 2007 near the source region of the k holo bay earthquake static stress drops range from 0 12 to 8 6 mpa and aftershocks around large slip patches of the main shock likely to have larger stress drops

effective diffusivity in the middle atmosphere based on general circulation model winds 1 the mixing of a passive tracer in the stratosphere and lower mesosphere is studied on the basis of the effective diffusivity which is obtained in the framework of the tracer based coordinate system this characteristic is proportional to the average diffusion flux over lagrangian contours and inversely proportional to the mean tracer gradient the tracer distribution used in the calculation of the effective diffusivity is obtained after integration of the advection diffusion equation using general circulation model winds and a new numerical advection scheme with small numerical diffusivity using some theoretical and experimental arguments it is shown that the interpretation of the seasonal variability of the effective diffusivity field cannot be done on the basis of the momentary wind field alone but some flow history should be taken into account the climatology of the effective diffusivity for different months is presented up to the lower mesosphere and compared with previous studies in the stratosphere some new features of the effective diffusivity distribution are obtained for example there is a local maximum of the effective diffusivity at midlatitudes of the northern hemisphere of the summer middle stratosphere the effective diffusivity fields in the lower mesosphere show a strong increase of the mean effective diffusivity from the upper stratosphere to the lower mesosphere and the existence of a complex latitudinal structure of the effective diffusivity at mesospheric heights in the lower mesosphere there is a marked interannual variability during the southern hemisphere easterly wind development a possible explanation for the obtained structure is discussed on the basis of in situ rossby wave generation and rossby wave breaking effects

satellite observations of thermospheric tides results from the wind imaging interferometer on uars thermospheric winds measured by the wind imaging interferometer windii on the upper atmosphere research satellite are analyzed for migrating solar tides the data cover a 2 year period commencing february 1992 and are obtained from the atomic oxygen o 1s 557 7 nm emission which provides observations of the 90 to 200 km altitude range during daytime and the 90 to 110 km range at night the subtropical lower thermosphere is dominated by the diurnal propagating tide which exhibits a vertical wavelength of approximately 22 km grows in amplitude up to 95 km and decays rapidly above where molecular diffusion greatly reduces the vertical shears although the phase remains fairly uniform throughout the year a pronounced semiannual oscillation is observed in the diurnal tide amplitude at both 20 n and 20 s the meridional and zonal wind components attain their maximum values at equinox of approximately 70 and 40 m s respectively while the solstitial minima are nearly a factor of 2 smaller at 35 n the diurnal tide semiannual amplitude oscillation is still present in the lower thermosphere but above 100 km it is replaced by an annual cycle with a maximum in july and august this contrasts with 35 s where the july august peak is absent and the semiannual oscillation extends to 110 km at midlatitudes the zonal and meridional winds are of similar magnitude and no significant hemispheric asymmetries in amplitudes are observed in the lower thermosphere the semidiurnal tide amplitude exhibits an annual oscillation with maximum values of 30 to 40 m s occurring in june july near 100 km at 35 n 35 s and the equator a bimodal structure in the seasonal variation of the semidiurnal phase is observed this feature is characterized by rapid equinoctial transitions and is particularly well defined at the equator examination of the equatorial middle thermosphere indicates that the semidiurnal tide attains its maximum amplitude at 140 km and exhibits a vertical wavelength of approximately 60 km these findings indicate the predominance of the antisymmetric 2 3 hough mode in the tropics

the extended canadian middle atmosphere model first results from the extended canadian middle atmosphere model are presented the model extends from the surface to the middle thermosphere and includes relevant physical parameterizations to simulate unresolved gravity waves a non orographic gravity wave drag parameterization has also been implemented results from a two year simulation are discussed and are shown to be in good qualitative agreement with upper atmosphere research satellite wind observations the realistic features which are reproduced include the reversal of the mesospheric jets and the semi annual variation of the migrating diurnal tide in the low latitude mesosphere

inversion for the rupture details of the 1987 east chiba earthquake japan using a fault model based on the distribution of relocated aftershocks the rupture of the 1987 east chiba earthquake was investigated using an empirical green s function method and a nonlinear inversion technique in order to discuss what happened on the fault plane in more detail we relocated the aftershocks as well as the main shock relocations were done using a joint hypocenter method in order to increase the accuracy of the relative locations based on this aftershock distribution we propose a double fault model the main fault strikes n11 w with nearly vertical dip a subsidiary fault is perpendicular to the main fault the results of a source inversion using strong motion waveforms are as follows 1 moment release was large both in the deep northern area of the main fault and on the subsidiary fault where aftershock activity was weak 2 rupture propagated on both main and subsidiary faults simultaneously and rupture velocity decreased at the depth of around 35 km near the upper boundary of the philippine sea plate and 3 stress drop was high on the lower northern corner of the main fault where the rupture terminated and on the deep southern area of main fault where the rupture initiated a negative correlation was found between moment release during the main shock and that during the aftershock sequence but they did not compensate each other quantitatively this suggests the existence of a heterogeneous strain distribution in the fault area either before the main shock or after its aftershock sequence

heterogeneity of the 1980 izu hanto toho oki earthquake rupture process summary the rupture process of the 1980 izu hanto toho oki japan earthquake is examined by inverting strong motion seismic waveform records synthetic seismograms for waveform inversion are calculated using the empirical green s function technique after selecting the optimum small events for the empirical green s function the rupture starting time moment release and rise time at each sub fault are estimated in this procedure the results are as follows in the northern part of the fault surface where the swarm activity occurred the rupture propagation is irregular less moment is released and the rise time is short the local stress drop distribution is also estimated using the results of the inversion it is found that the heterogeneities of the rupture propagation moment release distribution and stress drop distribution have some relation to the geological setting in this area or to the seismic swarm activity

quantification of the gravity wave forcing of the migrating diurnal tide in a gravity wave resolving general circulation model 1 the interaction of gravity waves gws and the migrating diurnal tide are studied in a gw resolving general circulation model gcm by calculating the tidal components of zonal wind accelerations and equivalent rayleigh friction due to tidal induced gw dissipation two 15 day periods for perpetual equinoctial and solstice simulations are analyzed which are performed with the japanese atmospheric general circulation model for upper atmosphere research jaguar high resolution gcm the model can directly simulate gws with horizontal wavelengths greater than about 190 km and thus reproduce the general features of the mean winds and temperatures from the surface to the mesosphere and lower thermosphere mlt the amplitudes of the migrating diurnal tide are successfully simulated during both seasons and the tidal winds affect the altitudes of gw dissipation in the low latitude mlt the tidal component of gw forcing has maximal values of about 15 m s 1 d 1 near the maximal vertical shears of the tidal winds and generally works to shorten the vertical wavelength of the migrating diurnal tide the phase relationship between the tidal winds and the tidal induced gw forcing is not exactly 90 out of phase causing amplification suppression of the tide the gw forcing amplifies the migrating diurnal tide during the equinox while during the solstice it suppresses the tidal winds in the upper mesosphere of both hemispheres this difference in behavior can be attributed to a seasonal variation of the mean zonal winds because combination of the mean and tidal winds affects the altitudes of gw dissipation

source model of the 1995 hyogo ken nanbu earthquake and simulation of near source ground motion the 1995 hyogo ken nanbu earthquake struck the heavily populated kobe and adjacent cities in western japan more than 6400 people were killed and more than 150 000 buildings were destroyed the characteristics of mainshock ground motions in the heavily damaged area are needed to understand how buildings and bridges performed and why they reached failure unfortunately very few strong ground motions were recorded in the heavily damaged area during the mainshock in this study we attempt to estimate mainshock ground motions by using the em pirical green s function method egf method first we assume an initial source model with the asperities based on the rupture process obtained by inversion of strong ground motion records for simplicity we consider each asperity as a sub event with uniform stress drop in a finite extent then the initial model was improved by matching the synthetic and observed ground motions using a trial and error pro cedure the final model consists of three subevents subevent 1 with stress drop of 163 bars under the akashi strait around the rupture starting point subevent 2 with stress drop of 86 bars under the nojima fault in awaji island and subevent 3 with stress drop of 86 bars under kobe finally we estimate strong ground motions using aftershock records at sites where the mainshock was not recorded the near source motions in kobe synthesized with the best fit model are characterized by two large pulses with a duration of 1 to 3 sec the pulses are caused by forward rupture directivity effects from subevents 1 and 3 peak horizontal acceleration and velocity of the synthesized motions at the heavily damaged sites are about 1000 cm sec 2 and 130 crrdsec respectively while those at a rock site in the near source region are about 300 crrdsec 2 and 60 cm sec

inversion of strong ground motion and teleseismic waveform data for the fault rupture history of the 1979 imperial valley california earthquake a least squares point by point inversion of strong ground motion and teleseismic body waves is used to infer the fault rupture history of the 1979 imperial valley california earthquake the imperial fault is represented by a plane embedded in a half space where the elastic properties vary with depth the inversion yields both the spatial and temporal variations in dislocation on the fault plane for both right lateral strike slip and normal dip slip components of motion inversions are run for different fault dips and for both constant and variable rupture velocity models effects of different data sets are also investigated inversions are compared which use the strong ground motions alone the teleseismic body waves alone and simultaneously the strong ground motion and teleseismic records the inversions are stabilized by adding both smoothing and positivity constraints the moment is estimated to be 5 0 10 25 dyne cm and the fault dip 90 5 dislocation in the hypocentral region south of the united states mexican border is relatively small and almost dies out near the border dislocation then increases sharply north of the border to a maximum of about 2 m under interstate 8 dipslip motion is minor compared to strike slip motion and is concentrated in the sediments the best fitting constant rupture velocity is 80 per cent of the local shear wave velocity however there is a suggestion that the rupture front accelerated from the hypocenter northward the 1979 imperial valley earthquake can be characterized as a magnitude 5 earthquake at the hypocenter which then grew into or triggered a magnitude 6 earthquake north of the border

observations of the diurnal tide from space abstract this study presents a climatology of mesospheric and lower thermospheric diurnal tidal winds obtained with the high resolution doppler imager hrdi on the upper atmosphere research satellite uars the observations reveal that although tidal structures are present at all times like the prevailing zonal winds they exhibit significant semiannual as well as other shorter term variations in amplitude results are presented for a period extending over more than one year from november 1991 to july 1993 the 1 1 diurnal tidal amplitude of the meridional component characterized by the value at an altitude of 90 km and a latitude of 22 ranges from a minimum at solstice of less than 20 m s 1 to an equinox maximum of over 70 m s 1 the vertical wavelength and phase of the tide show only slight variations throughout the year with a suggestion of semiannual variations in both

effects of horizontal geometrical spreading on the parameterization of orographic gravity wave drag part i numerical transform solutions abstractnumerical transform solutions for hydrostatic gravity waves generated by both uniform and sheared flow over elliptical obstacles are used to quantify effects of horizontal geometrical spreading on amplitude evolution with height both vertical displacement and steepness amplitudes are considered because of their close connections to drag parameterizations in weather and climate models novel diagnostics quantify the location and value of the largest wavefield amplitudes most likely to break at each altitude these horizontal locations do not stray far from the obstacle peak even at high altitudes resulting vertical profiles of wave amplitude are normalized to remove density and refraction effects thereby quantifying the horizontal geometrical spreading contribution currently absent from parameterizations horizontal geometrical spreading produces monotonic amplitude decreases with height through wave action conservation as waves propagate into progressively larger horizontal areas accumulated

comparison of l1 and l2 norms in a teleseismic waveform inversion for the slip history of the loma prieta california earthquake abstract broadband digital teleseismic body waves p and sh are inverted to obtain the history of slip for the 17 october 1989 loma prieta earthquake comparison of finite fault waveform inversion results using an l 1 norm and an l 2 norm reveal some basic uncertainties in the rupture history that are not usually appreciated in earthquake source studies however the more robust features of the source can be identified that are less dependent on the choice of minimization norm from a consideration of first motion data and forward modeling of body waves the strike and dip of our model fault plane are set at 126 and 67 respectively the waveform inversion results indicate that the majority of the moment release occurs over a fault length of 35 to 40 km and a depth range of 2 to 18 km the hypocenter 18 km depth experienced a relatively small amount of slip the rupture propagated upward and bilaterally with a peak slip of about 3 m occurring at a depth of 10 km the average velocity of the rupture propagation is 2 5 km sec the details of the slip distribution depend on the minimization norm but two main sources can be identified with the larger source above and to the southeast of the hypocenter the rake vector varies considerably over the fault the southern source is mostly strike slip while the northern source is approximately equally partitioned between dip slip and strike slip there is also a trend for more strike slip motion at shallower depths the majority of the moment release occurs in the first 7 5 sec with an integrated moment of 2 6 10 26 dyne cm including smaller peripheral sources that are less well resolved raises the moment estimate to 3 0 10 26 dyne cm

impact of a spectral gravity wave parameterization on the stratosphere in the met office unified model abstract the impact of a parameterized spectrum of gravity waves on the simulation of the stratosphere in the met office unified model um is investigated in the extratropical mesosphere the gravity wave forcing acts against the mean zonal wind and it dominates over the resolved wave forcing in the extratropical stratosphere the gravity wave forcing gives a small acceleration in the direction of the mean zonal wind both summer and winter stratospheric jets have improved maximum strength and tilt with height when the parameterized gravity wave forcing is included although the southern winter jet is still more vertically aligned than in observational analyses the timing of the seasonal breakdown of the southern winter vortex is also improved by the addition of gravity wave forcing in the tropics the most obvious impact is that the model reproduces the quasi biennial oscillation qbo with a realistic mean and range of periods it also reproduces most of the observed asymmetries between the easterl