Identification of changes in heavy rainfall events in Ontario,Canada

To assess whether changes in the frequency of heavy rainfall events are occurring over time, annual maximum records from 21 rainfall gauges in Ontario are examined using frequency analysis methods. Relative RMSE and related boxplots are used to characterize assessment for selecting distributions; the Gumbel distribution is verified as one of the most suitable distributions to provide accurate quantile estimates. Records were divided into two time periods, and tested using the Mann-Kendall test and lag-1 autocorrelations to ensure that data in each period are identically distributed. The confidence intervals of design rainfalls for each return period (2, 5, 10, and 25-year) are derived by using resampling method, and compared at 90 % confidence levels. The changes in heavy rainfall intensities are tested at gauges across the Province of Ontario. Several significant decreases in heavy rainfall intensities are identified in central and southern Ontario. Increases in heavy rainfall intensities are identified in gauges at Sioux Lookout and Belleville. The sensitivity analysis of changes identified with respect to the year of splitting indicates changes are occurring during the 1980s and 1990s.     

Deflection amplification factor for estimating seismic lateral deformations of RC frames

Different values have been assigned to the ratio of the defl ection amplifi cation factor(Cd) to the response modifi cation factor(R) for a specifi ed force-resisting system in the seismic design provisions while the same application is defi ned for it. An analytical study of the seismic responses of several reinforced concrete frames subjected to a suite of earthquake records performed in this research indicate that the stories’ overstrength and stiffness distribution along the structural height can affect local defl ections more than global ones. Therefore, the Cd/R ratio is calculated based on the ratio of both maximum inelastic to maximum elastic displacements and interstory drifts. Due to damage concentration in some specifi c stories, the defl ection amplifi cation factor calculated based on inelastic interstory drifts was larger than that of the inelastic displacements. Consequently, a minimum value of 1.0 is recommended for the Cd/R ratio in order to estimate maximum inelastic drifts. The ratio of inelastic to elastic displacement was generally found to increase slightly along the structural height for the studied RC models. In addition, it was detected that the story damage indices of the studied RC frames decrease when the inverted value of inelastic interstory drift ratios are increased through a(negative) power form.     

Full-scale tests of wind effects on a long span roof structure

Full-scale measurements are regarded as the most reliable method to evaluate wind effects on large buildings and structures. Some selected results are presented in this paper from the full-scale measurement of wind effects on a long-span steel roof structure during the passage of Typhoon Fanapi. Some fi eld data, including wind speed and direction, acceleration responses, etc., were continuously and simultaneously recorded during the passage of the typhoon. Comprehensive analysis of the measured data is conducted to evaluate the typhoon-generated wind characteristics and its effects on a long-span steel roof. The fi rst four natural frequencies and their vibration mode shapes of the Guangzhou International Sports Arena(GISA) roof are evaluated by the stochastic subspace identifi cation(SSI) method and comparisons with those from fi nite element(FE) analysis are made. Meanwhile, damping ratios of the roof are also identifi ed by the SSI method and compared with those identifi ed by the random decrement method; the amplitude-dependent damping behaviors are also discussed. The fullscale measurement results are further compared with the corresponding wind tunnel test results to evaluate its reliability. The results obtained from this study are valuable for academic and professional engineers involved in the design of large-span roof structures.     

Seismic performance assessment and base-isolated floor protection of statues exhibited in museum halls

A study concerning the evaluation of seismic response of statues exhibited in art museums, and a base-isolated floor strategy for their enhanced protection, are presented in this paper. Attention is particularly focused on statues made of small tensile strength materials, whose behaviour is simulated by a smeared-crack finite element approach. Seismic performance is assessed by referring to four levels specially postulated herein, and namely: (1) Rest conditions; (2) No rocking; (3) Damage control; and (4) Collapse prevention. The response is investigated via incremental dynamic analysis, by progressively increasing the amplitude of the ground motion histories adopted as inputs, and by relating output data to the limit conditions fixed for the above-mentioned performance levels. The assessment procedure is applied to a demonstrative case study, represented by a marble statue to be exhibited in the museum wing situated at the ground level of a medieval castle in Italy, according to an architectural hypothesis of partial rebuilding and reuse of the stronghold. The design solution for the base-isolated floor consists in a system of double-friction pendulum isolators. The finite element model constitutive laws and parameters, the dynamic analyses carried out in fixed-base and base-isolated floor conditions, and the practical implementation of the assumed performance assessment criteria are reported for the statue examined, along with a selection of technical details of the floor design.     

