Evaluation of the advanced artificial athlete and Hall effect sensors for measuring strain in multi-layer sports surfaces

Computer models are a useful means to explore the loading behaviour of third generation (3G) artificial turf sports surfaces; however, measuring the material stress–strain behaviour under realistic high loading rates is challenging. Therefore, the purpose of this study was two-fold: to evaluate the advanced artificial athlete (AAA) for measuring strain behaviour of polymeric sports surfaces under high loading rates typical of player interactions; and to evaluate Hall effect sensors (HES) for measuring strain behaviour of an individual layer within multi-layer sports surfaces. An independent optical measurement system (GOM) provided gold standard sample deformation and strain. Forty-five impacts onto three test samples were measured simultaneously using the three systems. Poor agreement was found between AAA and GOM peak sample deformations and strain (systematic bias 2.4 mm, 95% confidence intervals ± 1.3 mm, strain RMSD 29%), largely attributable to errors in the AAA time of initial contact. Using a regression equation to correct this time led to much better agreement in AAA peak deformations and strain (systematic bias 0.1 mm, 95% confidence intervals ± 0.7 mm, strain RMSD 8%). Good agreement was found between the HES and GOM (systematic bias 0.2 mm, 95% confidence intervals ± 0.4 mm, strain RMSD 11%). The corrected AAA and HES methods can measure deformation of polymeric sports surfaces under realistic loading rates to an accuracy of < 1 mm. In terms of strain, errors increase with decreasing peak sample deformation indicating that both systems should be used with caution for peak deformations < ~ 4–5 mm.


Introduction
Third generation (3G) artificial turf is a common alternative to natural turf for sports pitches in climates where extreme weather, such as freezing temperatures or lack or water, can make it difficult to maintain a high quality natural turf [1]. Third generation artificial turf is constructed from a number of layers with a range of materials and properties ( Fig. 1) [2,3]. The bottom layer consists of an elastomeric shockpad designed to absorb impacts. On top of this lies the carpet, consisting of tufts of polymer fibres stitched into a canvas backing. The carpet layer also contains two infills: at the base is a stabilising infill, commonly sand, which adds weight and helps to keep the fibres upright; and above this is a performance infill, commonly rubber crumb particles, to provide a suitable compressible interface for player and equipment interactions.
The detailed design of 3G turf surfaces focuses on meeting the performance and injury risk requirements of the target sport [2]. Thus, understanding the response of the individual layers, as well as the combined surface system, to loading is key. This is recognised by the FIFA Quality Concept for football turf which uses a variety of mechanical tests to assess the performance of 3G turf to ensure it is safe for use [4]. Of particular importance is the vertical loading response, i.e. hardness or shock absorption, as it has been linked to both performance and injury risk [5,6]. From a research and development perspective, computer simulation models of the surface system, incorporating the individual layers, are used to improve understanding of the vertical loading response [7]. These models require as inputs the stress-strain properties of the individual layers of the surface system.
The most common method to characterise the stress-strain response of a material is to apply a known load and measure the resulting displacement [8]. Typically, this data is obtained using a specialist mechanical test device, such as a tensometer, due to the ease of use, controllability and repeatability they provide [9]. However, these devices are generally limited to relatively low loading rates (~ 3-10 Hz) which are not representative of players interactions with sports surfaces, e.g. the heel strike of heel-toe running has a frequency of 30-50 Hz [10,11]. This is particularly relevant for 3G turf where the material response is highly strain rate dependent [2]. A potential solution is to capture stress-strain data using a non-standardised test, an approach that has been used successfully to capture the high strain rate response of materials during dynamic impacts [12,13].
The advanced artificial athlete (AAA) is a mechanical device used in the FIFA quality concept to measure and assess the stiffness response of a surface (Fig. 2) [14]. It is designed to replicate the vertical force and contact time of a heel strike during heel-toe running [10,11,14]. The basic principle of the device is to produce a controlled vertical impact on a surface and output the impact-rebound properties [15,16]. Alongside the certification outputs of force reduction, vertical deformation and energy restitution [14], further processing of the data allows the stress-strain material response throughout the contact phase to be calculated [7]. However, whilst the AAA outputs have been compared to those from its predecessor, the Artificial Athlete [15,16], they have not been validated against an independent measurement system of known accuracy. For the AAA to be used to quantify the relevant strain response of sports surface systems, the accuracy of the outputs needs to be assessed.
For evaluation of multi-layered surfaces, such as 3G turf with a shockpad, it is important to be able to understand the deformation contributions of the individual layers to the total deformation. Hall effect sensors (HES) represent a relevant non-intrusive measurement method with the potential to assess deformation of the individual layers in a multi-layer surface system. A magnet and sensor act to indirectly measure the distance between them by producing a voltage proportional to the magnetic field at the sensor location [17]. By embedding a sensor between the carpet and shockpad layers of a 3G surface system (Fig. 1) and a magnet in the AAA test foot (Fig. 2), it may be possible to assess the deformation of the individual layers.
Thus, the purpose of this study was twofold: firstly, to evaluate the AAA for measurement of strain behaviour of sports surfaces; and secondly, to evaluate the potential  The advanced artificial athlete (AAA): 1. Support frame; 2. Electric magnet; 3. Falling mass; 4. Accelerometer; 5. Spring; 6. Test foot; 7. Test surface; 8. Concrete floor [14] use of HES to measure strain behaviour of an individual layer in multi-layer sports surfaces. To achieve the first purpose, the sample deformation output from the AAA was compared to that from the gold standard reference measurement system GOM (GOM mbH, Germany) during drop tests onto three different polymeric sports surface test samples. To achieve the second purpose, a bespoke fixture was created to allow a magnet to be embedded within an extension to the AAA test foot and a sensor grid designed to sit beneath each of the three test samples. This set-up allowed all three measurement systems to operate concurrently.

