Geomorphic evolution of the Southern Alps,New Zealand

The Southern Alps are the topographic expression of late Cenozoic (<8 Ma ago) uplift of the crust of the leading edge of the Pacific plate in South Island, New Zealand. New fission track data on the basement exposed in the Southern Alps quantify the age, amount, and rate of rock uplift, and in combination with geomorphic parameters permit the construction of a new model of the geomorphic evolution of the Southern Alps. The model emphasizes the development over time and space of rock uplift, mean surface elevation, exhumation of crustal section, and relief. The earliest indications of mean surface uplift are between 4 and 5 Ma ago at the Alpine Fault. Mean surface uplift, which lagged the start of rock uplift, propagated southeastward from the Alpine Fault at a rate of 30 km/Ma. By about 4 Ma ago, exhumation had exposed greywacke basement adjacent to and east of the entire 300 km long central section of the Alpine Fault. At 3 Ma ago, greenschist was exposed in the southern parts of the Southern Alps near Lake Wanaka, and since then has become exhumed along a narrow strip east of the Alpine Fault. The model infers that amphibolite grade schist has been exhumed adjacent to the Alpine Fault only in the last 0·3 Ma. The age of the start of rock uplift and the amount and rate of rock uplift, all of which vary spatially, are considered to be the dominant influences on the development of the landscape in the Southern Alps. The Southern Alps have been studied in terms of domains of different rock uplift rate. At present the rate of rock uplift varies from up to 8–10 mm/a adjacent to the Alpine Fault to 0·8–1·0 mm/a along the southeastern margin of the Southern Alps. This spectrum can be divided into two domains, one northwest of the Main Divide where the present rock uplift rates are very high (up to 8–10 mm/a) and exceed the long-term value of 0·8–1·0 mm/a, and another to the southeast of the Main Divide where the long-term rate is 0·8–1·0 mm/a. A domain of no uplift lies immediately to the east of the Southern Alps, and is separated from them by a 1·0–1·5 km step in the basement topography. We argue that this spatial sequence of uplift rate domains represents a temporal one. The existing models of the geomorphic development of the Southern Alps—the dynamic cuesta model of J. Adams and the numerical model of P. Koons—are compared with the new data and evolutionary model. Particular constraints unrealized by these two earlier models include the following: the earlier timing of the start of rock uplift of the Southern Alps (8 Ma ago); the spatial variation in the timing of the start of rock uplift (8 Ma ago to 3 Ma ago); the lower long-term rock uplift rate (0·8–1·0 mm/a) of the Southern Alps for most of the late Cenozoic; the lag between the start of rock uplift and the start of mean surface uplift; and the patterns of the amounts of late Cenozoic rock uplift and erosion across the Southern Alps.     

Revised palaeoclimatic significance of Mueller Glacier moraines,Southern Alps,New Zealand

A sequence of Late Holocene moraines on the foreland of the Mueller Glacier, Southern Alps, New Zealand, forms part of a local moraine‐age database used to establish a regional glacier chronology and subsequently to investigate potential intra‐hemispheric and global climate forcing mechanisms. We present new sedimentological and geomorphological evidence that a set of these moraine ridges, previously considered to represent individual advances, constitutes a single moraine complex (the ‘Mueller Memorial Moraine’) formed by supraglacial transport of a large volume of landslide debris to the glacier terminus. Because a moraine formed in this way is not necessarily associated with an advance triggered by a climate event, we question the palaeoclimatic significance of the Mueller Memorial Moraine, as well as that of other moraines in comparable settings. Our findings suggest that the mode of formation and glacio‐dynamical context of moraines whose ages contribute to existing palaeoclimate reconstructions need to be re‐examined in order to assess the reliability of these reconstructions. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantitative relationships between uplift and relief parameters for the Southern Alps,New Zealand,as determined by fission track analysis

