REAPPRAISAL OF MAGNITUDE OF 20TH CENTURY EARTHQUAKES IN SWITZERLAND

We reappraise the material needed to assess the 20th-century seismicity of Switzerland and of the adjacent areas of the Alps in terms of magnitude. For this we make use of macroseismic reports and literature and by calculating the surface-wave magnitude from the Prague formula of all significant earthquakes in the region. No attempt is made to relocate earthquake positions; instead their reliability is ranked using existing solutions and macroseismic observations We find that for small earthquakes (M _s >4.5), which constitute the bulk of the events in the region, the calculation of M _s observed at relatively short distances, requires station and distance corrections which can be significant. Also we find that recomputed Ms estimates differ from those reported in other earthquake catalogues. From this reappraisal of M _s and from a uniform re-evaluation of the associated macroseismic data, we derive a stable correlations between M _s and felt areas that can be used to assess the magnitude of historical, pre-instrumental, events for which only isoseismal radii (r _i ) and the associated intensities (I _i ) are available. We examined the conversion of surface-wave magnitude into moment magnitude, a conversion that presents some interesting problems for relatively small events for which the M _s -log(Mo) scaling changes. We conclude that the rate of moment release derived from events of M _s >4.0 is small, and that it should be associated with horizontal and vertical velocity rates of less than 1 mm/yr, too small vo be confirmed by GPS measurements over short periods of time so that can be used to constrain hazard assessment.     

Repairs to an injured early Middle Cambrian Trilobite,Elkedra area,Northern Territory

Wang, Y., Shih, C., Szwedo, J. &; Ren, D. iFirst article. New fossil palaeontinids (Hemiptera, Cicadomorpha, Palaeontinidae) from the Middle Jurassic of Daohugou, China. Alcheringa, 1–12. ISSN 0311-5518.

A new genus and species assigned to the extinct family Palaeontinidae, Synapocossus sciacchitanoae Wang, Shih &; Ren, is described from the Middle Jurassic of Daohugou in Inner Mongolia, China. This new genus is established based on well-preserved fossil specimens with body and complete forewings and hind wings. It differs from other described genera by the following characters: small body size, RP and M₁ coalescing for an interval on the forewings and M_3 + 4 without bifurcation on the hind wings. The RP coalescence with M₁ in Synapocossus Wang, Shih &; Ren previously reported only in Turgaiella Becker-Migdisova &; Wootton, seems to be associated with strengthening of the anterior wing margin. The intra-specific and individual variations of Synapocossus and numerous other insect fossils of northeastern China probably indicate long-lasting ecological stresses and a competitive environment in the Middle Jurassic to Early Cretaceous ecosystems.     

A Note on the Extremities of Local Moran's Iis and Their Impact on Global Moran's I

By defining local Moran's I_i as a ratio of quadratic forms and making use of its overall additivity to match global Moran's I, we can identify spatial objects with a strong impact on global Moran's I. First, we concentrate on the spatial properties of local Moran's I_i expressed by the local linkage degree. Depending on whether we use the W- or C-coding of the spatial connectivity matrix, the variance of local Moran's I_i for a small local linkage degree will be either large or small. Note that spatial objects associated with a local Moran's I_i with a large variance affect the global statistic much more than spatial objects associated with a local Moran's I_i with a small variance. Counterintuitively, global Moran's I defined in the W-coding is most influenced by spatial objects with a small number of spatial neighbors. In contrast, spatial objects with a large number of spatial neighbors exert more impact on global Moran's I setup in the C-coding. Second, we investigate the impact of the empirical data on local Moran's I_i and show that local Moran's I_i will only be significant for extreme absolute residuals at and around the reference location. Clusters of average regression residuals cannot be detected by local Moran's I_i. Consequently, spatial cliques of extreme residuals contribute more to significance tests on global autocorrelation.     

Local Spatial Autocorrelation in a Biological Model