A cyclic UH model for sand

A simple cyclic elastoplastic constitutive model for sand is proposed based on the UH model for overconsolidated clay. The proposed model has the following features. First, in order to describe the stress-induced anisotropy in sand, a rotational hardening rule is introduced for the evolution of the yield surface axis with development of plastic deviatoric strain in the principal stress space. Second, the relationship between the rotational axis and stress-induced anisotropy is modeled by introducing the slope of rotational axis into the yield function. The fl atness of the yield surface can be determined by the slope of rotational axis. Finally, a revised unifi ed hardening parameter is proposed to incorporate the stress-induced anisotropy. The model capability in describing the cyclic response of sand is verifi ed by comparing the simulations with available test results.     

The structural characteristics,age of origin,and tectonic attribute of the Erguna Fault,NE China

The Erguna Fault runs along the east bank of the Erguna River in NE China and is a large-scale ductile shear zone comprising granitic mylonites. This paper reports on the geometry, kinematic indicators, and 40Ar/39 Ar biotite ages of the granitic mylonites, to constrain the structural characteristics, forming age, and tectonic attribute of the Erguna ductile shear zone. The zone strikes NE and records a top-to-the-NW sense of shear. A mylonitic foliation and stretching lineation are well developed in the mylonites, which are classified as S-L tectonites. Logarithmic flinn parameters(1.18–2.35) indicate elongate strain which approximates to plane strain. Kinematic vorticity numbers are 0.42–0.92 and 0.48–0.94, based on the polar Mohr diagram and the oblique foliation in quartz ribbons, respectively, suggesting that the ductile shear zone formed under general shear, or a combination of simple and pure shear. According to finite strain and kinematic vorticity analyses, the Erguna Fault is a lengthening-thinning ductile shear zone that formed by extension. The deformation behavior of minerals in the mylonites indicates that the fault was the site of three stages of deformation: an initial stage of middle- to deep-level, high-temperature shear, a post-stress recovery phase of high-temperature static recrystallization, and a final phase of low-temperature uplift and cooling. The 40Ar/39 Ar plateau ages of biotite from the granitic mylonites are 106.16 ± 0.79 and 111.55 ± 0.67 Ma, which constrain the timing of low-temperature uplift and cooling but are younger than the ages of metamorphic core complexes(MCCs) in the Transbaikalia-northeast Mongolia region. Using measured geological sections, microtectonics, estimates of finite strain and kinematic vorticity, and regional correlations and geochronology, we conclude that the Erguna Fault is an Early Cretaceous, NNE-trending, large-scale, sub-horizontal, and extensional ductile shear zone. It shares a similar tectonic background with the MCCs, volcanic fault basins, and large and super-large volcanic-hydrothermal deposits in Transbaikalia-northeast Mongolia and the western Great Khingan Mountains, all of which are the result of overthickened crust that gravitationally collapsed and extended in the Early Cretaceous after plate collision along the present-day Sino-Russia-Mongolia border tract.     

A Nd- and O-isotope study of the REE-rich peralkaline Strange Lake granite: implications for Mesoproterozoic A-type magmatism in the Core Zone (NE-Canada)

