Skip to main content

Dynamic changes of region-specific cortical features and scalp-to-cortex distance: implications for transcranial current stimulation modeling



Transcranial current stimulation in rehabilitation is a fast-growing field featured with computational and biophysical modeling. Cortical features and scalp-to-cortex distance (SCD) are key variables for determining the strength and distribution of the electric field, yet longitudinal studies able to capture these dynamic changes are missing. We sought to investigate and quantify the ageing effect on the morphometry and SCD of left primary motor cortex (M1) and dorsolateral prefrontal cortex (DLPFC) in normal ageing adults and mild cognitive impairment (MCI) converters.


Baseline, 1-year and 3-year follow-up structural magnetic resonance imaging scans from normal ageing adults (n = 32), and MCI converters (n = 22) were drawn from the Open Access Series of Imaging Studies. We quantified the changes of the cortical features and SCDs of left M1 and DLPFC, including grey matter volume, white matter volume, cortical thickness, and folding. Head model was developed to simulate the impact of SCD on the electric field induced by transcranial current stimulation.


Pronounced ageing effect was found on the SCD of left DLPFC in MCI converters. The SCD change of left DLPFC from baseline to 3-year follow-up demonstrated better performance to discriminate MCI converters from normal ageing adults than the other morphometric measures. The strength of electric field was consequently decreased with SCD in MCI converters.


Ageing has a prominent, but differential effect on the region-specific SCD and cortical features in older adults with cognitive impairments. Our findings suggest that SCD, cortical thickness, and folding of the targeted regions could be used as valuable imaging markers when conducting transcranial brain stimulation in individuals with brain atrophy.


Ageing is a complex nonlinear process that has negative impacts on the structures and functions of central neural system (CNS) [1]. Over the last decades, brain morphometry, derived from structural magnetic resonance imaging (sMRI), has emerged as a robust measure to quantify the cortical features [2, 3], which can enhance the accuracy of clinical diagnosis and prognosis, and further facilitate brain-based interventions in elderly population, such as in Alzheimer’s disease (AD) [4].

Advanced age is consistently associated with the reduction in total brain volume and cortical thickness [5, 6]; however, evidence regarding regional cortical changes measured using different scales is inconsistent [3]. For instance, the volumetric MRI-based measures of prefrontal cortex could discriminate the effect of age from the early morphometric changes resulting from neurodegeneration instead of global brain volume [7,8,9]. Meanwhile, studies that focused on surface-based measures, such as cortical thickness and folding, demonstrate more promising results in identifying the adults with mild cognitive impairment (MCI) from normal ageing adults and dementia patients [10]. Indeed, increasing evidence has confirmed that cortical regions exhibit ageing-related changes in morphometric features that are no larger than those measured across the whole cortex [11], implying a lack of anatomical specificity of progressive brain changes during normal ageing. Especially, the grey matter of lateral prefrontal cortex (PFC), as the vulnerable cortical region to AD pathology, has shown an inverted U-shaped across time [12], suggesting that global brain atrophy may not depict the regional brain changes of increased vulnerability due to pathological ageing.

Importantly, identifications of regional cortical features that are affected during ageing also emphasize on the key step toward developing personalized therapy for neurodegenerative diseases. Recently, transcranial current stimulation that delivers a small amount of electric current to the scalp and modulates the cortical activities is increasingly considered as a safe and effective intervention for individuals suffering from age-related brain diseases, such as late-life depression [13], MCI [14], and stroke [15]. Clinically, left primary motor cortex (M1) and left dorsolateral prefrontal cortex (DLPFC) are two commonly used therapeutic targets. The cortical features of the stimulation targets have been found to be related to the discrepancies in treatment outcome [16]. Except for cortical features, scalp-to-cortex distance (SCD), as an important parameter, could prominently influence the focality and strength of the electric field induced by brain stimulation and lead to heterogeneous treatment outcome [17,18,19]. For example, our previous studies have shown that older adults with dementia have an increased SCD of left M1 and DLPFC; therefore, same stimulation output may be insufficient to induce the desired therapeutic response in older adults with brain atrophy [20, 21]. Notably, when adjusting the stimulation output with the SCD of left M1, the SCD-adjusted output can induce effective stimulation intensity and improve the treatment response accordingly [22].

As mentioned, the links between region-specific cortical features, SCD, and electric field support the complex biophysical nature in ageing brains [20, 23, 24]. Considering the continuum of ageing to ageing-related neurodegenerative diseases [25], whether ageing has similar or differential effect on region-specific cortical features and SCD in individuals with increased risk of developing dementia is still unclear. Hence, the main purpose of this study was to determine the ageing effect on the cortical features and SCD of left M1 and DLPFC in normal ageing adults and mild cognitive impairment converters. A second objective was to examine and quantify the effect of SCD on the electric field through computational simulation model.



