Search

found 163 results

Research papers, University of Canterbury Library

he strong motion station at Heathcote Valley School (HVSC) recorded unusually high peak ground accelerations (2.21g vertical and 1.41g horizontal) during the February 2011 Christchurch earthquake. Ground motions recorded at HVSC in numerous other events also exhibited consistently higher intensities compared with nearby strong motion stations. We investigated the underlying causes of such high intensity ground motions at HVSC by means of 2D dynamic finite element analyses, using recorded ground motions during the 2010-2011 Canterbury earthquake sequence. The model takes advantage of a LiDAR-based digital elevation model (DEM) to account for the surface topography, while the geometry and dynamic properties of the surficial soils are characterized by seismic cone penetration tests (sCPT) and Multi-Channel Analyses of Surface Waves (MASW). Comparisons of simulated and recorded ground motions suggests that our model performs well for distant events, while for near-field events, ground motions recorded at the adopted reference station at Lyttelton Port are not reasonable input motions for the simulation. The simulations suggest that Rayleigh waves generated at the inclined interface of the surficial colluvium and underlying volcanic rock strongly affect the ground motions recorded at HVSC, in particular, being the dominant contributor to the recorded vertical motions.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

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.

Research papers, University of Canterbury Library

Pumice materials, which are problematic from an engineering viewpoint, are widespread in the central part of the North Island. Considering the impacts of the 2010-2011 Christchurch earthquakes, a clear understanding of their properties under earthquake loading is necessary. For example, the 1987 Edgecumbe earthquake showed evidence of localised liquefaction of sands of volcanic origin. To elucidate on this, research was undertaken to investigate whether existing empirical field-based methods to evaluate the liquefaction potential of sands, which were originally developed for hard-grained soils, are applicable to crushable pumice-rich deposits. For this purpose, two sites, one in Whakatane and another in Edgecumbe, were selected where the occurrence of liquefaction was reported following the Edgecumbe earthquake. Manifestations of soil liquefaction, such as sand boils and ejected materials, have been reported at both sites. Field tests, including cone penetration tests (CPT), shear-wave velocity profiling, and screw driving sounding (SDS) tests were performed at the sites. Then, considering estimated peak ground accelerations (PGAs) at the sites based on recorded motions and possible range of ground water table locations, liquefaction analysis was conducted at the sites using available empirical approaches. To clarify the results of the analysis, undisturbed soil samples were obtained at both sites to investigate the laboratory-derived cyclic resistance ratios and to compare with the field-estimated values. Research results clearly showed that these pumice-rich soils do not fit existing liquefaction assessment frameworks and alternate methods are necessary to characterise them.

Research papers, University of Canterbury Library

Liquefaction is a phenomenon that results in a loss of strength and stability of a saturated soil mass due to dynamic excitation such as that imposed by an earthquake. The granular nature of New Zealand soils and the location of many of our cities and towns on fluvial foundations are such that the effects of liquefaction can be very important. Research was undertaken to build on the past work undertaken at the University of Canterbury studying the effects of the 1929 Murchison earthquake, the 1968 Inangahua earthquake and the 1991 Hawks Crag earthquakes on the West Coast. Additional archival information has been gathered from newspapers and reports and from discussions with people who experienced one or all of these large earthquakes that occurred on the West Coast during the 20th Century. Further, some twenty Cone Penetrometer Tests were carried out, with varying success, in Greymouth and Karamea using the Department of Civil Engineering's Drilling Rig. These, combined with the basic site investigation information, consolidate and add to the liquefaction case history data bank at the University of Canterbury. Many of the sites have liquefied in some but not all of the three earthquakes and thus provide both upper and lower bounds for the calibration of empirical models. While a lack of knowledge of the 1929 source location reduces the value of information from that event, the data form a useful set of liquefaction case histories and will become more so as further earthquakes occur. A list of critical sites for checking of the future earthquakes is provided and recommendations are made for the installation of downhole arrays of accelerometers and pore water pressure transducers at a number of sites.

Research papers, Lincoln University

