The Coupled Influence of Relative Density,CSR, Plasticity and Content of Fines on Cyclic Liquefaction Resistance of Sands

In this study, stress-controlled cyclic simple shear tests were performed on sand specimens with up to 10% silt or clay contents, and coupled effects of plasticity, fines content (FC), relative density (D_R) and CSR (cyclic stress ratio) were investigated. The results demonstrated that for sands with low fines content, reasonable trends were obtained when the packing index control parameter was selected as D_R. Also, clean sand specimens demonstrated highest liquefaction strength compared with that of sands with fines up to 10% FC. The effect of fines’ plasticity became apparent as the FC increases and CSR reduces in relatively denser specimens.     

Editorial I

Communication strategies emphasize concerns for content (what is said) and process (the way things are said). Scientists have a responsibility to communicate the findings of their research, enhancing prospects that their insights can meaningfully inform management practice. When used effectively, principles from geomorphology provide critical guidance for environmental management. Three key geomorphic messages are outlined here: Respect diversity through communicating spatial and temporal controls upon landscape character and behaviour; Work with nature in conveying the range of behaviour of any given system; and Be proactive by determining the trajectory of landscape change. Uptake of geomorphic understanding reflects, among many factors, our ability to communicate findings clearly. In communicating these messages, different approaches are likely to prove more effective when addressing differing audiences, whether school/university classes, stakeholders, managers, politicians, or the broader community. Collective approaches to learning through dialogue are encouraged as they promote deeper learning, prospectively enhancing the uptake of geomorphic understanding, thereby promoting healthier environmental futures.     

Quantification of Linked Column Frame Seismic Performance Factors for Use in Seismic Design

Previous research has proposed the Linked Column Frame (LCF) as a lateral load-resisting system capable of providing rapid return to occupancy for buildings impacted by moderate earthquake events and collapse prevention in very large events. The LCF consists of flexible moment frames (MF) and linked columns (LC), which are closely spaced dual columns interconnected with bolted links. The linked columns (LC) are designed to limit seismic forces and provide energy dissipation through yielding of the links, while preventing damage to the moment frame under certain earthquake hazard levels. The proposed design procedure ensures the links of the linked column yield at a significantly lower story drift than the beams of the moment frame, enabling design of this system for two distinct performance states: rapid repair, where only link damage occurs and quick link replacement is possible; and collapse prevention, where both the linked column and moment frame may be damaged.

Here, the seismic performance factors for the LCF system, including the response modification factor, R, the system over-strength factor, Ω₀, and the deflection amplification factor, C_d, are established following the procedures described in FEMA P695 [2009]. These parameters are necessary for inclusion of the system in the building code. This work describes the development of archetype structures, numerical models of the LCF systems, incremental dynamic analyses, and interpretation of the results. From the results, it is recommended that R, Ω₀, and C_d values of 8, 3, and 5.5 be used for seismic design of the LCF system. A height limit of 35 m (115ft) is recommended at this time as taller LCFs are not considered in this study.     

An Investigation of Characteristic Periods of Seismic Ground Motions

The design seismic base shear was obtained from the spectral elastic acceleration S_a divided by a system behavior factor R, accounting for ductility and overstrength. The behavior factor is currently taken as a constant for a given type of structures in various codes regardless of structural periods. In fact, the behavior factor is also a spectrum varying with the natural periods of structures. In order to understand the relationship between the spectral values and the corresponding characteristic periods in these two spectra, S_a and R_μ, this article carries out an investigation into the characteristic periods of 370 seismic ground motions from 4 site types. It is found that the periods T_ga at which the peak values appear in the S_a spectra are much less than the periods T _gR at which the R_μ spectra take a maximum value. Two characteristic periods are necessary to determine the seismic action if a more elaborate procedure is required in practice. Statistical study on these two periods is carried out for the 370 records, and results are presented. For site types A–D, the ratio of these two periods has a statistically averaged value of 5.5–6.7.

The maximum input energy S _EI , relative velocity S _v , power density P _SD , and the Fourier amplitude F _S spectra were constructed to determine their characteristic periods, respectively. These four spectra predict similar characteristic periods to T _gR . T _gR is very close to the characteristic period T _gd of the elastic displacement spectra.

Analysis of SDOF systems under combined harmonic excitations shows that the S_a spectrum is more sensitive to high-frequency excitations, while the displacement spectrum is more sensitive to long period excitations. For the elastic-plastic S_a spectra, peak values tend to appear at shorter periods even the amplitudes of the longer periods are greater than that of the shorter period. This provides an explanation on different characteristic periods in the S_a and R_μ spectra.     

