We present a tensor-network approach for two-dimensional strong-coupling QCD with staggered quarks at nonzero chemical potential. After integrating out the gauge fields at infinite coupling, the partition func-tion can be written as a full contraction of a tensor network consisting of coupled local numeric and Grassmann tensors. To evaluate the partition function and to compute observables, we develop a Grassmann higher-order tensor renormalization group method, specifically tailored for this model. During the coarsen-ing procedure, the blocking of adjacent Grassmann tensors is performed analytically, and the total number of Grassmann variables in the tensor network is reduced by a factor of two at each coarsening step. The coarse-site numeric tensors are truncated using higher-order singular value decompositions. The method is validated by comparing the partition function, the chiral condensate and the baryon density computed with the tensor method with exact analytical results on small lattices up to volumes of 4 x 4. For larger volumes, we present first tensor results for the chiral condensate as a function of the mass and volume, and observe that the chiral symmetry is not broken dynamically in two dimensions. We also present tensor results for the number density as a function of the chemical potential, which hint at a first-order phase transition.(c) 2022 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/). Funded by SCOAP3.