The 2013 Seddon earthquake (Mw 6.5), the 2013 Lake Grassmere earthquake (Mw 6.6), and the 2016 Kaikōura earthquake (Mw 7.8) provided an opportunity to assemble the most extensive damage database to wine storage tanks ever compiled worldwide. An overview of this damage database is presented herein based on the in-field post-earthquake damage data collected for 2058 wine storage tanks (1512 legged tanks and 546 flat-based tanks) following the 2013 earthquakes and 1401 wine storage tanks (599 legged tanks and 802 flat-based tanks) following the 2016 earthquake. Critique of the earthquake damage database revealed that in 2013, 39% and 47% of the flat-based wine tanks sustained damage to their base shells and anchors respectively, while due to resilience measures implemented following the 2013 earthquakes, in the 2016 earthquake the damage to tank base shells and tank anchors of flat-based wine tanks was reduced to 32% and 23% respectively and instead damage to tank barrels (54%) and tank cones (43%) was identified as the two most frequently occurring damage modes for this type of tank. Analysis of damage data for legged wine tanks revealed that the frame-legs of legged wine tanks sustained the greatest damage percentage among different parts of legged tanks in both the 2013 earthquakes (40%) and in the 2016 earthquake (44%). Analysis of damage data and socio-economic findings highlight the need for industry-wide standards, which may have socio-economic implications for wineries.

Research papers, University of Canterbury Library

This dissertation addresses a diverse range of topics in the area of physics-based ground motion simulation with particular focus on the Canterbury, New Zealand region. The objectives achieved provide the means to perform hybrid broadband ground motion simulation and subsequently validates the simulation methodology employed. In particu- lar, the following topics are addressed: the development of a 3D seismic velocity model of the Canterbury region for broadband ground motion simulation; the development of a 3D geologic model of the interbedded Quaternary formations to provide insight on observed ground motions; and the investigation of systematic effects through ground motion sim- ulation of small-to-moderate magnitude earthquakes. The paragraphs below outline each contribution in more detail. As a means to perform hybrid broadband ground motion simulation, a 3D model of the geologic structure and associated seismic velocities in the Canterbury region is devel- oped utilising data from depth-converted 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 explicitly characterises five significant and regionally recognisable geologic surfaces that mark the boundaries between geologic units with distinct lithology and age, including the Banks Peninsula volcanics, which are noted to strongly influence seismic wave propagation. The Basement surface represents the base of the Canterbury sedimentary basin, where a large impedance contrast exists re- sulting in basin-generated waves. Seismic velocities for the lithological units between the geologic surfaces are derived from well logs, seismic reflection surveys, root mean square stacking velocities, empirical correlations, and benchmarked against a regional crustal model, thus providing the necessary information for a Canterbury velocity model for use in broadband seismic wave propagation. A 3D high-resolution model of the Quaternary geologic stratigraphic sequence in the Canterbury region is also developed utilising datasets of 527 high-quality water well logs, and 377 near-surface cone penetration test records. The model, developed using geostatistical Kriging, represents the complex interbedded regional Quaternary geology by characterising the boundaries between significant interbedded geologic formations as 3D surfaces including explicit modelling of the formation unconformities resulting from the Banks Peninsula volcanics. The stratigraphic layering present can result in complex wave propagation. The most prevalent trend observed in the surfaces was the downward dip from inland to the eastern coastline as a result of the dominant fluvial depositional environment of the terrestrial gravel formations. The developed model provides a benefi- cial contribution towards developing a comprehensive understanding of recorded ground motions in the region and also providing the necessary information for future site char- acterisation and site response analyses. To highlight the practicality of the model, an example illustrating the role of the model in constraining surface wave analysis-based shear wave velocity profiling is illustrated along with the calculation of transfer functions to quantify the effect of the interbedded geology on wave propagation. Lastly, an investigation of systematic biases in the (Graves and Pitarka, 2010, 2015) ground motion simulation methodology and the specific inputs used for the Canterbury region is presented considering 144 small-to-moderate magnitude earthquakes. In the simulation of these earthquakes, the 3D Canterbury Velocity Model, developed as a part of this dissertation, is used for the low-frequency simulation, and a regional 1D velocity model for the high-frequency simulation. Representative results for individual earthquake sources are first presented to highlight the characteristics of the small-to-moderate mag- nitude earthquake simulations through waveforms, intensity measure scaling with source- to-site distance, and spectral bias of the individual events. Subsequently, a residual de- composition is performed to examine the between- and within-event residuals between observed data, and simulated and empirical predictions. By decomposing the residuals into between- and within-event residuals, the biases in source, path and site effects, and their causes, can be inferred. The residuals are comprehensively examined considering their aggregated characteristics, dependence on predictor variables, spatial distribution, and site-specific effects. The results of the simulation are also benchmarked against empir- ical ground motion models, where their similarities manifest from common components in their prediction. Ultimately, suggestions to improve the predictive capability of the simulations are presented as a result of the analysis.

Research papers, Lincoln University