Equipment set up and alignment
The AAA was set up as detailed in the FIFA quality concept for football turf [14]. To replicate the human foot strike, a convex test foot is mounted to a linear spring above which a 20 kg mass is situated. A remote controlled electromagnet releases the mass which it is guided by a support frame throughout the drop and subsequent impact onto the surface. An accelerometer attached to the base of the falling mass records the acceleration throughout the drop (at a frequency of 9600 Hz) and is used to calculate the peak force, impact and rebound velocities and vertical displacement of the 70 mm diameter test foot into the surface [14].
In order to apply the added instrumentation of the GOM markers and the magnets for the HES an extended test foot was attached to the existing AAA test foot (Fig. 3). The extended test foot had an identical impact face to that of the AAA but was 50 mm longer allowing for additional room to embed magnets and apply GOM markers. The foot was created from solid nylon such that the additional mass was negligible but also stiff enough relative to the test samples not to affect the deformation results. It was also non-conductive and, therefore, did not interfere with the HES readings. The foot was attached to the AAA through an interference press fit and adhesive tape. A solid concrete floor was used as a base for all testing and AAA calibration drops were performed directly onto this concrete floor prior to sample testing to ensure that the AAA was operating within the required limits [14]. Pilot testing on a range of test samples confirmed that the addition of the extended test foot had a negligible effect on the AAA test measurements.
The GOM setup comprised of two Photron Fastcam SA1.1 (Photron, San Diego, CA) high speed video cameras with Titanar 50 mm lenses used to track 5 mm diameter passive markers positioned on the AAA falling mass and nylon impact foot (Figs. 3a, 4). The cameras were set up according to GOM recommendations [18] and placed on a raised bar just above floor level, parallel to the plane of impact and 915 mm away from the measurement volume. The cameras were gen-locked and calibrated using an independently certified GOM calibration board based on a capture volume centred on the AAA falling mass and extended test foot drop zone of 300 × 300 × 300 mm 3 (Fig. 4). This process gave a calibration deviation of less than 0.04 pixels, within the manufacturers recommendations [18]. A frame rate of 5400 Hz was selected to utilise the full camera resolution of 1024 × 1024 pixels. Shutter speed was set to 1/frame rate as no motion blur was apparent at the impact speeds of approximately 1 m/s based on a standard AAA drop height of 55 mm [14]. Lens aperture was set to f/5.6 again following GOM recommendations; however, this combined with the fast shutter speed meant additional lighting was needed to ensure adequate contrast to track the markers. A single flicker-free ARRI light was directed upon the AAA extended test foot (Fig. 4). The cameras were controlled manually to record from the approximate moment of release of the AAA falling mass through to after the first rebound.
Five linear sensors (Honeywell, USA) were used to estimate the vertical position of the extended test foot. A magnet was attached to the end of a socket screw that was embedded into the extended test foot and positioned so the magnet end face sat flush with the bottom of the foot (Fig. 3b). The sensors were arranged in a quincunx configuration with 15 mm spacing and embedded onto a flexible printed circuit board. The board was positioned under the AAA extended test foot and securely fixed directly to the concrete floor using double-sided carpet tape (Fig. 3a). The AAA was horizontally aligned directly above the central sensor manually using the live voltage output from this sensor; alignment was taken as the horizontal position that produced maximum voltage. Alignment was assessed after every five trials and realigned if appropriate. For each trial the voltage output from all sensors was recorded via SignalExpress (National Instruments, Austin, Texas, USA) sampling at 9600 Hz to match that of the AAA accelerometer.
Prior to testing the central sensor was calibrated using a tensometer (Instron Electropuls E3000). The sensor board was securely fastened to a flat aluminium plate directly below the tensometer crosshead where the magnet was attached. The magnet was then lowered towards the sensor from 30 mm until touching in 1 mm increments and at each distance the corresponding voltage from the HES recorded (Fig. 5). Whilst the field sensitivity (V/T) of the HES is linear with respect to magnetic field, the sharp increase in magnetic field with respect to distance from the pole resulted in an increase of the spatial sensitivity (V/ mm) as the magnet approaches the sensor. For the chosen probe current (6 mA) and the data acquisition system used (NI USB-6212, National Instruments, Austin, Texas) the sensor was effectively saturated at a separation of 6 mm. This was due to the voltage generated by the HES reaching the maximum deflection of the data acquisition card. As this saturation distance was not expected to be exceeded during sample testing, the choice of probe current ensured maximum sensitivity over the range of the measurements presented here. Cubic spline interpolation was used to create a transfer equation to convert a recorded voltage to a magnet-sensor distance (Fig. 5).
Three test samples were chosen for testing which differed in material and thickness to provide a range of loading response conditions (Table 1). All three samples were manufactured from polymeric materials for use as or within a sports surface. The first sample was designed for racket sports with a high stiffness to allow for fast play and high ball rebound (test sample RS). The remaining two samples were FIFA approved shockpads for 3G turf and were, therefore, designed to absorb impact forces (test samples SP1 and SP2). All samples were at least 300 × 300 mm in size. The thickness of each sample was obtained through calliper measurements taken at 10 positions on each sample while the density and material were obtained from the manufacturer datasheet.