The Southern Alps mountain chain, New Zealand, has formed as a consequence of late Cenozoic collision of the continental parts of the Pacific and Australia plates. Fission track analysis has yielded estimates of the amount, age of initiation, and rate of late Cenozoic rock uplift for 82 surface samples taken from transects across the Southern Alps. The mean surface, summit and valley elevations in the vicinity of each of the rock sample sites have also been measured. Regression of the geomorphic variables on the uplift variables has been used to establish quantitative relationships between uplift and geomorphology. There are strong and consistent linear associations between uplift and the elevations of the mean surface, summits and valleys. The preferred regression models have uniform slope but varying elevation response between transects. Substitution of space for time has allowed the evolution of landforms to be studied. To the east of the Main Divide, elevation and relief are proportional to, and closely related to, the age of initiation of rock uplift (‘uplift age’) and the amount of rock uplift (r² > 0·8). Mean surface uplift was delayed for ～2 Ma after the start of rock uplift, a result of the stripping of a soft cover rock succession that, prior to rock uplift, overlaid the harder greywacke basement. Inter-transect variations in regression response and x-intercept are inferred, therefore, to reflect the variable preuplift thickness of cover rocks. However, the regular regression slope for the transects reflects the consistent nature of the interaction between uplift and the erodibility of greywacke basement. Uplift of the mean surface proceeded at 0·4 km/km and 0·4 km/Ma of rock uplift, while the rock uplift rate was 0·8 km/Ma. Summit elevations have increased at a rate of 0·6 km/Ma and valley elevations have increased at 0·2 km/Ma. Regression lines relating mean surface, summit and valley elevations to rock uplift and uplift age diverge from common intercepts; it is concluded, therefore, that the mountains east of the Main Divide have continued to increase in elevation and relief and change in form over time since the start of mean surface uplift. Mountain elevation has little relationship with late Cenozoic mean rock uplift rates of 0·8–1·0 km/Ma or inferred contemporary rock uplift rates (r² ～ 0·3). In contrast, to the west of the Main Divide, elevation is shown to be closely related to rock uplift rate (r² > 0·3). In contrast, to the west of the Main Divide, elevation is shown to be closely related to rock uplift rate (r² > 0·8). Transects with higher rock uplift rates support higher topography. Landforms are therefore in a stable equilibrium with rock uplift rate, and the landscape contains no residual evidence of the total amount of rock uplift, or the age of uplift. Lithological variation appears to have no relationship with elevation.     

Thermal modeling of the Southern Alps,New Zealand

Finite-element modeling of the thermal regime across the Southern Alps of New Zealand has been carried out along two profiles situated near the Franz Josef and Haast valleys. The modeling involves viscous deformation beneath the Southern Alps, including both uplift and erosion, and crustal/lithospheric thickening, as a result of crustal shortening extending to 20 mm/y of a 25-km thick crust. Published uplift rates and crustal thickness variations along the two profiles are used to constrain the modeled advection of crustal material, and results are compared with the recent heat flow determinations, 190±50 mW/m² in the Franz Josef valley and 90±25 mW/m² in the Haast valley. Comparisons of the model with published K–Ar and fission track ages, show that the observed heat flow in the Franz Josef valley is consistent with observed zircon fission track ages of around 1 Ma, if the present-day uplift rate is close to 10 mm/y. Major thermal differences between the Franz Josef and Haast profiles appear to be due to different uplift and erosion rates. There is weak evidence that frictional heating close to the Alpine fault zone is not significant. The modeling provides explanations for the distribution of seismicity beneath the Southern Alps, and predicts a low surface heat flow over the eastern foothills due to the dominant thermal effect of crustal thickening beneath this region. Predicted temperatures at mid-crustal depth beneath the zone of maximum uplift rate are 50–100°C cooler than those indicated in previously published models, which implies that thermal weakening of the crust may not be the main factor causing the aseismicity of the central Southern Alps. The results of the modeling demonstrate that the different types of reset age data in the region within 25 km of the Alpine fault are critical for constraining models of the deformation and the thermal regime beneath the Southern Alps.     

