We consider the minimization of a sum of an expectation-valued coordinate-wise L i -smooth nonconvex function and a nonsmooth block-separable convex regularizer. Prior schemes are characterized by the following shortcomings: (a) Steplengths require global knowledge of Lipschitz constants; (b) Batchsizes of gradients are centrally updated and require knowledge of the global clock; (c) a.s. convergence guarantees are unavailable; (d) Rates are inferior compared to deterministic counterparts. Specifically, (a) and (b) require coordination across blocks and necessitate global information, leading to potentially larger constants in the rate and the oracle complexity bounds, impeding decentralized implementations, and resulting in relatively poor empirical behavior. We address these shortcomings by proposing an asynchronous variance-reduced algorithm, where in each iteration, a single block is randomly chosen to update its estimates by a proximal variable sample-size stochastic gradient scheme, while the remaining blocks are kept invariant. Notably, each block employs a steplength that is in accordance with its blockspecific Lipschitz constant while block-specific batch-sizes are random variables updated at a rate that grows either at a geometric or polynomial rate with the (random) number of times that block is selected. We show that every limit point for almost every sample path is a stationary point and establish the ergodic non-asymptotic rate O(1/K). Iteration and oracle complexity to obtain an -stationary point are shown to be O(1/ ) and O(1/ 2 ), respectively. Furthermore, under a µ-proximal Polyak-Łojasiewicz (PL) condition with the batch size increasing at a geometric rate, we prove that the suboptimality diminishes at a geometric rate, the optimal deterministic rate while iteration and oracle complexity to obtain an -optimal solution are proven to be O((L max /µ) ln(1/ )) and O (L ave /µ)(1/ ) 1+c with c ≥ 0, respectively. In the single block setting, we obtain the optimal oracle complexity bound O(1/ ). In pursuit of less aggressive sampling rates, when the batch sizes increase at a polynomial rate of degree v ≥ 1, suboptimality decays at a corresponding polynomial rate while the iteration and oracle complexity to obtain an −optimal solution are provably O(v(1/ ) 1/v ) and O e v v 2v+1 (1/ ) 1+1/v , respectively. Finally, preliminary numerics support our theoretical findings, displaying significant improvements over schemes where steplengths are based on global Lipschitz constants. * The authors are with the