DUCTILITY REDUCTION FACTOR OF MDOF SHEAR-BUILDING STRUCTURES

In this paper, the results of recent studies on inelastic seismic response of MDOF shear-building structures are presented. In the last few decades, the concept of response modification factor R has been introduced and developed to account for inelastic nonlinear behaviour of structures under earthquakes. In this paper, an attempt has been made to adjust and extend this concept through introducing a modifying factor R _T . This factor is used for dynamic analysis of MDOF structures, including the calculation of inelastic response spectra. Sensitivity analysis was carried out to identify the parameters that have influence on R _T . It has been demonstrated that R _T is predominantly a function of number of stories, and accordingly a relationship has been suggested. Finally, an approximate approach has been developed for evaluating the seismic strength and ductility demands of MDOF structures.     

CALIBRATION OF FORCE REDUCTION FACTORS OF RC BUILDINGS

A comprehensive study is undertaken to assess and calibrate the force reduction factors (R) adopted in modern seismic codes. Refined expressions are employed to calculate the R factors “supply” for 12 buildings of various characteristics represent a wide range of medium-rise RC buildings. The “supply” values are then compared with the “design” and “demand” recommended in the literature. A comprehensive range of response criteria at the member and storey levels, including shear as a failure criterion, alongside a detailed modelling approach and an extensively verified analytical tool are utilised. A rigorous technique is employed to evaluate R factors, including inelastic pushover and incremental dynamic collapse analyses employing eight natural and artificial records. In the light of the information obtained from more than 1500 inelastic analyses, it is concluded that including shear and vertical motion in assessment and calculations of R factors is necessary. Force reduction factors adopted by the design code (Eurocode 8) are over-conservative and can be safely increased, particularly for regular frame structures designed to lower PGA and higher ductility levels.     

Organic–inorganic rock–fluid interactions in stylolitic micro‐environments of carbonate rocks: a FIB‐TEM study combined with a hydrogeochemical modelling approach

Stylolites and the interfaces to the host limestone have been investigated by means of a multidisciplinary analytical approach (thin section microscopy, FIB‐TEM, organic geochemistry and petrography). Carbonate dissolution assuming different boundary conditions was simulated by applying a generic hydrogeochemical modelling approach. It is the conceptual approach to characterize and quantify traceable organic–inorganic interactions in stylolites dependent on organic matter type and its thermal maturity, and to follow stylolite formation in carbonates as result of organic matter reactivity rather than pressure solution as a main control. The investigated stylolite samples are of Upper Permian (Lopingian, Zechstein), Middle Triassic (Muschelkalk) and Late Cretaceous (Maastrichtian) age and always contain marine organic matter. The thermal maturity of the organic matter ranges from the pre‐oil generation zone (0.4–0.5% R_r) to the stage of dry gas generation (>1.3% R_r). The results of the generic hydrogeochemical modelling indicate a sharp increase of calcite dissolution and the beginning of stylolite formation at approximately 40°C, which is equivalent to a depth of less than 800 m under hydrostatic conditions considering a geothermal gradient of 30°C and a surface mean temperature of 20°C. This temperature corresponds to the pre‐oil window when kerogens release an aqueous fluid enriched in carbon dioxide and organic acids. This aqueous fluid may change the existing pore water pH or alkalinity and causes dissolution of carbonate, feldspar and quartz, and clay mineral precipitation along the stylolite. Dissolution of limestone and dolostone leads to reprecipitation of calcite or dolomite opposite of the dissolution side, which indicates only localized mass redistribution. All these integrated hydrogeochemical processes are coupled to the generation of water during organic matter maturation. In all of the calculated hydrogeochemical scenarios, H₂O is a reaction product and its formation supports the suggested hypothesis.     

Permeability of the Mercia Mudstone: suitability as caprock to carbon capture and storage sites

