Long-term seismotectonic influence on the hydrochemical composition of a spring located at Koryaksky-Volcano, Kamchatka: deduced from aggregated earthquake information

The analysis of earthquake-related signals in hydrochemical time series is still a challenging task. Mostly it is unclear how the geometrical and energetic distribution of earthquakes is influencing variation in the hydrochemical composition of monitoring sites, e.g. located close to volcanoes. Past research showed that linear stress-release models alone are not capable to explain sufficiently observed variation in hydrochemical time series due to earthquake activity. A spring located at the base of Koryaksky Volcano, which has shown hydrochemical variation close to 5, major earthquakes, was chosen to analyse relation patterns between hydrochemical variation and seismicity. A possible mechanism, explaining observed hydrochemical variation, that seismic waves trigger an underground water pumping caused by nucleation of gas bubbles in magma was proposed. Consequences are an increase of discharge, gas content in water and changes in the mixing ratios of waters of different genesis. Based on functions of aggregated earthquake information (E) it is herein shown that seismotectonic-triggered processes have a significant influence on the variation of the hydrochemistry of the spring, lasting longer than two decades. At least seven categories of relation patterns between hydrochemical variation and seismotectonic activity E can be identified. A conducted spectral analysis shows that earthquake activity and hydrochemistry share spikes in frequencies. Results prove that the use of functions of transformed aggregated seismic observations is useful to represent the seismotectonic activity for analysing earthquake-related signals in hydrochemical time series.     

The influence of seismotectonics on precursory changes in groundwater composition for the 1995 Kobe earthquake, Japan

It is shown that temporal variation in hydrochemistry may be attributed, in part, to continuous seismotectonic activity occurring before the onset of the 1995 “Kobe” earthquake, Japan, challenging the “one earthquake–one signal” hypothesis with respect to potential precursory signals to this devastating event. A possible continuous seismotectonic influence on chloride and sulphate ion-concentration is evaluated with aggregate earthquake-information by transforming a multivariate earthquake time series (including coordinates and magnitude) into a one-dimensional time series (considering geometric relationships between earthquakes and the well-site). A piecewise analysis of ion-concentration and seismotectonic-activity time series compares trends and change points between the two time series: a positive correlation (before the proposed onset of the preparation stage) is followed by a negative correlation (during the proposed preparation stage) which, in turn, is succeeded by a positive correlation (after the heaviest aftershock sequences). This suggests that seismotectonic processes occurring before the Kobe earthquake dynamically influenced aquifer characteristics, leading to temporal variations in the hydrochemical time series. Accordingly, a dynamic change in the mixing ratios of waters with different hydrochemical characteristics is proposed as a mechanism for explaining observed variation. The research can be extended by considering vectors in the stress field that lead to changes in the aquifer-well system.The influence of seismotectonics on precursory changes in groundwater composition for the 1995 Kobe earthquake, Japan.     

Seismic reversal pattern for the 1999 Chi-Chi, Taiwan, MW 7.6 earthquake

We in this study have calculated the standard normal deviate Z-value to investigate the variations in seismicity patterns in the Taiwan region before and after the Chi-Chi earthquake. We have found that the areas with relatively high seismicity in the eastern Taiwan became abnormally quiet before the Chi-Chi earthquake while the area in the central Taiwan with relatively low seismicity showed unusually active. Such a spatially changing pattern in seismicity strikingly demonstrates the phenomenon of “seismic reversal,” and we here thus present a complete, representative cycle of “seismic reversal” embedding in the changes of seismicity patterns before and after the Chi-Chi earthquake.     

Determination of spectral properties of earthquakes from their magnitudes