Nature has endowed New Zealand with unique geologic, climatic, and biotic conditions. Her volcanic cones and majestic Southern Alps and her verdant plains and rolling hills provide a landscape as rugged and beautiful as will be found anywhere. Her indigenous fauna and flora are often quite different from that of the rest of the world and consequently have been of widespread interest to biologists everywhere. Her geologic youth and structure and her island climate, in combination with the biological resources, have made a land which is ecologically on edge. These natural endowments along with the manner in which she has utilized her land, have given New Zealand some of the most spectacular and rapid erosion to be found. It is quite evident that geologic and climatic conditions combine to give unusually high rates of natural erosion. Present topographic features indicate the past occurrence of large-scale flooding as well. Prior to the arrival of the Maori, it is very likely that most of the land mass of New Zealand below present bush lines was covered with indigenous bush or forest. Forest fires of a catastrophic nature undoubtedly occurred as a result of lightning, and volcanic eruptions. The exposed soils left by these catastrophes contributed to natural deterioration. While vast areas of forest cover were destroyed, they probably were healed by nature with forest or with grass or herbaceous cover. Further, it is probable that large areas in the mountains were, as they are now, subject to landslides and slipping due to earthquakes and excessive local rainfall. Again, the healing process was probably rapid in most of such exposed areas.

Research papers, University of Canterbury Library

The Mw 6.2 February 22nd 2011 Christchurch earthquake (and others in the 2010-2011 Canterbury sequence) provided a unique opportunity to study the devastating effects of earthquakes first-hand and learn from them for future engineering applications. All major events in the Canterbury earthquake sequence caused widespread liquefaction throughout Christchurch’s eastern suburbs, particularly extensive and severe during the February 22nd event. Along large stretches of the Avon River banks (and to a lesser extent along the Heathcote) significant lateral spreading occurred, affecting bridges and the infrastructure they support. The first stage of this research involved conducting detailed field reconnaissance to document liquefaction and lateral spreading-induced damage to several case study bridges along the Avon River. The case study bridges cover a range of ages and construction types but all are reinforced concrete structures which have relatively short, stiff decks. These factors combined led to a characteristic deformation mechanism involving deck-pinning and abutment back-rotation with consequent damage to the abutment piles and slumping of the approaches. The second stage of the research involved using pseudo-static analysis, a simplified seismic modelling tool, to analyse two of the bridges. An advantage of pseudo-static analysis over more complicated modelling methods is that it uses conventional geotechnical data in its inputs, such as SPT blowcount and CPT cone resistance and local friction. Pseudo-static analysis can also be applied without excessive computational power or specialised knowledge, yet it has been shown to capture the basic mechanisms of pile behaviour. Single pile and whole bridge models were constructed for each bridge, and both cyclic and lateral spreading phases of loading were investigated. Parametric studies were carried out which varied the values of key parameters to identify their influence on pile response, and computed displacements and damages were compared with observations made in the field. It was shown that pseudo-static analysis was able to capture the characteristic damage mechanisms observed in the field, however the treatment of key parameters affecting pile response is of primary importance. Recommendations were made concerning the treatment of these governing parameters controlling pile response. In this way the future application of pseudo-static analysis as a tool for analysing and designing bridge pile foundations in liquefying and laterally spreading soils is enhanced.

Research papers, University of Canterbury Library

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

Research papers, University of Canterbury Library