The Upper Triassic Mercia Mudstone is the caprock to potential carbon capture and storage (CCS) sites in porous and permeable Lower Triassic Sherwood Sandstone reservoirs and aquifers in the UK (primarily offshore). This study presents direct measurements of vertical (k_v) and horizontal (k_h) permeability of core samples from the Mercia Mudstone across a range of effective stress conditions to test their caprock quality and to assess how they will respond to changing effective stress conditions that may occur during CO₂ injection and storage. The Mercia samples analysed were either clay‐rich (muddy) siltstones or relatively clean siltstones cemented by carbonate and gypsum. Porosity is fairly uniform (between 7.4 and 10.7%). Porosity is low either due to abundant depositional illite or abundant diagenetic carbonate and gypsum cements. Permeability values are as low as 10^?20 m² (10nD), and therefore, the Mercia has high sealing capacity. These rocks have similar horizontal and vertical permeabilities with the highest k_h/k_v ratio of 2.03 but an upscaled k_h/k_v ratio is 39, using the arithmetic mean of k_h and the harmonic mean of k_v. Permeability is inversely related to the illite clay content; the most clay‐rich (illite‐rich) samples represent very good caprock quality; the cleaner Mercia Mudstone samples, with pore‐filling carbonate and gypsum cements, represent fair to good caprock quality. Pressure sensitivity of permeability increases with increasing clay mineral content. As pore pressure increases during CO₂ injection, the permeability of the most clay‐rich rocks will increase more than carbonate‐ and gypsum‐rich rocks, thus decreasing permeability heterogeneity. The best quality Mercia Mudstone caprock is probably not geochemically sensitive to CO₂ injection as illite, the cause of the lowest permeability, is relatively stable in the presence of CO₂–water mixtures.     

Inelastic Displacement Ratios for Seismic Assessment of Structures Subjected to Forward-Directivity Near-Fault Ground Motions

This article presents results of a statistical study focused on evaluating inelastic displacement ratios (i.e., ratio of maximum inelastic displacement with respect to maximum elastic displacement demand) of degrading and non degrading single-degree-of-freedom (SDOF) systems subjected to forward-directivity near-fault ground motions. C_R spectra are computed for normalized periods of vibration with respect to the predominant period of the ground motion to provide a better ground motion characterization. This period normalization allows reducing the record-to-record variability in the estimation of C_R. An equation to obtain estimates of C_R for the seismic assessment of structures exposed to forward-directivity near-fault ground motions is proposed.     

ATTENUATION RELATIONSHIP FOR LOW MAGNITUDE EARTHQUAKES USING STANDARD SEISMOMETRIC RECORDS

Northwestern Italian weak-motion data were used to study attenuation characteristics of horizontal peak ground acceleration (PGA) and horizontal peak ground velocity (PGV) from earthquakes of local magnitudes (M _l ) up to 5.1. Data have been provided by the RSNI (Regional seismic network of Northwestern Italy) and RSLG (Regional seismic network of Lunigiana-Garfagnana) waveform database. The database consists of more than 14000 horizontal components recorded in the period 1999-2002 by both broadband and enlarged band seismometers. The accuracy of the procedure used to extract PGA values from the velocity recordings was verified comparing observed and derived PGA values at station STV2, which was equipped with both a temporary K2 Kinemctrtcs accelerometer and Guralp CMG40 broadband sensor. The attenuation of both peak ground acceleration and peak ground velocity was found to be logarithmically distributed with a strong attenuation for low distances (less than 50 km) and low M _l values (<3.0). The resulting equations are:

Log(PGA)=?3.19+0.87M?0.042M ²?1.92 Log(R)+0.249S,

Log(PGA)=?4.23+0.76M?0.018M²?1.56 Log(R)+0.230S,

where PGA is expressed in g, PGV is expressed in m/s, M is local magnitude, R is the hypocentral distance in kilometers and S is a dummy variable assuming values of 0 and 1 for rock and soil respectively. For increasing distance and magnitude, both PGA and PGV values show a linear distribution. The validity range of the obtained attenuation relationships is 0–200 km for distances and M _l up to 4.5. Sensitivity studies performed by analysis of residuals, showed that predicted PGA and PGV values are stable with respect to reasonable variations of the model and distances providing the data. Comparisons with attenuation relationships proposed for Italian region, derived from strong motion records, are also presented.     

MATHEMATICAL MODELLING OF THE RELATIONSHIP BETWEEN NEOLITHIC SITES AND THE RIVERS IN XI’AN (SHAANXI PROVINCE,CHINA)*

Based on analyses of the ecological environments of the Neolithic sites along the Chan and Ba Rivers in the Xi’an area (Shaanxi Province, China), it was found that the characteristic factors of the sites had a very close relationship with the nearby river. The characteristic factors of Neolithic sites have been analysed and defined as mathematical parameters by statistical methods. A mathematical model of the human–land relationship between Neolithic sites and rivers has been established using the stepwise regression procedure REGRESS in SPSS^®. The mathematical model is Y = 496.153 – 4.833X₁ + 113.692X3 + 725.097X₅ – 27.683X₃X₄. It shows that the distance (Y) from a site to a river has a marked relationship with the site area (X₁), the river flux (X₃), the site location (X₅) and river ratio drop (X₄). The model shows that the relationship between the sites and the river is in accordance with the conclusions of archaeological research. The mathematical model not only offers theoretical guidance for the archaeological excavation of Neolithic sites in Xi’an, but also contributes to the archaeological environmental research of Xi’an Neolithic sites.      