Protocol
A total of 15 impacts were completed on each sample, comprising five impacts at each of three drop heights of 25, 55 and 85 mm with the mid drop height (55 mm) representing that used in the FIFA Quality Concept for 3G turf [14]. A 2-min break was included between each drop to allow the sample to fully recover. The use of three drop heights resulted in a range of impact energies for evaluation of the AAA and HES.
Prior to any impacts on a sample the vertical position of the top level of the sample, i.e. the zero sample deformation position, was defined for the GOM and HES systems. For GOM this was achieved by slowly winding down the falling mass and extended test foot towards the sample until the initial point of contact between the sample and extended test foot. Visual identification of the point of contact was not possible, therefore a repeatable solution was adopted. A sheet of paper was placed between the sample and extended test foot and contact was defined at the first instance at which resistance was felt when trying to slide the paper out. Images of the position of the AAA falling mass and extended test foot were then recorded using the GOM camera set up. For the HES the sample thickness measurements were used to define the distance between the middle sensor and magnet at which initial contact with the sample occurred (Table 1).

Data processing
The AAA falling mass acceleration time series data for each trial was exported from the AAA software as a.csv file. The AAA software, by default, filtered the acceleration data using a 2nd order low-pass Butterworth filter at 600 Hz prior to export [14]. The corresponding AAA falling mass velocity and displacement data were calculated using numerical integration (trapezium rule) and force was also calculated using Newton's second law. Extended test foot displacement was calculated from the falling mass displacement and spring compression, the latter calculated using the force data and spring stiffness [14]. The extended test foot displacement was subsequently set to zero at the time of peak downward velocity, i.e. the criterion used to define the time of initial contact with the test sample [14] (Fig. 6). Sample deformation was assumed equal to extended test foot displacement during contact, with peak sample deformation corresponding to the lowest position of the extended test foot.
The GOM Correlate Professional software was used to output a.csv file containing the global positions of the markers on the falling mass and extended test foot from the video footage recorded throughout each drop and the images defining the top level of the sample. A Fourier transform analysis was performed to establish the frequency content for the marker trajectories in each trial. On this basis, a low-pass 4th order Butterworth filter was used to filter the marker trajectories before further processing, with the falling mass markers filtered at 250, 300 and 500 Hz and extended test foot markers at 60, 50 and 50 Hz respectively for ascending drop heights. The sample deformation was calculated as the vertical distance of the extended test foot markers below their position defining the top level of the sample. Velocity and acceleration of the falling mass were also calculated through numerical differentiation of the displacement data (first order finite difference method) for direct comparison to the AAA data.
HES output voltages were filtered using a low-pass 4th order Butterworth filter at 50 Hz again determined based on the frequency content of the signals. For each trial, the Recticel 3G shockpad Bound polyurethane foam 24.07 ± 0.14* output voltage was converted to a magnet (extended test foot)-sensor distance throughout each drop using the calibration data (Fig. 5). Sample deformation was calculated through subtracting this distance from the unloaded sample thickness over the contact phase. For each measurement system, sample strain was calculated as the sample deformation divided by sample thickness (Table 1). For the AAA and GOM sample stress was calculated using the falling mass acceleration and mass together with the extended test foot area. To time synchronise the data from the AAA, GOM and HES it was assumed that the peak sample deformation occurred concurrently for all three systems. All data processing was conducted in Matlab (The Mathworks, Natick, MA).

Statistics
To evaluate the accuracy of the AAA and HES measurements of sample deformation against the gold standard GOM, Bland-Altman limits of agreement analyses were performed between AAA-GOM and HES-GOM on the peak deformations [19]. This technique is commonly used to compare likewise measurements from two measurement systems, thus allowing for comparison of differences. The mean difference in the Bland-Altman analysis provides an indication of the systemic bias between the two measurement systems (closer to zero indicates lower systematic bias) whilst the 95% confidence intervals (± 1.96 standard deviations) gives an indication of the random differences between the two measurement systems (larger confidence intervals indicates larger random differences). To further assess the systematic bias between systems, repeated measures t-tests were also conducted on the peak sample deformations. Where relevant, regression analysis was conducted to explore the scope of correcting the AAA and/or HES measurements to remove the systematic bias and give improved absolute agreement with the GOM measurements. To compare the full sample time series data throughout the contact phase root mean square differences (RMSD) were calculated between AAA-GOM and HES-GOM in falling mass acceleration and velocity as well as sample strain. Absolute values of RMSD were reported for each variable, however, to allow for comparison across multiple variables, these values were also normalised by dividing by the peak values of each variable across the time series and reporting a normalised RMSD (expressed as a percentage of the peak value of the variable). All statistical analysis was conducted in Matlab with significance set at p < 0.05.
To assess whether the accuracy of the AAA and HES strain measures were sufficient for material characterisation, it was recognised that when recording material load-deformation data under high loading rates, oscillations (noise) are typically observed in the resulting stress-strain plots. However, provided these oscillations are relatively small, the (real) material stress-strain response can be discerned and used for material characterisation [20]. On this basis, a threshold root mean square deviation (RMSD) of 10% of the peak strain was selected as an appropriate criterion indicating that the real strain response had been captured with sufficient accuracy. Fig. 6 Typical a acceleration b velocity and c displacement data from the advanced artificial athlete (AAA) during an impact onto the SP1 shockpad. T1 represents the time of peak downward velocity and the time at which the extended test foot displacement was set to zero (assumed initial contact with the sample). Extended test foot displacement was calculated by subtracting the spring compression from the falling mass displacement

Results
All three measurement systems were successful in capturing data for all trials with the exception of the 85 mm drops onto the RS sample where only one successful capture was obtained for the HES due to the sensors breaking. Consequently, this trial was disregarded in the data analysis and the results presented below for the HES are based on the remaining 40 trials.

AAA evaluation
Outputs from AAA and GOM showed good qualitative agreement in the falling mass acceleration and velocity profiles (Fig. 7). Peak acceleration and downward velocity were within 5% and 1% respectively and the RMSD between acceleration curves was 4.0 ± 0.6 m/s 2 and between velocity curves was 0.014 ± 0.003 m/s ( Table 2). Key metrics from the contact phase showed poor agreement in falling mass displacement (AAA 2-3 mm lower), peak sample deformation (AAA 1-3 mm lower) and time to peak deformation (AAA 1.0-4.5 ms lower) but better agreement in spring compression (≤ 0.4 mm difference) ( Table 3).
The Bland-Altman analysis confirmed that the AAA under-estimated peak sample deformation for all conditions (Fig. 8a). The systematic bias increased as peak sample deformation increased; from approximately 1 mm for the lowest peak deformations to 3 mm at the highest peak deformations (giving a mean difference of 2.4 mm), with the limits of agreement (± 1.3 mm) largely defined by this trend rather than by random differences. Given the systematic nature of these errors in AAA sample deformation and good overall agreement between the fundamental measurements of acceleration and velocity, the main source of error can be identified as the definition of initial contact with the sample. For the AAA this was the time of peak downward velocity and for GOM it was based on the independently measured position of the top level of the sample. These times systematically differed ( Fig. 7; Table 3) leading to systematic differences in the absolute vertical position of zero sample deformation and, therefore, in peak sample deformation (Fig. 7c). Peak downward velocity occurred after contact had been made with the sample (Fig. 7b) (where the acceleration crossed 0 m/s 2 rather than first deviating from − 9.81 m/s 2 ) leading to the AAA systematically under-predicting sample deformation (Fig. 7c).
Given that the differences between the AAA and GOM sample deformation values were predominantly systematic, a regression analysis between the two was conducted to explore the feasibility of correcting the AAA data to give closer agreement with the gold standard. As highlighted above, the main source of error came from the incorrect identification of the time of initial contact with the sample (Fig. 7) and the magnitude of error was a function of the peak sample deformation (Fig. 8). Thus, the correction was based on altering the AAA time of initial contact to better match that defined by GOM as a function of peak acceleration. A plot of peak acceleration against time difference for all 45 trials supported an exponential regression equation to be fitted with a RMSD of only 0.27 ms (Fig. 9): where T1 is the time of peak downward velocity (in ms), T2 is the time of initial contact with the sample (in ms) and a PK is the peak acceleration (in m/s 2 ). Thus, from the AAA accelerometer values for T1 and a PK , the corrected time of initial contact with the sample (T2) can be calculated.
To further assess the applicability of this correction, the AAA outputs were re-calculated including this correction to the time of initial contact (AAAnew in Tables 2, 3; Figs. 8b, 10). In the Bland-Altman analysis this reduced the systematic bias (mean difference) from 2.4 mm down to 0.1 mm and the random differences (limits of agreement) from ± 1.3 mm down to ± 0.7 mm (Fig. 8). The improved agreement was also reflected in the strain data where the strain RMSD compared to GOM decreased from 0.11 ± 0.02 to 0.03 ± 0.01 (Table 2; Fig. 10). Upon breakdown, both shockpad samples met the criterion for material characterisation as their respective RMSDs were within 10% of the peak GOM strain (SP1 5 ± 2% and SP2 4 ± 1%). Despite an improvement in RMSD, the RS sample did not fall within the required criterion for material characterisation (15 ± 5%; Table 2).