Controls on slope‐to‐channel fine sediment connectivity in a largely ice‐free valley,Hoophorn Stream,Southern Alps,New Zealand

There has been little work to date into the controls on slope‐to‐channel fine sediment connectivity in alpine environments largely ice‐free for most of the Holocene. Characterization of these controls can be expected to result in better understanding of how landscapes “relax” from such perturbations as climate shock. We monitored fine sediment mobilization on a slope segment hydrologically connected to a stream in the largely ice‐free 8·3 km² Hoophorn Valley, New Zealand. Gerlach traps were installed in ephemeral slope channels to trap surficial material mobilized during rainfall events. Channel sediment flux was measured using turbidimeters above and below the connected slope, and hysteresis patterns in discharge‐suspended sediment concentrations were used to determine sediment sources. Over the 96 day measurement period, sediment mobilization from the slope segment was limited to rainfall events, with increasingly larger particles trapped as event magnitude increased. Less than 1% of the mass of particles collected during these events was fine sediment. During this period, 714 t of suspended sediment was transported through the lower gauging station, 60% of it during rainfall events. Channel sediment transfer patterns during these events were dominated by clockwise hysteresis, interpreted as remobilization of nearby in‐channel sources, further suggesting limited input of fine sediment from slopes in the lower valley. Strong counterclockwise hysteresis, representing input of fine sediment from slope segments, was restricted to the largest storm event (JD2 2009) when surfaces in the upper basin were activated. The results indicate that the slopes of the lower Hoophorn catchment are no longer functioning as sources of fine sediment, but rather as sources of coarse material, with flux rates controlled by the intensity and duration of rainfall events. Although speculative, these findings suggest a shift to a coarse sediment dominated slope‐to‐channel transfer system as the influence of pre‐Holocene glacial erosion declines. Copyright © 2011 John Wiley & Sons, Ltd.     

The lunar barometric tide in New Zealand

The lunar barometric tide has been determined with reasonable accuracy, on an annual and seasonal basis, at five stations on the mainland of New Zealand and at three of the outlying islands. The determinations show that in the New Zealand region the lunar tide has a larger amplitude and smaller phase than might have been expected from previously available southern-hemisphere results. In general, the seasonal variation of phase in the New Zealand region conforms to the currently recognised global pattern, with the J-season phase greater than that of the D-season. Similarly, the amplitude variation tends to support the suggestion that, south of latitude 30°S, the D-season amplitude is greater than that of the J-season. Approximate tests are introduced and used to assess the statistical significance of some of the apparent differences in amplitude and phase made evident by the analysis. These assessments indicate that although many of the apparent differences may be attributed to sampling fluctuations, the main regional anomalies in amplitude and phase are likely to be real. It is suggested that these anomalies may indicate a significant regional input of tidal energy to the atmosphere (at the lower boundary) from the Pacific oceanic tide.     

Solution of marble in the karst of the Pikikiruna range,Northwest Nelson,New Zealand

The paper reports the first research on karst solution processes and rates in New Zealand. The study area is an IHD representative basin in the northwest corner of the South Island in a mountain range consisting principally of Ordovician marble. The climate is sunny and warm (17°C) in summer and wet and cool (7°C) in winter. Average precipitation is 2,158 mm of which 525 mm evapotranspires, yielding a discharge of 51–75 l/s/km² in the river basin studied. Almost half of the catchment of 45.1 km² consists of karst which occurs mainly as a doline covered plateau at 600–900 m within which most drainage is subterranean. Water tracing is with fluorescein defined drainage patterns. Marble solution was established by estimating inputs, throughputs and outputs of water and dissolved calcium and magnesium in both autogenic and allogenic karst drainage systems. Particular attention was paid to estimating errors. Water samples for chemical analysis were taken irregularly for approximately one year, and a rating curve relating chemical load to discharge was established. The best estimate of solution loss from the basin yields a mean rate of 100 ±M24 m³/km²/a. Of this 80 per cent is derived from solution of marble by autogenic waters, mostly in the top 10–30 m of the marble outcrop. The remaining 20 per cent is accomplished by allogenic stream solution. Approximately 9.9 per cent of the dissolved calcium and magnesium load leaving the basin originates from non-karst rocks and 4.6 per cent is initially introduced by rainfall. River flows that are exceeded only 5 per cent of the time transport approximately 44 per cent of the annual dissolved load, while mean to low flows that occur for 75 per cent of the time transport 35 per cent of the annual solute load. This confirms the importance of low frequency-high magnitude events, but indicates also that in corrosion systems high frequency events of moderate to low magnitude can also accomplish significant work.     