In the period between September 2010 and December 2011, Christchurch was shaken by a series of strong earthquakes including the MW7.1 4 September 2010, Mw 6.2 22 February 2011, MW6.2 13 June 2011 and MW6.0 23 December 2011 earthquakes. These earthquakes produced very strong ground motions throughout the city and surrounding areas that resulted in soil liquefaction and lateral spreading causing substantial damage to buildings, infrastructure and the community. The stopbank network along the Kaiapoi and Avon River suffered extensive damage with repairs projected to take several years to complete. This presented an opportunity to undertake a case-study on a regional scale of the effects of liquefaction on a stopbank system. Ultimately, this information can be used to determine simple performance-based concepts that can be applied in practice to improve the resilience of river protection works. The research presented in this thesis draws from data collected following the 4th September 2010 and 22nd February 2011 earthquakes. The stopbank damage is categorised into seven key deformation modes that were interpreted from aerial photographs, consultant reports, damage photographs and site visits. Each deformation mode provides an assessment of the observed mechanism of failure behind liquefaction-induced stopbank damage and the factors that influence a particular style of deformation. The deformation modes have been used to create a severity classification for the whole stopbank system, being ‘no or low damage’ and ‘major or severe damage’, in order to discriminate the indicators and factors that contribute to ‘major to severe damage’ from the factors that contribute to all levels of damage a number of calculated, land damage, stopbank damage and geomorphological parameters were analysed and compared at 178 locations along the Kaiapoi and Avon River stopbank systems. A critical liquefiable layer was present at every location with relatively consistent geotechnical parameters (cone resistance (qc), soil behaviour type (Ic) and Factor of Safety (FoS)) across the study site. In 95% of the cases the critical layer occurred within two times the Height of the Free Face (HFF,). A statistical analysis of the geotechnical factors relating to the critical layer was undertaken in order to find correlations between specific deformation modes and geotechnical factors. It was found that each individual deformation mode involves a complex interplay of factors that are difficult to represent through correlative analysis. There was, however, sufficient data to derive the key factors that have affected the severity of deformation. It was concluded that stopbank damage is directly related to the presence of liquefaction in the ground materials beneath the stopbanks, but is not critical in determining the type or severity of damage, instead it is merely the triggering mechanism. Once liquefaction is triggered it is the gravity-induced deformation that causes the damage rather than the shaking duration. Lateral spreading and specifically the depositional setting was found to be the key aspect in determining the severity and type of deformation along the stopbank system. The presence or absence of abandoned or old river channels and point bar deposits was found to significantly influence the severity and type of deformation. A review of digital elevation models and old maps along the Kaiapoi River found that all of the ‘major to severe’ damage observed occurred within or directly adjacent to an abandoned river channel. Whilst a review of the geomorphology along the Avon River showed that every location within a point bar deposit suffered some form of damage, due to the depositional environment creating a deposit highly susceptible to liquefaction.

Research papers, University of Canterbury Library

Previous earthquakes demonstrated destructive effects of soil-structure interaction on structural response. For example, in the 1970 Gediz earthquake in Turkey, part of a factory was demolished in a town 135 km from the epicentre, while no other buildings in the town were damaged. Subsequent investigations revealed that the fundamental period of vibration of the factory was approximately equal to that of the underlying soil. This alignment provided a resonance effect and led to collapse of the structure. Another dramatic example took place in Adapazari, during the 1999 Kocaeli earthquake where several foundations failed due to either bearing capacity exceedance or foundation uplifting, consequently, damaging the structure. Finally, the Christchurch 2012 earthquakes have shown that significant nonlinear action in the soil and soil-foundation interface can be expected due to high levels of seismic excitation and spectral acceleration. This nonlinearity, in turn, significantly influenced the response of the structure interacting with the soil-foundation underneath. Extensive research over more than 35 years has focused on the subject of seismic soil-structure interaction. However, since the response of soil-structure systems to seismic forces is extremely complex, burdened by uncertainties in system parameters and variability in ground motions, the role of soil-structure interaction on the structural response is still controversial. Conventional design procedures suggest that soil-structure interaction effects on the structural response can be conservatively ignored. However, more recent studies show that soil-structure interaction can be either beneficial or detrimental, depending on the soil-structure-earthquake scenarios considered. In view of the above mentioned issues, this research aims to utilise a comprehensive and systematic probabilistic methodology, as the most rational way, to quantify the effects of soil-structure interaction on the structural response considering both aleatory and epistemic uncertainties. The goal is achieved by examining the response of established rheological single-degree-of-freedom systems located on shallow-foundation and excited by ground motions with different spectral characteristics. In this regard, four main phases are followed. First, the effects of seismic soil-structure interaction on the response of structures with linear behaviour are investigated using a robust stochastic approach. Herein, the soil-foundation interface is modelled by an equivalent linear cone model. This phase is mainly considered to examine the influence of soil-structure interaction on the approach that has been adopted in the building codes for developing design spectrum and defining the seismic forces acting on the structure. Second, the effects of structural nonlinearity on the role of soil-structure interaction in modifying seismic structural response are studied. The same stochastic approach as phase 1 is followed, while three different types of structural force-deflection behaviour are examined. Third, a systematic fashion is carried out to look for any possible correlation between soil, structural, and system parameters and the degree of soil-structure interaction effects on the structural response. An attempt is made to identify the key parameters whose variation significantly affects the structural response. In addition, it is tried to define the critical range of variation of parameters of consequent. Finally, the impact of soil-foundation interface nonlinearity on the soil-structure interaction analysis is examined. In this regard, a newly developed macro-element covering both material and geometrical soil-foundation interface nonlinearity is implemented in a finite-element program Raumoko 3D. This model is then used in an extensive probabilistic simulation to compare the effects of linear and nonlinear soil-structure interaction on the structural response. This research is concluded by reviewing the current design guidelines incorporating soil-structure interaction effects in their design procedures. A discussion is then followed on the inadequacies of current procedures based on the outcomes of this study.