It is well established that A-type granites enriched in high field strength elements, such as Zr, Nb and the REE, form in anorogenic tectonic settings. The sources of these elements and the processes controlling their unusual enrichment, however, are still debated. They are addressed here using neodymium and oxygen isotope analyses of samples from the 1.24 Ga Strange Lake pluton in the Paleoproterozoic Core Zone of Québec-Labrador, an A-type granitic body characterized by hyper-enrichment in the REE, Zr, and Nb. Age-corrected εNd values for bulk rock samples and sodic amphiboles (mainly arfvedsonite) from the pluton range from ?0.6 to ?5.7, and ?0.3 to ?5.3, respectively. The εNd values for the Napeu Kainiut quartz monzonite, which hosts the pluton, range from ?4.8 to ?8.1. The ¹⁴⁷Sm/¹⁴⁴Nd ratios of the suite and the host quartz monzonite range from 0.0967 to 0.1659, large variations that can be explained by in situ fractionation of early LREE-minerals (Strange Lake), and late hydrothermal HREE remobilization. Oxygen isotope analyses of quartz of both Strange Lake and the host yielded δ¹⁸O values between +8.2 and +9.1, which are considerably higher than the mantle value of 5.7 ± 0.2‰. Bulk rock oxygen isotope analyses of biotite-gneisses in the vicinity of the Strange Lake pluton yielded δ¹⁸O values of 6.3, 8.6 and 9.6‰. The negative εNd values and positive δ¹⁸O values of the Strange Lake and Napeu Kainiut samples indicate that both magmas experienced considerable crustal contamination. The extent of this contamination was estimated, assuming that the contaminants were sedimentary-derived rocks from the underlying Archean Mistinibi (para-) gneiss complex, which is characterized by low εNd and high δ¹⁸O values. Mixing of 5–15% of a gneiss, having an εNd value of ?15 and a δ¹⁸O value of +11, with a moderately enriched mantle source (εNd = +0.9, δ¹⁸O = +6.3) would produce values similar to those obtained for the Strange Lake granites. Based on analogies between the Nain Plutonic Suite and the Gardar alkaline igneous province (SW-Greenland), we conclude that the Strange Lake pluton and associated REE-mineralized anorogenic bodies formed from a combination of subduction-induced fertilization of the sublithospheric mantle, crustal extension and in situ magma evolution.     

The geochemical and petrological characteristics of prenatal caldera volcano: a case of the newly formed small dacitic caldera,Hijiori, Northeast Japan

Evaluating the magma depth and its physical properties is critical to conduct a better geophysical assessment of magma chambers of caldera volcanoes that may potentially cause future volcanic hazards. To understand pre-eruptive conditions of a magma chamber before its first appearance at the surface, this paper describes the case of Hijiori caldera volcano in northeastern Japan, which emerged approximately 12,000 years ago at a place where no volcano ever existed. We estimated the depth, density, bulk modulus, vesicularity, crystal content, and bulk H\(_2\)O content of the magma chamber using petrographic interpretations, bulk and microchemical compositions, and thermodynamic calculations. The chemical mass balance calculations and thermodynamic modeling of the erupted magmas indicate that the upper portion of the Hijiori magmatic plumbing system was located at depths between 2 and 4 km, and had the following characteristics: (1) pre-eruptive temperature: about 780 \(^{\circ }\)C; (2) bulk magma composition: 66 ± 1.5 wt% SiO\(_{2}\); (3) bulk magmatic H\(_2\)O: approximately 2.5 wt%, and variable characteristics that depend on depth; (4) crystal content: \(\le\)57 vol%; (5) bulk modulus of magma: 0.1–0.8 GPa; (6) magma density: 1.8–2.3 g/cm³; and (7) amount of excess magmatic H\(_2\)O: 11–32 vol% or 48–81 mol%. The range of melt water contents found in quartz-hosted melt inclusions (2–9 wt%) suggests the range of depth phenocrysts growth to be wide (2\(\sim\)13 km). Our data suggest the presence of a vertically elongated magma chamber whose top is nearly solidified but highly vesiculated; this chamber has probably grown and re-mobilized by repeated injections of a small amount of hot dacitic magma originated from the depth.     