The geomorphology and evolution of a large barrier spit: Farewell Spit,New Zealand

Farewell Spit is a 25 km long barrier spit that marks the end of a littoral drift system, almost 1000 km in length that runs along South Island, New Zealand. The spit is composed of barchan dunes over 20 m high, sand sheets over 1 km wide and vegetated linear dunes. Analysis of aerial photography indicates a rapid colonization of the spit by vegetation which has expanded in area by 75% since 1950. Vegetation colonization preferentially occurs on the southern side of the spit, with its northern margin characterized by barchan dunes which migrate at rates of up to 64 m/yr. Sand sourced from longshore drift appears to be the primary source of beach sediment, which is then transported into the dune field by the persistent westerly winds of the Roaring 40s. While there has been significant dune roll‐over on the surface of the spit, its overall area has remained much the same for the past 54 years. Occasional cyclone events cause erosion, but this is balanced by aeolian sediment transport. It would appear that extension of the subaerial portion of the spit is related to the development of shells banks at its downdrift end which are periodically welded to the main spit by dune extension. Farewell Spit therefore provides an ideal example of a barrier environment where longshore sediment supply and aeolian transport dominates geomorphic evolution. This differentiates the study site from other barrier environments where overwash or tidal inlet development often characterizes recent landform evolution. Copyright © 2010 John Wiley & Sons, Ltd.     

Semi-determinate hydraulic geometry of river channels,South Island,New Zealand

Data to describe the morphologic, hydrologic and sedimentologic characteristics of 72 South Island, New Zealand, rivers were collected and analysed. Nearly 70 per cent of variation in channel morphology is accounted for by differences in cross-sectional area, slope, and cross-section shape; only 53 per cent of the morphologic variability could be statistically ‘explained’ by the hydrologic and sediment variables used. The level of explanation varied for different morphologic variables; nearly 90 per cent of the variability in cross-sectional area could be explained, but aspect ratio (maximum depth divided by hydraulic radius) was completely independent. Apart from the inadequacy of the measured variables as indices of the true underlying controlling factors, and the imperfect measurement and sampling procedures, the low level of explanation is probably due to the influence of factors such as floodplain vegetation, high quasi-random variability in bark sediment character, boundary effects imposed by bedrock bluffs, and the precise sequence of flood events, none of which are easily quantified. In addition, observations indicate that there is a large random variation in channel form which cannot be related to any factor. An attempt to relate channel morphology to flow variability, using simple indices of the latter, was unsuccessful.     

Temporal variations in coda duration of local earthquakes in the Wellington region,New Zealand