Following the recommendation of the “IASPEI Commission on Magnitudes” from 1967, P-wave magnitudes of distant earthquakes are being computed with the aid of the calibration functions of Gutenberg and Richter (1956). With the development of the instrumental seismology, especially in view of the availability of broad-band recordings, the question arises whether the functions are still adequate.In this investigation new calibration functions for P- and S-waves are presented, which are not only dependent on epicentral distance and focal depth, but also on the period of the spectral component of the wave. The determination of calibration functions is based on the currently most plausible global models for the velocity and anelasticity distribution. Taking into account the effective bandwidth of the instruments employed, eventually leads to the so-called spectral magnitudes.From the spectral magnitudes m(T) , the energy density spectrum E(T) of the respective wave can readily be computed. The corresponding formula is E(T) = 10^{2m(T)− 1.4}, with E(T) in joules per hertz.     

Hydrochemical baseline condition of groundwater at the Mizunami underground research laboratory (MIU)

Hydrochemical conditions up to depths of 1000 m below ground level around the Mizunami Underground Research Laboratory were investigated to construct a “baseline condition model” describing the undisturbed hydrochemical environment prior to excavation of the underground facilities at Mizunami, Gifu, Japan. Groundwater chemistry in this area was classified into a Na–Ca–HCO₃ type of groundwater in the upper part of sedimentary rock sequence and a Na–(Ca)–Cl type of groundwater in the deeper part of the sedimentary rock sequence and basement granite. The residence time of the groundwaters was estimated from their ¹⁴C contents to be approximately 9.3 ka in the middle part of the sedimentary rock and older than 50 ka in the deep part of the granite. The evolution processes of these groundwaters were inferred to be water–rock interactions such as weathering of plagioclase, dissolution of marine sulphate/sulphide minerals and carbonate minerals in the Na–Ca–HCO₃ type of groundwater, and mixing between “low-salinity water” in the shallow part and “higher-salinity water” in the deeper part of the granite in the Na–(Ca)–Cl type of groundwater. The source of salinity in the deeper part of the granite was possibly a palaeo-hydrothermal water or a fossil seawater that recharged in the Miocene, subsequently being modified by long-term water–rock interaction. The Cl-depth trend in granitic groundwater changes at a depth of −400 m below sea level. The hydrogeological properties controlling the groundwater flow and/or mixing processes such as advection and diffusion were inferred to be different at this depth in the granite. This hydrochemical conceptual model is indispensable not only when constructing the numerical model for evaluating the hydrochemical disturbance during construction and operation of the MIU facility, but also when confirming a hydrogeological model.     

Relative stress variations as determined by b-values from earthquakes in circum-pacific subduction zones

Assuming a relation of “b” to stress state, the possibility of globe-wide stress variation and transmission was investigated. The NOAA earthquake data file served to determine the temporal change in “b” of log N = a − bM from 1963 to 1975.Periods of six to eight years are observed in the b-values (stress pattern) for most circum-Pacific areas (South America, Tonga, Kermadec, New Hebrides, Kamchatka and Eastern Aleutians).In the Kurils, fore- and aftershock sequences of large earthquakes seem to mask any characteristic global pattern that might exist. These sequences exhibit low b-values (high stress) through the time of foreshocks and early stages of aftershocks, followed by rapid increase in b-values (decrease in stress).Use of a worldwide earthquake data file clearly yields less resolution of the temporal “b” variation than the use of local network studies published by other authors.Incidental to the study, 1124 earthquakes of the NOAA data file yield the M_s − m_b relations: M_s = 1.16m_b − 0.835 for 4.5 m_b 6 and: log₁₀M_s = 0.1432m_b − 0.0629 form_b > 6 with correlative coefficients of 0.994 and 0.992 respectively.     

Seismogeological effects on rocks during the 12 April 1998 upper So a Territory earthquake (NW Slovenia)

On 12 April 1998, the strongest earthquake (M_L 5.7), with an epicentre in Slovenia, in the last 100 years shook the upper So a Territory (NW Slovenia). Its maximum intensity was between VII and VIII according to the European Macroseismic Scale (EMS-98). Its epicentre is located in the area SE of the town of Bovec in the Krn mountain range (Julian Alps). Apart from substantial material damage to nearby settlements, the earthquake caused considerable changes to the landscape through many rockfalls and landslides. In this paper, we try to assess the effects of this event on the natural surroundings. The presented results are useful to understand the relation between “seismogeological” effects and the intensity scale EMS-98 for this particular event, but might also be a first step towards a future statistically meaningful assessment of the intensity scale on the basis of “seismogeological” effects.     

