In a simultaneous equations model, applying Ordinary Least Squares (D) directly to an equation with an endogenous regressor (here, Yt appears on the right side of equation (2), but Yt is itself jointly determined with Zt) produces biased and inconsistent estimates due to simultaneity bias — (D) is inappropriate.
Indirect Least Squares (B) is a valid method, but only works for exactly identified equations — it recovers structural parameters by solving the reduced-form OLS estimates back out, which requires a unique (one-to-one) mapping that only exists under exact identification. Since equation (2) is overidentified, ILS does not apply directly here.
Weighted Least Squares (C) addresses heteroscedasticity, not simultaneity — irrelevant to this specific problem.
Two-Stage Least Squares (2SLS) (A) is the standard, appropriate estimator for both exactly identified and overidentified equations in a simultaneous system — it uses the exogenous variables as instruments in a first stage to purge the endogeneity, then estimates the structural equation in a second stage.
Therefore the answer is (A) Two-Stage Least Squares.
The trap: (B) Indirect Least Squares
Students who recall that 'Indirect Least Squares is used for simultaneous equations' without the specific caveat that ILS is restricted to exactly-identified cases may pick (B) as a generically correct answer for 'simultaneous equations estimation,' missing the overidentification detail that specifically rules ILS out and points to 2SLS instead.