For detecting the cortical changes, study participants were drawn from the longitudinal brain dataset of the Open Access Series of Imaging Studies (OASIS) ( [26]. The cases who were dementia patients or had a history of major neurologic or psychiatric disorders or serious cerebrovascular conditions were excluded from this study. Overall, we included 32 normal ageing (NA) adults and 22 mild cognitive impairment (MCI) converters who had valid baseline and follow-up assessments. A three-wave longitudinal investigation of global cognitive function and structural magnetic resonance imaging was conducted at baseline, 1 year, and 3 years.

Global cognitive function measured by Mini-Mental State Examination (MMSE) was assessed by the physicians at the Washington University Alzheimer Disease Research Center (ADRC). All cases participated in accordance with the guidelines of the Washington University Human Studies Committee. Approval for public sharing of the anonymized data was also specifically obtained at the study site. The demographics, in terms of age, sex, and years of education and MMSE score of the participants, were directly obtained from the OASIS dataset.

MRI data acquisition

Structural MRI (sMRI) data were acquired on a 1.5 T Vision scanner (Siemens, Erlangen, Germany) within a single session during which cushioning and a thermoplastic face mask were employed to minimize head movement. The T1-weighted magnetization prepared rapid gradient echo (MPRAGE) sequence was empirically optimized for grey–white matter contrast, with repetition time (TR) = 9.7 ms, echo time (TE) = 4.0 ms, inversion time = 20 ms, delay time = 200 ms, flip angle = 10°, resolution = 256 × 256 matrix (1 mm × 1 mm), slices = 128, and thickness = 1.25 mm.

Surface-based morphometry

BrainSuite 16a ( was employed for surface-based morphometry (SBM) analysis of cortical features [27]. BrainSuite is an automatic cortical surface identification integrated package with the optimized version of Brain Surface Extraction (BSE), which is widely used in dementia and ageing research [28, 29]. Following the standardized pipeline (Fig. 1), the pre-processing of SBM includes: (1) motion correction; (2) intensity normalization; (3) removal of non-brain voxels; (4) segmentation into grey matter (GM), white matter (WM), and cerebrospinal fluid (CSF) images; (5) tessellation of the GM/WM boundary, and automated topology correction. At each step, we visually checked the outputs and manually corrected when there are segmentation errors (i.e., non-brain tissue).

Fig. 1
figure 1

The pipeline of structural magnetic resonance imaging (sMRI) analysis: Path A represents the region-of-interest-based cortical mapping; Path B represents the measurement of scalp-to-cortex (SCD) and the construction of geometric simulation model

Volumetric measures

Based on the Automated Anatomical Labeling (AAL) template, the grey matter volume (GMV) of left M1 and left DLPFC was calculated individually. To adjust the heterogeneity of head size, Cendes method was used to correct the individual variance with total intracranial volume (TIV) through the formula [30]: corrected SV = (MBV × SV)/IBV. MBV is the mean brain volume in the group (a constant), SV is the regional grey matter volume, and IBV is the individual brain volume.

Surface-based measures

Cortical thickness was calculated by computing the average distance (in millimetres) between pial surface and GM/WM boundary at each vertex on cortical surface. Cortical folding was measured by gyrification index (GI), a ratio of inner surface area to an outer surface area that smoothly encloses the cerebral cortex [31]. As part of SBM analysis, advanced processing steps for calculating cortical thickness and surface area, including atlas registration, spherical surface map, and parcellation, were performed based on AAL template. The cerebral cortex was parcellated into 68 anatomical regions for determining the region-specific mean cortical thickness and surface area. Each cortical map was visually inspected for proper registration and normalization prior to further analysis.

Geometric measures

Region-specific scalp-to-cortex distance (SCD) was measured using the Brainsight neuronavigation system ( (Rogue Research, Montreal, Canada). Based on individual sMRI data, we first conducted the 3D curvilinear reconstruction of scalp and cortex, and then adjusted the MRI-to-head co-registration using the anterior commissure-posterior commissure (AC-PC) line in the Montreal Neurological Institute (MNI) space. After co-registration, the locations of the targets on cortex were labeled with MNI coordinates (x, y, z) (Fig. 1).

Based on reconstructed cortex, we identified and pinpointed the locations of left M1 and left DLPFC individually (Fig. 2a). The hand representation of left M1 was determined according to the MNI coordinates as [x = − 42, y = − 16, z = 68] [20, 21, 32], representing as the “hook sign” on sagittal plane. The location of left M1 was verified within the grey matter on the top of paracentral gyrus. Regarding the conservative landmark approximating the cytoarchitectonic definition of prefrontal junction, we targeted left DLPFC with the MNI coordinates as [x = − 46, y = 45, z = 38] [16, 20]. The location of left DLPFC was verified within the grey matter on the top of middle frontal gyrus (MFG).

Fig. 2
figure 2

Trajectories demonstrating the region-specific morphometric measures listed in old adults with different cognitive statuses. In each plot, blue line represents the results from normal ageing (NA) adults, and red line represents the results from MCI converters. Error bars represent the standard error (SEM). * indicates significant between-group differences (p < 0.05)

To better mimic the realistic brain stimulation, the corresponding locations of left M1 and DLPFC on the scalp were targeted in neuronavigation system by pointing the cursor to the scalp and then adjusting the orientation of the coil or electrode from the midline at 45°. Euclidean distance (Di), as a geometric index, was used to measure the distance between the two points locating on the scalp (xs, ys, zs) and the cortex (xc, yc, zc) with the following formula [20, 21, 33]:

Di = \(\sqrt{{\left({x}_{s}-{x}_{c}\right)}^{2}+{\left({y}_{s}-{y}_{c}\right)}^{2}+{\left({z}_{s}-{z}_{c}\right)}^{2}}\).

Simulation of electric field

To examine the effect of SCD on the E-field in different therapeutic protocols, 3D head models were created by SPHERES ( and the realistic volumetric approach to simulate transcranial electric simulation (ROAST) toolbox ( [34]. SPHERES is a stand-alone application that allows the considerations of arbitrary montages and the adjustments of brain parameters on a concentric sphere model by leveraging an analytical solution [35]. Regarding the most commonly used protocol in clinical practice and rehabilitation, we simulated the effect of SCD on the E-field induced by anodal tDCS over left DLPFC (i.e., F3 in international 10–20 system) (Fig. 3a).

Fig. 3
figure 3

Demonstration of region-specific scalp-to-cortex distance (SCD) in normal ageing (NA) adults and mild cognitive impairment (MCI) converters. a Importing individual structural MRI data to the neuronavigation system and constructing the scalp and cerebral cortex; b localizing the targets and measuring the SCD; c comparing the ageing effect on the SCD of left primary motor cortex (M1) in NA adults and MCI converters; d comparing the ageing effect on the SCD of left dorsolateral prefrontal cortex (DLPFC) in NA adults and MCI converters

To further quantify the effect of SCD on E-field, we calculated the entropy of the simulated E-field as described in information theory [36]. In brief, the entropy of the E-field induced by tDCS is quantified by calculating the pixel values within a two-dimensional centered region. Three geometric parameters of the grey-scale E-field include the following (Fig. 3c): represents the center of the E-field; represents the magnitude of the E-field; represents the edge of the E-field. The values of the three parameters were calculated and plotted using the Image Processing Toolbox embedded in MATLAB (

Statistical analysis

The differences of baseline cognitive performance and demographics in terms of age and sex were tested either with Chi-square (χ2) test for categorical variable or with independent sample t test for continuous variables. The comparisons of morphometric features were conducted using the code embedded in the BrainSuite 16a ( Multiple comparison correction was used by the above code using false discovery rate (FDR) estimation. Multilevel linear models were used to uncover the correlations from two levels: within group and between groups across different time points. Receiver-operating characteristic (ROC) analysis was used to evaluate the power of region-specific SCD in differentiating the seniors with different cognitive statuses. Pearson correlation coefficients were used to detect the relationship between morphometric measures, SCD and cognitive performance. The Chi-square test, pair-t test, multilevel linear models, and ROC analyses were performed using SPSS Statistics 24.0 (IBM, Armonk, NY).


Ageing effect on cortical features

As shown in Table 1, age, sex, years of education, global cognitive function, mean cortical thickness, and total intracranial volume (TIV) were comparable between the two groups at baseline. MCI converters had worse performance on MMSE than NA adults at 1-year follow-up (t = 2.91, p = 0.009) and 3-year follow-up (t = 4.23, p < 0.001) (baseline vs. 3-year follow-up).

Table 1 Baseline demographics, cognitive function, and global cortical features

Within each group, no significant ageing effect was found in region-specific morphometric measures between the time points. Furthermore, we compared the cortical thickness, GMV, WMV, and cortical folding of left M1 and left DLPFC between NA adults and MCI converters (Table 2). Figure 2 representatively illustrates the trajectories of global and region-specific morphometric measures, of which cortical changes, particularly volumetric atrophy, were found in NA adults, but inverse changes of GMV and WMV were found in MCI converters. No cortical thinning was found in the left M1 in NA adults, but noticeable thinning was observed in the left M1 and left DLPFC in MCI converters. Meanwhile, nonlinear changes of the cortical folding of left M1 were also captured between the two groups over time. Compared to NA adults, MCI converters showed a prominent increase in the GI of left M1 at 1-year follow-up (t = 2.344, p = 0.024).

Table 2 Longitudinal changes of cortical features of left M1 and DLPFC

Ageing effect on geometric features

Table 3 demonstrates the changes of region-specific SCD across the follow-up observations. Overall, similar patterns of SCD changes were detected in both groups. No ageing effect was observed on the SCDs of left M1 and left DLPFC in NA adults, but significant ageing effect was found on the SCD of left DLPFC in MCI converters (Baseline vs. 3-year follow-up: t = − 2.54, p = 0.02) (Fig. 3b). Moreover, the SCD change of left DLPFC was correlated with the decreased score of MMSE (r = − 0.571, p = 0.011), and significantly contributed to the global cognitive decline (R2 = 0.326, p = 0.011, 95% CI: − 1.539, − 0.234).

Table 3 Longitudinal changes of scalp-to-cortex distance of left M1 and DLPFC

ROC analysis

To classify the groups with different cognitive statuses, the value of the area under the ROC curve (AUC) was used to test the discriminant power of the changes of cortical features and SCD. Neither MMSE nor volumetric measures showed a significant discriminative power; while, surface-based measures, the change of cortical folding (i.e., GI) of left M1 from baseline to 3-year follow-up had better performance to differentiate MCI converters from NA adults (AUC = 0.71, p = 0.014). Moreover, geometric measures, the SCD change of left DLPFC from baseline to 3-year follow-up also could discriminate MCI converters from NA adults (AUC = 0.687, p = 0.028).

Simulation of SCD-associated E-field

Given the prominent effect of ageing on the SCD of left DLPFC in MCI converters, head models based on sMRI data were prepared for the E-field simulation of anodal tDCS by placing a rectangular 5 cm × 5 cm electrode on the scalp centered over left DLPFC (i.e., F3) (Fig. 3a). The SCD of left DLPFC was highlighted as a parameter of interest in SPHERES. In the head models, the isotropic conductivity values of the brain/non-brain tissues were adopted with default setting, including: 0.3 S/m for scalp, 0.03 S/m for skull, 2 S/m for CSF, and 0.03 S/m for cortical tissue [19,20,21].

Using the SCD of left DLPFC as the “depth toward scalp” in SPHERES, the spatial distribution of the SCD-associated E-field was markedly decreased in MCI converters at 3-year follow-up. To quantify the effect SCD change on the E-field induced by tDCS, the information entropy of left DLPFC was 6.69 × 103 at baseline (Fig. 4d) and 6.65 × 103 at 3-year follow-up (Fig. 4e) (t = 4.17, p < 0.001). The magnitude of the E-field was 19.82 × 103 at baseline and 14.02 × 103 at 3-year follow-up (t = 13.6, p < 0.001). The entropy of the edge of the E-field was 32.84 × 103 at baseline and 32.14 × 103 at 3-year follow-up (t = 6.57, p < 0.001).

Fig. 4
figure 4

Three-dimensional geometric simulation model of anodal transcranial current stimulation. a We first reconstructed the scalp and placed anodal electrode over left DLPFC (i.e., F3 in the international 10/20 system); b second, completed the simulation of SCD-associated E-field using SPHERES and ROAST toolbox, and c transferred the E-field to grey-scale images for entropy analysis. Finally, we calculated the entropy changes of simulated E-field in MCI converters at baseline (d) and 3-year follow-up (e)


Using the MRI scans collected over a 3-year observational period, we examined and quantified the ageing effect on region-specific morphometric and geometric measures in older adults with different cognitive statuses. Generally, different from global measures, longitudinal volumetric, geometric, and surface-based measures showed nonlinear trajectories as represented by region-specific structural changes in normal ageing adults and MCI converters over time. In the context of global brain atrophy, ageing has a pronounced, but differential effect on regional cortical features, especially on cortical thickness, folding, and scalp-to-cortex distance (SCD). Moreover, the SCD change of left DLPFC has shown a significant impact on the electric field induced by transcranial current stimulation.

For decades, a majority of the cross-sectional studies has reported the structural brain changes with a linear pattern of age-related global decline [37]. Given the distributed nature of structural changes across the adulthood [38, 39], the marked atrophy in parenchymal volumes (i.e., grey matter and white matter) was mostly reported in older adults relative to younger adults [40]. Of note, even in late adulthood, a steeper and more rapid cortical atrophy was found in old–old adults (+ 90 years) compared to young–old adults (60–66 years) [41, 42]. Subsequently, longitudinal studies appeared to confirm these structural brain changes occurring at different rates for different cortical regions and tissue types [3], of which a greater annualized decline was found in total white matter volume rather than grey matter volume during normal ageing [43]. Indeed, a consistent pattern of brain atrophy in parenchymal volumes was observed in normal ageing adults, but an inversed pattern was detected in MCI converters during a 3-year period. Additionally, our results lead directly to the hypothesis that not all cortical regions demonstrate the similar patterns of volumetric brain changes [44]. In several targeted investigations, the most prominent regional reduction was observed in the fronto-parietal neocortex [45], particularly in the grey matter of prefrontal cortex [7, 46]. In accordance with these findings, regional grey matter loss over time was identical in left M1 and DLPFC among normal ageing adults compared to MCI converters, indicating that accelerated regional volume loss is not simply an effect of ageing on brain.

Whereas earlier studies measured cortical changes almost exclusively in terms of volumetric measures, more recent work adopted with surface-based measure reveals a more complex picture of ageing brain, with region-specific cortical differences demonstrating exclusive and often nonlinear trajectories [3]. Similar to regional volumetric changes, ageing-related cortical thinning was also mainly observed in the fronto-parietal cortex, including prefrontal cortex and motor cortex, in normal ageing adults and dementia patients [3, 7]. The nonlinear trajectories of cortical atrophy are normally represented with two morphometric features: (1) cortical thinning and (2) sulcal widening. In our observations, an overall accelerated cortical thinning was found in both left M1 and DLPFC. However, a region-specific pattern of cortical thinning was only detected in MCI converters during a mid-term follow-up period, of which the reduced cortical thickness is disproportionately greater in left M1 than left DLPFC at 1-year follow-up, but the cortical thickness was decreased rapidly in left DLPFC than left M1 at 3-year follow-up. During ageing, neuron loss and cortical atrophy may lead to cortical thinning [7], while decreased cortical folding might mainly reflect the underlying white matter change [47, 48]. Recent evidence highlighted that the age-related differences in the sulcal width of the frontal and central sulci were significantly associated with differences in the adjacent parenchymal volumes in older adults [49] and AD patients [50]. It seems that the sulcal morphometry measured by gyrification index (i.e., folding) is closely related to the cortical features in frontal and parietal cortex. Indeed, the anti-correlated trajectories of gyrification index and regional white matter volume observed in MCI converters also support this assumption. Moreover, the change of cortical folding, not white matter, showed the highest discriminative value in differentiating MCI converters from normal controls, indicating that the differences in cortical folding are progressive and could be detected and monitored during the long-term management of neurodegenerative diseases.

Combined with surface-based measures, scalp-to-cortex distance, as a geometric measure in the field of transcranial brain stimulation, is first studied to detect the accelerated structural brain changes in the region with a high degree of folding. Interestingly, the trajectories of SCD and cortical folding were nearly paralleled in normal ageing adults, but anti-correlated in MCI converters. The results may reflect the complicated dynamic changes in cortical features during pathological ageing. Different from other morphometric measures, significant ageing effect was only found in the SCD of left DLPFC in MCI converters, of which the SCD change of left DLPFC from baseline to 3-year follow-up also showed a modest power to discriminate MCI converters from normal ageing adults. In addition, with SCD as the key parameter of interest, we constructed the head model of anodal transcranial direct current stimulation (tDCS) over left DLPFC. Both strength and spatial distribution of the SCD-associated E-field were prominently decreased in MCI converters at 3-year follow-up. Therefore, the geometric simulation model demonstrates a novel and useful way to describe and quantify the complex biophysical interactions between SCD and E-field.

Moreover, SCD may essentially “scale-up” the collection of MRI-based morphometric measures (i.e., scalar variables) to serve as a vector in transcranial brain stimulation therapies. For example, despite cortical thickness and SCD share the same measurement scale (i.e., mm), SCD can provide additional geometric information from 3D space, even four-dimensional space (i.e., time dimension). The dynamic changes of region-specific SCD and the simulated E-field featured with quantitative measures may open the way for optimizing the technical parameters of brain stimulation, such as tDCS electrode montages, and TMS coil configurations, when conducting the brain-based interventions in age- or disease-specific populations.

Implications for intervention

Understanding the trajectories of region-specific brain morphometry may allow us to harness the dynamic cortical changes to combat the cognitive decline during pathological ageing. Facilitating the scale-dependent cortical features, such as cortical thickness and scalp-to-cortex distance, shown to subserve brain-based interventions may expand the functional range of individuals contending with brain atrophy. Another potential application of MRI-based measures is to serve as an imaging marker to guide and inform the personalized rehabilitation and intervention.

Limitations and future directions

The findings in this study should be interpreted with caution due to its limitations, of which the major one is the low number of participants. Nevertheless, it is somewhat remarkable that we found relatively robust ageing effect on cortical features, which further alludes to region-specific brain changes in the context of global brain atrophy. Furthermore, we were unable to examine the associations between morphometric measures and the domain-specific cognitive function during normal and pathological ageing in the present study due to the limitations of OASIS dataset.

Future research will aim to examine the cortical features and the functioning of core domains of cognition concurrently, so as to obtain a specific picture on the pattern of brain morphometry in different types of age-related neurodegenerative diseases, such as frontotemporal dementia (FTD), Alzheimer’s disease (AD), and Parkinson’s disease (PD). This will help to construct disease-specific brain atlas and further expand the applications of MRI-based measures in tracking the progression of cognitive decline and providing valuable information facilitated the personalized rehabilitation. Also, methodologies such as closed-loop transcranial brain stimulation with combined imaging modalities (e.g., TMS-EEG, tDCS-EEG) may represent the next frontier of neurotherapeutics in decoding the complex interlinks between ageing, brain morphometry, and disease progression.


Structural MRI is a powerful non-invasive imaging modality that has enormous potential in studies of brain-based interventions. The longitudinal analyses of the region-specific cortical features in MCI converters help to shift towards greater awareness of nonlinear dynamic changes, indicating that the differential effect of ageing on cortical folding and scalp-to-cortex distance most likely occur in the course of pathological ageing. Furthermore, scalp-to-cortex distance demonstrates a promising and valuable geometric imaging marker for monitoring the cortical changes and optimizing the brain stimulation therapies for the individuals with brain atrophy.

Availability of data and materials

The dataset supporting the conclusions of this article can be found at the official website of the OASIS project:



Automated Anatomical Labeling


Alzheimer’s disease


Area under the ROC curve


Brain surface extraction


Central neural system


Dorsolateral prefrontal cortex


Gyrification index


Grey matter volume


Primary motor cortex


Mild cognitive impairment


Mini-Mental State Examination


Magnetic resonance imaging


Normal ageing


Open Access Series of Imaging Studies


Prefrontal cortex


Receiver-operating characteristic


Surface-based morphometry


Scalp-to-cortex distance


Transcranial direct current stimulation


Total intracranial volume


Transcranial magnetic stimulation


White matter volume


  1. Cole JH, Franke K. Predicting age using neuroimaging: innovative brain ageing biomarkers. Trends Neurosci. 2017;40:681–90.

    Article  CAS  PubMed  Google Scholar 

  2. Chandra A, Dervenoulas G, Politis M, Alzheimer’s Disease Neuroimaging Initiative, Magnetic resonance imaging in Alzheimer’s disease and mild cognitive impairment. J Neurol. 2019;266:1293–302.

    Article  PubMed  Google Scholar 

  3. Fjell AM, Westlye LT, Grydeland H, Amlien I, Espeseth T, Reinvang I, Initiative ADN. Accelerating cortical thinning: unique to dementia or universal in aging? Cereb Cortex. 2014;24:919–34.

    Article  PubMed  Google Scholar 

  4. Madan CR, Kensinger EA. Cortical complexity as a measure of age-related brain atrophy. NeuroImage. 2016;134:617–29.

    Article  PubMed  Google Scholar 

  5. Fjell AM, Walhovd KB, Fennema-Notestine C, McEvoy LK, Hagler DJ, Holland D, Dale AM. One-year brain atrophy evident in healthy aging. J Neurol. 2009;29:15223–31.

    CAS  Google Scholar 

  6. Hogstrom LJ, Westlye LT, Walhovd KB, Fjell AM. The structure of the cerebral cortex across adult life: age-related patterns of surface area, thickness, and gyrification. Cereb Cortex. 2013;23:2521–30.

    Article  PubMed  Google Scholar 

  7. Salat DH, Kaye JA, Janowsky JS. Prefrontal gray and white matter volumes in healthy aging and Alzheimer disease. Arch Neurol. 1999;56:338–44.

    Article  CAS  PubMed  Google Scholar 

  8. Thompson PM, Hayashi KM, De Zubicaray G, Janke AL, Rose SE, Semple J, Toga AW. Dynamics of gray matter loss in Alzheimer's disease. J Neurol. 2003;23:994–1005.

    CAS  Google Scholar 

  9. Dickstein DL, Kabaso D, Rocher AB, Luebke JI, Wearne SL, Hof PR. Changes in the structural complexity of the aged brain. Aging Cell. 2007;6:275–84.

    Article  CAS  PubMed  Google Scholar 

  10. Ruiz de Miras J, Costumero V, Belloch V, Escudero J, Ávila C, Sepulcre J. Complexity analysis of cortical surface detects changes in future Alzheimer's disease converters. Hum Brain Mapp. 2017;38:5905–18.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Cox SR, Bastin ME, Ritchie SJ, Dickie DA, Liewald DC, Maniega SM, Corley J. Brain cortical characteristics of lifetime cognitive ageing. Brain Struct Funct. 2018;223:509–18.

    Article  PubMed  Google Scholar 

  12. Douaud G, Groves AR, Tamnes CK, Westlye LT, Duff EP, Engvig A, Matthews PM. A common brain network links development, aging, and vulnerability to disease. Proc Natl Acad Sci U S A. 2014;111:17648–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Wong HL, Chan WC, Wong YL, Wong SN, Yung HY, Wong SMC, Cheng PWC. High-definition transcranial direct current stimulation-An open-label pilot intervention in alleviating depressive symptoms and cognitive deficits in late-life depression. CNS Neurosci Ther. 2019;25:1244–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Lu H, Chan SSM, Chan WC, Lin C, Cheng CPW, Lam LCW. Randomized controlled trial of TDCS on cognition in 201 seniors with mild neurocognitive disorder. Ann Clin Transl Neurol. 2019;6:1938–48.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Coppens MJ, Staring WH, Nonnekes J, Geurts AC, Weerdesteyn V. Offline effects of transcranial direct current stimulation on reaction times of lower extremity movements in people after stroke: a pilot cross-over study. J Neuroeng Rehabil. 2019;16:1–10.

    Article  Google Scholar 

  16. Fox MD, Buckner RL, White MP, Greicius MD, Pascual-Leone A. Efficacy of transcranial magnetic stimulation targets for depression is related to intrinsic functional connectivity with the subgenual cingulate. Biol Psychiatry. 2012;72:595–603.

    Article  PubMed  PubMed Central  Google Scholar 

  17. McConnell KA, Nahas Z, Shastri A, Lorberbaum JP, Kozel FA, Bohning DE, George MS. The transcranial magnetic stimulation motor threshold depends on the distance from coil to underlying cortex: a replication in healthy adults comparing two methods of assessing the distance to cortex. Biol Psychiatry. 2001;49:454–9.

    Article  CAS  PubMed  Google Scholar 

  18. Lee EG, Duffy W, Hadimani RL, Waris M, Siddiqui W, Islam F, Jiles DC. Investigational effect of brain-scalp distance on the efficacy of transcranial magnetic stimulation treatment in depression. IEEE Trans Magn. 2016;52:1–4.

    Google Scholar 

  19. Huang Y, Liu AA, Lafon B, Friedman D, Dayan M, Wang X, Parra LC. Measurements and models of electric fields in the in vivo human brain during transcranial electric stimulation. Elife. 2017;6:e18834.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Lu H, Chan SS, Lam LC. Localized analysis of normalized distance from scalp to cortex and personalized evaluation (LANDSCAPE): focusing on age-and dementia-specific changes. J Alzheimers Dis. 2019;67:1331–41.

    Article  PubMed  Google Scholar 

  21. Lu H, Lam LC, Ning Y. Scalp-to-cortex distance of left primary motor cortex and its computational head model: implications for personalized neuromodulation. CNS Neurosci Ther. 2019;25:1270–6.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Herbsman T, Forster L, Molnar C, Dougherty R, Christie D, Koola J, Nahas Z. Motor threshold in transcranial magnetic stimulation: the impact of white matter fiber orientation and skull-to-cortex distance. Hum Brain Mapp. 2009;30:2044–55.

    Article  PubMed  Google Scholar 

  23. Stokes MG, Chambers CD, Gould IC, Henderson TR, Janko NE, Allen NB, Mattingley JB. Simple metric for scaling motor threshold based on scalp-cortex distance: application to studies using transcranial magnetic stimulation. J Neurophysiol. 2005;94:4520–7.

    Article  PubMed  Google Scholar 

  24. Hanlon CA, Dowdle LT, Moss H, Canterberry M, George MS. Mobilization of medial and lateral frontal-striatal circuits in cocaine users and controls: an interleaved TMS/BOLD functional connectivity study. Neuropsychopharmacology. 2016;41:3032–41.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Aisen PS, Cummings J, Jack CR, Morris JC, Sperling R, Frölich L, Scheltens P. On the path to 2025: understanding the Alzheimer’s disease continuum. Alzheimers Res Ther. 2017;9:60.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  26. Marcus DS, Wang TH, Parker J, Csernansky JG, Morris JC, Buckner RL. Open Access Series of Imaging Studies (OASIS): cross-sectional MRI data in young, middle aged, nondemented, and demented older adults. J Cogn Neurosci. 2007;19:1498–507.

    Article  PubMed  Google Scholar 

  27. Shattuck DW, Mirza M, Adisetiyo V, Hojatkashani C, Salamon G, Narr KL, Toga AW. Construction of a 3D probabilistic atlas of human cortical structures. Neuroimage. 2008;39:1064–80.

    Article  PubMed  Google Scholar 

  28. Apostolova LG, Thompson PM, Rogers SA, Dinov ID, Zoumalan C, Steiner CA, Cummings JL. Surface feature-guided mapping of cerebral metabolic changes in cognitively normal and mildly impaired elderly. Mol Imaging Biol. 2010;12:218–24.

    Article  PubMed  Google Scholar 

  29. Lu H, Ma SL, Chan SSM, Lam LCW. The effects of apolipoprotein ε 4 on aging brain in cognitively normal Chinese elderly: a surface-based morphometry study. Int Psychogeriatr. 2016;28:1503–11.

    Article  PubMed  Google Scholar 

  30. Cendes F, Andermann F, Gloor P, Evans A, Jones-Gotman M, Watson C, Leroux G. MRI volumetric measurement of amygdala and hippocampus in temporal lobe epilepsy. Neurology. 1993;43:719–25.

    Article  CAS  PubMed  Google Scholar 

  31. Cao B, Mwangi B, Passos IC, Wu MJ, Keser Z, Zunta-Soares GB, Soares JC. Lifespan gyrification trajectories of human brain in healthy individuals and patients with major psychiatric disorders. Sci Rep. 2017;7:1–8.

    Article  CAS  Google Scholar 

  32. Mylius V, Ayache SS, Ahdab R, Farhat WH, Zouari HG, Belke M, Schmidt S. Definition of DLPFC and M1 according to anatomical landmarks for navigated brain stimulation: inter-rater reliability, accuracy, and influence of gender and age. Neuroimage. 2013;78:224–32.

    Article  CAS  PubMed  Google Scholar 

  33. Randazzo MJ, Kondylis ED, Alhourani A, Wozny TA, Lipski WJ, Crammond DJ, Richardson RM. Three-dimensional localization of cortical electrodes in deep brain stimulation surgery from intraoperative fluoroscopy. Neuroimage. 2016;125:515–21.

    Article  PubMed  Google Scholar 

  34. Huang Y, Datta A, Bikson M, Parra LC. Realistic volumetric-approach to simulate transcranial electric stimulation—ROAST—a fully automated open-source pipeline. J Neural Eng. 2019;16:056006.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Truong DQ, Hüber M, Xie X, Datta A, Rahman A, Parra LC, Bikson M. Clinician accessible tools for GUI computational models of transcranial electrical stimulation: BONSAI and SPHERES. Brain Stimul. 2014;7:521–4.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Allen B, Stacey BC, Bar-Yam Y. Multiscale information theory and the marginal utility of information. Entropy. 2017;19:273.

    Article  Google Scholar 

  37. Meunier D, Stamatakis EA, Tyler LK. Age-related functional reorganization, structural changes, and preserved cognition. Neurobiol Aging. 2014;35:42–544.

    Article  PubMed  Google Scholar 

  38. Fjell AM, Westlye LT, Amlien I, Espeseth T, Reinvang I, Raz N, Dale AM. High consistency of regional cortical thinning in aging across multiple samples. Cereb Cortex. 2009;19:2001–122.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Giorgio A, Santelli L, Tomassini V, Bosnell R, Smith S, De Stefano N, Johansen-Berg H. Age-related changes in grey and white matter structure throughout adulthood. Neuroimage. 2010;51:943–51.

    Article  PubMed  Google Scholar 

  40. Dotson VM, Szymkowicz SM, Sozda CN, Kirton JW, Green ML, O’Shea A, Woods AJ. Age differences in prefrontal surface area and thickness in middle aged to older adults. Front Aging Neurosci. 2016;7:250.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Shaw ME, Sachdev PS, Anstey KJ, Cherbuin N. Age-related cortical thinning in cognitively healthy individuals in their 60s: the PATH Through Life study. Neurobiol Aging. 2016;39:202–9.

    Article  PubMed  Google Scholar 

  42. Yang Z, Wen W, Jiang J, Crawford JD, Reppermund S, Levitan C, Trollor JN. Age-associated differences on structural brain MRI in nondemented individuals from 71 to 103 years. Neurobiol Aging. 2016;40:86–97.

    Article  PubMed  Google Scholar 

  43. Thambisetty M, Wan J, Carass A, An Y, Prince JL, Resnick SM. Longitudinal changes in cortical thickness associated with normal aging. Neuroimage. 2010;52:1215–23.

    Article  PubMed  Google Scholar 

  44. Driscoll I, Davatzikos C, An Y, Wu X, Shen D, Kraut M, Resnick S. Longitudinal pattern of regional brain volume change differentiates normal aging from MCI. Neurology. 2009;72:1906–13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Terribilli D, Schaufelberger MS, Duran FL, Zanetti MV, Curiati PK, Menezes PR, Busatto GF. Age-related gray matter volume changes in the brain during non-elderly adulthood. Neurobiol Aging. 2011;32:354–68.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Raz N, Gunning FM, Head D, Dupuis JH, McQuain J, Briggs SD, Acker JD. Selective aging of the human cerebral cortex observed in vivo: differential vulnerability of the prefrontal gray matter. Cereb Cortex. 1997;7:268–82.

    Article  CAS  PubMed  Google Scholar 

  47. Van Essen DC. A tension-based theory of morphogenesis and compact wiring in the central nervous system. Nature. 1997;385:313–8.

    Article  PubMed  Google Scholar 

  48. Lu H. Quantifying age-associated cortical complexity of left dorsolateral prefrontal cortex with multiscale measurements. J Alzheimers Dis. 2020;76:505–12.

    Article  PubMed  Google Scholar 

  49. Jin K, Zhang T, Shaw M, Sachdev P, Cherbuin N. Relationship between sulcal characteristics and brain aging. Front Aging Neurosci. 2018;10:339.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Bertoux M, Lagarde J, Corlier F, Hamelin L, Mangin JF, Colliot O, Sarazin M. Sulcal morphology in Alzheimer's disease: an effective marker of diagnosis and cognition. Neurobiol Aging. 2019;84:41–9.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors would like to thank the Principal Investigators of OASIS (Longitudinal): D. Marcus, R, Buckner, J. Csernansky, and J. Morris for providing clinical and MRI data. The OASIS (Longitudinal) was supported by NIH grants P50 AG05681, P01 AG03991, P01 AG026276, R01 AG021910, P20 MH071616, and U24 RR021382. The funders had no involvement in the study design, data collection and analysis, or the preparation of the article. The authors would like to express their gratitude to all the staffs from the Chan Wei Wei Therapeutic Physical Mental Exercise Centre. The authors would thank the Brain X ( for providing the platform to share the processed MRI data with open science community. The authors also thank all the reviewers and editors for their valuable comments to improve the scientific value of this paper.

Data used in preparation of this article were obtained from the Open Access Series of Imaging Studies (OASIS) ( database. The investigators at OASIS contributed to the design and implementation of OASIS and/or provided data, but did not participate in analysis or writing of this paper.

PI of the Open Access Series of Imaging Studies: Daniel Marcus:


This work was supported by the Direct Grant of the Chinese University of Hong Kong (Project Number: 2018.077).

Author information

Authors and Affiliations




HL was responsible for designing the study and preparing the manuscript. HL and JL performed data analyses. All authors were involved in data interpretation. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hanna Lu.

Ethics declarations

Ethics approval and consent to participate

For the purpose of this study, we used the structural images of the OASIS that were previously collected under several study protocols at the Washington University. This study was carried out in accordance with the recommendations of the University’s Institutional Review Board (IRB). The protocol was approved by the University’s IRB. All participants gave written informed consent in accordance with the Declaration of Helsinki. All participants were given written informed consent at the time of study participation. The University’s IRB also provided explicit approval for open sharing of the anonymized data.

Consent for publication

All participants from the OASIS have given written informed consent to publish these data. No individual details, images, or videos are included in the manuscript.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lu, H., Li, J., Zhang, L. et al. Dynamic changes of region-specific cortical features and scalp-to-cortex distance: implications for transcranial current stimulation modeling. J NeuroEngineering Rehabil 18, 2 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: