Understanding the spatial scale of local adaptation and the factors associated with adaptive diversity are important objectives for ecology and evolutionary biology, and have significant implications for effective conservation and management of wild populations and natural resources. In this study, we used an environmental association analysis to identify important bioclimatic variables correlated with putatively adaptive genetic variation in a benthic marine invertebrate-the giant California sea cucumber (Parastichopus californicus)-spanning coastal British Columbia and southeastern Alaska. We used a redundancy analysis (RDA) with 3,699 single nucleotide polymorphisms (SNPs) obtained using RAD sequencing to detect candidate markers associated with 11 bioclimatic variables, including sea bottom and surface conditions, across two spatial scales (entire study area and within subregions). At the broadest scale, RDA revealed 59 candidate SNPs, 86% of which were associated with mean bottom temperature. Similar patterns were identified when population structure was accounted for. Additive polygenic scores, which provide a measure of the cumulative signal across all candidate SNPs, were strongly correlated with mean bottom temperature, consistent with spatially varying selection across a thermal gradient. At a finer scale, 23 candidate SNPs were detected, primarily associated with surface salinity (26%) and bottom current velocity (17%). Our findings suggest that environmental variables may play a role as drivers of spatially varying selection for P. californicus. These results provide context for future studies to evaluate the genetic basis of local adaptation in P. californicus and help inform the relevant scales and environmental variables for in situ field studies of putative adaptive variation in marine invertebrates. K E Y W O R D S climate data, environmental association analysis, marine invertebrate, RAD-sequencing, redundancy analysis, seascape genomics