On the swelling potential of compacted high plasticity clays

We conducted a study of the spatial distributions of seismicity and earthquake hazard parameters for Turkey and the adjacent areas, applying the maximum likelihood method. The procedure allows for the use of either historical or instrumental data, or even a combination of the two. By using this method, we can estimate the earthquake hazard parameters, which include the maximum regional magnitude Mˆ_max, the activity rate of seismic events and the well-known bˆ value, which is the slope of the frequency-magnitude Gutenberg-Richter relationship. These three parameters are determined simultaneously using an iterative scheme. The uncertainty in the determination of the magnitudes was also taken into consideration. The return periods (RP) of earthquakes with a magnitude M ≥ m are also evaluated. The whole examined area is divided into 24 seismic regions based on their seismotectonic regime. The homogeneity of the magnitudes is an essential factor in such studies. In order to achieve homogeneity of the magnitudes, formulas that convert any magnitude to an M_S-surface scale are developed. New completeness cutoffs and their corresponding time intervals are also assessed for each of the 24 seismic regions. Each of the obtained parameters is distributed into its respective seismic region, allowing for an analysis of the localized seismicity parameters and a representation of their regional variation on a map. The earthquake hazard level is also calculated as a function of the form Θ = (Mˆ_max,RP_6.0), and a relative hazard scale (defined as the index K) is defined for each seismic region. The investigated regions are then classified into five groups using these parameters. This classification is useful for theoretical and practical reasons and provides a picture of quantitative seismicity. An attempt is then made to relate these values to the local tectonics.     

A comparison of evidence for late Holocene summer temperature variations in the northern hemisphere

Data on glacier, tree-line, tree-ring, pollen, and ice-core variations in North America, Greenland, and Europe during the last 2000 yr (up to A.D. 1800) are compared in detail on the century time scale. Only data that may be indicative of summer temperature changes are included, since these comprise most of the available paleoclimatic information, although some variations (especially of glaciers) may have been in response to precipitation changes instead. Radiocarbon dates and ¹⁴C-dated records are converted to calendar (dendrochronological) years using the calibration of M. Stuiver (1982, Radiocarbon 24, 1–26). Despite the basic uncertainties in dating, interpretation, response times, and “noise level” of proxy climatic data, there is at times good agreement among different kinds of evidence from within a region to suggest an episode of generally warmer or cooler summers. Three previously identified episodes find expression in records from all of the regions considered: the “Little Ice Age” of the last few centuries, a “Medieval Warm Period” around the 12th century A.D., and an earlier cold period some time between the 8th and 10th centuries. However, the timing of minima and maxima within these episodes seems to have varied from region to region (although the evidence is consistent within regions). In the 15th century, summers were warm in the eastern Canadian Arctic and southern Greenland while there was a cold episode in Europe and western North America.     

Methane-induced dolomite “chimneys” on the Kuroshima Knoll, Ryukyu islands, Japan

The Kuroshima Knoll is about 26 km south of Ishigaki Island in the southern part of the Ryukyu Arc. The area is considered to be the source area of “The 1771 Yaeyama Earthquake Tsunami”, which was due to the submarine landslide caused by an earthquake. It has been cleared from some investigations using “Dolphin 3K” and “Shinkai 2000” that there are large-scale dead Calyptogena colonies, many gravels of fallen dolomite chimneys and carbonates on the top of the Knoll [Matsumoto, T., Uechi, C., Kimura, M., 1997; Machiyama, H., Matsumoto, T., Matsumoto, R., Hattori, M., Okano, M., Iwase, R., Tomaru, H., 2001b.]. Carbonates of Kuroshima Knoll have various shapes and macroscopic textures. These have been classified into 4 types; shell crust (pavement), chimney, burrow, and nodule. It is clear that all chimney and burrow carbonates are composed of dolomite, while shell curst and nodule are composed of calcite, sometimes both calcite and dolomite. These carbonates are considered to have been formed by cold seep, because they are characterized by the light carbon isotopic ratio (semi-biogenic) and the heavy oxygen isotopic ratio. This suggests that methane hydrate layers develop under this survey area and the water that has the heavy oxygen and the light carbon isotopic ratio is derived from the dissociation of methane hydrate.     