Tectonic control over large-scale diffuse degassing in eastern Sicily (Italy)

Eastern Sicily (southern Italy) is characterised by the presence of many natural gas emissions (mofettes, mud volcanoes). These gases are mostly carbon dioxide and methane, with minor amounts of helium, hydrogen, carbon monoxide and hydrocarbons. In this study, the extent and orientation of soil gas anomalies (He and CO₂) were investigated on a wide area (approximately 110 km²) located just SW of Mt. Etna. From a structural point of view, this area lays on a typical foredeep–foreland system that marks the boundary between the southern part of the Eurasian plate and the northern part of the African plate in the central Mediterranean. No tectonic structure was revealed in this area by surface geological surveys. Very high soil emissions were found, and their spatial pattern reveals the existence of some active faults all directed about N50°E. This direction coincides with that of two major fault systems that cut eastern Sicily and are evident, respectively, NE and SW of the study area. Soil gas data suggest that these fault systems are the expression of a single continuous structural line which is probably responsible for the past and present magma uprise in eastern Sicily. Isotopic values of carbon of CO₂ suggest a minor contribution of organic carbon. Moreover, in the highest degassing sites the isotopic values of He found in association with CO₂ (He abundance = 11–70 p.p.m.; R/R_a between 6.0 and 6.2) suggest that both gases are mantle derived. The extent of the areas affected by high gas emissions and the amounts of deep CO₂ emitted in the investigated area (several hundred tonnes per day) may provide additional supporting evidence of a mantle upwelling taking place beneath this region.     

A palynological investigation of plesiosaur-bearing rocks from the Upper Cretaceous Tahora Formation,Mangahouanga, New Zealand

High-palaeolatitude plesiosaur, mosasaur and, more rarely, dinosaur fossils are well known from the Maungataniwha Sandstone Member of the Tahora Formation in Mangahouanga Stream, Hawkes Bay, New Zealand. A palynological investigation of strata exposed along Mangahouanga Stream and of transported boulders hosting vertebrate fossils reveals well-preserved assemblages dominated by terrestrial pollen and spores but also including marine dinoflagellate cysts in some samples. The palynofacies are strongly dominated by wood fragments including charcoal; one outcrop sample and the sample taken from a boulder hosting plesiosaur vertebrae contain entirely terrestrially derived palynoassemblages, suggesting a freshwater habitat for at least some of the plesiosaurs. The host unit spans the Santonian to lowermost Maastrichtian, while the key pollen taxa Nothofagidites senectus and Tricolpites lilliei, together with the dinocyst Isabelidinium pellucidum and the megaspore Grapnelispora evansii, indicate a late Campanian to early Maastrichtian age for the fossiliferous boulders. The palynoflora indicates a mixed local vegetation dominated by podocarp conifers and angiosperms with a significant tree-fern subcanopy. The presence of taxa with modern temperate distributions, such as Nothofagus (southern beech), Proteaceae and Cyatheaceae (tree-ferns), indicates a mild-temperate climate and lack of severe winter freezing during the latest Cretaceous, providing an ecosystem that most probably made it possible for polar dinosaurs to overwinter in this part of the world.     

Correlations Between Frequency Content Indicators of Strong Ground Motions and PGV

The frequency content of ground motions seems to be one of the most important parameters to explain the structural damage experienced during worldwide strong earthquakes. The frequency content of ground motions can be characterized by various stochastic and/or deterministic indicators: the frequency bandwidth indicator ? (Cartwright & Longuet-Higgins) related to the power spectral density function and, respectively, the control (corner) period T_c of the structural response spectra or the mean period T_M . Peak ground velocity (PGV) and the ratio PGA/PGV can be used as either damage potential parameters or frequency content indicators. A comparative analysis of stochastic and deterministic frequency content indicators and of PGV is applied to a set of 30 strong ground motion records having peak ground acceleration (PGA) from 0.2–0.8 g and recorded on 4 continents during the last 70 years.     

Darfur and the failure of the responsibility to protect

