The magnitude Mw 6.2 earthquake of February 22nd 2011 that struck beneath the city of Christchurch, New Zealand, caused widespread damage and was particularly destructive to the Central Business District (CBD). The shaking caused major damage, including collapses of structures, and initiated ground failure in the form of soil liquefaction and consequent effects such as sand boils, surface flooding, large differential settlements of buildings and lateral spreading of ground towards rivers were observed. A research project underway at the University of Canterbury to characterise the engineering behaviour of the soils in the region was influenced by this event to focus on the performance of the highly variable ground conditions in the CBD. This paper outlines the methodology of this research to characterise the key soil horizons that underlie the CBD that influenced the performance of important structures during the recent earthquakes, and will influence the performance of the rebuilt city centre under future events. The methodology follows post-earthquake reconnaissance in the central city, a desk study on ground conditions, site selection, mobilisation of a post-earthquake ground investigation incorporating the cone penetration test (CPT), borehole drilling, shear wave velocity profiling and Gel-push sampling followed by a programme of laboratory testing including monotonic and cyclic testing of the soils obtained in the investigation. The research is timely and aims to inform the impending rebuild, with appropriate information on the soils response to dynamic loading, and the influence this has on the performance of structures with various foundation forms.
This report summarizes the development of a region-wide surficial soil shear wave velocity (Vs ) model based on the unique combination of a large high-spatial-density database of cone penetration test (CPT) logs in the greater Christchurch urban area (> 15, 000 logs as of 1 February 2014) and the Christchurch-specific empirical correlation between soil Vs and CPT data developed by McGann et al. [1, 2]. This model has applications for site characterization efforts via maps of time-averaged Vs over specific depths (e.g. Vs30, Vs10), and for numerical modeling efforts via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin [3] currently being developed for the purpose of broadband ground motion simulation. The general development of these region-wide near-surface Vs models includes the following general phases, with each discussed in separate chapters of this report. • An evaluation of the available CPT dataset for suitability, and the definition of other datasets and assumptions necessary to characterize the surficial sediments of the region to 30 m depth. • The development of time-averaged shear wave velocity (Vsz) surfaces for the Christchurch area from the adopted CPT dataset (and supplementary data/assumptions) using spatial interpolation. The Vsz surfaces are used to explore the characteristics of the near-surface soils in the regions and are shown to correspond well with known features of the local geology, the historical ecosystems of the area, and observations made following the 2010- 2011 Canterbury earthquakes. • A detailed analysis of the Vs profiles in eight subregions of Christchurch is performed to assess the variablity in the soil profiles for regions with similar Vsz values and to assess Vsz as a predictive metric for local site response. It is shown that the distrubution of soil shear wave velocity in the Christchurch regions is highly variable both spatially (horizontally) and with depth (vertically) due to the varied geological histories for different parts of the area, and the highly stratified nature of the nearsurface deposits. This variability is not considered to be greatly significant in terms of current simplified site classification systems; based on computed Vs30 values, all considered regions can be categorized as NEHRP sites class D (180 < Vs < 360 m/s) or E (Vs < 180 m/s), however, detailed analysis of the shear wave velocity profiles in different subregions of Christchurch show that the expected surficial site response can vary quite a bit across the region despite the relative similarity in Vs30
Structural members made of laminated veneer lumber (LVL) in combination with unbonded post-tensioning have recently been proposed, which makes it possible to design moment-resisting frames with longer spans for multi-storey timber buildings. It has been shown that prefabricated and prestressed timber structures can be designed to have excellent seismic resistance, with enhanced re-centring and energy dissipation characteristics. The post-tensioning provides re-centring capacity while energy is dissipated through yielding of mild steel dissipating devices. This paper summarizes an experimental investigation into the seismic response of LVL columns to bi-directional seismic loading, performed as part of a research programme on timber structures at the University of Canterbury. The experimental investigation includes testing under both quasi-static cyclic and pseudo-dynamic protocols. The results show excellent seismic performance, characterized by negligible damage of the structural members and small residual deformations, even under the combined effect of loading in two directions. Energy is dissipated mostly through yielding of external dissipators connecting the column and the foundation, which can be easily removed and replaced after an earthquake. Since post-tensioning can be economically performed on site, the system can be easily implemented in multi-storey timber buildings
A three dimensional approximately half scale experimental subassemblage is currently being tested at the University of Canterbury to investigate the effect of precast-prestressed floor units, which do not span past the internal columns, on the seismic performance of reinforced concrete moment resisting frames. This paper reports the preliminary results from the test, with the focus on elongation within the plastic hinges and strength enhancement in the frames. The preliminary results have shown that elongation between the external and internal plastic hinges varies by more than two fold. With the addition of the prestressed floor units, the strength of the moment resisting frame used in the test was found to be 25% higher than the current code specified value. In other situations, particularly where there are more than 2 bays in a moment resisting frame, greater strength enhancement may be expected. Any underestimation of beam strength is undesirable as it may result in the development of nonductile failure modes in a major earthquake.
During the past two decades, the focus has been on the need to provide communities with structures that undergo minimal damage after an earthquake event while still being cost competitive. This has led to the development of high performance seismic resisting systems, and advances in design methodologies, in order respect this demand efficiently. This paper presents the experimental response of four pre-cast, post-tensioned rocking wall systems tested on the shake-table at the University of Canterbury. The wall systems were designed as a retrofit solution for an existing frame building, but are equally applicable for use in new design. Design of the wall followed a performance-based retrofit strategy in which structural limit states appropriate to both the post-tensioned wall and the existing building were considered. Dissipation for each of the four post-tensioned walls was provided via externally mounted devices, located in parallel to post-tensioned tendons for re-centring. This allowed the dissipation devices to be easily replaced or inspected following a major earthquake. Each wall was installed with viscous fluid dampers, tension-compression yielding steel dampers, a combination of both or no devices at all – thus relying on contact damping alone. The effectiveness of both velocity and displacement dependant dissipation are investigated for protection against far-field and velocity-pulse ground motion characteristics. The experimental results validate the behaviour of ‘Advanced Flag-Shape’ rocking, dissipating solutions which have been recently proposed and numerically tested. Maximum displacements and material strains were well controlled and within acceptable bounds, and residual deformations were minimal due to the re-centring contribution from the post-tensioned tendons. Damage was confined to inelastic yielding (or fluid damping) of the external dampers.
Timber-based hybrid structures provide a prospective solution for utilizing environmentally friendly timber material in the construction of mid-rise or high-rise structures. This study mainly focuses on structural damage evaluation for a type of timber-steel hybrid structures, which incorporate prefabricated light wood frame shear walls into steel moment-resisting frames (SMRFs). The structural damage of such a hybrid structure was evaluated through shake table tests on a four-story large-scale timber-steel hybrid structure. Four ground motion records (i.e., Wenchuan earthquake, Canterbury earthquake, El-Centro earthquake, and Kobe earthquake) were chosen for the tests, with the consideration of three different probability levels (i.e., minor, moderate and major earthquakes) for each record. During the shake table tests, the hybrid structure performed quite well with visual damage only to wood shear walls. No visual damage in SMRF and the frame-to-wall connections was observed. The correlation of visual damage to seismic intensity, modal-based damage index and inter-story drift was discussed. The reported work provided a basis of knowledge for performance-based seismic design (PBSD) for such timber-based hybrid structures.
A member of the Los Angeles County Fire Department Search and Rescue Team looking at the maps at the entrance to 'Camp Hollywood', the headquarters of the American Civil Defence Forces in Latimer Square. After the 22 February 2011 earthquake, emergency service agencies set up their headquarters in Latimer Square.
The entrance to 'Camp Hollywood', the American Civil Defence Force headquarters in Latimer Square. After the 22 February 2011 earthquake, emergency service agencies set up their headquarters in Latimer Square. The entrance way has a wooden sign reading, "Camp Hollywood" and a series of maps of Christchurch and New Zealand.
A member of the Los Angeles County Fire Department Search and Rescue Team looking at the maps at the entrance to 'Camp Hollywood', the headquarters of the American Civil Defence Forces in Latimer Square. After the 22 February 2011 earthquake, emergency service agencies set up their headquarters in Latimer Square.
Photograph captioned by Fairfax, "Local relief: Waimakariri District Mayor Ron Keating (left) receives a cheque for $100,000 from MainPower managing director Allan Berge at the Civil Defence headquarters set up in council's chambers. The money is targeted for Waimakariri families who have been hard hit by the earthquake".
A photograph of members of the Wellington Emergency Management Office lining up outside a Maggi van which is distributing out soup. The van is parked outside Cowles Stadium on Pages Road, which was set up by Civil Defence as temporary accommodation for those displaced by the 4 September 2010 earthquake.
A member of the Los Angeles County Fire Department Search and Rescue Team with an Australian Police Officer at the entrance to 'Camp Hollywood', the headquarters of the American Civil Defence Forces in Latimer Square. After the 22 February 2011 earthquake, emergency service agencies set up their headquarters in Latimer Square.
University of Canterbury staff members are escorted by Civil Defence members in order to retrieve essential items from their offices. The photographer comments, "E-learning team after clearing their offices. Brendon Stillwell (ICTS tech helping with PCs), Antoine Monti, Susan Tull, Herbert Thomas, Paul Nicholls, Gregor Ronald, Lei Zhang, Jess Hollis".
University of Canterbury staff members are escorted by Civil Defence members in order to retrieve essential items from their offices. The photographer comments, "E-learning team after clearing their offices. Brendon Stillwell (ICTS tech helping with PCs), Antoine Monti, Susan Tull, Herbert Thomas, Paul Nicholls, Gregor Ronald, Lei Zhang, Jess Hollis".
A press release from the United States of America Embassy New Zealand, about their response to the 22 February 2011 earthquake in Christchurch.
A story submitted by Tracy to the QuakeStories website.
Summary of oral history interview with Alia Afzali about her experiences of the Canterbury earthquakes.
Shows two visitors to Christchurch who can enjoy the excitement of the danger of earthquakes without car bombings. Context: refers to the Christchurch earthquakes and continuing aftershocks. Quantity: 1 digital cartoon(s).
This paper investigates the effects of variability in source rupture parameters on site-specific physics-based simulated ground motions, ascertained through the systematic analysis of ground motion intensity measures. As a preliminary study, we consider simulations of the 22 February 2011 Christchurch earthquake using the Graves and Pitarka (2015) methodology. The effects of source variability are considered via a sensitivity study in which parameters (hypocentre location, earthquake magnitude, average rupture velocity, fault geometry and the Brune stress parameter) are individually varied by one standard deviation. The sensitivity of simulated ground motion intensity measures are subsequently compared against observational data. The preliminary results from this study indicate that uncertainty in the stress parameter and the rupture velocity have the most significant effect on the high frequency amplitudes. Conversely, magnitude uncertainty was found to be most influential on the spectral acceleration amplitudes at low frequencies. Further work is required to extend this preliminary study to exhaustively consider more events and to include parameter covariance. The ultimate results of this research will assist in the validation of the overall simulation method’s accuracy in capturing various rupture parameters, which is essential for the use of simulated ground motion models in probabilistic seismic hazard analysis.
A 3D high-resolution model of the geologic structure and associated seismic velocities in the Canterbury, New Zealand region is developed utilising data from depthconverted seismic reflection lines, petroleum and water well logs, cone penetration tests, and implicitly guided by existing contour maps and geologic cross sections in data sparse subregions. The model, developed using geostatistical Kriging, explicitly represents the significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age. The model is examined in the form of both geologic surface elevation contour maps as well as vertical cross sections of shear wave velocity, with the most prominent features being the Banks Peninsula Miocene-Pliocene volcanic edifice, and the Pegasus and Rakaia late Mesozoic-Neogene sedimentary basins. The adequacy of the modelled geologic surfaces is assessed through a residual analysis of point constraints used in the Kriging and qualitative comparisons with previous geologic models of subsets of the region. Seismic velocities for the lithological units between the geologic surfaces have also been derived, thus providing the necessary information for a Canterbury velocity model (CantVM) for use in physics-based seismic wave propagation. The developed model also has application for the determination of depths to specified shear wave velocities for use in empirical ground motion modelling, which is explicitly discussed via an example.
On February 22, 2011, a magnitude Mw 6.2 earthquake affected the Canterbury region, New Zealand, resulting in many fatalities. Liquefaction occurred across many areas, visible on the surface as ‘‘sand volcanoes’’, blisters and subsidence, causing significant damage to buildings, land and infrastructure. Liquefaction occurred at a number of sites across the Christchurch Boys High School sports grounds; one area in particular contained a piston ground failure and an adjacent silt volcano. Here, as part of a class project, we apply near-surface geophysics to image these two liquefaction features and determine whether they share a subsurface connection. Hand auger results enable correlation of the geophysical responses with the subsurface stratigraphy. The survey results suggest that there is a subsurface link, likely via a paleo-stream channel. The anomalous responses of the horizontal loop electromagnetic survey and electrical resistivity imaging highlight the disruption of the subsurface electrical properties beneath and between the two liquefaction features. The vertical magnetic gradient may also show a subtle anomalous response in this area, however the results are inconclusive. The ground penetrating radar survey shows disruption of the subsurface stratigraphy beneath the liquefaction features, in particular sediment mounding beneath the silt ejection (‘‘silt volcano’’) and stratigraphic disruption beneath the piston failure. The results indicate how near-surface geophysics allow the characteristics of liquefaction in the subsurface to be better understood, which could aid remediation work following liquefaction-induced land damage and guide interpretation of geophysical surveys of paleoliquefaction features.
his poster presents the ongoing development of a 3D Canterbury seismic velocity model which will be used in physics-based hybrid broadband ground motion simulation of the 2010-2011 Canterbury earthquakes. Velocity models must sufficiently represent critical aspects of the crustal structure over multiple length scales which will influence the results of the simulations. As a result, numerous sources of data are utilized in order to provide adequate resolution where necessary. Figure 2: (a) Seismic reflection line showing P-wave velocities and significant geologic horizons (Barnes et al. 2011), and (b) Shear wave profiles at 10 locations (Stokoe et al. 2013). Figure 4: Cross sections of the current version of the Canterbury velocity model to depths of 10km as shown in Figure 1: (a) at a constant latitude value of -43.6˚, and (b) at a constant longitude value of 172.64˚. 3. Ground Surface and Geologic Horizon Models Figure 3: (a) Ground surface model derived from numerous available digital elevation models, and (b) Base of the Quaternary sediments derived from structural contours and seismic reflection line elevations. The Canterbury region has a unique and complex geology which likely has a significant impact on strong ground motions, in particular the deep and loose deposits of the Canterbury basin. The Canterbury basin has several implications on seismic wave phenomena such as long period ground motion amplification and wave guide effects. Using a realistic 3D seismic velocity model in physics-based ground motion simulation will implicitly account for such effects and the resultant simulated ground motions can be studied to gain a fundamental understanding of the salient ground motion phenomena which occurred during the Canterbury earthquakes, and the potential for repeat occurrences in the Canterbury region. Figure 1 shows the current model domain as a rectangular area between Lat=[-43.2˚,-44.0˚], and Lon=[171.5˚,173.0˚]. This essentially spans the area between the foot of the Southern Alps in the North West to Banks Peninsula in the East. Currently the model extends to a depth of 50km below sea level.
This poster provides a comparison between the strong ground motions observed in the 22 February 2011 Mw6.3 Christchurch earthquake with those observed in Tokyo during the 11 March 2011 Mw9.0 Tohoku earthquake. The destuction resulting from both of these events has been well documented, although tsunami was the principal cause of damage in the latter event, and less attention has been devoted to the impact of earthquake-induced ground motions. Despite Tokyo being located over 100km from the nearest part of the causative rupture, the ground motions observed from the Tohoku earthquake were significant enough to cause structural damage and also significant liquefaction to loose reclaimed soils in Tokyo Bay. The author was fortunate enough (from the perspective of an earthquake engineer) to experience first-hand both of these events. Following the Tohoku event, the athor conducted various ground motion analyses and reconniassance of the Urayasu region in Tokyo Bay affected by liquefaction in collaboration with Prof. Kenji Ishihara. This conference is therefore a fitting opportunity in which to discuss some of authors insights obtained as a result of this first hand knowledge. Figure 1 illustrates the ground motions recorded in the Christchurch CBD in the 22 February 2011 and 4 September 2010 earthquakes, with that recorded in Tokyo Bay in the 11 March 2011 Tohoku earthquake. It is evident that these three ground motions vary widely in their amplitude and duration. The CBGS ground motion from the 22 February 2011 event has a very large amplitude (nearly 0.6g) and short duration (approx. 10s of intense shaking), as a result of the causal Mw6.3 rupture at short distance (Rrup=4km). The CBGS ground motion from the 4 September 2010 earthquake has a longer duration (approx. 30s of intense shaking), but reduced acceleration amplitude, as a result of the causal Mw7.1 rupture at a short-to-moderate distance (Rrup=14km). Finally, the Urayasu ground motion in Tokyo bay during the 11 March 2011 Tohoku earthquake exhibits an acceleration amplitude similar to the 4 September 2010 CBGS ground motion, but a significantly larger duration (approx 150s of intense shaking). Clearly, these three different ground motions will affect structures and soils in different ways depending on the vibration characteristics of the structures/soil, and the potential for strength and stiffness degradation due to cumulative effects. Figure 2 provides a comparison between the arias intensities of the several ground motion records from the three different events. It can be seen that the arias intensities of the ground motions in the Christchurch CBD from the 22 February 2011 earthquake (which is on average AI=2.5m/s) is approximately twice that from the 4 September 2010 earthquake (average AI≈1.25). This is consistent with a factor of approximately 1.6 obtained by Cubrinovski et al. (2011) using the stress-based (i.e.PGA-MSF) approach of liquefaction triggering. It can also be seen that the arias intensity of the ground motions recorded in Tokyo during the 2011 Tohoku earthquake are larger than ground motions in the Christchurch CBD from the 4 September 2011 earthquake, but smaller than those of the 22 February 2011 earthquake. Based on the arias intensity liquefaction triggering approach it can therefore be concluded that the ground motion severity, in terms of liquefaction potential, for the Tokyo ground motions is between those ground motions in Christchurch CBD from the 4 September 2010 and 22 February 2011 events.
This poster provides a summary of the development of a 3D shallow (z<40m) shear wave velocity (Vs) model for the urban Christchurch, New Zealand region. The model is based on a recently developed Christchurch-specific empirical correlation between Vs and cone penetration test (CPT) data (McGann et al. 2014a,b) and the large high-density database of CPT logs in the greater Christchurch urban area (> 15,000 logs as of 01/01/2014). In particular, the 3D model provides shear wave velocities for the surficial Springston Formation, Christchurch Formation, and Riccarton gravel layers which generally comprise the upper 40m in the Christchurch urban area. Point-estimates are provided on a 200m-by- 200m grid from which interpolation to other locations can be performed. This model has applications for future site characterization and numerical modeling efforts via maps of timeaveraged Vs over specific depths (e.g. Vs30, Vs10) and via the identification of typical Vs profiles for different regions and soil behaviour types within Christchurch. In addition, the Vs model can be used to constrain the near-surface velocities for the 3D seismic velocity model of the Canterbury basin (Lee et al. 2014) currently being developed for the purpose of broadband ground motion simulation.
Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
Generalized conditional intensity measure (GCIM) method is extended to ground motion selection for scenario ruptures. Using different rupture scenarios and site conditions, various aspects of the GCIM methodology are scrutinized, including: (i) implementation of different weight vectors and the composition of the IM vector; (ii) quantifying the importance of replicate selections for different number of desired ground motions; and (iii) the effect of considering bounds on the implicit causal parameters of the prospective ground motions. Using the extended methodology, representative ground motion ensembles for several major earthquake scenarios in New Zealand are developed. Cases considered include representative ground motions for the occurrence of Alpine, Hope, and Porters Pass earthquakes in Christchurch city, and the occurrence of Wellington, Wairarapa, and Ohariu fault ruptures in Wellington city. Challenges in the development of ground motion ensembles for subduction zone earthquakes are also highlighted. The selected scenario-based ground motion sets can be used to complement ground motions which are often selected in conjunction with probabilistic seismic hazard analysis, in order to understand the performance of structures for the question “what if this fault ruptures?”
Glazing systems are non-structural elements in a building that, more often than not, appear to be given little consideration in seismic design. Recent experimental work into glazing systems at the University of Canterbury, however, has shown that glazing systems can be very susceptible to serviceability damage, defined as loss of water-tightness. The focus of this paper is to highlight the difference in vulnerability of standard and seismic glazing systems and consider the implications of this for future repair costs and losses. The paper first describes the damage states chosen for glazing units according to the repair strategies required and expected repair costs. This includes three damage states: DS1: Water Leakage, DS2: Gasket Failure and DS3: Frame/Glass Failure. Implementing modern performance-based earthquake engineering, the paper proceeds to highlight a case study comparing costs and expected losses of a standard glazing unit and a seismic glazing unit installed on a case study building. It is shown that the use of seismic glazing units is generally beneficial over time, due to the early onset of serviceability damage in standard glazing units. Finally, the paper provides suggestions for designers aimed at reducing costs related to earthquake induced repairs of glazing.
Disasters can create the equivalent of 20 years of waste in only a few days. Disaster waste can have direct impacts on public health and safety, and on the environment. The management of such waste has a great direct cost to society in terms of labor, equipment, processing, transport and disposal. Disaster waste management also has indirect costs, in the sense that slow management can slow down a recovery, greatly affecting the ability of commerce and industry to re-start. In addition, a disaster can lead to the disruption of normal solid waste management systems, or result in inappropriate management that leads to expensive environmental remediation. Finally, there are social impacts implicit in disaster waste management decisions because of psychological impact we expect when waste is not cleared quickly or is cleared too quickly. The paper gives an overview of the challenge of disaster waste management, examining issues of waste quantity and composition; waste treatment; environmental, economic, and social impacts; health and safety matters; and planning. Christchurch, New Zealand, and the broader region of Canterbury were impacted during this research by a series of shallow earthquakes. This has led to the largest natural disaster emergency in New Zealand’s history, and the management of approximately 8 million tons of building and infrastructure debris has become a major issue. The paper provides an overview of the status of disaster waste management in Christchurch as a case study. A key conclusion is the vital role of planning in effective disaster waste management. In spite of the frequency of disasters, in most countries the ratio of time spent on planning for disaster waste management to the time spent on normal waste management is extremely low. Disaster waste management also requires improved education or training of those involved in response efforts. All solid waste professionals have a role to play to respond to the challenges of disaster waste management.
This paper examines the consistency of seismicity and ground motion models, used for seismic hazard analysis in New Zealand, with the observations in the Canterbury earthquakes. An overview is first given of seismicity and ground motion modelling as inputs of probabilistic seismic hazard analysis, whose results form the basis for elastic response spectra in NZS1170.5:2004. The magnitude of earthquakes in the Canterbury earthquake sequence are adequately allowed for in the current NZ seismicity model, however the consideration of ‘background’ earthquakes as point sources at a minimum depth of 10km results in up to a 60% underestimation of the ground motions that such events produce. The ground motion model used in conventional NZ seismic hazard analysis is shown to provide biased predictions of response spectra (over-prediction near T=0.2s , and under-predictions at moderate-to-large vibration periods). Improved ground motion prediction can be achieved using more recent NZ-specific models.
The September 2010 Canterbury and February 2011 Christchurch earthquakes and associated aftershocks have shown that the isolator displacement in Christchurch Women's Hospital (Christchurch City's only base-isolated structure) was significantly less than expected. Occupant accounts of the events have also indicated that the accelerations within the hospital superstructure were larger than would usually be expected within a base-isolated structure and that residual low-level shaking lasts for a longer period of time following the strong-motion of an event than for non-isolated structures.