We review the recently developed local spatial autocorrelation statistics I_i, c_i, G_i, and G_i*. We discuss two alternative randomization assumptions, total and conditional, and then newly derive expectations and variances under conditional randomization for I_i and c_i, as well as under total randomization for c_i. The four statistics are tested by a biological simulation model from population genetics in which a population lives on a 21 × 21 lattice of stepping stones (sixty-four individuals per stone) and reproduces and disperses over a number of generations. Some designs model global spatial autocorrelation, others spatially random surfaces. We find that spatially random designs give reliable test results by permutational methods of testing significance. Globally autocorrelated designs do not fit expectations by any of the three tests we employed. Asymptotic methods of testing significance failed consistently, regardless of design. Because most biological data sets are autocorrelated, significance testing for local spatial autocorrelation is problematic. However, the statistics are informative when employed in an exploratory manner. We found that hotspots (positive local autocorrelation) and coldspots (negative local autocorrelation) are successfully distinguished in spatially autocorrelated, biologically plausible data sets.     

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.     

Measurements of gradients in ionospheric parameters with a new nine-position experiment at Millstone Hill

A new nine-position experiment is now routinely carried out with the Millstone Hill incoherent scatter radars which allows estimation of spatial gradients in the measured ionospheric scalar parameters N_e, T_e, and T_i, and in the components of the ion velocity vector v_i. Use of this technique results in improved estimates of basic and derived parameters from incoherent scatter data at times of significant gradients. We detail the data analysis method and present the first results from this new experiment. The gradients in N_e and in the components of vi are used to compute the motion term in the ionospheric F region continuity equation ▿ · (Nv), which is then combined with ∂N/∂t to estimate the O⁺ recombination rate β at night. Meridional neutral winds U_mer are computed from the field-aligned ion velocity v_∥ and a calculation of the O⁺ diffusion velocity v_d, and it is found that horizontal gradients in the ion velocity field at times significantly affect the calculation of the neutral winds.     

The Analysis of Spatial Association by Use of Distance Statistics

Introduced in this paper is a family of statistics, G, that can be used as a measure of spatial association in a number of circumstances. The basic statistic is derived, its properties are identified, and its advantages explained. Several of the G statistics make it possible to evaluate the spatial association of a variable within a specified distance of a single point. A comparison is made between a general G statistic and Moran's I for similar hypothetical and empirical conditions. The empirical work includes studies of sudden infant death syndrome by county in North Carolina and dwelling unit prices in metropolitan San Diego by zip-code districts. Results indicate that G statistics should be used in conjunction with I in order to identify characteristics of patterns not revealed by the I statistic alone and, specifically, the G_i and G_i* statistics enable us to detect local “pockets” of dependence that may not show up when using global statistics.     

Local Spatial Autocorrelation Statistics: Distributional Issues and an Application

The statistics G_i(d) and G_i*(d), introduced in Getis and Ord (1992) for the study of local pattern in spatial data, are extended and their properties further explored. In particular, nonbinary weights are allowed and the statistics are related to Moran's autocorrelation statistic, I. The correlations between nearby values of the statistics are derived and verified by simulation. A Bonferroni criterion is used to approximate significance levels when testing extreme values from the set of statistics. An example of the use of the statistics is given using spatial-temporal data on the AIDS epidemic centering on San Francisco. Results indicate that in recent years the disease is intensifying in the counties surrounding the city.     

A new genus and species of dasyurid from the Pliocene Chinchilla Local Fauna of south-eastern Queensland