Temporal variations of coda duration, relative to event magnitude, for local earthquakes near Wellington have been investigated. The region is one of plate convergence and subduction. The data consist of routinely made observations for events from 1978 through 1985 (1552 events), magnitude 1.6 to 5.2, depth 0.5 to 90.0 km. The observed average (over time) correspondence between duration (as measured from the origin time) and magnitude is reasonably well predicted by the single backscattering theory of coda formation for events of magnitude 4 or less; for larger events the observed durations are longer than predicted. This theory predicts that a temporal increase in scattering attenuation will reduce the coda duration relative to magnitude. Thus, any temporal changes in the duration-magnitude relation can be interpreted in terms of changes inQ. However, it is necessary to consider spatial biases since the observed durations are relatively long for shallow events and, for events of all depths, at stations situated in the south-eastern half of the region, usually believed part of the accretionary border. For both these situations other evidence would suggest that the scattering coefficient would be relatively high (lowQ). These observations may be due to a relatively high component of surface wave scattering and the importance of multiple scattering in the later part of the coda in regions of lowQ as suggested by finite element studies of coda formation. Despite the lack of any significant earthquakes during the 1978–1985 period there is nevertheless a significant temporal change in duration observed in the Wellington region: a change from relatively long to relatively short codas occurred in mid-1981. This change correlates well with changes in the rate of activity,b-value, radon emission, and ground tilt as derived from lake levels. It is not yet clear how all these parameters are related physically, but an episode of aseismic slip, or creep, along the plate interface below the region may have been the cause.     

Seismicity at White Island volcano, New Zealand: a revised classification and inferences about source mechanism

The classification of earthquakes at White Island volcano, New Zealand, has been revised to address problems in existing classification schemes, to better reflect new data and to try to focus more on source processes. Seismicity generated by the direct involvement of magmatic or hydrothermal fluids are referred to as volcanic, and that generated by fault movement in response to stresses caused by those fluids, regional stresses, thermal effects and so on are referred to as volcano-tectonic. Spasmodic bursts form a separate category, as we have insufficient information to classify them as volcanic or volcano-tectonic. Volcanic seismicity is divided into short-duration, long-period volcanic earthquakes, long-duration volcanic earthquakes, and harmonic- and non-harmonic volcanic tremor, while volcano-tectonic seismicity is divided into shallow and deep volcano-tectonic earthquakes. Harmonic volcanic tremor is related to sub-surface intrusive processes, while non-harmonic volcanic tremor originates close to active craters at shallow depth, and usually occurs during eruptive activity. Short-duration, long-period volcanic earthquakes come from a single source close to the active craters, but originate deeper than non-harmonic volcanic tremor, and are not related to eruptive activity. Long-duration volcanic earthquakes often accompany larger discrete eruptions. The waveform of these events consists of an initial low-frequency part from a deep source, and a later cigar-shaped part of mixed frequencies from a shallow crater source.     

Timing of volcanic, plutonic and geothermal activity at Ngatamariki, New Zealand

Four ⁴⁰Ar/³⁹Ar dates on mineral separates from fresh and hydrothermally altered volcanic and plutonic rocks from the Ngatamariki geothermal field indicate that andesitic volcanism took place in the eastern portion of the Taupo Volcanic Zone (TVZ) prior to 1.2 Ma and probably considerably earlier. These data significantly extend the onset and duration of andesitic volcanism in the east-central TVZ over previous estimates. Intrusive activity is represented at Ngatamariki by a dioritic pluton, the only such pluton yet recognized in the entire TVZ. Hornblende from the pluton yields a crystallization age of near 550 ka. Hydrothermal alteration spatially associated with the pluton produced sericite of a similar age. Overlying and postdating the most intense hydrothermal alteration zone is the Whakamaru Ignimbrite (or its equivalent) which was emplaced at 330 ka. Two distinct geothermal systems may have been active at nearly the same site from 550 ka to present. The most intense activity occurred before 330 ka and was associated with emplacement of the Ngatamariki diorite. This was followed by the less intense system that is currently active. The geothermal regime at Ngatamariki has, therefore, probably been active intermittently for at least 550 ka.     

Morphometric controls and geomorphic responses on fans in the Southern Alps,New Zealand