HES evaluation
Outputs from HES and GOM showed good overall agreement in peak sample deformation ( Table 3). The Bland-Altman analysis revealed that the peak sample deformation from the HES had no significant systematic bias (mean difference 0.2 mm) and small random differences of ± 0.3 mm (Fig. 11). There was a trend for the HES to over-predict peak sample deformation for the stiffer response conditions (smallest peak deformations) but under-predict peak sample deformation for the least stiff response conditions (largest peak deformations). Similarly, time to peak sample deformation was over-estimated for the RS and SP1 samples (by 1.5-2 ms) but was well matched by the SP2 sample (< 0.5 ms difference; Table 3; Fig. 12). This was also reflected in the sample strain data where the SP2 sample produced a better fit (0.017 ± 0.003) compared to the SP1 (0.058 ± 0.007) and RS (0.058 ± 0.004) samples (Table 2; Fig. 12). Only SP2 met the criterion for material characterisation with a strain RMSD of 4 ± 1% of peak strain with the remaining samples ≥ 12%. Linked to this, despite the good agreement in peak sample deformation, at initial contact (based on GOM) the differences were much bigger with sample deformation overpredicted by at least 1 mm for both the RS and SP1 samples, while agreement was much better for the SP2 sample. RS and SP1 samples both deformed in the range 7-15 mm from the sensor. In contrast the SP2 sample deformed 11-24 mm from the sensors and, therefore, operated in a lower sensitivity region (~ 0.16 V/mm) compared to the other two samples (~ 0.04 V/mm) (Fig. 5).