Research papers, University of Canterbury Library

In 2010 and 2011 Christchurch, New Zealand experienced a series of earthquakes that caused extensive damage across the city, but primarily to the Central Business District (CBD) and eastern suburbs. A major feature of the observed damage was extensive and severe soil liquefaction and associated ground damage, affecting buildings and infrastructure. The behaviour of soil during earthquake loading is a complex phenomena that can be most comprehensively analysed through advanced numerical simulations to aid engineers in the design of important buildings and critical facilities. These numerical simulations are highly dependent on the capabilities of the constitutive soil model to replicate the salient features of sand behaviour during cyclic loading, including liquefaction and cyclic mobility, such as the Stress-Density model. For robust analyses advanced soil models require extensive testing to derive engineering parameters under varying loading conditions for calibration. Prior to this research project little testing on Christchurch sands had been completed, and none from natural samples containing important features such as fabric and structure of the sand that may be influenced by the unique stress-history of the deposit. This research programme is focussed on the characterisation of Christchurch sands, as typically found in the CBD, to facilitate advanced soil modelling in both res earch and engineering practice - to simulate earthquake loading on proposed foundation design solutions including expensive ground improvement treatments. This has involved the use of a new Gel Push (GP) sampler to obtain undisturbed samples from below the ground-water table. Due to the variable nature of fluvial deposition, samples with a wide range of soil gradations, and accordingly soil index properties, were obtained from the sampling sites. The quality of the samples is comprehensively examined using available data from the ground investigation and laboratory testing. A meta-quality assessment was considered whereby a each method of evaluation contributed to the final quality index assigned to the specimen. The sampling sites were characterised with available geotechnical field-based test data, primarily the Cone Penetrometer Test (CPT), supported by borehole sampling and shear-wave velocity testing. This characterisation provides a geo- logical context to the sampling sites and samples obtained for element testing. It also facilitated the evaluation of sample quality. The sampling sites were evaluated for liquefaction hazard using the industry standard empirical procedures, and showed good correlation to observations made following the 22 February 2011 earthquake. However, the empirical method over-predicted liquefaction occurrence during the preceding 4 September 2010 event, and under-predicted for the subsequent 13 June 2011 event. The reasons for these discrepancies are discussed. The response of the GP samples to monotonic and cyclic loading was measured in the laboratory through triaxial testing at the University of Canterbury geomechanics laboratory. The undisturbed samples were compared to reconstituted specimens formed in the lab in an attempt to quantify the effect of fabric and structure in the Christchurch sands. Further testing of moist tamped re- constituted specimens (MT) was conducted to define important state parameters and state-dependent properties including the Critical State Line (CSL), and the stress-strain curve for varying state index. To account for the wide-ranging soil gradations, selected representative specimens were used to define four distinct CSL. The input parameters for the Stress-Density Model (S-D) were derived from a suite of tests performed on each representative soil, and with reference to available GP sample data. The results of testing were scrutinised by comparing the data against expected trends. The influence of fabric and structure of the GP samples was observed to result in similar cyclic strength curves at 5 % Double Amplitude (DA) strain criteria, however on close inspection of the test data, clear differences emerged. The natural samples exhibited higher compressibility during initial loading cycles, but thereafter typically exhibited steady growth of plastic strain and excess pore water pressure towards and beyond the strain criteria and initial liquefaction, and no flow was observed. By contrast the reconstituted specimens exhibited a stiffer response during initial loading cycles, but exponential growth in strains and associated excess pore water pressure beyond phase-transformation, and particularly after initial liquefaction where large strains were mobilised in subsequent cycles. These behavioural differences were not well characterised by the cyclic strength curve at 5 % DA strain level, which showed a similar strength for both GP samples and MT specimens. A preliminary calibration of the S-D model for a range of soil gradations is derived from the suite of laboratory test data. Issues encountered include the influence of natural structure on the peak-strength–state index relationship, resulting in much higher peak strengths than typically observed for sands in the literature. For the S-D model this resulted in excessive stiffness to be modelled during cyclic mobility, when the state index becomes large momentarily, causing strain development to halt. This behaviour prevented modelling the observed re- sponse of silty sands to large strains, synonymous with “liquefaction”. Efforts to reduce this effect within the current formulation are proposed as well as future research to address this issue.