To empirically corroborate vendor-provided gradient nonlinearity (GNL) characteristics and demonstrate efficient GNL bias correction for human brain apparent diffusion coefficient (ADC) across 3T MR systems and spatial locations. Spatial distortion vector fields (DVF) were mapped in 3D using a surface fiducial array phantom for individual gradient channels on three 3T MR platforms from different vendors. Measured DVF were converted into empirical 3D GNL tensors and compared with their theoretical counterparts derived from vendor-provided spherical harmonic (SPH) coefficients. To illustrate spatial impact of GNL on ADC, diffusion weighted imaging using three orthogonal gradient directions was performed on a volunteer brain positioned at isocenter (as a reference) and offset superiorly by 10-17cm (>10% predicted GNL bias). The SPH tensor-based GNL correction was applied to individual DWI gradient directions, and derived ADC was compared with low-bias reference for human brain white matter (WM) ROIs. Empiric and predicted GNL errors were comparable for all three studied 3T MR systems, with <1.0% differences in the median and width of spatial histograms for individual GNL tensor elements. Median (±width) of ADC (10-3mm2/s) histograms measured at isocenter in WM reference ROIs from three MR systems were: 0.73±0.11, 0.71±0.14, 0.74±0.17, and at off-isocenters (before versus after GNL correction) were respectively 0.63±0.14 versus 0.72±0.11, 0.53±0.16 versus 0.74±0.18, and 0.65±0.16 versus 0.76±0.18. The phantom-based spatial distortion measurements validated vendor-provided gradient fields, and accurate WM ADC was recovered regardless of spatial locations and clinical MR platforms using system-specific tensor-based GNL correction for routine DWI.