Adaptation of reach action to a novel force-field is not predicted by acuity of dynamic proprioception in either older or younger adults

Healthy ageing involves degeneration of the neuromuscular system which impacts movement control and proprioception. Yet the relationship between these sensory and motor deficits in upper limb reaching has not been examined in detail. Recently, we reported that age-related proprioceptive deficits were unrelated to accuracy in rapid arm movements, but whether this applied in motor tasks more heavily dependent on proprioceptive feedback was not clear. To address this, we have tested groups of younger and older adults on a force-field adaptation task under either full or limited visual feedback conditions and examined how performance was related to dynamic proprioceptive acuity. Adaptive performance was similar between the age groups, regardless of visual feedback condition, although older adults showed increased after-effects. Physically inactive individuals made larger systematic (but not variable) proprioceptive errors, irrespective of age. However, dynamic proprioceptive acuity was unrelated to adaptation and there was no consistent evidence of proprioceptive recalibration with adaptation to the force-field for any group. Finally, in spite of clear age-dependent loss of spatial working memory capacity, we found no relationship between memory capacity and adaptive performance or proprioceptive acuity. Thus, non-clinical levels of deficit in dynamic proprioception, due to age or physical inactivity, do not affect force-field adaptation, even under conditions of limited visual feedback that might require greater proprioceptive control. Supplementary Information The online version contains supplementary material available at 10.1007/s00221-020-05997-3.