Reaction-induced grain boundary cracking and anisotropic fluid flow during prograde devolatilization reactions within subduction zones

Devolatilization reactions during prograde metamorphism are a key control on the fluid distribution within subduction zones. Garnets in Mn-rich quartz schist within the Sanbagawa metamorphic belt of Japan are characterized by skeletal structures containing abundant quartz inclusions. Each quartz inclusion was angular-shaped, and showed random crystallographic orientations, suggesting that these quartz inclusions were trapped via grain boundary cracking during garnet growth. Such skeletal garnet within the quartz schist formed related to decarbonation reactions with a positive total volume change (?V _t > 0), whereas the euhedral garnet within the pelitic schists formed as a result of dehydration reaction with negative ?V _t values. Coupled hydrological–chemical–mechanical processes during metamorphic devolatilization reactions were investigated by a distinct element method (DEM) numerical simulation on a foliated rock that contained reactive minerals and non-reactive matrix minerals. Negative ?V _t reactions cause a decrease in fluid pressure and do not produce fractures within the matrix. In contrast, a fluid pressure increase by positive ?V _t reactions results in hydrofracturing of the matrix. This fracturing preferentially occurs along grain boundaries and causes episodic fluid pulses associated with the development of the fracture network. The precipitation of garnet within grain boundary fractures could explain the formation of the skeletal garnet. Our DEM model also suggests a strong influence of reaction-induced fracturing on anisotropic fluid flow, meaning that dominant fluid flow directions could easily change in response to changes in stress configuration and the magnitude of differential stress during prograde metamorphism within a subduction zone.     

How do granitoid magmas mix with each other? Insights from textures,trace element and Sr–Nd isotopic composition of apatite and titanite from the Matok pluton (South Africa)

In plutonic systems, magma mixing is often modelled by mass balance based on whole-rock geochemistry. However, magma mixing is a chaotic process and chemical equilibration is controlled by non-linear diffusive–advective processes unresolved by the study of bulk samples. Here we present textural observations, LA-(MC-)ICP-MS trace element and Sr–Nd isotopic data of accessory apatites and titanites from a hybrid granodiorite of the Neoarchean Matok pluton (South Africa), collected in a zone of conspicuous mixing between mafic and felsic magmas. Apatite grains mostly show a pronounced zoning in CL images, corresponding to abrupt changes in REE and HFSE concentrations recording their transfer through compositionally different melt domains during mixing. These grains crystallized early, at temperatures of 950–1000 °C. Titanite grains crystallized at temperatures of 820–900 °C (Zr-in-sphene thermometry). They show limited intra-grain chemical variations but huge inter-grain compositional scatter in REE and HFSE, pinpointing crystallization within a crystal mush, from isolated melt pockets having different composition from one another owing to incomplete chemical homogenization and variable Rayleigh fractionation. These chemical–textural characteristics, in combination with partitioning models and Polytopic Vector Analysis, point to “self-mixing” between co-genetic dioritic and granodioritic/granitic magmas. Both resulted from differentiation of mantle-derived mafic melts, showing that mixing does not necessarily involve magmas from contrasted (crust vs. mantle) sources. Systematic variations in εNd _t (?4.5 to ?2.5) and ⁸⁷Sr/⁸⁶Sr_(i) (0.703–0.707) of titanite and apatite grains/domains crystallized from the two magmas point to an isotopically inhomogeneous mantle source, which is not resolved by bulk-rock isotopic data. Interaction between the two magmas must have occurred at relatively high temperatures (ca. 900°C) so that their viscosity contrast remained low, allowing efficient mechanical mixing. Despite this, chemical homogenization was incomplete, as recorded by diffusive fractionation between REE–HFSE and Sr. Modelling thereof reveals that chemical exchange between the liquid phases of the two mixed magmas did not last more than a few tens to hundreds of years. The chemical equilibration between mixed magmas thus strongly depends on the considered elements and observational length scales.