A new genus and species of dasyurid is described on the basis of P³, M¹–M³ and M_2–3 collected from the middle Pliocene Chinchilla Local Fauna of south-eastern Queensland. Character analysis supports the placement of this new taxon within the dasyurid subfamily Dasyurinae. Derived dasyuromorphian features common to A. chinchillaensis and most dasyurines include: reduction of P³, reduction of stylar cusp B on M¹; reduction of metaconids; approximation of stylar cusp B and the metacone on M¹; alignment of the postmetacrista on M¹ with the long axis of the maxilla and an oblique orientation of the hypocristids with respect to the long axis of the teeth. However in most instances, A. chinchillaensis is less specialised than other members of the subfamily for characters considered. One derived dental feature, reduction of the entoconid, may indicate a special relationship with ‘parantechinines’, Dasyuroides and/or Dasycercus. Additionally, the presence of this feature is suggestive of an open woodland/seasonal habitat. Autapomorphic P³ morphology compounds difficulties with the phylogenetic placement of this taxon.     

Comparison of incoherent scatter observations of electron density,and electron and ion temperature at Millstone Hill with the International Reference Ionosphere

Millstone Hill incoherent scatter (IS) observations of electron density (N_e, electron temperature (T_e) and ion temperature (T_i) are compared with the International Reference Ionosphere (IRI-86) for both noon and midnight, for summer, equinox and winter, at both solar maximum (1979–1980) and solar minimum (1985–1986). The largest difference inN_e is found in the topside, where values of N_e given by IRI-86 are generally larger than those obtained from IS measurements, by a factor which increases with increasing height, and which has a mean value near two at 600 km. Apart from the bottom of the profile, which is tied to the CIRA neutral temperature, the IRI-86 T_e model has no solar cycle variation. However, the IS measurements during the summer reveal larger T_e at solar maximum than at solar minimum. At other seasons higher T_e at solar maximum occurs only during the daytime at the greater heights. Nighttime T_e is shown by the IS radar to be generally larger in winter than in summer, an effect not included in the IRI. This is apparently due to photoelectron heating during winter from the sunlit ionosphere conjugate to Millstone Hill. The day-night difference in T_i given by IRI-86 above 600km is not as large in the IS measurements.     

Gas breakthrough experiments on pelitic rocks: comparative study with N2, CO2 and CH4

The capillary‐sealing efficiency of intermediate‐ to low‐permeable sedimentary rocks has been investigated by N₂, CO₂ and CH₄ breakthrough experiments on initially fully water‐saturated rocks of different lithological compositions. Differential gas pressures up to 20 MPa were imposed across samples of 10–20 mm thickness, and the decline of the differential pressures was monitored over time. Absolute (single‐phase) permeability coefficients (k_abs), determined by steady‐state fluid flow tests, ranged between 10^?22 and 10^?15 m². Maximum effective permeabilities to the gas phase k_eff(max), measured after gas breakthrough at maximum gas saturation, extended from 10^?26 to 10^?18 m². Because of re‐imbibition of water into the interconnected gas‐conducting pore system, the effective permeability to the gas phase decreases with decreasing differential (capillary) pressure. At the end of the breakthrough experiments, a residual pressure difference persists, indicating the shut‐off of the gas‐conducting pore system. These pressures, referred to as the ‘minimum capillary displacement pressures’ (P_d), ranged from 0.1 up to 6.7 MPa. Correlations were established between (i) absolute and effective permeability coefficients and (ii) effective or absolute permeability and capillary displacement pressure. Results indicate systematic differences in gas breakthrough behaviour of N₂, CO₂ and CH₄, reflecting differences in wettability and interfacial tension. Additionally, a simple dynamic model for gas leakage through a capillary seal is presented, taking into account the variation of effective permeability as a function of buoyancy pressure exerted by a gas column underneath the seal.     

FIRING ANALYSIS OF SOUTH-EASTERN MISSOURI INDIAN POTTERY USING IRON MÖSSBAUER SPECTROSCOPY*

The Mössbauer spectra of prehistoric pottery and clay from the Malden Plain of south-eastern Missouri, USA, were examined. The pottery dates from AD 500 to 1400. The earlier sherds were tempered with sand and the later ones with shell. The sherds and clay contain a mixture of illite, smectite, and kaolinite. Most samples of clay and sherds contained both ferric and ferrous species. The spectral parameters for the ferric species were %DLE_Q～0. 7—1.1 mms^-1 and δ～0.0–0.3 mms^-1 for the ferrous species, ΔE_Q～2.1–2.7mms^-1 and δ～0.7–1.0mms^-1. The ferric to ferrous ratio is more accurately determined from the spectra than from consideration of the colour or the firing core of the sherds, which are not always related to the iron species ratio. Differences among sherds in the ferrous ΔE_Q can be related to original temperature of heating but with low precision.     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

Transient fluid flow in and around a fault

We present the results of simple numerical experiments in which we study the evolution with time of fluid flow around and within a permeable fault embedded in a less permeable porous medium. Fluid movement is driven by an imposed vertical pressure gradient. The results show that fluid flow is controlled by two timescales: τ_f = Sl²/κ_F and τ_F = Sl²/κ_M, where S is the specific storage of the porous material, l the length of the fault, and κ_M and κ_F are the hydraulic conductivities of the porous material and the fault, respectively. Fluid flow and the associated fluid pressure field evolve through three temporal stages: an early phase [t < τ_f] during which the initial fluid pressure gradient within the fault is relaxed; a second transient stage [τ_f < t < τ_F] when fluid is rapidly expelled at one end of the fault and extracted from the surrounding rocks at the other end leading to a reduction in the pressure gradient in the intact rock; a third phase [t < τ_F] characterized by a steady‐state flow. From the numerical experiments we derive an expression for the steady‐state maximum fluid velocity in the fault and the values of the two timescales, τ_f and τ_F. A comparison indicates excellent agreement of our results with existing asymptotic solutions. For km‐scale faults, the model results suggest that steady‐state is unlikely to be reached over geological timescales. Thus, the current use of parameters such as the focusing ratio defined under the assumption of steady‐state conditions should be reconsidered.     

Studies of high latitude mesospheric turbulence by radar and rocket 1: energy deposition and wave structure

During early spring, 1985, the MAE-3 (Middle Atmospheric Electrodynamics) Program was conducted at Poker Flat Research Range, Alaska to study the origin of wintertime mesospheric echoes observed with the Poker Flat MST radar there, by probing the mesosphere with in situ rocket measurements when such echoes occurred. Pre-launch criteria required the appearance of echoes exhibiting some wave structure on the MST radar display; these could be met even under weak precipitation conditions with riometer absorption near or above 1.0 dB. Two morning rockets were launched under such conditions, the first (31.048) on 29 March 1985, at 1703 UT and the second (31.047) on 1 April 1985, at 1657 UT. Both payloads were deployed on a high altitude parachute near a 95 km apogee to provide a stable platform for data acquisition within the mesosphere (below 80 km). Each payload carried a solid state detector to measure energetic electrons between 0.1 and 1.0 MeV and an NaI crystal detector to measure x-rays from >5 to >80 keV. Payload 31.048 also carried a positive ion ‘turbulence’ probe which measured ion density changes (ΔN_i/N_i) during payload descent, whereas 31.047 carried a nose tip ‘turbulence’ probe designed to measure electron density changes (ΔN_e/N_e) during upleg ram conditions plus a Gerdien condenser for the measurement of bulk ion properties during downleg. The energy deposition curves for each event exhibited peak deposition rates between 75 and 80 km with a half width of 16–18 km, almost exclusively induced by precipitating relativistic electrons. They also showed a maximum bottomside gradient between 65 and 75 km. Radar echoes and atmospheric turbulence were observed in the same altitude domain, consistent with the anticipated need for adequate free thermal electron gradients to make such phenomena visible on the radar. The vertical wave structure from radar echoes was found to be consistent with that observed in horizontal wind and temperature profiles measured by Datasondes flown shortly after each large rocket. An analysis of the wave structure from radar data has shown that although large scale waves (λ_z ~ 7 km) were found to be present, a higher frequency shorter wavelength (∼ 1–3 km) component probably played a more significant role in modulating the signal-to-noise structure of the radar echoes.     

Effects of the ionospheric conductivity on the longitudinal structure of the F-region equatorial anomaly

A longitudinal comparison has been undertaken of the vertical ion velocities in terms of ion mobilities and the eastward electrojet currents in terms of the electric conductivities overhead the magnetic equator at African and West Asian regions to search for longitudinal differences in the developments of the equatorial anamoly. It is shown that the vertical ion velocities and electrojet currents in both regions are different. It is proposed that longitudinal differences in (k²_1i+k²_2i)/k_2i and K_1i/K_2i ratios give rise to a different vertical ion velocities in these regions. This is likely to cause different equatorial F-layer plasma fountain magnitudes resulting in different developments of the equatorial anomaly.     

Invasion of a karst aquifer by hydrothermal fluids: evidence from stable isotopic compositions of cave mineralization

Mineral deposits in the Cupp‐Coutunn/Promeszutochnaya cave system (Turkmenia, central Asia) record a phase of hydrothermal activity within a pre‐existing karstic groundwater conduit system. Hydrothermal fluids entered the caves through fault zones and deposited sulphate, sulphide and carbonate minerals under phreatic conditions. Locally, intense alteration of limestone wall rocks also occurred at this stage. Elsewhere in the region, similar faults contain economic quantities of galena and elemental sulphur mineralization. Comparisons between the Pb and S isotope compositions of minerals found in cave and ore deposits confirm the link between economic mineralization and hydrothermal activity at Cupp‐Coutunn. The predominance of sulphate mineralization in Cupp‐Coutunn implies that the fluids were more oxidized in the higher permeability zone associated with the karst aquifer. A slight increase in the δ³⁴S of sulphate minerals and a corresponding δ³⁴S decrease in sulphides suggest that partial isotopic equilibration occurred during oxidation. Carbonate minerals indicate that the hydrothermal fluid was enriched in ¹⁸O (δ¹⁸O_SMOW ～ + 10‰) relative to meteoric groundwater and seawater. Estimated values for δ¹³C_DIC (δ¹³C_PDB ～ ? 13‰) are consistent with compositions expected for dissolved inorganic carbon (DIC) derived from the products of thermal decomposition of organic matter and dissolution of marine carbonate. Values derived for δ¹³C_DIC and δ¹⁸O_water indicate that the hydrothermal fluid was of basinal brine origin, generated by extensive water–rock interaction. Following the hydrothermal phase, speleothemic minerals were precipitated under vadose conditions. Speleothemic sulphates show a bimodal sulphur isotope distribution. One group has compositions similar to the hydrothermal sulphates, whilst the second group is characterized by higher δ³⁴S values. This latter group may either record the effects of microbial sulphate reduction, or reflect the introduction of sulphate‐rich groundwater generated by the dissolution of overlying evaporites. Oxygen isotope compositions show that calcite speleothems were precipitated from nonthermal groundwater of meteoric origin. Carbonate speleothems are relatively enriched in ¹³C compared to most cave deposits, but can be explained by normal speleothem‐forming processes under thin, arid‐zone soils dominated by C4 vegetation. However, the presence of sulphate speleothems, with isotopic compositions indicative of the oxidation of hydrothermal sulphide, implies that CO₂ derived by reaction of limestone with sulphuric acid (‘condensation corrosion’) contributed to the formation of ¹³C‐enriched speleothem deposits.     

On the geochemistry of gases and noble gas isotopes (including 222Rn) in deep crustal fluids: the 4000 m KTB-pilot hole fluid production test 2002–03

The concentrations of H₂, O₂, CO₂, and concentrations and isotopic composition of the noble gases (including ²²²Rn), N₂, CH₄, and higher hydrocarbons dissolved in 4000 m deep‐seated fluids from a 12‐month fluid production test in the KTB pilot hole were analyzed. This determination of the gas geochemistry during the test in combination with the knowledge of the hydraulic data provides relevant information about the fluid hydraulics of the deep system. All gas concentrations and isotopic signatures, except for ²²²Rn, showed constancy during the course of the test. This, in combination with large fluid flow rates at a moderate water table drawdown, imply an almost infinite fluid reservoir in 4000 m depth. From the change in ²²²Rn‐activity as a function of pump rate, the contribution of smaller and wider pores to the overall fluid flow in an aquifer can be deduced. This ²²²Rn‐activity monitoring proved therefore to be a valuable instrument for the qualitative observation of the scavenging of pore and fracture surfaces, a hydraulic feature invisible to standard hydraulic testing tools. The observance of this scavenging effect is due to (i) the continuous on‐line geochemical monitoring, (ii) the durability of the test, (iii) a change in pump rate during the course of the test, and (iv) due to the short half‐life of ²²²Rn. The fluids have a 5.9% mantle He component, and a δ²¹Ne excess of 14%, and a noble gas model age of about (5.5–6.2) ± 2.0 Myr. The mean N₂/Ar‐ratio of 516 and δ¹⁵N‐data of about +1.5‰ indicates sedimentary or metamorphic origin of N₂. The hydrocarbons, amounting to 33 vol.% in the gas phase, are derived from thermal decomposition of marine organic matter of low maturity. But a key question, the identification of the potential source region of the fluids and the migration pathway, is still unidentified.     

MODELLING OF THE GROUND MOTION AT RUSSE SITE (NE BULGARIA) DUE TO THE VRANCEA EARTHQUAKES

An approach, capable of synthesising strong ground motion from a basic understanding of fault mechanism and of seismic wave propagation in the Earth, is applied to model the seismic input at a set of 25 sites along a chosen profile at Russe, NE Bulgaria, due to two intermediate-depth Vrancea events (August 30, 1986, M _ω=7.2, and May 30, 1990, M _ω=6.9). Accordingly to our results, once a strong ground motion parameter has been selected to characterise the ground motion, it is necessary to investigate the relationships between its values and the features of the earthquake source, the path to the site and the nature of the site. Therefore, a proper seismic hazard assessment requires an appro-priate parametric study to define the different ground shaking scenarios corresponding to the relevant seismogenic zones affecting the given site. Site response assessment is provided simultaneously in frequency and space domains, and thus the applied procedure differs from the traditional engineering approach that discusses the site as a single point. The applied procedure can be efficiently used to estimate the ground motion for different purposes like microzonation, urban planning, retrofitting or insurance of the built environment.     

Increasing returns to scale in the towns of early Tudor England

Abstract

Urban agglomeration economies make cities central to theories of modern economic growth. There is historical evidence for the presence of Smithian growth and agglomeration effects in English towns c.1450-1670, but seminal assessments deny the presence of agglomeration effects and productivity gains to Early Modern English towns. This study evaluates the presence of increasing returns to scale (IRS) in aggregate urban economic outputs—the empirical signature of feedbacks between Smithian growth and agglomeration effects—among the towns of 16th century England. To do so, we test a model from settlement scaling theory against the 1524/5 Lay Subsidy returns. Analysis of these data indicates that Tudor towns exhibited IRS—a finding that is robust to alternative interpretations of the data. IRS holds even for the smallest towns in our sample, suggesting the absence of town size thresholds for the emergence of agglomeration effects. Spatial patterning of scaling residuals further suggests regional demand-side interactions with Smithian-agglomeration feedbacks. These findings suggest the presence of agglomeration effects and Smithian growth in pre-industrial English towns. This begs us to reconsider the economic performance of Early Modern English towns, and suggests that the qualitative economic dynamics of contemporary cities may be applicable to premodern settlements in general.