Small-scale spatial variation of the stress field in the back-arc Aegean area: Results from the seismotectonic study of the broader area of Mygdonia basin (N. Greece)

In the present work a detailed seismotectonic study of the broader area of the Mygdonia basin (N. Greece) is performed. Digital data for earthquakes which occurred in the broader Mygdonia basin and were recorded by the permanent telemetric network of the Geophysical Laboratory of the Aristotle University of Thessaloniki during the period 1989–1999 were collected and fault plane solutions for 50 earthquakes which occurred in the study area were calculated with a modified first motions approach which incorporates amplitude and radiation pattern information. Fault plane solutions for the 3 main shocks of Volvi (23/05/78, M_W = 5.8 and 20/06/78, M_W = 6.5) and Arnaia (04/05/95, M_W = 5.8) events and the 1978 aftershock sequence were additionally used. Moreover, data from two local networks established in the Mygdonia basin were also incorporated in the final dataset.Determination of the stress field was realized by the use of the method of Gephart and Forsyth [Gephart, J.W., Forsyth, D.W., 1984. An improved method for determining the regional stress tensor using earthquake focal mechanism data: application to the San Fernando earthquake sequence: Jour. Geophys. Res., v.89, no. B11, p. 9305–9320] for the stress tensor inversion and the results were compared with independent estimates based on the calculation of the average moment tensor [Papazachos, C.B.,Kiratzi, A.A., 1992. A formulation for reliable estimation of active crustal deformation and its application to central Greece. Geophys. J. Int. 111, 424–432]. The obtained stress results show a relatively good agreement between the two approaches, with differences in the azimuth of the dominant extension axis of the order of 10°. Furthermore, comparison with independent information for the mean stress axes provided by the study of kinematics on neotectonic faults [Mountrakis, D., Kilias, A., Tranos, M., Thomaidou, E., Papazachos, C., Karakaisis, G., Scordilis, E., Chatzidimitriou, P., Papadimitriou, E., Vargemezis, G., Aidona, E., Karagianni, E., Vamvakaris, D. Skarlatoudis, A. 2003. Determination of the settings and the seismotectonic behavior of the main seismic-active faults of Northern Greece area using neotectonic and seismological data. Earthquake Planning and Protection Organisation (OASP) (in Greek)] shows a similar agreement with typical misfit of the order 10°. The stress inversion method was modified in order to select one or both nodal planes of the focal mechanism which corresponds to the “true” fault plane of the occurred earthquakes and was able to select a single fault plane in the majority of examined cases. Using this approach, the obtained fault plane rose diagrams are in agreement with results from various neotectonic studies. Moreover, several secondary active fault branches were identified, which are still not clearly observed in the field.     

Seismicity patterns before the M=5.8 2002, Palermo (Italy) earthquake: seismic quiescence and accelerating seismicity

Seismic quiescence and accelerating seismic energy release are considered as possible spatio-temporal patterns of the preparation process of the 6 September 2002 Palermo, Italy, earthquake (M 5.8). The detailed properties of the quiescence are analyzed applying the RTL algorithm. The RTL algorithm is based on the analysis of the RTL prognostic parameter, which is designed in such a way that it has a negative value if, in comparison with long-term background, there is a deficiency of events in the time–space vicinity of the tested point. The RTL parameter increases if activation of seismicity takes place. The RTL algorithm identified that a seismic quiescence started from the beginning of November 2001 and reached its minimum at the end of May 2002. The Palermo 2002 earthquake occurred 2 months after the RTL parameter restored its long-term background level. The application of a log-periodic time-to-failure model gives a “predicted” (in retrospect) magnitude M=6.2 main shock on 5 May 2002.     

Stress transfer in earthquakes, hazard estimation and ensemble forecasting: Inferences from numerical simulations

Observations indicate that earthquake faults occur in topologically complex, multi-scale networks driven by plate tectonic forces. We present realistic numerical simulations, involving data-mining, pattern recognition, theoretical analyses and ensemble forecasting techniques, to understand how the observable space–time earthquake patterns are related to the fundamentally inaccessible and unobservable dynamics. Numerical simulations can also help us to understand how the different scales involved in earthquake physics interact and influence the resulting dynamics. Our simulations indicate that elastic interactions (stress transfer) combined with the nonlinearity in the frictional failure threshold law lead to the self-organization of the statistical dynamics, producing 1) statistical distributions for magnitudes and frequencies of earthquakes that have characteristics similar to those possessed by the Gutenberg–Richter magnitude–frequency distributions observed in nature; and 2) clear examples of stress transfer among fault activity described by stress shadows, in which an earthquake on one group of faults reduces the Coulomb failure stress on other faults, thereby delaying activity on those faults. In this paper, we describe the current state of modeling and simulation efforts for Virtual California, a model for all the major active strike slip faults in California. Noting that the Working Group on California Earthquake Probabilities (WGCEP) uses statistical distributions to produce earthquake forecast probabilities, we demonstrate that Virtual California provides a powerful tool for testing the applicability and reliability of the WGCEP statistical methods. Furthermore, we show how the simulations can be used to develop statistical earthquake forecasting techniques that are complementary to the methods used by the WGCEP, but improve upon those methods in a number of important ways. In doing so, we distinguish between the “official” forecasts of the WGCEP, and the “research-quality” forecasts that we discuss here. Finally, we provide a brief discussion of future problems and issues related to the development of ensemble earthquake hazard estimation and forecasting techniques.     

A study of the kinetic parameters of individual macerals from Upper Permian coals in South China via open-system pyrolysis

A unique Upper Permian coal, Leping coal, is widely distributed in South China. The coal samples studied in the paper were collected from two mines in the Shuicheng coalfield of Guizhou Province, southwest China. The geochemical works including coal petrography, maceral content, Rock–Eval pyrolysis, and kinetic modelling of hydrocarbon-generating have been carried out on whole coal and individual macerals. The higher contents of volatile matter, elemental hydrogen, and tar yield, and the high hydrocarbon generation potential of the Leping coals are attributed to their high content of “barkinite”, a special liptinite maceral.The hydrocarbon generation potential of “barkinite” (S₂=287 mg/g, hydrogen index (HI)=491 mg/g TOC) is greater than that of vitrinite (S₂=180 mg/g, HI=249 mg/g TOC), and much higher than that of fusinite (S₂=24 mg/g, HI=35 mg/g TOC). At the same experimental conditions, “barkinite” has a higher threshold and a narrower “oil window” than those of vitrinite and fusinite, and consequently, can generate more hydrocarbons in higher coalification temperature and shorter geological duration. Data from the activation energy distributions indicate that “barkinite” has a more homogenous chemical structure than that of vitrinite and fusinite. The above-mentioned characteristics are extremely important for exploring hydrocarbon derived from the Leping coals in South China.     

Evidence for cyclic reactions between andalusite, “sericite” and sillimanite, mount franks area, willyama complex N.S.W

In the Mt. Franks area of the Willyama Complex, microfabric evidence suggests that the alteration of andalusite to sillimanite has taken place by a process similar to that suggested by Carmichael (1969). Andalusite is pre- to syn-S₂ in age. Alteration to “sericite” has resulted in the formation of “sericite” laths, some of which are crenulated about S₂, and some which are syn- and post-S₂. “Fibrolite” occurs in these andalusite—“sericite” aggregates within the sillimanite zone and is wholly embedded in “sericite”. “Fibrolite” is pre- to syn-S₂ in age. This evidence is interpreted as suggesting that the formation of sillimanite from andalusite took place via a “sericite” phase.Further microfabric observations are interpreted to imply constant volume for the reaction aluminosilicate → “sericite”. This suggests a situation in which Al³⁺ is relatively mobile but Al⁴⁺ is relatively immobile. This suggestion differs from Carmichael's (1969) idea of Al³⁺ immobility.     