When official representatives of more than 170 countries adopted the principle of the ‘responsibility to protect’ (R₂P) at the September 2005 World Summit, Darfur was quickly identified as the test case for this new doctrine. The general verdict is that the international community has failed the test due to lack of political will. This article argues that the failure is real but that it is more fundamentally located within the doctrine of R₂P itself. Fulfilling the aspiration of R₂P demands an international protection capability that does not exist now and cannot be realistically expected. The critical weakness in R₂P is that the ‘responsibility to react’ has been framed as coercive protection, which attempts to be a middle way between classic peacekeeping and outright military intervention that can be undertaken without the consent of the host government. Thus far, theoretical and practical attempts to create this intermediate space for coercive protection have failed to resolve basic strategic and operational issues. In addition, the very act of raising the prospect of external military intervention for human protection purposes changes and distorts the political process and can in fact make a resolution more difficult. Following an introductory section that provides background to the war in Darfur and international engagement, this article examines the debates over the R₂P that swirled around the Darfur crisis and operational concepts developed for the African Union Mission in Sudan (AMIS) and its hybrid successor, the UN–African Union Mission in Darfur (UNAMID), especially during the Abuja peace negotiations. Three operational concepts are examined: ceasefire, disarmament and civilian protection. Unfortunately, the international policy priority o bringing UN troops to Darfur had an adverse impact on the Darfur peace talks without grappling with the central question of what international forces would do to resolve the crisis. Advocacy for the R₂P set an unrealistic ideal which became the enemy of achievable goals.     

SEISMIC DESIGN REGULATION CODES: CONTRIBUTION OF K-NET DATA TO SITE EFFECT EVALUATION

The purpose of this study is to compare the site effect section of building codes (EC8 and UBC97) with the set of data provided by the Kyoshin network. In order to obtain a set of site coefficients and spectral shapes, we have first deduced an attenuation law for both horizontal and vertical motion. Site conditions are represented by the shear velocity averaged over the upper 30 m (V _s ³⁰). Our site classification (4 categories similar to those proposed in the new ECS and the UBC97) is based on borehole investigations at every station. This classification has permitted to distinguish clearly four response spectra which demonstrates the efficiency of V _s ³⁰ as characterising site conditions. Our law is then used to test site coefficients and spectral shapes of building codes ECS and UBC97. Concerning spectral shapes and site coefficients, our results are found to be in good agreement with EC8 and UBC97 only if category B (400<V _s ³⁰<800 m/s) is taken as reference. We also conclude that a site which is characterised as “rock” on geological criteria can not generally be classified in category A (V _s ³⁰>800 m/s). This suggests that classification in category A should be based only on field measurements. Concerning vertical motion, our analysis of the K-NET data shows that the ratio av/ah (vertical peak ground acceleration over horizontal peak ground acceleration) is between 0.50 and 0.68.     

Fracture‐focused fluid flow in an acid and redox‐influenced system: diagenetic controls on cement mineralogy and geomorphology in the Navajo Sandstone

Differential cement mineralogy is influenced by depositional textures, structural deformation, pore fluid chemistry, and ultimately influences landscape evolution by introducing heterogeneities in erodibility. In Southern Utah, the region West of the Kaibab uplift known as Mollies Nipple (Mollies) in Grand Staircase‐Escalante National Monument exhibits a complex history of fluid–sediment interactions, which has resulted in a localized zone of anomalous diagenetic iron sulfate (jarosite) mineralogy in a well‐exposed dune–interdune deposit within the Navajo Sandstone. Mineralogy and geochemistry of cements within this region are examined using reflectance and imaging spectroscopy, field investigations, microscopy, and whole‐rock geochemical analyses. These data show that the in‐situ jarosite cement is localized to a plane along the highest ridge of the butte, providing an armor along with other secondary cements, which controls the butte's geomorphic evolution. The jarosite cement is associated with other mineralogies suggesting that the sulfate was one of the latest fluid‐related precipitates in the paragenetic sequence. It was preceded by a regional bleaching event, precipitation of clay cements, some localized concretionary iron oxide precipitation, and formation of deformation bands. At least one generation of dense iron oxide mineralization is associated with cataclastic brittle deformation predating the sulfate precipitation. Trace element geochemistry of cements shows certain metal oxide populations associated with extremely high (>2000 ppm) arsenic values. We interpret the combination of spatial mineral distribution, observed paragenetic sequence, and trace element geochemistry to suggest this region experienced acid sulfate diagenesis along fracture‐controlled fluid conduits related to weathering of proximal, unidentified, sulfides, or H₂S associated with deep source beds. Jarosite is highly soluble, and its presence suggests that abundant fluid flow has not occurred in this region since its formation. This terminal cement‐forming event must have occurred prior to sandstone exhumation and erosion to form the current extreme landscape at Mollies. This site exhibits the influence that fluid geochemistry, sedimentary mineralogy, and structural fabric have on geomorphic evolution.     