Discussion
Computer simulation models are a useful tool to improve understanding related to the vertical loading response, and hence cushioning effects, of sports surface systems  [7]. However, characterising the component material stress-strain properties under the high loading rates typical of player and equipment interactions with sports surfaces using traditional test methods is challenging. The AAA is a mechanical drop test device specifically designed to replicate certain features of player-surface interactions, including the high loading rates. However, its ability to provide accurate stress-strain characteristics throughout surface contact has yet to be investigated. Hence, the first objective was to evaluate the use of the AAA to measure the strain behaviour of sport surfaces. Whilst the AAA provides data on the overall deformation of a sport surface system, it cannot provide information on the contribution of the individual layers to this overall deformation in a multi-layer system (as often used in sport). Measurement of the individual layer deformations is relevant to surface system research and development through either physical testing or virtual modelling. Hall effect sensors provide a non-intrusive distance measurement method with potential to be applied to intra-layer deformation measurement. Thus, the second objective was to evaluate the potential of HES to support measurement of the strain behaviour of individual layers in multi-layer surface systems. Three test samples and three AAA drop heights were used to provide a wide range of sample responses to dynamic loading (confirmed by the wide range in peak falling mass decelerations from 60 to 330 m/s 2 and in peak sample deformations from 2.8 to 12.0 mm) while simultaneously providing good control over the test conditions to support the evaluation of the AAA and HES measurement systems. The optical measurement system GOM was used as the gold standard reference for surface sample deformation and strain. It has a manufacturer reported accuracy of 25 μm/m field of view which was also demonstrated in a recent independent study [21]. Measurement volume, camera set-up and calibration procedures similar to Leach [21] were applied in this study to support its use for the evaluation of the deformation and strain outputs from the AAA and HES measurement systems. It should be noted that while there was high confidence in the accuracy of the GOM positional data as the gold standard, numerical differentiation was required to obtain GOM velocity and acceleration. Consequently, confidence in these values was lower and, on this basis, GOM velocity and acceleration were considered for comparative purposes only.
Strain, calculated from the AAA/HES sample deformation data and sample thickness, was deemed suitable for material characterisation if the RMSD (compared to the GOM strain data) was less than 10% of the peak strain. Previous studies investigating impacts under high loading rate conditions found that material characterisation could successfully be achieved if the main features of the material response could be discerned from the stress-strain Table 3 Contact phase peak values Falling mass displacement and spring compression for the AAA and GOM. Sample deformation and time from initial contact to peak sample deformation for all three measurement system Sample  [20]. However, this assessment was not quantified and, therefore, an estimated criterion of strain RMSD less than 10% of the peak strain was proposed here, which visually appeared to meet the qualitative criterion detailed above. Expressing the RMSD relative to the peak strain provided a relevant means of accounting for the differences in sample thickness. Calculating the strain RMSD over the whole stress-strain curve was considered to be important since the samples used in this study, along with many of the materials used in 3G surfaces, are hyperelastic suggesting that the rate of change in strain as well as energy loss in recovery are important features to capture. Poor agreement was obtained between AAA and GOM for peak sample deformation and strain. The AAA underpredicted peak deformation by between 25 and 47% and the overall average strain RMSD was 29% of peak strain (Tables 2, 3; Figs. 8a, 10). With much better agreement in acceleration and velocity data, the poor agreement was identified as resulting from the method used by the AAA to determine the time of initial contact of the extended test foot with the sample, i.e. the zero sample deformation conditions. The AAA used the time of peak downward velocity of the falling mass as initial contact corresponding to the acceleration crossing 0 m/s 2 ; however, initial sample contact occurs at an earlier time point where the acceleration first deviates from gravitational acceleration with the GOM time of initial contact providing a much better estimate of this instant (Fig. 7). As expected based on the foregoing, the error in AAA time of initial contact increased for the more compliant test sample, characterised by lower peak impact decelerations (Fig. 9); however, given its systematic nature correction methods were explored.
A correction method for the AAA time of initial contact was developed based on fitting an exponential regression equation relating the AAA falling mass peak deceleration to the time difference between GOM initial contact and the original AAA initial contact based on peak downward velocity (Fig. 9). The applicability of this regression equation was explored by re-calculating the AAA sample deformation data (AAAnew) and repeating the comparison to the GOM sample deformation and strain. The correction resulted in much better agreement; for peak sample deformation there was no significant systematic bias (mean difference of < 2% or < 0.1 mm) and random differences of ± 0.7 mm (95% confidence intervals). The strain RMSD was reduced to 8% of peak strain (Tables 2, 3; Figs. 8b, 10), thus meeting the criterion for material characterisation. The correction appeared to work equally across all test samples and drop heights for peak sample deformation; however, for strain RMSD, tests onto sample RS gave higher RMSDs (15 ± 5% of peak strain) than those onto the remaining two samples (≤ 5%), largely due to the smaller (unloaded) thickness of this sample and smaller peak deformations. Given the empirical nature of the correction method, the recommendation is that it should only be applied to impacts with peak falling mass decelerations in the range over which it was developed (60-330 m/s 2 ). With this caveat in place, the results support the use of the AAA drop test, with the inclusion of the correction to the time of initial contact, as a suitable method for measuring the deformation of sports surface systems under high loading rate conditions, giving sub 1 mm accuracy. In terms of strain measurement, the errors naturally scale with the surface or layer thickness indicating that it should be used with caution for thin, stiff surfaces experiencing peak deformations of < 4-5 mm. Very good agreement was obtained between HES and GOM for peak sample deformation but poorer agreement for strain RMSD (Tables 2, 3; Figs. 11,12). Peak sample deformation systematic bias and random differences were 0.2 mm and ± 0.3 mm (95% confidence intervals) respectively, while the overall average strain RMSD was 11% of peak strain, just outside the criterion for material characterisation. There were, however, clear inter-sample differences in the level of agreement; notably, while SP2 sample had the lowest strain RMSD (4%, i.e. within the criterion for material characterisation), it had the highest mean difference in peak sample deformation (underestimated by 0.5 mm).
To explore the reasoning behind these observations, initially the HES calibration curve was considered with respect to measurement sensitivity (Fig. 5). Differences in the three sample thicknesses and peak deformations resulted in each sample deformation range falling within a different range on the HES calibration curve and, therefore, within a different sensitivity range. SP1 and RS samples were thin and stiff, relative to SP2, thus the magnet was closer to the sensor board during deformation, resulting in a higher average sensitivity (0.15 and 0.17 V/mm respectively). In contrast, SP2 was double the thickness and had a sample deformation range that sat on the steeper part of the calibration curve corresponding to a much lower average sensitivity (0.04 V/mm). On this basis, it might be expected that the HES errors would increase with increasing magnet-sensor distance; however, this wasn't observed in the test results indicating that sensor sensitivity wasn't the main factor limiting the performance of the HES measurement system. Given the relatively limited range of magnet to sensor distance over which a HES operates (particularly with high sensitivity), these results reinforce the importance of careful selection of magnet properties for the specific application.
To further explore the errors in the HES measurements, horizontal movement of the AAA device during testing leading to poor vertical alignment between the magnet and sensor was also considered. The HES system relied on an accurate manual vertical alignment between the sensor and magnet to precisely maintain them on the same vertical axis. Whilst the alignment protocol was completed at the start of every set of five trials (sample and drop height), vibration of the AAA device during each drop may have caused some horizontal movement of the device relative to the sensor which was hard to perceive visually. Analysis of the GOM markers on the extended test foot in the horizontal plane across each set of five trials revealed a small permanent horizontal displacement of the extended test foot of up to 0.4 mm. Although the magnitude of this horizontal displacement is too small to substantially explain the HES test sample results, it reinforces the importance of careful vertical alignment of the magnet over the sensor to maximise the accuracy of the HES deformation and strain measurements.
The AAA is used globally for certification of a sport surface system in comfort and regard safety, for example in the FIFA Quality Concept for football turf [14]. The three variables derived from the AAA test used to define surface performance are; force reduction, vertical deformation and energy restitution. It is relevant to consider the implications of the results of this study to these outputs. The calculation of force reduction relies on the peak falling mass deceleration and, therefore, is unaffected by this study findings. Vertical deformation is the equivalent of the peak sample deformation reported in this study and; the AAA has been demonstrated to systematically underpredict this variable due to a systematic error in the time of initial contact with the surface/sample. Given that the acceptable limits of force reduction in the FIFA standard correspond to peak falling mass decelerations in the range 90-145 m/s 2 , the corresponding error in time of initial surface contact can then be estimated as 3-4 ms which may result in an under-estimation of vertical deformation by as much as 2-3 mm. It should be noted the FIFA limits of acceptability for vertical deformation of a football turf surface is 4-11 mm. Although the potential error of 2-3 mm may appear large, the football turf standards are based on benchmarking to good quality natural grass such that the current limits are empirical rather than based on any fundamental absolute values. Hence systematic bias may not be an issue in regard to surface safety implications. Notably, the Quality Concept requires vertical deformation to be reported to the nearest 0.5 mm which corresponds to the standard deviation of the random errors observed in this study. Energy restitution is the ratio of the falling mass energy immediately post-impact to pre-impact. Its calculation is based on the peak upward velocity (postimpact) and peak downward velocity (pre-impact) of the falling mass. This study has demonstrated that the peak downward velocity occurs after impact (by 3-4 ms for surfaces tested within the Quality Concept) and by a similar reasoning, the peak upward velocity will occur before the extended test foot loses contact with the surface. For the three surface samples and three drop heights used in this study, applying the erroneous peak velocity criterion to estimate pre-and post-impact velocities resulted in energy restitution values 4.0 ± 1.8%, 2.0 ± 3.0% and 8.4 ± 4.5% smaller (for RS, SP1 and SP2 respectively) than using the correct pre-and post-impact velocity method. At present, however, the FIFA standard requires energy restitution to be reported, but does not stipulate limits of acceptability. In contrast the World Rugby standard [22], based on the same AAA test method and empirical basis, stipulates an acceptable range of energy restitution as 20-50%.