Introduction
Typical healthy ageing brings about a number of notable changes to the neuromuscular system which influences the ability to control simple movements. This includes loss and remodelling of musculature and motor units (Lexell 1995;Morley et al. 2001;Slack et al. 1979), as well as degeneration of the peripheral nerves and neuromuscular junction (Jacobs and Love 1985;Valdez et al. 2010). Collectively, these changes contribute to the characteristic increases in variability, duration and online corrections that are seen in basic upper limb movements of older adults (Contreras-Vidal et al. 1998;Darling et al. 1989;Helsen et al. 2016;Ketcham et al. 2002;Kitchen and Miall 2019;Seidler et al. 2002;Yan et al. 2000). Some of these kinematic features are also thought to play a compensatory role to preserve endpoint accuracy for older adults during targeted reaching tasks (Helsen et al. 2016;Kitchen and Miall 2019;Lee et al. 2007;. Ageing also involves degeneration of proprioceptive sensory organs in the neuromuscular system responsible for sensations of limb position and motion (Proske and Gandevia 2012). Human and animal studies have shown age-related degeneration of the muscle spindle, including a reduced number of intrafusal fibres (Kararizou et al. 2005;Swash and Fox 1972), increased capsular thickness (Swash and Fox 1972) and degraded morphology of primary spindle endings (Kim et al. 2007) that place limits on the perceptual acuity of proprioceptive sensations in advanced age. Indeed, a broad range of proprioceptive assessment techniques-both 1 3 passive and active in nature-have shown age-dependent loss of proprioceptive acuity (Adamo et al. 2007;Cressman et al. 2010;Herter et al. 2014;Lei and Wang 2018;Schaap et al. 2015; for review see Goble et al. 2009), and acceleration of loss with physical inactivity (Adamo et al. 2009;Helsen et al. 2016;Kitchen and Miall 2019;Wright et al. 2011). Given that complete loss of proprioceptive sensation has profoundly detrimental effects on upper limb movements (Miall et al. 2018;Sarlegna et al. 2010;Yousif et al. 2015), it is clear that controlled motor performance is dependent on these sensory signals. Yet in spite of reports showing age-dependent proprioceptive decline and age-dependent reduction in motor function in the lower limb (Hurley et al. 1998;Lord et al. 1991;Sorock and Labiner 1992;Wingert et al. 2014), only recently has there been increased interest in this issue for the upper limb. Helsen et al. (2016) investigated this sensorimotor relationship in the wrist using two proprioceptive assessment tasks (passive movement detection and ipsilateral position matching) and a rapid movement task to visual targets. In line with other reports (Adamo et al. 2009;Kitchen and Miall 2019;Wright et al. 2011), physically inactive older adults were found to have reduced proprioceptive acuity, with increased movement detection thresholds and position matching errors. They also made more online corrective adjustments during movements, though endpoint accuracy was relatively unaffected. Critically, no relationships were found between proprioceptive acuity and either aiming errors or movement kinematics. However, in order to understand this relationship in a more naturalistic or realworld setting, the proprioceptive tests should be active, to incorporate sense of effort and motor efference which can influence position sense ('t Hart and Henriques 2016;Smith et al. 2009). Active testing can even reduce position matching errors in both older and younger adults (Erickson and Karduna 2012;Langan 2014;Lönn et al. 2000). We therefore tested groups of younger and older adults on two tasks that separately assessed active proprioceptive acuity and rapid target reaching performance (Kitchen and Miall 2019). Like Helsen et al. (2016), we found age effects on motor performance. We also found that physically inactive older adults had increased systematic (but not variable) proprioceptive errors, but there were no associations between proprioceptive acuity and motor errors for either age group (Kitchen and Miall 2019).
However, those tasks emphasised speed, where it seems likely that sensory feedback control is minimised in favour of predictive, feedforward processes (Miall and Wolpert 1996;Wolpert et al. 1995). Thus, the question remains as to how sensory impairments might influence performance on tasks which require greater sensory feedback control (for review see Shadmehr et al. 2010). To perturb proprioceptive feedback specifically, external forces can be applied during discrete targeted movements (Krakauer et al. 1999;Sarlegna et al. 2010;Shadmehr and Mussa-Ivaldi 1994). Unlike tasks with visual feedback perturbation Bock 2005;Buch et al. 2003;Contreras-Vidal et al. 2002;Hegele and Heuer 2010;Seidler 2006; Vandevoorde and Orban de Xivry 2019), a number of studies have indicated a minimal role of ageing in the ability to adapt movements to novel force environments (Cesqui et al. 2008;Huang and Ahmed 2014;Rajeshkumar and Trewartha 2019;Reuter et al. 2018;Trewartha et al. 2014). Nonetheless, several age-dependent performance predictors have been highlighted, including increased muscle co-contraction (Huang and Ahmed 2014), reduced central processing of kinematic errors (Reuter et al. 2018) and reduced spatial working memory capacity (Trewartha et al. 2014). Furthermore, younger adults demonstrate a tight relationship between proprioceptive perception and motor performance when adapting to novel force-fields (Haith et al. 2008;Mattar et al. 2013;Ohashi et al. 2019;Ostry et al. 2010). Moreover, whilst adaptation to forces is possible in younger adults without visual feedback of hand position (Franklin et al. 2007;Lefumat et al. 2015;Scheidt et al. 2005), all previous studies with older adults have included either full (Cesqui et al. 2008;Huang and Ahmed 2014;Rajeshkumar and Trewartha 2019;Trewartha et al. 2014) or early visual feedback (first half of the movement; Reuter et al. 2018). Hence, the interactions of age and proprioceptive acuity on dynamic motor adaptation remains poorly understood, and the contributions of visual and proprioceptive feedback to motor adaptation have not been dissociated.
We have therefore investigated the relationship between dynamic proprioception and force-field adaptation in older and younger adults. Participants were given either full or limited visual feedback during force-field adaptation, and dynamic proprioceptive assessments were made before and at intervals throughout the adaptation task. Spatial working memory capacity and physical activity status were also examined, collectively providing a comprehensive overview of reach adaptation to novel field dynamics and its interaction with proprioceptive acuity across the lifespan.

Participants
A total of 34 older (13 males, 75.6 ± 7.2 years) and 35 younger (3 males, 19.1 ± 0.9 years) adults participated in the experiment, all of which were right-handed (30 or higher on the 10-item Edinburgh Handedness Inventory; Oldfield 1971). Participants were excluded from the study for history of neurological illness, carpal tunnel syndrome, arthritis or similar movement pains or limitations in the arm, 1 3 wrist or fingers. Participants were also screened to ensure no prior experience with adaptation to novel field dynamics. Older adults completed the Montreal Cognitive Assessment (MoCA) and were only included in the analysis if they scored 26 or above out of 30, which is indicative of normal cognitive functioning (Nasreddine et al. 2005). This study was approved by the University of Birmingham ethics panel and all participants gave written consent prior to experiment performance.

Physical activity
Self-reported physical activity levels were measured using the Physical Activity Scale for the Elderly (PASE; Washburn et al. 1993) and the IPAQ-Short questionnaire (Craig et al. 2003) for older and younger participants, respectively. Based on a median score threshold of 155.4 from the PASE questionnaire, older adults were assigned to either a physically active (n = 17, 219.4 ± 61.9) or inactive (n = 17, 108.7 ± 40.7) group. Likewise, a median score of 1988 MET-min per week on the IPAQ was used to sub-divide younger adults into physically active (n = 17, 3763.9 ± 1253.2 MET-min per week) and inactive (n = 18, 1082.7 ± 531.1 MET-min per week) groups. These sub-groups were used to examine interactions of physical activity and proprioception.

Spatial working memory assessment
Prior to completing the experiment, spatial working memory capacity was assessed using a modified search task previously used with older adults and cognitively impaired patient populations (Kessels et al. 2010;Van Asselen et al. 2005). Participants were presented with sets of 4-10 blue tiles on a computer monitor. They were instructed to search for one hidden target, by selecting tiles one at a time. They continued this until they chose the target tile, briefly revealed as a green tick symbol before reverting back to blue (a sequence we term a "search"). Participants were required to remember this tile location and then continue searching for a new target, which might appear on any of the remaining tiles. They had to complete searching until all possible targets had been found (i.e. as many searches as there were tiles). In any given search, if the participant chose a tile which they had previously selected in that search, a red cross symbol briefly appeared; this was defined as a within-search error. Similarly, if participants selected a tile which had been the target in an earlier search, they also saw a red cross; this was a between-search error. The purpose of the task was, therefore, to find all of the green ticks whilst minimising the number of search errors. All participants completed 2 trials at each of 4, 6, 8 and 10 tiles after an initial 2 trials with 4 tiles as a familiarisation. To minimise distractions, participants also wore headphones that played distinct audio tones associated with selection of a blue tile, finding a target, or making an error. Between-and within-search errors were averaged across the two attempts at each number of tiles and then compared between age groups. After reviewing the data, we found that only a single error was made across all participants in the 4-tile condition, so it was not included in the main analysis; this left the 6-, 8-and 10-tile conditions. Total search errors (the sum of between-and within-search errors across these three tile conditions) were also used as a correlate for performance in the other tasks.

Main experimental set-up
The main tasks were performed on a low-inertia and nearfrictionless 2D-planar manipulandum (vBOT; Howard et al. 2009b) which made it possible to record reaching movements in a 40 × 64 cm workspace (Fig. 1a). Participants placed their forehead on a padded headrest and looked down onto a large, horizontally mounted, mirrored surface that reflected targets and visual feedback images from an LCD screen directly above. Participants grasped the manipulandum handle with their right hand underneath the mirrored surface, which blocked any direct vision of the arm or hand. The start position and target were displayed as white and grey 1 cm radius markers, respectively, with the handle position displayed as a white 0.5 cm radius marker, when feedback was available. Start position was located in the midline, 8 cm into the workspace (roughly 28 cm from participant's chest) with the target located 20 cm directly ahead along the mid-sagittal axis. The start position and target remained in the same position and were kept visible at all times during task performance.

Task A. Force-field adaptation
Procedure A summary of this task is displayed in Fig. 1b. To begin, participants first moved to the start position, which turned orange, and then made a reaching movement towards the target. Once they had moved the required 20 cm ahead, they intersected a soft virtual wall which ran orthogonally through the target. Regardless of whether the target was hit or not, it then displayed an 'explosion' graphic that provided participants with colour-coded feedback regarding the speed of their movement compared to a desired movement duration of 350-450 ms (Fig. 1b). The target then turned grey again and participants were actively guided back to the start position by a spring force (500 N m −1 , 1 N m s −1 damping) for the next trial.
Movements during null-field trials were unconstrained. However, force-field trials were performed in the presence of a velocity-dependent force-field applied to the vBOT handle as described by Eq. 1.
where f x and f y refer to the imposed forces on the vBOT handle, FS refers to the field strength which was held constant at 15 N m −1 s −1 , and v x and v y refer to the velocity of the vBOT handle. This "curl-field" (Shadmehr and Mussa-Ivaldi 1994) perturbed movements in a clockwise direction that deviated movements from their normally linear path towards the target (Fig. 1b).
(1) Fig. 1 Experimental set-up, tasks and design. a vBOT handle underneath reflective surface that displayed images during task performance. b Force-field task and visual feedback conditions. The Vision− group (left column) were shown a horizontal white bar as a hand position cursor (which provided movement feedback only along y-axis regarding distance to target), as well as white cross terminal position feedback. The Vision+ group (right column) were provided with full visual feedback (white circular cursor) throughout. Purple arrows show forces imposed by velocity-dependent field (lower row), with furthest right panel illustrating feedback provided at the target, after trial termination, on movement speed. c Dynamic propriocep-tion task. Participants actively moved through a minimum jerk pathway constrained by stiff virtual walls before verbally indicating which side of the target they felt they had been guided to ("Square" or "Circle"). d Experimental design. Participants were randomly allocated to either the Vision+ or Vision− group before completing the spatial working memory task (SWM) and brief familiarisation task. Blocks (P1-P4) of channel trials (Chan) for proprioceptive assessment were interleaved between target reaching blocks of null-field (red outlined boxes) and force-field (green outlined boxes) trials. Data were averaged across early and late portions of blocks for statistical analyses Catch trials were presented randomly within the final three of every set of five null-and force-field trials. In these trials, hand position was constrained to a linear path between start position and target using stiff virtual walls (stiffness: 2000 N m −1 with 10 N m s −1 damping imposed by vBOT motors). Forces against the channel wall were used as a measure of participants' compensation for the force-field during the course of adaptation. Visual feedback of distance without direction was presented as a white semi-circular arc whose radius increased proportionally to the linear distance travelled towards the target (see Lago-Rodriguez and Miall 2016). Coloured speed feedback and active guidance back to start were still provided for these trials.
Participants were informed that they may experience some forces at the vBOT handle during target reaching movements, but that their objective was to try and keep movements as straight and as accurate as possible (keeping the white cross on target for Vision− and hand position cursor on target for Vision+), whilst also trying to maintain appropriate movement speed (green explosions).

Visual feedback conditions
Participants were given either full visual feedback (Vision+) or limited visual feedback (Vision−) during the experiment. For the Vision+ sub-group, the white circular hand position cursor was available at all times. The Vision− group was shown a moving horizontal white bar as a hand position cursor, which only provided visual feedback regarding movement distance to the target (Fig. 1b, left column). This meant that any lateral movements (and hence features of the perturbation) were hidden from view at all times. To reduce chances of proprioceptive drift occurring (Brown et al. 2003a, b), terminal position feedback (the location at which participant intersected the endpoint orthogonal soft wall) was also provided at the end of each movement (1 cm width white cross). Movement speed feedback (coloured explosion) and active guidance back to the start position after trial termination were identical for Vision+ and Vision− subgroups. Visual feedback during catch trials was identical, regardless of sub-group (see Task A-Procedure above).

Outcome measures and analysis
Lateral hand deviation from the linear path connecting start position and target was recorded at peak hand velocity (peakvelocity lateral deviation; PVLD). For catch trials, the force generated against the channel wall was compared with the ideal force necessary to perfectly compensate for the forcefield given the velocity profile of the movement. Specifically, we fitted a least-squared linear regression model (without an intercept) to the observed versus optimal orthogonal forces for each trial and used the slope as an adaptation index (Huang and Ahmed 2014;Trewartha et al. 2014); a value of +1 indicated complete compensation for the force-field.
To examine kinematic performance, we recorded peak hand velocity, total movement duration and the time to peak velocity (expressed as a percentage of total movement duration). Movement duration was defined as the time taken to move from 1 cm beyond start position to 1 cm before target. Trials where total movement duration was greater than 1.5 s were not included in the analysis (this equated to < 1% of trials).
To monitor performance change across the experiment, adaptation and kinematic measures were averaged across the first and last 15 trials (3 catch trials, and 12 null/force field) in a given block which we term "early" or "late", respectively (this excluded the additional 20 null and 20 forcefield trials at the start of the adaptation and washout blocks, respectively; see Fig. 1d for details). Extent of adaptation was then defined as the positive difference between early and late performance in the adaptation block: reduction in lateral error (early minus late PVLD) and increase in adaptation index (late minus early).

Procedure
The details of this task have been described previously (Kitchen and Miall 2019) and are illustrated in Fig. 1c. Participants made reaching movements towards the target that were tightly constrained to a pre-determined trajectory using stiff virtual walls (stiffness: 2000 N m −1 with 10 N m s −1 damping imposed by vBOT motors). This channel deviated the hand path laterally from the target through a minimum jerk profile. No forces were applied in the forward direction; no visual feedback was provided. At the end of the movement, a white square and circle appeared at constant positions on the left-and right-hand sides of the target, respectively. The participant then made a verbal 2AFC response ("Square" or "Circle") to indicate the side of the target that they felt they had been guided to. Following their verbal response, participants were actively guided back to the start position by a spring force (500 N m −1 , 1 N m s −1 damping) before starting the next trial. The size and direction of the lateral deviation was manipulated on a trial-by-trial basis by two adaptive staircase sequences.

PEST sequences
The lateral deviation imposed by the virtual channel walls was defined by two randomly interleaved PEST sequences (Taylor and Creelman 1967), with one starting from the left side of the target ("Square") and the other from the right ("Circle"). This was designed to minimise awareness 1 3 of sequence progression and improve convergence of lateral deviations towards participants true perceptual thresholds. Each sequence started at a deviation magnitude of 3 cm (± 0.05 cm added noise) with an initial step size of ± 1 cm. This step either increased or decreased the deviation magnitude (or "level") according to the cumulative accuracy of the verbal responses from the two repeats which were performed at each level of the sequence. If both responses were correct, then the deviation magnitude would decrease, and if both were incorrect it would increase. When there was a tie (1 correct and 1 incorrect response), a third repeat at that level was performed, with the subsequent response determining the overall accuracy of the responses to that level. Whenever the sequence reversed direction (a successive step increase and decrease in deviation magnitude or vice versa), the new step size became half of the previous one, i.e. from 1 to 0.5 cm at the first reversal. The interleaved PEST sequences continued in this manner for 50 trials, which constituted one proprioceptive assessment block (Fig. 1d).

Outcome measures and analysis
The verbal responses indicating perceived limb position were converted to binary values ("Circle" = 1, "Square" = 0) and expressed as a proportional response probability for each level of the deviation before being fitted with a logistic function using the Matlab function glmfit. To remove the effects of outlying response values on curve fitting, data points with a Pearson residual that was greater than 2.5 standard deviations away from residual mean were excluded from analysis (this equated to < 1% of data). From the logistic, we then estimated the bias and the uncertainty range as indices of proprioceptive acuity. The bias represents a systematic error in limb perception and corresponds to the 50th percentile of the logistic function; a positive bias here represents a perception of hand position to be further to the right ("Circle") of the target and negative bias indicates a perceptual error to the left ("Square"). The uncertainty range gives a variable error of perceived hand position and is given as the interval between the 25th and 75th percentiles of the logistic function.
In addition to verbal perceptual responses, the average movement velocity was calculated for the portion of the movement from 1 cm beyond the start position to 1 cm before the target. The average orthogonal forces exerted on the virtual channel walls were also recorded in the middle of the final, straight portion of the movement (16-19 cm of the 0-20 cm movement). These measures were used as correlates for proprioceptive measures to ensure minimal confounding effects of effort on perceptual judgements of limb position (Smith et al. 2009).

Experimental design
The full design is shown in Fig. 1d. The experiment started with a familiarisation of 20 null-field trials (self-selected speed with full hand position cursor feedback regardless of visual feedback condition) and 8 proprioceptive channel trials. Participants then performed the first block of 50 proprioception trials (P1), followed by a block of 60 target reaching movements trials (48 null-field and 12 catch; "Null"), before another proprioception block (P2). The main block of target reaching trials was then performed, consisting of 16 null-field trials (interleaved with 4 catch trials) and 120 force-field trials (interleaved with 30 catch trials; "Adaptation"). This was immediately followed by the third proprioceptive block (P3) and a final block of target reaching trials consisting of 16 force-field trials (interleaved with 4 catch trials) and 48 null-field trials (interleaved with 12 catch trials; "Washout"). The fourth and final proprioceptive block (P4) was then performed to finish.

Statistical analysis
Both proprioceptive and adaptation data were normalised to baseline performance (block P1 for proprioception, late Null data for adaptation; Fig. 1d) before analysis with mixeddesign ANOVAs that included between-subjects factors of age group (older or younger) and visual feedback condition (Vision+ or Vision− ) and a within-subjects factors of time-point (early or late in the block for adaptation data; test blocks P2, P3, P4 for proprioceptive data). Adaptation and washout blocks were analysed separately to examine finer aspects of group performance for these respective stages of the experiment. To probe the relationship between proprioception and adaptation, adaptation extent was separately correlated with baseline proprioceptive acuity and with change in bias after the adaptation (P2 subtracted from P3).
Since the baseline (P1) proprioceptive assessment was made prior to exposure to the visual feedback condition, these data were separately analysed in two-way ANOVAs that included only age group and physical activity status (active or inactive) as between-subjects factors. Physical activity indices were converted to z-scores within the two age groups, and used as a predictor of adaptation extent in multiple linear regression models that controlled for age and visual feedback condition grouping. Finally, spatial working memory scores were analysed in mixed ANOVAs with a between-subjects factor of age group and within-subjects factor of memory load (6, 8 or 10 tiles).
In an attempt to meet parametric assumptions whilst minimising data exclusion, we examined several approaches to deal with outlying data. Namely, we compared the number of significant Shapiro-Wilk tests reported for all subsets of data included in our statistical analyses, testing the raw data, log-transformed data, and data with outliers removed according to two different thresholds (3 and 2.5 standard deviations [SD] from the group mean at each repeated measure level). Whilst the log transformation and the 3 SD threshold outlier removal did reduce the number of significant Shapiro-Wilk tests compared to raw data, there were still a concerningly high number of non-normal distributions from these approaches. Conversely, the 2.5 SD threshold removal minimised the number of non-normal distributions with relatively low data exclusion overall (4.0% of all data), so we selected this approach for our analyses. All values are presented as group means ± standard deviation unless otherwise stated, and in all cases where the sphericity assumption was violated a Greenhouse-Geisser correction was used, with significance assessed at the p < 0.05 level. In cases of more than two multiple comparisons (post-hoc tests, correlations and regression models), we report p values which have been adjusted using a false discovery rate analysis (Benjamini and Hochberg 1995;Yekutieli and Benjamini 1999; also see use in Vandevoorde and Orban de Xivry 2019), denoted as p adj . The Bayes factor (BF 01 ) associated with correlations between measures from different tasks was also calculated using the program JASP (JASP Team 2020) to provide further evidence for (or against) non-significant Pearson's correlations (Wagenmakers et al. 2018b). All correlations were calculated as two-tailed tests using a uniform prior distribution. Values for BF 01 ranged from 0.51 (anecdotal evidence for the alternative hypothesis) to 4.49 (moderate evidence for the null hypothesis), though some were also close to 1 indicating equivocal evidence for both (Wagenmakers et al. 2018a).

Channel trials
Adaptation index values across the experiment (before outlier exclusion) are shown in Fig. 3a (Vision+ groups) and 3B (Vision− groups), with group adaptation index (normalised to baseline) for early and late portions of the adaptation and washout block shown in Fig. 3c and d, respectively. Channel trial data showed that adaptation index increased from early to late adaptation as expected (F[1, 61] = 75.09, p < 0.001, η 2 p = 0.55); however, there were no two-or three-way interactions involving timepoint (all p ≥ 0.356). The main effects and interactions of age and visual feedback on adaptation index were also non-significant (all p ≥ 0.456).

Physical activity and adaptation
We used physical activity z-scores alongside the category variables of age and visual feedback group to predict adaptation and washout extent in a series of multiple linear regression models. We focus mainly on the strength of the physical activity coefficient as a predictor in these models, with the categorical group variables included to control for some of the interaction effects reported in the previous section where we outline a more detailed analysis of these factors. The models showed that physical activity did not predict adaptation extent (deviation error modelβ = 0.12, t(62) = 0.97, p adj = 0.505; adaptation index modelβ = 0.12, t(62) = 0.89, p adj = 0.505) or washout extent (deviation error modelβ = 0.02, t(63) = 0.18, p adj = 0.859; adaptation index modelβ = − 0.20, t(64) = − 1.60, p adj = 0.455), with none of the models being significant overall (all R 2 ≤ 0.12, p adj ≥ 0.181). This indicates that physical activity was not strongly associated with adaptation.

Movement kinematics
Peak velocity of movements made over course of the experiment are shown in Fig. 4a (Vision+) and B (Vision−); peak Group peak velocity lateral deviation (PVLD) data across the experiment for older (purple) and younger (green) adults in a Vision+ condition (circles) and b Vision− condition (triangles). For visualisation purposes only, epochs of 4 trials were created for each participant and then averaged within groups to generate the group means (± 1SE) shown in panels a and b, with grey shaded regions representing peri-ods where the force-field was on. Lower panels show PVLD data (normalised to baseline) for the early and late portions of the adaptation block (c) and the washout block (d). In these panels, the means of individual participants are shown as symbols, with group averages (± 1SE) displayed in neighbouring bars. Significant post-hoc comparisons for group interactions also shown above respective panels velocity averaged across the adaptation block (when the force-field was switched on) is shown in Fig. 4c.
A two-way (age × visual feedback group) ANOVA conducted on average peak velocity in the adaptation block showed that movement speed in the force-field was comparable between all groups (age group -F[1, 65] = 1.47, p = 0.230; visual feedback group -F[1, 65] = 0.06, p = 0.804; age × visual feedback group -F[1, 65] = 0.17, p = 0.686). This indicates that force exposure was similar between groups during adaptation to the velocity-dependent field. Detailed statistical analyses of baseline-normalised peak velocity, movement duration and time to peak velocity across the different phases of the experiment can be found in the Supplementary 1. As a brief summary, participants dropped peak velocity in early adaptation, but it was otherwise similar across experimental phases and between groups (Supplementary Figure S1-A). Movements became longer in duration when the perturbation was introduced, then reduced back to baseline levels and remained stable, before finally becoming shorter in duration by the end of the experiment (Supplementary Figure  S1-B). Finally, time to peak velocity (expressed as percentage of movement duration) became shorter in early adaptation, recovering slightly by late adaptation, before returning to near baseline values for the remainder of the experiment. In addition, older adults tended to have shorter time to peak velocity relative to baseline than younger adults (Supplementary Figure S1-C).

Fig. 3
Change in adaptation index across the experiment for older (purple) and younger (green) adults in the a Vision+ condition (circles) and b Vision− condition (triangles). Lower panels show adapta-tion index data (normalised to baseline) for the early and late portions of the adaptation block (c) and the washout block (d). The format for each panel is as described for Fig. 2

Proprioceptive acuity during adaptation
The change in bias relative to baseline (P1) across the repeated proprioceptive assessments (P2-P4) is shown for all four groups in Fig. 6, with raw baseline performance violin plots inset for reference. A three-way (age group × visual feedback × repetition) ANOVA showed no effect of repetition on the baseline-normalised bias (F[2, 122] = 1.43, Fig. 4 Peak velocity data shown across the experiment for older (purple) and younger (green) adults in the a Vision+ condition (circles) and b Vision− condition (triangles). For visualisation, epochs of 4 trials were created for each participant and then averaged within groups to generate data points above (± 1SE), with grey shaded regions rep-resenting periods where the force-field was on. c Violin plots of peak velocity averaged across the adaptation block when the force-field was switched on (indicated by grey dashed box outline in a and b) for all four sub-groups  For normalised uncertainty range, there was no main effect of repetition (F[2, 124] = 0.22, p = 0.804) nor any two-or three-way interactions of repetition with either age or visual feedback group (all p ≥ 0.311). Furthermore, there were no main effects or interactions of age or visual feedback group (all p ≥ 0.532).
Thus, there were no consistent changes in proprioceptive bias or uncertainty after successful adaptation to the force-field.

Kinematics and proprioceptive acuity
After combining the raw, non-normalised data from all four proprioceptive assessments, we found that movement speed was not correlated with either uncertainty range (all |r| ≤ 0.17, p adj ≥ 0.685) or absolute bias (all |r| ≤ 0.15, p adj ≥ 0.717) for any of the four separate groups. This ruled out any effect of differences in proprioception due to selfselection of movement speed. Furthermore, whilst we did find that forces exerted against the channel walls were correlated with the bias for the younger vision-group, this effect did not survive adjustment for multiple comparisons (r = −0.27, p = 0.026, p adj = 0.105) and the correlations were not significant for the remaining three groups (all |r| ≤ 0.20, p adj ≥ 0.233). Thus, our estimates of proprioceptive acuity were largely independent of movement kinematics.

Proprioceptive recalibration and adaptation extent
Based on previous work (Ohashi et al. 2019;Ostry et al. 2010), we expected that larger perceptual shifts in the direction of the applied perturbing force field (indicated here by positive shift in bias) would be positively associated with greater adaptation extent. In fact, the correlations were of mixed directions and non-significant for all four groups (Supplementary Table S3; PVLD extent -all |r| ≤ 0.39, p adj ≥ 0.505, 1.13 ≤ BF 01 ≤ 3.34; adaptation index extent -all |r| ≤ 0.46, p adj ≥ 0.285, 0.72 ≤ BF 01 ≤ 3.05). Accordingly, we found no indication that proprioceptive recalibration was related to extent of adaptation.  Figure S3-B) and were significantly larger in older (5.12 ± 2.48) than younger adults (2.44 ± 2.10; F[1, 63] = 22.0, p < 0.001, η 2 p = 0.26). Post-hoc comparisons of the significant memory load × age group interaction (F[1.3, 82.4] = 8.26, p = 0.002, η 2 p = 0.12) showed that older adults made more BSEs than younger adults in the higher load conditions of 8-(t[49.5] = 4.72, p adj < 0.001) and 10 tiles (t[63] = 3.54, p adj = 0.001) specifically, whilst performance was similar between ages at the lowest memory load condition of 6 tiles (t[53.5] = 1.41, p adj = 0.165). Thus, older adults showed reduced spatial working memory capacity that became more pronounced with task complexity.

Working memory and baseline proprioceptive acuity
Total search errors were not correlated with absolute baseline bias for either older (r = −0.054, p = 0.762, BF 01 = 4.49) or younger adults (r = −0.182, p = 0.303, BF 01 = 2.82). However, there was a trend towards an association with baseline uncertainty range in older adults (r = 0.311, p = 0.08, BF 01 = 1.05) which was not seen in the younger group (r = 0.168, p = 0.343, BF 01 = 3.04). This suggests that SWM capacity was largely independent of dynamic proprioceptive errors, although, there was a tendency for older participants with lower SWM capacity to have lower perceptual acuity.

Discussion
We aimed to investigate the relationship of dynamic proprioception with adaptation to novel dynamic forces, in older and younger adults. We found that the level of adaptation was similar between the age groups, regardless of visual feedback limitations, although older adults showed increased after-effects of force production. Surprisingly, we found systematic (but not variable) proprioceptive errors were larger in physically inactive participants, regardless of age. However, there was no association between baseline proprioceptive acuity and adaptation, nor any consistent evidence of proprioceptive recalibration. Finally, despite clear agerelated deficits in spatial working memory capacity there was no relation of this cognitive measure with adaptation. Taken together, we find little evidence to support a relationship between dynamic upper limb proprioception and adaptation to novel field dynamics in either older or younger adults.
This study adds to the increasing, but still limited, number of reports indicating a minimal effect of ageing on reach adaptation to novel forces (Cesqui et al. 2008;Rajeshkumar and Trewartha 2019;Reuter et al. 2018;Trewartha et al. 2014). In addition, we demonstrate this finding holds true under conditions of limited visual feedback (providing only distance information during movement, and terminal error), conditions which might emphasise proprioceptive control. Whereas studies of force-field adaptation in the ageing population are still scarce, there is a wealth of research which has focused on visuomotor transformations in this group (usually a rotated or altered gain hand position feedback; Anguera et al. 2011;Bock 2005;Buch et al. 2003;Contreras-Vidal et al. 2002;Hegele and Heuer 2010;Seidler 2006;Vandevoorde and Orban de Xivry 2019). Although the magnitudes of age-specific deficits reported in these studies are mixed, a recurring observation is that performance is most impaired for older adults when the perturbation is salient and requires greater reliance on explicit strategies for adaptation (Buch et al. 2003;Cressman et al. 2010;Heuer 2010, 2013;McNay and Willingham 1998). Converging evidence suggests this is directly related to age-dependent cognitive decline and in particular spatial working memory capacity Uresti-Cabrera et al. 2015;Vandevoorde and Orban de Xivry 2019;Wolpe et al. 2020; for review see Seidler et al. 2010), which has been directly associated with the explicit component of adaptation (Christou et al. 2016). Despite clear evidence of reduced spatial working memory capacity in our sample of older adults, we did not find any association of this cognitive measure with adaptation. One interpretation of these data might, therefore, be that force-field adaptation relies less on explicit strategies than do visuomotor rotation tasks. Specifically, this could explain why older adults appear to rely more on implicit adaptation processes (Vandevoorde and Orban de Xivry 2019; Wolpe et al. 2020) but do not present notable deficits in force-field adaptation, which also appears to be unrelated to spatial working memory impairments. This is not to say that force-field adaptation relies only on implicit processes. Indeed, whilst better documented for visuomotor tasks (Benson et al. 2011;Mazzoni and Krakauer 2006;Neville and Cressman 2018;Taylor et al. 2014;Werner et al. 2015), recent work has sought to measure the explicit component of force-field adaptation (Schween et al. 2020) and the fast and slow processes of motor adaptation distinguished within force-field tasks (Smith et al. 2006) have been shown to be closely related to explicit and implicit learning on visuomotor tasks respectively (McDougle et al. 2015). However, the lack of any relationship with age or with age-related cognitive decline suggests explicit strategies may be reduced in force adaptation compared with visuomotor paradigms. Nevertheless, this interpretation of results was beyond the immediate scope of this study; further research will be necessary to test this hypothesis directly.
We did not observe any consistent, direction-specific shifts in perceived hand position following adaptation, regardless of age or visual feedback conditions. Sensory perceptual recalibration with adaptation has been reported elsewhere (Cressman and Henriques 2009;Modchalingam et al. 2019;Rossi et al. 2019;Sexton et al. 2019;Shiller et al. 2009), including proprioceptive shifts related to force-field adaptation (Haith et al. 2008;Mattar et al. 2013;Ohashi et al. 2019;Ostry et al. 2010). In the latter case, the magnitude of recalibration is typically in the range of 1-4 mm (cf. Haith et al. 2008 where this is closer to 20 mm), which means that factors like measurement sensitivity might contribute to difficulties in capturing these shifts empirically. Although different proprioceptive assessment methods can often provide distinct estimates of acuity (Elangovan et al. 2014;Hoseini et al. 2015), we specifically chose a task which emulates several previous studies that have reported perceptual recalibration after adaptation (Mattar et al. 2013;Ohashi et al. 2019;Ostry et al. 2010). Furthermore, we are confident in the sensitivity of our measures since we were able to detect an average 4 mm bias difference between physically inactive and active participants (Fig. 5b). We also note that proprioceptive recalibration has been observed despite variations of channel trajectory, movement type (passive instead of active) and staircase procedure (Ohashi et al. 2019), as well as completely different tasks altogether (Haith et al. 2008). Moreover, whilst we employed a slightly weaker force-field (15 N m −1 s −1 ) than other studies reporting recalibration effects (≥ 18 N m −1 s −1 ; Mattar et al. 2013;Ohashi et al. 2019;Ostry et al. 2010;Haith et al. 2008), the magnitude of movement error was in a similar range and stronger fields do not always evoke proprioceptive shifts (Sexton et al. 2018). Collectively, these results suggest that our choice of perceptual task parameters is not the primary basis for the absence of adaptation-dependent shifts in this experiment. Consequently, factors contributing to the presence (or absence) of proprioception recalibration with forcefield motor adaptation are still unclear and are likely made difficult to assess by the small scale on which recalibration occurs.
Previously, we reported that physically inactive older adults had larger systematic (but not variable) proprioceptive errors compared with a single group of younger adults (n = 20; Kitchen and Miall 2019), which we suggested may have resulted from use-dependent intrafusal fibre sparing in advanced age that biased the tuning of sensed limb position (Bergenheim et al. 2000;Jones et al. 2001). In this experiment, an increased sample size allowed us to assign both older and younger adults into physical activity subgroups, which revealed physically inactive participants have larger proprioceptive biases, regardless of age. Since selective loss of intrafusal fibres in low activity younger adults seems unlikely, an alternative explanation (as first suggested by Wilson et al. 2010) is that sedentary behaviour biases the experience of limb movements to a more laterally constrained region than is typically seen in naturalistic, active settings (Howard et al. 2009a). If so, perception may be drawn towards the location of highest or most frequent sensory experience as a prior (Gritsenko et al. 2007;Körding and Wolpert 2006), increasing systematic proprioceptive biases during periods of sensory uncertainty (as in the perceptual test). Further study with methods that permit measurements of everyday movements over extended periods of time will be necessary to test this directly.
We found that younger adults made larger lateral errors than older adults when the perturbation was first introduced (Fig. 2c). Since any systematic difference in movement velocity between the age groups disappeared early in the baseline phase ( Fig. 4; Supplementary Figure S1-A), this raises questions as to the basis for this age-dependent difference in movement error. One explanation could be that older adults increase muscle co-contraction during movement to reduce age-dependent motor variability (Huang and Ahmed 2014;. Increased co-contraction has been reported for older adults during force-field adaptation, but interestingly, was negatively correlated with adaptation extent (Huang and Ahmed 2014). Unlike Huang and Ahmed (2014), we did not observe age differences in adaptation. Accordingly, it could be that our older adults reduced their initially high co-contraction over time, as observed in younger adults (Heald et al. 2018;Huang et al. 2012;Huang and Ahmed 2014), thus minimising maladaptive consequences on their overall performance. Ability to modulate co-contraction and subsequent limb stiffness may, therefore, be a factor determining adaptation in advanced age.
An unexpected finding was that older adults showed increased after-effects of force production profiles during washout (Fig. 3d). This is particularly surprising given the absence of any major age effects on the adaptation phase of the experiment. Although sometimes overlooked, the washout phase or "de-adaptation" of sensorimotor adaptation tasks is considered an active process and not a simple switching-off or gradual forgetting of the behaviour (Davidson and Wolpert 2004) that can be modified with explicit knowledge and instruction (Benson et al. 2011;Mazzoni and Krakauer 2006). Prior work with older adults has indicated reduced retention of force-field adaptation with advanced age (Rajeshkumar and Trewartha 2019;Trewartha et al. 2014), however, retention was measured in these studies using error-clamped channel trials without vision, which differs from our null-field washout phase. Elsewhere, retention differences were either equivocal or not studied in detail (Cesqui et al. 2008;Huang and Ahmed 2014;Reuter et al. 2018). In visual perturbation tasks, after-effects are generally unaffected by age (Bock 2005; Hegele and Heuer 2010; see also Table 1 in Buch et al. 2003) although some studies have indicated age-specific increases in magnitude of after-effects (Fernández-Ruiz et al. 2000;Wolpe et al. 2020). This result has been suggested to depend on over-reliance on non-strategic adaptation that persists into washout (Fernández-Ruiz et al. 2000), but this notion is difficult to resolve with our data. If over-reliance on non-strategic processes was strong enough to present as increased after-effects for older adults, we would expect to see some evidence of it during the adaptation phase (as observed by Fernández-Ruiz et al. 2000), yet this wasn't the case. Furthermore, the absence of an association with working memory at least suggests that an age-dependent imbalance in the use of cognitive strategies was unlikely to have contributed to the effect. Accordingly, further study on factors which influence this aspect of performance in advanced age are necessary to better understand the basis of this novel finding.
In conclusion, we found minimal age-dependent differences in force-field adaptation and only limited evidence for proprioceptive recalibration with adaptation, regardless of visual feedback extent. We did find differences in proprioceptive acuity between high-and low-physical activity individuals, but this was across the age range. Dynamic proprioceptive acuity deficits at the non-clinical level (with healthy ageing or physical inactivity) do not, therefore, appear to influence force-field adaptation. Age-dependent cognitive decline has been closely linked with adaptation deficits in visuomotor tasks, but we found no association between spatial working memory capacity and adaptation performance. Hence, we suggest that force-field adaptation may be sufficiently weighted towards implicit processes to abolish these age-dependent performance impairments.