Culture and Legal Policy Punctuation in the Supreme Court's Gender Discrimination Cases

For the most part, punctuated equilibrium scholarship has ignored the legal policy change generated by the Supreme Court. In this study, I address this gap though an examination of the Court's equal protection and gender cases from the 1970s. My case study here has two aims. First, I offer an adaptation of the jurisprudential regimes framework as a device for framing and identifying legal policy punctuations. After identifying Reed v. Reed (1971) as the cut point of such a regime, I then use Reed and its progeny to illustrate the promise of culture in explaining stasis and change, specifically focusing on the concepts of cultural cognition and cultural surprise.     

Geometry‐coupled reactive fluid transport at the fracture scale: application to CO2 geologic storage

Water acidification follows CO₂ injection and leads to reactive fluid transport through pores and rock fractures, with potential implications to reservoirs and wells in CO₂ geologic storage and enhanced oil recovery. Kinetic rate laws for dissolution reactions in calcite and anorthite are combined with the Navier‐Stokes law and advection–diffusion transport to perform geometry‐coupled numerical simulations in order to study the evolution of chemical reactions, species concentration, and fracture morphology. Results are summarized as a function of two dimensionless parameters: the Damköhler number Da which is the ratio between advection and reaction times, and the transverse Peclet number Pe defined as the ratio between the time for diffusion across the fracture and the time for advection along the fracture. Reactant species are readily consumed near the inlet in a carbonate reservoir when the flow velocity is low (low transverse Peclet number and Da > 10^?1). At high flow velocities, diffusion fails to homogenize the concentration field across the fracture (high transverse Peclet number Pe > 10^?1). When the reaction rate is low as in anorthite reservoirs (Da < 10^?1), reactant species are more readily transported toward the outlet. At a given Peclet number, a lower Damköhler number causes the flow channel to experience a more uniform aperture enlargement along the length of the fracture. When the length‐to‐aperture ratio is sufficiently large, say l/d > 30, the system response resembles the solution for 1D reactive fluid transport. A decreased length‐to‐aperture ratio slows the diffusive transport of reactant species to the mineral fracture surface, and analyses of fracture networks must take into consideration both the length and slenderness of individual fractures in addition to Pe and Da numbers.     

Reconstructing Diet of the Early Qin (ca. 700–400 BC) at Xishan,Gansu Province,China

Here we report the bone collagen carbon and nitrogen isotopic results of humans (n = 33) and animals (n = 58) to reconstruct the dietary practices of an early Qin population dating to the Zhou Dynasty (Late Western–Early Eastern period ca. 700–400 BC ) at the Xishan site in Gansu Province, China. The humans have a very large range of δ¹³C (−23.3‰ to −7.1‰) and δ¹⁵N (4.3‰ to 10.9‰) values which reflects extraordinarily diverse diets and included individuals with predominately C₃ as well as those with exclusive C₄ diets. This wide span of isotopic results produced a subtle linear trend (R² = 0.62) in the human data, which paralleled the animals across the C₃ and C₄ environmental gradient. However, the majority of the individuals had a predominately C₄ diet based on millet with δ¹⁵N results only slightly elevated above the animals, except for the pigs and cattle. This is evidence that many of the animals were likely used for their secondary products, labour or as sacrificial offerings and that pork and beef were the main sources of animal protein for the population. High status individuals had elevated δ¹⁵N values (10.2 ± 0.6‰) compared to medium (8.9 ± 0.3‰) and lower status (8.8 ± 0.8‰) individuals, possibly related to increased animal protein in the diet. Differences related to gender were also found with females having elevated δ¹³C (−11.2 ± 1.9‰) and δ¹⁵N (9.4 ± 0.8‰) values compared to the males (δ¹³C = −14.1 ± 4.2‰; δ¹⁵N = 7.9 ± 1.9‰), but these results necessitate caution given the large number of individuals that could not be sexed. The results of this study support the view that the early Qin people were a more sedentary society focused on millet agriculture and animal husbandry, and that they were influenced by the pre‐existing populations of the central Gansu region. Copyright © 2015 John Wiley & Sons, Ltd.