Let n, b, and c be integers with \(n\ge 0\) . The generalized central trinomial coefficient \(T_n(b,c)\) is defined as the coefficient of \(x^n\) in the expansion of \((x^2+bx+c)^n\) . Let \(p\ge 5\) be a prime. In this paper, by means of several combinatorial identities, we establish new congruences of the form \( \sum _{k=0}^{p-1}(2k+1)^a \epsilon ^k H_k^j \frac{T_k(b,c)^2}{d^k} \) modulo \(p^3\) and \(p^2\) , where \(a=1,3\) , \(j=1,2\) , and \(\epsilon \in \{-1,1\}\) . Here \(H_n\) denotes the nth harmonic number and \(d:=b^2-4c\) satisfies \(p\not \mid d\) . As an illustration, in the case \(b=c=1\) , we obtain \( \sum _{k=0}^{p-1}(2k+1)H_k^2(-1)^k\frac{T_k^2}{3^k} \equiv -\frac{p}{2}\left( 3-\frac{3}{2}q_p(3)+\left( \frac{p}{3}\right) \right) \pmod {p^2}, \) where \(T_n\) denotes the central trinomial coefficient, \(\left( \frac{\cdot }{\cdot }\right) \) is the Legendre symbol, and \(q_p(3)\) is the Fermat quotient.