Morphometric variables associated with 41 debris‐?ow and 18 ?uvial fans and their basins in the Southern Alps of New Zealand are examined. The results show statistically signi?cant differences in the area, maximum elevation, relief and ruggedness (Melton's R) of the basin and the area, gradient, and apex and toe elevations of the fan between debris‐?ow and ?uvial sites. Concavity of the fan longitudinal pro?le also differs between the two fan types, although this could not be tested statistically. Most of these morphometric differences re?ect differences in processes and environmental controls on them. Discriminant analysis indicates that basin area and fan gradient best differentiate the two fan types by process. Moderately strong correlations exist, on both debris‐?ow and ?uvial fans, between basin area or Melton's R and fan area. Correlations between basin area or Melton's R and fan gradient are generally weaker. The results of this study also indicate that on debris‐?ow‐prone fans the fan gradient and basin Melton's R have lower thresholds which overlap little with upper thresholds associated with basins where only stream?ow reaches the fan. These thresholds may therefore have value in preliminary morphometric assessments of debris‐?ow hazard on fans in the Southern Alps. Copyright © 2004 John Wiley & Sons, Ltd.     

Sediment variability and bed material sampling in gravel-bed rivers

In a 350 m wide reach of the braided, gravel-bed Ashley River, the surface layer of the bed material was sampled in 141 areas of homogeneous graded sediment along seven cross-sections, and 30 kg bulk samples were collected at 86 randomly selected locations along the cross-sections. At one location, a single 854 kg sample composed of 28 subsamples was also collected. Analysis of the single large sample indicate that accurate determination of mean grain size D at that site requires, desirably, a sample of ～ 100 kg, but that samples in which the weight of the largest stone is less than 5 per cent of the total weight have unbiased estimates of D. Spatial variability of bulk material is such that 228 and 50 samples are needed to estimate D to ± 10 and ± 20 per cent respectively of the true value; requirements for estimating inclusive graphic standard duration are only 11 and 3 respectively. The grain-size distribution of the surface layer is only weakly related to the bulk material beneath. The results of ‘Wolman sampling’ along 12 cross-sections at two pace intervals (average 120 stones per cross-section) indicate that estimation of overall surface D to ±10 and ±20 per cent would require sampling along 64 and 14 cross-sections respectively. It is concluded that accurate characterization of bed sediment in gravel-bed rivers is very demanding of labour and resources, and that careful planning is needed to ensure efficiency and meaningful results.     

New predictive equations for Arias intensity from crustal earthquakes in New Zealand

Arias Intensity (Arias, MIT Press, Cambridge MA, pp 438–483, 1970) is an important measure of the strength of a ground motion, as it is able to simultaneously reflect multiple characteristics of the motion in question. Recently, the effectiveness of Arias Intensity as a predictor of the likelihood of damage to short-period structures has been demonstrated, reinforcing the utility of Arias Intensity for use in both structural and geotechnical applications. In light of this utility, Arias Intensity has begun to be considered as a ground-motion measure suitable for use in probabilistic seismic hazard analysis (PSHA) and earthquake loss estimation. It is therefore timely to develop predictive equations for this ground-motion measure. In this study, a suite of four predictive equations, each using a different functional form, is derived for the prediction of Arias Intensity from crustal earthquakes in New Zealand. The provision of a suite of models is included to allow for epistemic uncertainty to be considered within a PSHA framework. Coefficients are presented for four different horizontal-component definitions for each of the four models. The ground-motion dataset for which the equations are derived include records from New Zealand crustal earthquakes as well as near-field records from worldwide crustal earthquakes. The predictive equations may be used to estimate Arias Intensity for moment magnitudes between 5.1 and 7.5 and for distances (both r_jb and r_rup) up to 300 km.     

Tilting of active folds and drainage asymmetry on the Manawatu anticlines,New Zealand: a preliminary investigation

The apparent southward tilting of five anticlines in the Manawatu, North Island, New Zealand, is investigated using two drainage basin analysis techniques: the transverse topographic symmetry factor (T) and the asymmetry factor (AF). The techniques are used to determine whether westward‐flowing trunk streams on the gentle western limbs of five anticlines are migrating southwards across their drainage basins, concomitant with down‐tilting of the fold axes. While results show potential for the techniques to be applied in such a tectonic setting, our findings are inconclusive. This could be due to the incipient nature of the folds (c. 0·5 Ma in age), and it would seem more appropriate to apply the techniques in concert with detailed field analyses of river migration patterns. Copyright © 2008 John Wiley & Sons, Ltd.     

Estimation of Tsunami Hazard in New Zealand due to South American Earthquakes

We develop a probabilistic model for estimating the tsunami hazard along the coast of New Zealand due to plate-interface earthquakes along the South American subduction zone. To do this we develop statistical and physical models for several stages in the process of tsunami generation and propagation, and develop a method for combining these models to produce hazard estimates using a Monte-Carlo technique. This process is largely analogous to that used for seismic hazard modelling, but is distinguished from it by the use of a physical model to represent the tsunami propagation, as opposed to the use of empirical attenuation models for probabilistic seismic hazard analysis.     

类新西兰地震波形分析方法

对泰安地震台记录的2010年9月4日新西兰地震波形,进行震相、震源深度等地震参数辨析和分析,以具有代表性的2009年7月15日新西兰南岛西海岸远海地震波形和2010年4月5日墨西哥地震波形分析为辅,提出震中距大于83°、小于105°的类新西兰地震分析方法.根据两个辅震的波形特性,结合中国地震台网中心定位结果,对震级误差...     

In situ cosmogenic 10Be production-rate calibration from the Southern Alps,New Zealand

We present a ¹⁰Be production-rate calibration derived from an early Holocene debris-flow deposit at about 1000 m above sea level in the central Southern Alps, New Zealand, in the mid-latitude Southern Hemisphere. Ten radiocarbon ages on macrofossils from a soil horizon buried by the deposit date the deposit to 9690 ± 50 calendar years before AD2008. Surface ¹⁰Be concentrations of seven large boulders partially embedded in the stable surface of the deposit are tightly distributed, yielding a standard deviation of ～2%. Conversion of the ¹⁰Be measurements to sea level/high-latitude values using each of five standard scaling methods indicates ¹⁰Be production rates of 3.84 ± 0.08, 3.87 ± 0.08, 3.83 ± 0.08, 4.15 ± 0.09, and 3.74 ± 0.08 atoms g^?1 a^?1, relative to the ‘07KNSTD’ ¹⁰Be AMS standard, and including only the local time-integrated production-rate uncertainties. When including a sea level high-latitude scaling uncertainty the overall error is ～2.5% (1σ) for each rate. To test the regional applicability of this production-rate calibration, we measured ¹⁰Be concentrations in a set of nearby moraines deposited before 18 060 ± 200 years before AD2008. The ¹⁰Be ages are only consistent with minimum-limiting ¹⁴C age data when calculated using the new production rates. This also suggests that terrestrial in situ cosmogenic-nuclide production did not change significantly from Last Glacial Maximum to Holocene time in New Zealand. Our production rates agree well with those of a recent calibration study from northeastern North America, but are 12–14% lower than other commonly adopted values. The production-rate values presented here can be used elsewhere in New Zealand for rock surfaces exposed during or since the last glacial period.     

Research on the multifractal characteristics of the temporal—spatial distribution of earthquakes over New Zealand area

Using the Hill estimator,general multifractal characteristics of events in the New Zealand area have been dis-cussed.Results show that the spatial distribution of shallow events has apparent clustering characteristics,inde-pendent of the threshold magnitude;but for deep events these characteristics are not clear.While the time interval distribution has obvious clustering characteristics both for deep and shallow events,although with a different scal-ing range,the Hill estimates tend to indicate that the time interval distribution has a unifractal rather than a multi-fractal nature.All above reveal that the seismicity nature for shallow and deepevents is apparently different.