FURTHER EDIT: Here's a sketch of a proof of the first result. The space of Fourier seriesof weight 2 cusp forms for gamma_0 (9) has a basis of Eisenstein elements F,G, and Hlying in Z[[x^3]], xZ[[x^3]], and x^2 Z[[x^3]] respectively. In Z[[x]], F is congruent to1 mod 12 x^3. Furthermore the coefficient of x^n in G is sigma_1(n) when n is 1 mod 3, while the coefficient of x^n in H is (1/3)(sigma_1 (n)) when n is 2 mod 3.
Let C=x-8x^4+20x^7+.. be the Fourier expansion of the weight 4 form (eta(3z))^8 for gamma_0(9). A comparison of the coefficients of x^n for small n gives the identities C=FG-27H^2, and G^2=FH. So mod 3, C^2=G^2=H, and the coefficient of x^p in C^2, when p is a prime congruent to 2 mod 3 is, modulo 3, equal to (1/3)(sigma_1(p))=(p+1)/3. Now the cube of C^2 is the square of f(x^3), where f is the Fourier expansion of delta. It follows that mod 3,the coefficient of x^p in f^2 is (p+1)/3 when p is a prime congruent to 2 mod 3. This isprecisely 1. above.