Senescence- and death-related alteration of chlorophylls and carotenoids in marine phytoplankton

This report extends previous work ([Louda et al., 1998a] and [Louda et al., 1998b]. Chlorophyll degradation during senescence and death. Organic Geochemistry 29, 1233–1251.) in which we detailed type-I (alteration) and -II (destruction) degradation of chlorophyll with representative fresh water phytoplankton. The present study covers similar experiments with marine phytoplankton, namely, a cyanobacterium (“ANA” Anacystis sp), a coccolithophore (“COC” Coccolithophora sp.), a dinoflagellate (“GYM” Gymnodinium sp.) and two diatoms (“CYC” Cyclotella meneghiniana and “THAL” Thalassiosira sp.). Mg loss (‘pheophytinization') was rapid and continuous in all species under room-oxic conditions and slow or sporadic under anoxic conditions. The proportion of dephytylated pigments (pheophorbides-a, chlorophyllides-a), relative to the phytylated forms (chlorophyll-a, pheophytins-a), increased over the first year under room-oxic conditions and in room-anoxic conditions only in “CYC”. Pheophorbide-a was converted to pyropheophorbide-a within 15 months only in “THAL” and “ANA”, and slightly in “COC”. After 9–15 months of oxic incubation, “COC” was found to contain traces of purpurin-18 phytyl ester. Consideration of carotenoid pigments is also included herein. All fucoxanthin containing species, except “THAL”, exhibited conversion of fucoxanthin to fucoxanthinol in room-oxic conditions. Diadinoxanthin was rapidly de-epoxidized to give diatoxanthin within the first 2–4 weeks. Diatoxanthin then disappeared from all species by 15 months with a concurrent increase in a pigment which we tentatively identify as a cis-zeaxanthin. Incubations of pure cultures are found to be an effective way by which to model the early type-I reactions for both chlorophylls and carotenoids. The influence of oxygen during senescence-death and the onset of early diagenesis is of paramount importance. The absence of oxygen and, by inference, aerobic microbiota, retards the breakdown of these pigments dramatically.     

The Cenomanian/Turonian Boundary Event (CTBE) at Wunstorf, north-west Germany, as reflected by marine palynology

The Cenomanian/Turonian Boundary Event (CTBE) at Wunstorf, north-west Germany, has been analysed palynologically by high resolution sampling to reconstruct changes in relative sea-level and water mass character within photic zone waters. Based on changes in the ratio of terrigenous sporomorphs to marine palynomorphs (t/m index), the distribution of the organic-walled algal taxa as well as of selected dinocyst taxa and groups the section can largely be subdivided into pre-“plenus-bed” and post-“plenus-bed” intervals, reflecting different stages of third-order relative sea-level cycles and/or changes in water mass influence in the photic zone. Accordingly, the pre-“plenus-bed” interval is placed in a transgressive systems tract starting at the “facies change” event (C. guerangeri/M. geslinianum ammonite Zone boundary) with the maximum flooding surface at the top of the “Chondrites II” bed (top of R. cushmani Biozone). A highstand systems tract is suggested from the base of the “plenus-bed” up the base of the “fish-shale” event. Within the “fish-shale” event interval, a transgressive systems tract is suggested to start at the base of the thin, grey-green marly interbed. The Cenomanian/Turonian boundary proper, as defined by the first occurrence of Mytiloides spp., as well as the lowermost Turonian are located within the initial phase of a transgressive systems tract. With respect to water mass characteristics within photic-zone waters, the pre-“plenus-bed” interval is predominantly characterized by warm water masses that changed gradually towards the deposition of the “Chondrites II” bed, where a strong influence of cool and/or salinity-reduced waters is indicated by various palynological proxies. Within the post-“plenus-bed” interval a mixture and/or alternation of warmer and cooler waters is indicated, with the warmer water influence increasing gradually towards and within the Lower Turonian stage. The increased proportions of prasinophytes within the “Chondrites II” bed and parts of the “fish-shale” interval may indicate availability of reduced nitrogen chemospecies, especially ammonium, within photic-zone waters as a function of a vertical expansion of the oceanic O₂-minimum zone.     

Magnitudes of great shallow earthquakes from 1953 to 1977

The surface-wavemagnitudes M_s are determined for 30 great shallow earthquakes that occurred during the period from 1953 to 1977. The determination is based on the amplitude and period data from all available station bulletins, and the same procedure as that employed in Gutenberg and Richter's “Seismicity of the Earth” is used. During this period, the Chilean earthquake of 1960 has the largest M_s, 8.5. The surface-wave magnitudes listed in “Earthquake Data Reports” are found to be higher than M_s on the average. By using the same method as that used by Gutenberg, the broad-band body-wave magnitudes m_B are determined for great shallow shocks for the period from 1953 to 1974. m_B is based on the amplitudes of P, PP and S waves which are measured on broadband instruments at periods of about 4–20 s. The 1-s body-wave magnitudes listed in “Bulletin of International Seismological Center” and “Earthquake Data Reports” are found to be much smaller than m_B on the average. Through the examination of Gutenberg and Richter's original worksheets, the relation between m_B and M_sis revised to m_B = 0.65 M_s+ 2.5 which well satisfies the mg and M_sdata for M_sbetween 5.2 and     

A Late Jurassic fossil assemblage in Gondwana: Biostratigraphy and correlations of the Tacuarembó Formation, Parana Basin, Uruguay

The Tacuarembó Formation has yielded a fossil assemblage that includes the best known body fossils, consisting of isolated scales, teeth, spines, and molds of bones, recovered from thin and patchy bonebeds, from the Botucatu Desert, Parana Basin, South America. The remains are preserved in the sandstones widespread around the city of Tacuarembó. We propose a new formalized nomenclature for the Tacuarembó Formation, naming its “Lower” and “Upper” members as the Batoví (new name) and Rivera (new rank) members, respectively. An assemblage zone is defined for the Batoví Member (fluviolacustrine and aeolian deposits). In this unit, the freshwater hybodontid shark Priohybodus arambourgi D’Erasmo is well represented. This species was previously recorded in Late Jurassic–Early Cretaceous units of the Sahara and the southern Arabian Peninsula. Globally considered, the fossil assemblage of this member (P. arambourgi, dipnoan fishes, Ceratosaurus-like theropods, and conchostracans) is indicative of a Kimmeridgian–Tithonian age, which in combination with the stratigraphic relationships of the Tacuarembó Formation with the overlying basalts of the Arapey Formation (132 My average absolute age) implies that the latter was deposited during the Kimmeridgian–Hauterivian interval.     

Gravitational potential as a source of earthquake energy

Some degree of tectonic stress within the earth originates from gravity acting upon density structures. The work performed by this “gravitational tectonics stress” must have formerly existed as gravitational potential energy contained in the stress-causing density structure.According to the elastic rebound theory (Reid, 1910), the energy of earthquakes comes from an elastic strain field built up by fairly continuous elastic deformation in the period between events. For earthquakes resulting from gravitational tectonic stress, the elastic rebound theory requires the transfer of energy from the gravitational potential of the density structures into an elastic strain field prior to the event.An alternate theory involves partial gravitational collapse of the stress-causing density structures. The earthquake energy comes directly from a net decrease in gravitational potential energy. The gravitational potential energy released at the time of the earthquake is split between the energy released by the earthquake, including work done in the fault zone and an increase in stored elastic strain energy. The stress associated with this elastic strain field should oppose further fault slip.